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 [10–12], 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.
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.
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.
where t is the time duration of the signal, in our case, 20 s.
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].
where is the angular velocity (rad/s) at x-axis.
where f could be any factors mentioned above to evaluate the postural stability, subscripts 1 and 2 indicate the device number.
Condition/Device | JERK | DIST | RMS | PATH | RANGE | MV | MF | AREA | PWR | F50 | F95 | CFREQ | FREQD | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.195 | 0.008 | 0.009 | 2.322 | 0.081 | 0.017 | 2.432 | 0.005 | 0.001 | 0.250 | 4.560 | 1.922 | 0.896 |
Std | 0.011 | 0.000 | 0.000 | 0.086 | 0.044 | 0.006 | 0.064 | 0.000 | 0.000 | 0.000 | 0.263 | 0.046 | 0.006 | ||
Relative Std | 5.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 2 | Mean | 0.150 | 0.004 | 0.006 | 1.830 | 0.133 | 0.002 | 3.494 | 0.002 | 0.000 | 0.260 | 8.305 | 3.138 | 0.884 | |
Std | 0.030 | 0.000 | 0.001 | 0.042 | 0.082 | 0.001 | 0.040 | 0.000 | 0.000 | 0.020 | 0.129 | 0.103 | 0.011 | ||
Relative Std | 19.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 difference | 26.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 closed | Device 1 | Mean | 0.143 | 0.007 | 0.008 | 1.973 | 0.089 | 0.013 | 2.306 | 0.004 | 0.001 | 0.250 | 4.379 | 1.817 | 0.890 |
Std | 0.019 | 0.001 | 0.001 | 0.085 | 0.051 | 0.007 | 0.082 | 0.000 | 0.000 | 0.000 | 0.278 | 0.081 | 0.006 | ||
Relative Std | 13.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 2 | Mean | 0.083 | 0.004 | 0.004 | 1.487 | 0.055 | 0.002 | 3.199 | 0.002 | 0.000 | 0.250 | 6.986 | 2.637 | 0.891 | |
Std | 0.003 | 0.000 | 0.000 | 0.018 | 0.023 | 0.000 | 0.036 | 0.000 | 0.000 | 0.000 | 0.247 | 0.066 | 0.002 | ||
Relative Std | 3.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 difference | 53.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 leg | Device 1 | Mean | 0.588 | 0.009 | 0.011 | 3.631 | 0.117 | 0.010 | 3.139 | 0.009 | 0.002 | 0.260 | 8.150 | 2.931 | 0.890 |
Std | 0.052 | 0.000 | 0.001 | 0.158 | 0.008 | 0.005 | 0.081 | 0.001 | 0.000 | 0.020 | 0.226 | 0.109 | 0.004 | ||
Relative Std | 8.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 2 | Mean | 0.512 | 0.007 | 0.010 | 3.149 | 0.166 | 0.002 | 3.556 | 0.004 | 0.002 | 0.330 | 8.355 | 3.219 | 0.867 | |
Std | 0.036 | 0.000 | 0.001 | 0.043 | 0.051 | 0.001 | 0.083 | 0.001 | 0.001 | 0.024 | 0.631 | 0.233 | 0.010 | ||
Relative Std | 6.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 difference | 13.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/Device | JERK | DIST | RMS | PATH | RANGE | MV | MF | AREA | PWR | F50 | F95 | CFREQ | FREQD | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.195 | 0.008 | 0.009 | 2.322 | 0.081 | 0.017 | 2.432 | 0.005 | 0.001 | 0.250 | 4.560 | 1.922 | 0.896 |
Std | 0.011 | 0.000 | 0.000 | 0.086 | 0.044 | 0.006 | 0.064 | 0.000 | 0.000 | 0.000 | 0.263 | 0.046 | 0.006 | ||
Relative Std | 5.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 2 | Mean | 0.150 | 0.004 | 0.006 | 1.830 | 0.133 | 0.002 | 3.494 | 0.002 | 0.000 | 0.260 | 8.305 | 3.138 | 0.884 | |
Std | 0.030 | 0.000 | 0.001 | 0.042 | 0.082 | 0.001 | 0.040 | 0.000 | 0.000 | 0.020 | 0.129 | 0.103 | 0.011 | ||
Relative Std | 19.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 difference | 26.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 closed | Device 1 | Mean | 0.143 | 0.007 | 0.008 | 1.973 | 0.089 | 0.013 | 2.306 | 0.004 | 0.001 | 0.250 | 4.379 | 1.817 | 0.890 |
Std | 0.019 | 0.001 | 0.001 | 0.085 | 0.051 | 0.007 | 0.082 | 0.000 | 0.000 | 0.000 | 0.278 | 0.081 | 0.006 | ||
Relative Std | 13.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 2 | Mean | 0.083 | 0.004 | 0.004 | 1.487 | 0.055 | 0.002 | 3.199 | 0.002 | 0.000 | 0.250 | 6.986 | 2.637 | 0.891 | |
Std | 0.003 | 0.000 | 0.000 | 0.018 | 0.023 | 0.000 | 0.036 | 0.000 | 0.000 | 0.000 | 0.247 | 0.066 | 0.002 | ||
Relative Std | 3.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 difference | 53.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 leg | Device 1 | Mean | 0.588 | 0.009 | 0.011 | 3.631 | 0.117 | 0.010 | 3.139 | 0.009 | 0.002 | 0.260 | 8.150 | 2.931 | 0.890 |
Std | 0.052 | 0.000 | 0.001 | 0.158 | 0.008 | 0.005 | 0.081 | 0.001 | 0.000 | 0.020 | 0.226 | 0.109 | 0.004 | ||
Relative Std | 8.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 2 | Mean | 0.512 | 0.007 | 0.010 | 3.149 | 0.166 | 0.002 | 3.556 | 0.004 | 0.002 | 0.330 | 8.355 | 3.219 | 0.867 | |
Std | 0.036 | 0.000 | 0.001 | 0.043 | 0.051 | 0.001 | 0.083 | 0.001 | 0.001 | 0.024 | 0.631 | 0.233 | 0.010 | ||
Relative Std | 6.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 difference | 13.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/Device | JERK | DIST | RMS | PATH | RANGE | MV | MF | AREA | PWR | F50 | F95 | CFREQ | FREQD | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.558 | 0.009 | 0.011 | 3.628 | 0.168 | 0.008 | 3.340 | 0.006 | 0.002 | 0.290 | 7.737 | 3.078 | 0.873 |
Std | 0.035 | 0.000 | 0.000 | 0.091 | 0.061 | 0.001 | 0.071 | 0.001 | 0.000 | 0.020 | 0.244 | 0.090 | 0.005 | ||
Relative Std | 6.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 2 | Mean | 0.511 | 0.007 | 0.010 | 3.458 | 0.100 | 0.004 | 3.733 | 0.005 | 0.002 | 0.300 | 8.786 | 3.342 | 0.878 | |
Std | 0.037 | 0.000 | 0.001 | 0.091 | 0.020 | 0.001 | 0.128 | 0.001 | 0.000 | 0.000 | 0.462 | 0.194 | 0.006 | ||
Relative Std | 7.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 difference | 8.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 closed | Device 1 | Mean | 0.501 | 0.008 | 0.012 | 3.229 | 0.210 | 0.013 | 3.092 | 0.007 | 0.002 | 0.280 | 7.186 | 2.751 | 0.882 |
Std | 0.081 | 0.001 | 0.002 | 0.182 | 0.072 | 0.007 | 0.202 | 0.002 | 0.000 | 0.024 | 0.560 | 0.216 | 0.014 | ||
Relative Std | 16.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 2 | Mean | 0.319 | 0.006 | 0.008 | 2.761 | 0.097 | 0.002 | 3.763 | 0.003 | 0.001 | 0.300 | 8.606 | 3.420 | 0.872 | |
Std | 0.021 | 0.000 | 0.000 | 0.059 | 0.046 | 0.000 | 0.067 | 0.000 | 0.000 | 0.000 | 0.080 | 0.061 | 0.004 | ||
Relative Std | 6.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 difference | 44.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 leg | Device 1 | Mean | 2.062 | 0.014 | 0.020 | 6.282 | 0.342 | 0.007 | 3.506 | 0.013 | 0.007 | 0.300 | 8.788 | 3.465 | 0.865 |
Std | 0.159 | 0.001 | 0.001 | 0.201 | 0.082 | 0.003 | 0.020 | 0.001 | 0.001 | 0.000 | 0.175 | 0.089 | 0.003 | ||
Relative Std | 7.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 2 | Mean | 2.888 | 0.018 | 0.027 | 7.427 | 0.349 | 0.010 | 3.268 | 0.015 | 0.013 | 0.340 | 7.849 | 3.145 | 0.852 | |
Std | 0.193 | 0.002 | 0.003 | 0.142 | 0.081 | 0.004 | 0.208 | 0.003 | 0.003 | 0.020 | 0.308 | 0.196 | 0.012 | ||
Relative Std | 6.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 difference | 33.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/Device | JERK | DIST | RMS | PATH | RANGE | MV | MF | AREA | PWR | F50 | F95 | CFREQ | FREQD | ||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.558 | 0.009 | 0.011 | 3.628 | 0.168 | 0.008 | 3.340 | 0.006 | 0.002 | 0.290 | 7.737 | 3.078 | 0.873 |
Std | 0.035 | 0.000 | 0.000 | 0.091 | 0.061 | 0.001 | 0.071 | 0.001 | 0.000 | 0.020 | 0.244 | 0.090 | 0.005 | ||
Relative Std | 6.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 2 | Mean | 0.511 | 0.007 | 0.010 | 3.458 | 0.100 | 0.004 | 3.733 | 0.005 | 0.002 | 0.300 | 8.786 | 3.342 | 0.878 | |
Std | 0.037 | 0.000 | 0.001 | 0.091 | 0.020 | 0.001 | 0.128 | 0.001 | 0.000 | 0.000 | 0.462 | 0.194 | 0.006 | ||
Relative Std | 7.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 difference | 8.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 closed | Device 1 | Mean | 0.501 | 0.008 | 0.012 | 3.229 | 0.210 | 0.013 | 3.092 | 0.007 | 0.002 | 0.280 | 7.186 | 2.751 | 0.882 |
Std | 0.081 | 0.001 | 0.002 | 0.182 | 0.072 | 0.007 | 0.202 | 0.002 | 0.000 | 0.024 | 0.560 | 0.216 | 0.014 | ||
Relative Std | 16.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 2 | Mean | 0.319 | 0.006 | 0.008 | 2.761 | 0.097 | 0.002 | 3.763 | 0.003 | 0.001 | 0.300 | 8.606 | 3.420 | 0.872 | |
Std | 0.021 | 0.000 | 0.000 | 0.059 | 0.046 | 0.000 | 0.067 | 0.000 | 0.000 | 0.000 | 0.080 | 0.061 | 0.004 | ||
Relative Std | 6.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 difference | 44.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 leg | Device 1 | Mean | 2.062 | 0.014 | 0.020 | 6.282 | 0.342 | 0.007 | 3.506 | 0.013 | 0.007 | 0.300 | 8.788 | 3.465 | 0.865 |
Std | 0.159 | 0.001 | 0.001 | 0.201 | 0.082 | 0.003 | 0.020 | 0.001 | 0.001 | 0.000 | 0.175 | 0.089 | 0.003 | ||
Relative Std | 7.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 2 | Mean | 2.888 | 0.018 | 0.027 | 7.427 | 0.349 | 0.010 | 3.268 | 0.015 | 0.013 | 0.340 | 7.849 | 3.145 | 0.852 | |
Std | 0.193 | 0.002 | 0.003 | 0.142 | 0.081 | 0.004 | 0.208 | 0.003 | 0.003 | 0.020 | 0.308 | 0.196 | 0.012 | ||
Relative Std | 6.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 difference | 33.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/Device | gx_max | gx_min | gx_rms | gy_max | gy_min | gy_rms | gz_max | gz_min | gz_rms | ||
---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.036 | −0.038 | 0.010 | 0.003 | −0.005 | 0.001 | 0.009 | −0.009 | 0.003 |
Std | 0.001 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.002 | 0.002 | 0.000 | ||
Relative Std | 3.0% | 1.9% | 0.2% | 9.1% | 6.6% | 1.6% | 17.4% | 18.4% | 7.9% | ||
Device 2 | Mean | 0.037 | −0.041 | 0.011 | 0.002 | −0.002 | 0.001 | 0.008 | −0.008 | 0.002 | |
Std | 0.001 | 0.001 | 0.000 | 0.001 | 0.000 | 0.000 | 0.002 | 0.002 | 0.000 | ||
Relative Std | 2.2% | 3.1% | 0.5% | 33.9% | 3.5% | 5.2% | 22.4% | 25.2% | 4.4% | ||
Devices relative difference | 3.3% | 8.0% | 11.4% | 69.9% | 64.3% | 59.7% | 14.7% | 12.2% | 17.2% | ||
Eyes closed | Device 1 | Mean | 0.026 | −0.035 | 0.008 | 0.007 | −0.008 | 0.001 | 0.010 | −0.011 | 0.003 |
Std | 0.001 | 0.000 | 0.000 | 0.002 | 0.001 | 0.000 | 0.003 | 0.003 | 0.000 | ||
Relative Std | 3.1% | 1.4% | 0.2% | 23.6% | 14.8% | 3.7% | 30.6% | 31.9% | 11.1% | ||
Device 2 | Mean | 0.029 | −0.037 | 0.008 | 0.001 | −0.001 | 0.001 | 0.007 | −0.006 | 0.002 | |
Std | 0.001 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 2.2% | 1.6% | 0.2% | 7.6% | 6.5% | 2.3% | 14.3% | 12.7% | 6.6% | ||
Devices relative difference | 9.0% | 5.6% | 8.5% | 143.4% | 146.0% | 94.7% | 40.8% | 48.5% | 14.0% | ||
One leg | Device 1 | Mean | 0.128 | −0.110 | 0.018 | 0.018 | −0.016 | 0.003 | 0.014 | −0.013 | 0.003 |
Std | 0.001 | 0.001 | 0.000 | 0.002 | 0.001 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 1.1% | 1.0% | 0.2% | 8.6% | 3.9% | 2.7% | 7.8% | 8.6% | 3.4% | ||
Device 2 | Mean | 0.146 | −0.128 | 0.022 | 0.011 | −0.015 | 0.002 | 0.011 | −0.011 | 0.003 | |
Std | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 2.4% | 1.1% | 0.9% | 4.1% | 3.1% | 0.4% | 7.9% | 9.5% | 4.4% | ||
Devices relative difference | 12.9% | 15.6% | 17.3% | 43.8% | 8.9% | 8.7% | 19.1% | 19.4% | 2.7% |
Condition/Device | gx_max | gx_min | gx_rms | gy_max | gy_min | gy_rms | gz_max | gz_min | gz_rms | ||
---|---|---|---|---|---|---|---|---|---|---|---|
Eyes open | Device 1 | Mean | 0.036 | −0.038 | 0.010 | 0.003 | −0.005 | 0.001 | 0.009 | −0.009 | 0.003 |
Std | 0.001 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.002 | 0.002 | 0.000 | ||
Relative Std | 3.0% | 1.9% | 0.2% | 9.1% | 6.6% | 1.6% | 17.4% | 18.4% | 7.9% | ||
Device 2 | Mean | 0.037 | −0.041 | 0.011 | 0.002 | −0.002 | 0.001 | 0.008 | −0.008 | 0.002 | |
Std | 0.001 | 0.001 | 0.000 | 0.001 | 0.000 | 0.000 | 0.002 | 0.002 | 0.000 | ||
Relative Std | 2.2% | 3.1% | 0.5% | 33.9% | 3.5% | 5.2% | 22.4% | 25.2% | 4.4% | ||
Devices relative difference | 3.3% | 8.0% | 11.4% | 69.9% | 64.3% | 59.7% | 14.7% | 12.2% | 17.2% | ||
Eyes closed | Device 1 | Mean | 0.026 | −0.035 | 0.008 | 0.007 | −0.008 | 0.001 | 0.010 | −0.011 | 0.003 |
Std | 0.001 | 0.000 | 0.000 | 0.002 | 0.001 | 0.000 | 0.003 | 0.003 | 0.000 | ||
Relative Std | 3.1% | 1.4% | 0.2% | 23.6% | 14.8% | 3.7% | 30.6% | 31.9% | 11.1% | ||
Device 2 | Mean | 0.029 | −0.037 | 0.008 | 0.001 | −0.001 | 0.001 | 0.007 | −0.006 | 0.002 | |
Std | 0.001 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 2.2% | 1.6% | 0.2% | 7.6% | 6.5% | 2.3% | 14.3% | 12.7% | 6.6% | ||
Devices relative difference | 9.0% | 5.6% | 8.5% | 143.4% | 146.0% | 94.7% | 40.8% | 48.5% | 14.0% | ||
One leg | Device 1 | Mean | 0.128 | −0.110 | 0.018 | 0.018 | −0.016 | 0.003 | 0.014 | −0.013 | 0.003 |
Std | 0.001 | 0.001 | 0.000 | 0.002 | 0.001 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 1.1% | 1.0% | 0.2% | 8.6% | 3.9% | 2.7% | 7.8% | 8.6% | 3.4% | ||
Device 2 | Mean | 0.146 | −0.128 | 0.022 | 0.011 | −0.015 | 0.002 | 0.011 | −0.011 | 0.003 | |
Std | 0.003 | 0.001 | 0.000 | 0.000 | 0.000 | 0.000 | 0.001 | 0.001 | 0.000 | ||
Relative Std | 2.4% | 1.1% | 0.9% | 4.1% | 3.1% | 0.4% | 7.9% | 9.5% | 4.4% | ||
Devices relative difference | 12.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.
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 6–9 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.
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).
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.
3.2.3 Consistency Between Different Trials on Devices—Test With Postural Stability Input Signal.
Figures 12–14 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.
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
- =
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