Abstract

Parkinson's disease (PD) is the second most common neurodegenerative disease in the United States. The cardinal symptoms of PD are tremor, rigidity, slowed movement, and impaired balance. These symptoms often interfere with the daily activities of people with Parkinson's (PwPD) and negatively affect quality of life (QoL). Therefore, monitoring PD symptoms is essential for clinical evaluations and adjusting medication to help maintain QoL for PwPD. We are developing a mobile app to conduct at-home PD symptom monitoring to provide more timely, frequent, and accurate measurements of PD symptoms. While the tremor and finger-tapping results collected in the mobile app have been discussed in previous publications, this paper focuses on the design and evaluation of postural stability tests in the app and validating the reliability of the embedded accelerometers and gyroscopes in smartphones. During the test, a shaker was employed to provide vibration in amplitude and frequency ranges similar to human postural stability signals, and both the accelerometer and gyroscope measurements were evaluated. We used signal processing algorithms to extract postural stability factors, such as the root-mean-square (RMS) value, the derivative of acceleration, and frequency factors for the accelerations, and the ranges and RMS for the angular velocity. Our findings show that smartphone devices have good consistency over multiple trials and between devices, and motion patterns achieved from multiple data points are reliable for postural stability analysis.

1 Introduction

Parkinson's disease (PD) is a common neurodegenerative disease that, typically, is diagnosed in older adulthood (>65 years of age) [1]. PD affects over a million people in the United States [2]. It is predicted that the number of people with PD (PwPD) will nearly double by 2030 [2].

Parkinson's disease is associated with a wide range of heterogeneous symptoms. Symptoms include nonmotor symptoms that affect a person's memory and ability to focus, plan, and follow directions (executive functioning) [3], as well as a range of motor symptoms, including difficulty initiating movements (akinesia), slowness of movement (bradykinesia), tremors, stiffness, and muscle rigidity [4]. These nonmotor and motor symptoms manifest themselves in a unique range of severity in PwPD and significantly interfere with many daily activities involving fine motor skills, such as difficulty in writing, typing, eating, drinking, and carrying out tasks of daily living. A loss of gross motor skill functioning can also lead to increased difficulty with walking (gait) and balance, which we focus on in this paper. Balance impairments are particularly challenging for PwPD. Both motor and nonmotor symptoms are present because PD disrupts a person's ability to coordinate information between the nervous system and its corresponding muscular effectors [5]. The result is that voluntary movement is inhibited in PwPD and they experience a loss of control in the precision and accuracy of their movements as well as irregular movements in response to attempting to recover their balance [6].

Problems with balance can have devastating effects on PwPD, with the most devastating effect being falls. Falls are especially dangerous for PwPD who are typically older adults [1]. Injuries, such as a broken hip, can be life-altering. Statistics show that if a person over 75 years of age breaks their hip, they have a 50% chance of dying within 6 months of their injury [7]. Quality of life (QoL) has also been shown to diminish with falls in older adults and PwPD, limiting their ability to function independently [8].

Postural stability is the ability of a person to control their bodily position in space for the purpose of balance and movement. Postural stability can be quantified with motion measures as the postural sway that occurs in a person in a medial-lateral (ML) plane and an anterior–posterior (AP) plane. This measure is important to track for PwPD because an increased measure of postural instability is correlated with an increased risk of falls [6]. Also, for PwPD who display increasing postural instability, that measure could be used to indicate worsening symptoms [9]. Meanwhile, wearable devices and smartphones have become popular and researchers started investigating using smartphone/wearable devices in measuring motions for medical-related information [1012], including those for PwPD [13,14].

For this reason, we are developing a mobile application to empower PwPD with the ability to track their symptoms through a smartphone which, in turn, can potentially help with fall prevention for PwPD. People who present with increased postural instability would be able to track their progression and discuss it with their healthcare providers in a timely manner to proactively avoid harmful consequences. Additionally, our mobile symptom tracking application has great potential in meeting the urgent needs of the increasing number of PwPD among the aging population [1,15].

Our current research employs the use of a shaker to test the sensitivity and accuracy of sensors embedded in smartphones in measuring postural instability. We selected sample data of postural stability from a dataset that we previously collected from one healthy, male adult and used it as input signals for the shaker, simulating different conditions of a participant. These postural stability tests, such as the Romberg test, were chosen based on clinical tests for neurological diseases. The Romberg test is a proven clinical method to assess proprioception in cases of neurological diseases. It is designed to identify abnormalities in postural stability with the participant assuming certain positions with their eyes open or closed while the physician observes their ability to control their movement to balance [16]. We employed a modified version of this test in our research—the participant assumed positions with (1) eyes open while standing on both legs (“eyes open”), (2) eyes closed while standing on both legs (“eyes closed”), or (3) eyes open while standing on one leg (“one leg”). In this way, a test commonly used to assess balance will also aid in the testing of our technology for PwPD.

The input dataset was run multiple times using the shaker with two different smartphones. The time–domain signals collected were analyzed to evaluate the accuracy and precision of the measurements captured by the smartphones' embedded accelerometers and gyroscopes for predicting postural instability.

2 Materials and Methods

2.1 Experiment Setup.

In the experiment, two smartphones were selected to conduct the measurement. We chose the smartphones based on affordability, as we thought these phones would reflect the types of phones used by a majority of PwPD. In the smartphones ranging between $150 and 200 available in the market, we chose Samsung Galaxy A31 and Motorola Mogo G Power (2021).

The reliability test of smartphone accelerometers and gyroscopes was conducted by measuring mechanical vibration that mimicked postural stability tests. The vibration stimulus was introduced by one-dimensional input signals from the shaker and two different setups (Figs. 1 and 2) introduced vibration in three-dimensional (3D), in translation and rotation dominate motions. The amplitude of vibration provided at the touch point to the smartphone was adjusted to the range of postural stability at 0.1–0.2 m/s2. The equipment used included the Modal Shop SmartAmp 2100E21 and shaker amplifier.

Fig. 1
Translational motion setup
Fig. 1
Translational motion setup
Close modal
Fig. 2
Rotational motion setup
Fig. 2
Rotational motion setup
Close modal

While the shaker amplifier could be controlled by a knob, the vibration was adjusted at the input waveform and/or amplitude to the desired value. Two experimental setups that introduced vibration in 3D were designed to test the translational and rotational motions collected by the accelerometers and gyroscopes, respectively. The translational motion setup included other physical components such as a smartphone holder connected to the plate of the shaker (Fig. 1). The design was inspired by both the range of height and mobility of the shaker and the consistent efficiency required to accurately conduct the research. The phone stand was fastened to the shaker head expander with a screw and a washer, reinforced with tape. As shown in Fig. 1, while the shaker could only provide displacement in one direction, the phone stand was set at a small angle to the shaker plate; thus, the smartphone would receive vibrational motion on different axes, which can help in reliability tests and combat the common problem one may occur when collecting data from multiple axes. The mobile app was built in such a way that the smartphone is to be vertically attached to the back of the user; therefore x, y, and z-axis corresponded to the motion at Medial-Lateral (ML), Resultant Distance (RD), and AP in postural stability. During the experiment, the shaker provided excitation primarily in the AP plane, and the overall phone stand setup system also introduced movement in the ML and RD planes. Therefore, the motion patterns in 3D could be tested. Once the desired vibration amplitude was met, each smartphone was put on the stand for tests in separate trials, and the mobile app was open with a button-press to start recording measurements.

Different signal types were used during the experiment. First, a sweep signal was sent to the shaker to test the smartphone measurement ability in the targeted frequency range of 0.01–5 Hz, results shown in Fig. 3. Second, three signal samples previously recorded by a healthy male adult volunteer were used as input to the shaker, in the conditions of balancing with “eyes open,” “eyes closed,” and “one leg.” The “eyes open” and “eyes closed” signals had amplitude up to 0.12 m/s2 along time. The “one leg” signal was more dynamic as shown in Fig. 4. It started with ∼ 0.2 m/s2 amplitude, followed by a peak up to 0.5 m/s2, then back to lower amplitude, indicating that there was a perturbation in balance that the human body corrected this perturbation. All input signals were at a 50-Hz sampling rate, each lasted 20 s and was played repeatedly by the shaker. The shaker provided primarily the AP motions by vibration in z direction, yet the setup of the smartphone holder with angles also introduced motions in RD and ML in x and y directions, respectively.

Fig. 3
Measurement of sweep signal at two devices—translational motion setup
Fig. 3
Measurement of sweep signal at two devices—translational motion setup
Close modal
Fig. 4
Input signal to the shaker and measurements at two devices
Fig. 4
Input signal to the shaker and measurements at two devices
Close modal

The rotational motion setup was designed to provide mainly rotational motion and allow the mobile application to measure both the acceleration by accelerometers and rotational information by gyroscopes. As shown in Fig. 2, a curved aluminum plate was attached between the smartphone and the holder. The bottom of the plate was screwed to the stand, and its top was connected to the shaker with a stinger. The shaker provided the same input linear motion to the touchpoint as in the translational motion setup. This allowed the smartphone to move with the aluminum plate and rotate around the bottom of the stand, with a major translational motion along the z-axis (AP direction) as well as rotation around the x-axis. All the input signals at the shaker were the same as the translational motion setup in Fig. 1, while the measurement was collected by both the accelerometers and gyroscopes.

2.2 Signal Processing Methods.

The mobile app recorded 60 s in each trial with a 50 Hz sampling rate built in the app. In the results, device 1 (one smartphone) achieved a 50.05 Hz sampling rate, and device 2 (another smartphone) achieved a 52.08 Hz sampling rate. During the experiment, in order to minimize any influence caused by the finger touch on the smartphone screen to start the app and the time difference between starting the app and the shaker, we first started the app and waited for a few seconds, then started the shaker. We decided to use the first 20-s signal once the shaker started. The target signal was extracted with a rectangular window and the DC component was removed as the first step in signal processing.

The x, y, and z directions were used to describe ML, RD, and AP, respectively. Besides the direct comparison of time and frequency domain data, we also processed the data based upon the methods of Mancini et al. [17] and Prieto et al. [18]. We used their algorithm and applied to our acceleration measures and the signal processing algorithms are as followed.

A 3.5 Hz low-pass Butterworth filter was applied to extract signals representing postural stability motions. The acceleration a is the vector that combines the ax and az, amag is the magnitude (absolute value) of a. Time–domain factors, including JERK, DIST, root-mean-square (RMS), PATH, RANGE, mean velocity (MV), mean frequency (MF), and AREA, are calculated based on Mancini et al. [17], and adjusted for acceleration evaluation. JERK described the jerkiness of the motion, DIST the distance by counting the cumulative value of a, and RMS the root-mean-square
(1)
(2)
(3)
PATH was the total path along with the acceleration and calculated the cumulative value of magnitudes. RANGE was derived from the maximum value of amag
(4)
(5)
Mean velocity estimated the mean velocity based on current acceleration information
(6)
Mean frequency was the mean frequency
(7)

where t is the time duration of the signal, in our case, 20 s.

AREA was the total spanned by acceleration
(8)

Multitaper power spectral density estimate was used to obtain frequency information. 8 taper is used. PWR is the overall power in the frequency domain, F50 is the frequency at which 50% of power is covered below the frequency value, and F95 is where it covers 95%. CFREQ is the centroidal frequency and FREQD is frequency dispersion [18].

For the gyroscope measurement, the maximum and the minimum values within the 20-s signal were found to identify the range of angular motion. A RMS value was also calculated
(9)

where gx is the angular velocity (rad/s) at x-axis.

During the reliability test for mobile app measurement, each condition was measured five times by each device. Factors of each trial were recorded, and the mean, standard deviation, and relative standard deviation in percentage were calculated and listed in Tables 13. To compare the difference between the two devices, we introduced the parameter of relative difference dr
(10)

where f could be any factors mentioned above to evaluate the postural stability, subscripts 1 and 2 indicate the device number.

Table 1

Summary of test results in translational motion setup

Condition/DeviceJERKDISTRMSPATHRANGEMVMFAREAPWRF50F95CFREQFREQD
Eyes openDevice 1Mean0.1950.0080.0092.3220.0810.0172.4320.0050.0010.2504.5601.9220.896
Std0.0110.0000.0000.0860.0440.0060.0640.0000.0000.0000.2630.0460.006
Relative Std5.6%4.8%4.4%3.7%54.2%36.0%2.6%8.9%10.4%0.0%5.8%2.4%0.7%
Device 2Mean0.1500.0040.0061.8300.1330.0023.4940.0020.0000.2608.3053.1380.884
Std0.0300.0000.0010.0420.0820.0010.0400.0000.0000.0200.1290.1030.011
Relative Std19.9%3.2%16.5%2.3%61.2%54.8%1.1%7.2%7.4%7.7%1.6%3.3%1.3%
Devices relative difference26.0%58.4%43.8%23.7%49.1%159.5%35.8%95.4%102.5%3.9%58.2%48.1%1.4%
Eyes closedDevice 1Mean0.1430.0070.0081.9730.0890.0132.3060.0040.0010.2504.3791.8170.890
Std0.0190.0010.0010.0850.0510.0070.0820.0000.0000.0000.2780.0810.006
Relative Std13.4%7.6%10.1%4.3%57.4%56.0%3.6%11.4%14.7%0.0%6.4%4.4%0.6%
Device 2Mean0.0830.0040.0041.4870.0550.0023.1990.0020.0000.2506.9862.6370.891
Std0.0030.0000.0000.0180.0230.0000.0360.0000.0000.0000.2470.0660.002
Relative Std3.7%2.1%3.8%1.2%43.0%24.6%1.1%5.9%2.3%0.0%3.5%2.5%0.3%
Devices relative difference53.4%59.4%57.7%28.1%48.1%158.3%32.4%94.6%105.2%0.0%45.9%36.8%0.1%
One legDevice 1Mean0.5880.0090.0113.6310.1170.0103.1390.0090.0020.2608.1502.9310.890
Std0.0520.0000.0010.1580.0080.0050.0810.0010.0000.0200.2260.1090.004
Relative Std8.8%4.2%4.6%4.3%6.7%45.0%2.6%14.8%9.2%7.7%2.8%3.7%0.4%
Device 2Mean0.5120.0070.0103.1490.1660.0023.5560.0040.0020.3308.3553.2190.867
Std0.0360.0000.0010.0430.0510.0010.0830.0010.0010.0240.6310.2330.010
Relative Std6.9%3.7%11.7%1.4%30.6%32.5%2.3%12.1%23.9%7.4%7.6%7.2%1.2%
Devices relative difference13.7%26.6%8.1%14.2%34.9%124.7%12.5%63.4%13.5%23.7%2.5%9.3%2.6%
Condition/DeviceJERKDISTRMSPATHRANGEMVMFAREAPWRF50F95CFREQFREQD
Eyes openDevice 1Mean0.1950.0080.0092.3220.0810.0172.4320.0050.0010.2504.5601.9220.896
Std0.0110.0000.0000.0860.0440.0060.0640.0000.0000.0000.2630.0460.006
Relative Std5.6%4.8%4.4%3.7%54.2%36.0%2.6%8.9%10.4%0.0%5.8%2.4%0.7%
Device 2Mean0.1500.0040.0061.8300.1330.0023.4940.0020.0000.2608.3053.1380.884
Std0.0300.0000.0010.0420.0820.0010.0400.0000.0000.0200.1290.1030.011
Relative Std19.9%3.2%16.5%2.3%61.2%54.8%1.1%7.2%7.4%7.7%1.6%3.3%1.3%
Devices relative difference26.0%58.4%43.8%23.7%49.1%159.5%35.8%95.4%102.5%3.9%58.2%48.1%1.4%
Eyes closedDevice 1Mean0.1430.0070.0081.9730.0890.0132.3060.0040.0010.2504.3791.8170.890
Std0.0190.0010.0010.0850.0510.0070.0820.0000.0000.0000.2780.0810.006
Relative Std13.4%7.6%10.1%4.3%57.4%56.0%3.6%11.4%14.7%0.0%6.4%4.4%0.6%
Device 2Mean0.0830.0040.0041.4870.0550.0023.1990.0020.0000.2506.9862.6370.891
Std0.0030.0000.0000.0180.0230.0000.0360.0000.0000.0000.2470.0660.002
Relative Std3.7%2.1%3.8%1.2%43.0%24.6%1.1%5.9%2.3%0.0%3.5%2.5%0.3%
Devices relative difference53.4%59.4%57.7%28.1%48.1%158.3%32.4%94.6%105.2%0.0%45.9%36.8%0.1%
One legDevice 1Mean0.5880.0090.0113.6310.1170.0103.1390.0090.0020.2608.1502.9310.890
Std0.0520.0000.0010.1580.0080.0050.0810.0010.0000.0200.2260.1090.004
Relative Std8.8%4.2%4.6%4.3%6.7%45.0%2.6%14.8%9.2%7.7%2.8%3.7%0.4%
Device 2Mean0.5120.0070.0103.1490.1660.0023.5560.0040.0020.3308.3553.2190.867
Std0.0360.0000.0010.0430.0510.0010.0830.0010.0010.0240.6310.2330.010
Relative Std6.9%3.7%11.7%1.4%30.6%32.5%2.3%12.1%23.9%7.4%7.6%7.2%1.2%
Devices relative difference13.7%26.6%8.1%14.2%34.9%124.7%12.5%63.4%13.5%23.7%2.5%9.3%2.6%
Table 2

Summary of test results in rotational motion setup

Condition/DeviceJERKDISTRMSPATHRANGEMVMFAREAPWRF50F95CFREQFREQD
Eyes openDevice 1Mean0.5580.0090.0113.6280.1680.0083.3400.0060.0020.2907.7373.0780.873
Std0.0350.0000.0000.0910.0610.0010.0710.0010.0000.0200.2440.0900.005
Relative Std6.3%3.1%2.1%2.5%36.0%19.9%2.1%10.4%4.3%6.9%3.2%2.9%0.6%
Device 2Mean0.5110.0070.0103.4580.1000.0043.7330.0050.0020.3008.7863.3420.878
Std0.0370.0000.0010.0910.0200.0010.1280.0010.0000.0000.4620.1940.006
Relative Std7.3%5.8%8.0%2.6%19.8%22.0%3.4%14.0%13.2%0.0%5.3%5.8%0.7%
Devices relative difference8.8%15.7%12.0%4.8%51.1%62.8%11.1%31.4%17.3%3.4%12.7%8.2%0.6%
Eyes closedDevice 1Mean0.5010.0080.0123.2290.2100.0133.0920.0070.0020.2807.1862.7510.882
Std0.0810.0010.0020.1820.0720.0070.2020.0020.0000.0240.5600.2160.014
Relative Std16.1%11.1%13.6%5.6%34.1%50.4%6.5%28.6%26.0%8.7%7.8%7.8%1.5%
Device 2Mean0.3190.0060.0082.7610.0970.0023.7630.0030.0010.3008.6063.4200.872
Std0.0210.0000.0000.0590.0460.0000.0670.0000.0000.0000.0800.0610.004
Relative Std6.6%2.7%5.5%2.2%47.0%16.1%1.8%5.7%3.3%0.0%0.9%1.8%0.5%
Devices relative difference44.3%35.5%42.1%15.6%73.7%139.4%19.6%87.5%62.2%6.9%18.0%21.7%1.1%
One legDevice 1Mean2.0620.0140.0206.2820.3420.0073.5060.0130.0070.3008.7883.4650.865
Std0.1590.0010.0010.2010.0820.0030.0200.0010.0010.0000.1750.0890.003
Relative Std7.7%3.6%5.2%3.2%23.8%40.8%0.6%6.6%7.8%0.0%2.0%2.6%0.3%
Device 2Mean2.8880.0180.0277.4270.3490.0103.2680.0150.0130.3407.8493.1450.852
Std0.1930.0020.0030.1420.0810.0040.2080.0030.0030.0200.3080.1960.012
Relative Std6.7%8.6%9.9%1.9%23.1%43.3%6.4%21.3%20.0%5.9%3.9%6.2%1.4%
Devices relative difference33.4%24.0%28.2%16.7%2.0%34.3%7.0%15.6%62.8%12.4%11.3%9.7%1.5%
Condition/DeviceJERKDISTRMSPATHRANGEMVMFAREAPWRF50F95CFREQFREQD
Eyes openDevice 1Mean0.5580.0090.0113.6280.1680.0083.3400.0060.0020.2907.7373.0780.873
Std0.0350.0000.0000.0910.0610.0010.0710.0010.0000.0200.2440.0900.005
Relative Std6.3%3.1%2.1%2.5%36.0%19.9%2.1%10.4%4.3%6.9%3.2%2.9%0.6%
Device 2Mean0.5110.0070.0103.4580.1000.0043.7330.0050.0020.3008.7863.3420.878
Std0.0370.0000.0010.0910.0200.0010.1280.0010.0000.0000.4620.1940.006
Relative Std7.3%5.8%8.0%2.6%19.8%22.0%3.4%14.0%13.2%0.0%5.3%5.8%0.7%
Devices relative difference8.8%15.7%12.0%4.8%51.1%62.8%11.1%31.4%17.3%3.4%12.7%8.2%0.6%
Eyes closedDevice 1Mean0.5010.0080.0123.2290.2100.0133.0920.0070.0020.2807.1862.7510.882
Std0.0810.0010.0020.1820.0720.0070.2020.0020.0000.0240.5600.2160.014
Relative Std16.1%11.1%13.6%5.6%34.1%50.4%6.5%28.6%26.0%8.7%7.8%7.8%1.5%
Device 2Mean0.3190.0060.0082.7610.0970.0023.7630.0030.0010.3008.6063.4200.872
Std0.0210.0000.0000.0590.0460.0000.0670.0000.0000.0000.0800.0610.004
Relative Std6.6%2.7%5.5%2.2%47.0%16.1%1.8%5.7%3.3%0.0%0.9%1.8%0.5%
Devices relative difference44.3%35.5%42.1%15.6%73.7%139.4%19.6%87.5%62.2%6.9%18.0%21.7%1.1%
One legDevice 1Mean2.0620.0140.0206.2820.3420.0073.5060.0130.0070.3008.7883.4650.865
Std0.1590.0010.0010.2010.0820.0030.0200.0010.0010.0000.1750.0890.003
Relative Std7.7%3.6%5.2%3.2%23.8%40.8%0.6%6.6%7.8%0.0%2.0%2.6%0.3%
Device 2Mean2.8880.0180.0277.4270.3490.0103.2680.0150.0130.3407.8493.1450.852
Std0.1930.0020.0030.1420.0810.0040.2080.0030.0030.0200.3080.1960.012
Relative Std6.7%8.6%9.9%1.9%23.1%43.3%6.4%21.3%20.0%5.9%3.9%6.2%1.4%
Devices relative difference33.4%24.0%28.2%16.7%2.0%34.3%7.0%15.6%62.8%12.4%11.3%9.7%1.5%
Table 3

Summary of gyroscope results in rotational motion setup

Condition/Devicegx_maxgx_mingx_rmsgy_maxgy_mingy_rmsgz_maxgz_mingz_rms
Eyes openDevice 1Mean0.036−0.0380.0100.003−0.0050.0010.009−0.0090.003
Std0.0010.0010.0000.0000.0000.0000.0020.0020.000
Relative Std3.0%1.9%0.2%9.1%6.6%1.6%17.4%18.4%7.9%
Device 2Mean0.037−0.0410.0110.002−0.0020.0010.008−0.0080.002
Std0.0010.0010.0000.0010.0000.0000.0020.0020.000
Relative Std2.2%3.1%0.5%33.9%3.5%5.2%22.4%25.2%4.4%
Devices relative difference3.3%8.0%11.4%69.9%64.3%59.7%14.7%12.2%17.2%
Eyes closedDevice 1Mean0.026−0.0350.0080.007−0.0080.0010.010−0.0110.003
Std0.0010.0000.0000.0020.0010.0000.0030.0030.000
Relative Std3.1%1.4%0.2%23.6%14.8%3.7%30.6%31.9%11.1%
Device 2Mean0.029−0.0370.0080.001−0.0010.0010.007−0.0060.002
Std0.0010.0010.0000.0000.0000.0000.0010.0010.000
Relative Std2.2%1.6%0.2%7.6%6.5%2.3%14.3%12.7%6.6%
Devices relative difference9.0%5.6%8.5%143.4%146.0%94.7%40.8%48.5%14.0%
One legDevice 1Mean0.128−0.1100.0180.018−0.0160.0030.014−0.0130.003
Std0.0010.0010.0000.0020.0010.0000.0010.0010.000
Relative Std1.1%1.0%0.2%8.6%3.9%2.7%7.8%8.6%3.4%
Device 2Mean0.146−0.1280.0220.011−0.0150.0020.011−0.0110.003
Std0.0030.0010.0000.0000.0000.0000.0010.0010.000
Relative Std2.4%1.1%0.9%4.1%3.1%0.4%7.9%9.5%4.4%
Devices relative difference12.9%15.6%17.3%43.8%8.9%8.7%19.1%19.4%2.7%
Condition/Devicegx_maxgx_mingx_rmsgy_maxgy_mingy_rmsgz_maxgz_mingz_rms
Eyes openDevice 1Mean0.036−0.0380.0100.003−0.0050.0010.009−0.0090.003
Std0.0010.0010.0000.0000.0000.0000.0020.0020.000
Relative Std3.0%1.9%0.2%9.1%6.6%1.6%17.4%18.4%7.9%
Device 2Mean0.037−0.0410.0110.002−0.0020.0010.008−0.0080.002
Std0.0010.0010.0000.0010.0000.0000.0020.0020.000
Relative Std2.2%3.1%0.5%33.9%3.5%5.2%22.4%25.2%4.4%
Devices relative difference3.3%8.0%11.4%69.9%64.3%59.7%14.7%12.2%17.2%
Eyes closedDevice 1Mean0.026−0.0350.0080.007−0.0080.0010.010−0.0110.003
Std0.0010.0000.0000.0020.0010.0000.0030.0030.000
Relative Std3.1%1.4%0.2%23.6%14.8%3.7%30.6%31.9%11.1%
Device 2Mean0.029−0.0370.0080.001−0.0010.0010.007−0.0060.002
Std0.0010.0010.0000.0000.0000.0000.0010.0010.000
Relative Std2.2%1.6%0.2%7.6%6.5%2.3%14.3%12.7%6.6%
Devices relative difference9.0%5.6%8.5%143.4%146.0%94.7%40.8%48.5%14.0%
One legDevice 1Mean0.128−0.1100.0180.018−0.0160.0030.014−0.0130.003
Std0.0010.0010.0000.0020.0010.0000.0010.0010.000
Relative Std1.1%1.0%0.2%8.6%3.9%2.7%7.8%8.6%3.4%
Device 2Mean0.146−0.1280.0220.011−0.0150.0020.011−0.0110.003
Std0.0030.0010.0000.0000.0000.0000.0010.0010.000
Relative Std2.4%1.1%0.9%4.1%3.1%0.4%7.9%9.5%4.4%
Devices relative difference12.9%15.6%17.3%43.8%8.9%8.7%19.1%19.4%2.7%

3 Results and Discussion

3.1 Translational Motion Setup

3.1.1 Consistency Between Devices - Test With Sweep Signal.

Here, the sweep signal was used to measure the consistency between the two smartphones. The sweep signal was administered continuously between 0.01 Hz and 5 Hz for 2 s and repeated for up to 60 s until the mobile app completed the measurement. As shown in Fig. 3, the time–domain signals overlapped with each other with small differences in the amplitude, indicating acceptable precision between the two devices. There were differences in the sensitivity of the devices depending on the frequencies, with the largest offset being observed at low frequencies ∼3 Hz.

3.1.2 Consistency Between Devices - Test With Postural Stability Input Signal.

The input signal of “one leg” was used to further evaluate the consistency of the device in the application of measuring postural stability. Figure 4 shows the input signal to the shaker and the measurement at two devices for comparison. It can be seen that a peak happened in the middle of the duration, which indicated a perturbation in balance and then balance recovery. The overall tendency of measured signals agreed with the input signal. The differences in measurement can be explained by the stand and the smartphone, both introducing more vibration to the signal measurement compared to the input signal to the shaker.

Figure 5 shows the measurements in two different devices. To test the consistency between the smartphones, results from the two different phones measuring postural stability at the same time were overlaid. As a comparison, the two smartphone measurements matched well in the time domain as shown in Fig. 5, but had slightly different frequency distribution ∼3 Hz, indicating different sensitivities along with frequencies between devices as the embedded sensor model and configuration could be different. This result also agreed with frequency distribution results in Fig. 3, when the devices were tested by sweep signal 0.01 Hz–5 Hz.

Fig. 5
Acceleration measurement in Z-direction at two devices—translational motion setup
Fig. 5
Acceleration measurement in Z-direction at two devices—translational motion setup
Close modal

3.1.3 Consistency Between Different Trials on Devices - Test With Postural Stability Input Signal.

Multiple trials with each device were conducted to test the consistency of the mobile app measurement if an individual uses one particular device for multiple times. For each device and each condition, five trials were conducted. The postural stability factors introduced in the section “Signal Processing Methods” were calculated for each trial and the mean, standard deviation, and relative standard deviation in percentage are listed in Table 1, showing the consistency among different trials when using the same device. As highlighted in light gray, the difference between trials was acceptable, with relative standard deviation values varying from 0% to 61.2%, yet mostly below 10%. Some of the factors have shown better consistency than others. For example, relative standard deviation of FREDQ for trials in two devices and three conditions were 0.7%, 1.3%, 0.6%, 0.3%, 0.4%, 1.2%; JERK were 5.6%, 19.9%, 13.4%, 3.7%, 8.8%, 6.9%, respectively; and RANGE were 54.2%, 61.2%, 57.4%, 43.0%, 6.7%, 30.6%, respectively. This could be explained by the nature of factors. Factors by frequency calculation showed good consistency in general. For the factors from the time-domain information, statistical calculation, such as JERK, cumulatively evaluated the whole 20-s data. Thus, they are more reliable than those that only depend on one or several particular data point(s) along time, such as RANGE.

Differences between devices were evaluated with dr values, highlighted in darker gray in Table 1. The dr values are influenced by different types of postural stability factors and also vary with different devices and test conditions. JERK and PATH have shown acceptable results in time domain signal measurements, and F50 and FREQD for frequency factors. Figures 69 show the comparison between devices of JERK, PATH, F50, and FREQD as examples. The error bars represented the standard deviation. The measurement factors were obtained from acceleration vector a, including the motions in ML and AP planes, and calculated through algorithms introduced in the section “Signal Processing Methods.” These four measures have shown relatively lower dr than others (Table 1), indicating better reliability, thus were recommended by our research results. It can be seen that the two devices show the same tendency on measuring multiple conditions.

Fig. 6
JERK measurement comparison between two devices—translational motion setup
Fig. 6
JERK measurement comparison between two devices—translational motion setup
Close modal
Fig. 7
PATH measurement comparison between two devices—translational motion setup
Fig. 7
PATH measurement comparison between two devices—translational motion setup
Close modal
Fig. 8
F50 measurement comparison between two devices—translational motion setup
Fig. 8
F50 measurement comparison between two devices—translational motion setup
Close modal
Fig. 9
FREQD measurement comparison between two devices—translational motion setup
Fig. 9
FREQD measurement comparison between two devices—translational motion setup
Close modal

3.2 Rotational Motion Setup (Test on Both Translational and Rotational Motion).

The rotational motion setup shown in Fig. 2 provided both translation and rotation motions, that the translational motion captured by accelerometers and the rotational motion by gyroscopes have produced promising results.

3.2.1 Consistency Between Devices—Test With Sweep Signal.

The sweep signal was tested in the rotational motion setup shown in Fig. 2. The time and frequency domain by the accelerometer showed good consistency, as shown in Fig. 10(a). The time domain signals collected by the gyroscopes by two devices were also compared and well-matched, as shown in Fig. 10(b).

Fig. 10
Measurement of sweep signal at two devices—rotational motion setup: (a) acceleration measures at two devices and (b)gyroscope measures at two devices
Fig. 10
Measurement of sweep signal at two devices—rotational motion setup: (a) acceleration measures at two devices and (b)gyroscope measures at two devices
Close modal

3.2.2 Consistency Between Devices—Test With Postural Stability Input Signal.

The postural stability signals were also tested in the same manner. As expected, both accelerometers and gyroscopes measurements were matched well, as shown in Fig. 11, with light offset in accelerometers at low frequency ∼ 3 Hz.

Fig. 11
Measurement of “on one leg” at two devices—rotational motion setup: (a) acceleration measures at two devices and (b)gyroscope measures at two devices
Fig. 11
Measurement of “on one leg” at two devices—rotational motion setup: (a) acceleration measures at two devices and (b)gyroscope measures at two devices
Close modal

3.2.3 Consistency Between Different Trials on Devices—Test With Postural Stability Input Signal.

Figures 1214 show the acceleration results during the rotational motion setup. Detailed results are shown in Table 2. Similar to the test results in translational motion setup, some measures proved to be more precise than others in replicating the input signal across the devices. Some of the factors that performed the best in the translational setup also performed strongly in the rotational motion setup. An example of this being the FREQD, which demonstrated relative differences of 0.6%, 1.1%, and 1.5% across devices for the different conditions of “eyes open,” “eyes closed,” and “one leg,” respectively.

Fig. 12
PATH measurement comparison between two devices—rotational motion setup
Fig. 12
PATH measurement comparison between two devices—rotational motion setup
Close modal
Fig. 13
F50 measurement comparison between two devices—rotational motion setup
Fig. 13
F50 measurement comparison between two devices—rotational motion setup
Close modal
Fig. 14
FREQD measurement comparison between two devices—rotational motion setup
Fig. 14
FREQD measurement comparison between two devices—rotational motion setup
Close modal

Figure 15 shows the range of angular velocity gx, which indicates similarity between the gyroscope measurements by the two devices. The detailed results are shown in Table 3.

Fig. 15
gx_min and gx_max measurements between two devices—rotational motion setup
Fig. 15
gx_min and gx_max measurements between two devices—rotational motion setup
Close modal

4 Conclusion

The measured signals from the different smartphones were able to capture the input signals with precision relative to each other. They were also able to replicate the input signals with accuracy. Differences between the two devices' signal measurements could be accounted for by the difference in embedded sensors and also the vibrations introduced by the use of the stand and smartphones themselves, which could also happen when PwPD uses the mobile app. Sensor sensitivities varied along with frequencies and have shown differences between devices at ∼ 3 Hz, but the variation did not result in major impacts on the results of postural stability factors such as JERK, PATH, F50, and FREQD.

Multiple trials with the same devices have shown good consistency. Based on our test results, some of the postural stability factors were more reliable than others when using mobile devices. Frequency factors were better than time-domain signal calculations in general, and those factors involving the whole signal information performed better than those focusing on one or several particular data points. This information helped us understand the reliability of measuring postural stability using mobile devices, and will help in selecting data processing algorithms when we improve the mobile app.

In the future, we plan to test more variety of smartphones and compare with clinical grade sensors. We also plan to conduct the test with PwPD with variation in age, gender, and different stages of PD, and develop a numerical model to match the sensor measures to the clinical scales, for a better interpretation of the participants' postural stability levels.

Acknowledgment

The research work is supported by the University of Michigan M-Cubed Fund, Undergraduate Research Opportunity Program (UROP), and Summer Undergraduate Research Experience (SURE) at the University of Michigan—Flint.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request.

Nomenclature

a =

acceleration in vector

amag =

magnitude of acceleration, m/s2

ax =

acceleration in x direction, m/s2

ay =

acceleration in y direction, m/s2

az =

acceleration in z direction, m/s2

AREA =

total spanned by acceleration, m2/s4

CFREQ =

centroidal frequency, Hz

dr =

relative difference

DIST =

distance, m/s2

FREQD =

frequency dispersion

F50 =

frequency at 50% of power covered, Hz

F95 =

frequency at 95% of power covered, Hz

gx =

angular velocity at x-axis, rad/s

JERK =

jerkiness, m2/s5

MF =

mean frequency, Hz

MV =

mean velocity, m/s

PATH =

total path along with the acceleration, m/s2

RANGE =

range of acceleration, m/s2

RMS =

root-mean-square, m/s2

RMS_gx =

root-mean-square of angular velocity at x-axis, rad/s

t =

time of signal, seconds

References

1.
Wong
,
S. L.
,
Gilmour
,
H.
, and
Ramage-Morin
,
P. L.
,
2014
, “
Parkinson's Disease: Prevalence, Diagnosis, and Impact
,”
Health Rep.
,
25
(
11
), pp.
10
14
.https://pubmed.ncbi.nlm.nih.gov/25408491/
2.
Marras
,
C.
,
Beck
,
J. C.
,
Bower
,
J. H.
,
Roberts
,
E.
,
Ritz
,
B.
,
Ross
,
G. W.
,
Abbott
,
R. D.
, et al.,
2018
, “
Prevalence of Parkinson's Disease Across North America
,”
Npj Parkinson's Dis.
,
4
(
1
), p.
21
.10.1038/s41531-018-0058-0
3.
Sammer
,
G.
,
Reuter
,
I.
,
Hullmann
,
K.
,
Kaps
,
M.
, and
Vaitl
,
D.
,
2006
, “
Training of Executive Functions in Parkinson's Disease
,”
J. Neurol. Sci.
,
248
(
1–2
), pp.
115
119
.10.1016/j.jns.2006.05.028
4.
Xia
,
R.
, and
Mao
,
Z.-H.
,
2012
, “
Progression of Motor Symptoms in Parkinson's Disease
,”
Neurosci. Bull.
,
28
(
1
), pp.
39
48
.10.1007/s12264-012-1050-z
5.
Park
,
J.-H.
,
Kang
,
Y.-J.
, and
Horak
,
F. B.
,
2015
, “
What is Wrong With Balance in Parkinson's Disease?
,”
J. Mov. Disord.
,
8
(
3
), pp.
109
114
.10.14802/jmd.15018
6.
Bloem
,
B. R.
,
1992
, “
Postural Instability in Parkinson's Disease
,”
Clin. Neurol. Neurosurg.
,
94
, pp.
41
45
.10.1016/0303-8467(92)90018-X
7.
Anpalahan
,
M.
,
Morrison
,
S. G.
, and
Gibson
,
S. J.
,
2013
, “
Hip Fracture Risk Factors and the Discriminability of Hip Fracture Risk Vary by Age: A Case-Control Study
,”
Geriatr. Gerontol. Int.
,
14
(
2
), pp.
413
419
.10.1111/ggi.12117
8.
Vaapio
,
S. S.
,
Salminen
,
M. J.
,
Ojanlatva
,
A.
, and
Kivela
,
S.-L.
,
2008
, “
Quality of Life as an Outcome of Fall Prevention Interventions Among the Aged: A Systematic Review
,”
Eur. J. Public Health
,
19
(
1
), pp.
7
15
.10.1093/eurpub/ckn099
9.
Kim
,
S. D.
,
Allen
,
N. E.
,
Canning
,
C. G.
, and
Fung
,
V. S.
,
2012
, “
Postural Instability in Patients With Parkinson's Disease
,”
CNS Drugs
,
27
(
2
), pp.
97
112
.10.1007/s40263-012-0012-3
10.
Saraiva
,
M.
,
Fernandes
,
O. J.
,
Vilas-Boas
,
J. P.
, and
Castro
,
M. A.
,
2022
, “
Standing Posture in Motor and Cognitive Dual-Tasks During Smartphone Use: Linear and Nonlinear Analysis of Postural Control
,”
Eur. J. Invest. Health, Psychol. Educ.
,
12
(
8
), pp.
1021
1033
.10.3390/ejihpe12080073
11.
Promsri
,
A.
,
Bangkomdet
,
K.
,
Jindatham
,
I.
, and
Jenchang
,
T.
,
2023
, “
Leg Dominance—Surface Stability Interaction: Effects on Postural Control Assessed by Smartphone-Based Accelerometry
,”
Sports
,
11
(
4
), p.
75
.10.3390/sports11040075
12.
Albarello
,
J. C. D. S.
,
Laett
,
C. T.
,
da Silva
,
S. C.
, and
Lemos
,
T.
,
2024
, “
Reliability and Sensitivity of Smartphone Accelerometer to Assess Postural Balance in Healthy Individuals
,”
Muscles, Ligaments Tendons J.
,
14
(
1
), pp.
197
202
.10.32098/mltj.01.2024.17
13.
Borzì
,
L.
,
Olmo
,
G.
,
Artusi
,
C. A.
,
Fabbri
,
M.
,
Rizzone
,
M. G.
,
Romagnolo
,
A.
,
Zibetti
,
M.
, and
Lopiano
,
L.
,
2020
, “
A New Index to Assess Turning Quality and Postural Stability in Patients With Parkinson's Disease
,”
Biomed. Signal Process. Control
,
62
, p.
102059
.10.1016/j.bspc.2020.102059
14.
Borzì
,
L.
,
Fornara
,
S.
,
Amato
,
F.
,
Olmo
,
G.
,
Artusi
,
C. A.
, and
Lopiano
,
L.
,
2020
, “
Smartphone-Based Evaluation of Postural Stability in Parkinson's Disease Patients During Quiet Stance
,”
Electronics
,
9
(
6
), p.
919
.10.3390/electronics9060919
15.
Tang
,
C.
,
Shuva
,
I. K.
,
Thelen
,
M.
,
Zhu
,
L.
, and
Miller
,
N. S.
,
2024
, “
Exploring the Strategies People With Parkinson's Disease Use to Self-Track Symptoms and Medications
,”
ACM Trans. Accessible Comput.
,
17
(
1
), pp.
1
21
.10.1145/3649454
16.
Ickenstein
,
G. W.
,
Ambach
,
H.
,
Klöditz
,
A.
,
Koch
,
H.
,
Isenmann
,
S.
,
Reichmann
,
H.
, and
Ziemssen
,
T.
,
2012
, “
Static Posturography in Aging and Parkinson's Disease
,”
Front. Aging Neurosci.
,
4
, p.
20
.10.3389/fnagi.2012.00020
17.
Mancini
,
M.
,
Salarian
,
A.
,
Carlson-Kuhta
,
P.
,
Zampieri
,
C.
,
King
,
L.
,
Chiari
,
L.
, and
Horak
,
F. B.
,
2012
, “
ISway: A Sensitive, Valid and Reliable Measure of Postural Control
,”
J. Neuroeng. Rehabil.
,
9
(
1
), p.
59
.10.1186/1743-0003-9-59
18.
Prieto
,
T. E.
,
Myklebust
,
J. B.
,
Hoffmann
,
R. G.
,
Lovett
,
E. G.
, and
Myklebust
,
B. M.
,
1996
, “
Measures of Postural Steadiness: Differences Between Healthy Young and Elderly Adults
,”
IEEE Trans. Biomed. Eng.
,
43
(
9
), pp.
956
966
.10.1109/10.532130