 Research
 Open Access
 Published:
Assessing postural stability via the correlation patterns of vertical ground reaction force components
BioMedical Engineering OnLinevolume 15, Article number: 90 (2016)
Abstract
Background
Many methods have been proposed to assess the stability of human postural balance by using a force plate. While most of these approaches characterize postural stability by extracting features from the trajectory of the center of pressure (COP), this work develops stability measures derived from components of the ground reaction force (GRF).
Methods
In comparison with previous GRFbased approaches that extract stability features from the GRF resultant force, this study proposes three feature sets derived from the correlation patterns among the vertical GRF (VGRF) components. The first and second feature sets quantitatively assess the strength and changing speed of the correlation patterns, respectively. The third feature set is used to quantify the stabilizing effect of the GRF coordination patterns on the COP.
Results
In addition to experimentally demonstrating the reliability of the proposed features, the efficacy of the proposed features has also been tested by using them to classify two age groups (18–24 and 65–73 years) in quiet standing. The experimental results show that the proposed features are considerably more sensitive to aging than one of the most effective conventional COP features and two recently proposed COM features.
Conclusions
By extracting information from the correlation patterns of the VGRF components, this study proposes three sets of features to assess human postural stability during quiet standing. As demonstrated by the experimental results, the proposed features are not only robust to intertrial variability but also more accurate than the tested COP and COM features in classifying the older and younger age groups. An additional advantage of the proposed approach is that it reduces the force sensing requirement from 3D to 1D, substantially reducing the cost of the force plate measurement system.
Background
Postural control, the foundation of our ability to stand and walk, is a complex motor skill derived from the interaction of multiple sensorimotor processes. Deterioration of postural control can lead to balance impairment and has been found to be associated with the risk of falling [1–5]. Since fall is a major cause of morbidity and mortality for older people [6, 7], there is a crucial need for a simpletoimplement method to assess postural stability in order to predict fall risk.
Conventionally, postural stability has been assessed by static and/or dynamic posturography [8–11]. Dynamic posturography evaluates postural stability by disrupting a stable stance and measuring the postural response to the external perturbation. In contrast, static posturography measures the postural steadiness of a subject standing quietly (typically on a force plate) without any external disturbance. Numerous experimental studies suggest that subjects who exhibit larger sway during quiet standing have poorer postural stability [5, 12–16].
Mechanically, balance is achieved when the external forces and torques applied to a body are in static equilibrium. During quiet standing, there are two external forces: the weight and the ground reaction force (GRF). The weight acts at the center of mass (COM), which is the point where all of the mass of the human body can be considered to be concentrated. The GRF acts the center of pressure (COP), which is the point of application of the resultant GRF below the body. With these two external forces, any difference between the vertical projection of the COM and the COP results in a moment that destabilizes postural balance [17, 18]. To avoid postural instability, the COP regulates the COM by adjusting the location of the vertical projection of COM to maintain it within a region between the two feet (the polygon of support).
By modelling the standing body be as an inverted pendulum pivoting about ankle, the postural balance system has often been viewed as a feedback control system with the COM and COP as the output and control variables, respectively [17]. The time responses of the COP and COM have been used to provide insights into the control mechanism of human balance. Determining the COM is much more difficult than measuring the COP, which can be found easily by a force plate. As such, many studies have employed COP features to characterize the postural stability during quiet standing [19–24].
Since force plates measure the GRF, a number of studies have also tried to use GRFbased features to assess postural stability. For example, it was found that there was a significant correlation between the standard deviation of the vertical component of the GRF (VGRF) and the Berg Balance Scale (BBS) score [25], which is a 14item validated scale that evaluates balance abilities during sitting, standing and positional changes. A previous study found that young adults tend to have larger GRF variability than older adults during the first 5 s of single leg stance [15]. For poststroke patients, it was shown that the standard deviation of the anterior–posterior GRF component is the most prominent feature that correlates to the BBS score [26]. Statistically significant differences in the VGRF were also found between Parkinson patients and healthy controls [27]. After reviewing thirteen articles that studied the relationship between the history of lowerextremity stress fractures and GRF, the loading rate of the VGRF was concluded to be different for patient and control groups [28]. In summary, the results of previous studies demonstrate the importance of the GRF and its utility in assessing postural stability.
In interacting with the environment and objects, humans typically make multiple adjustments of forces that require agile and subtle coordination. For example, in multifingered grasping of an object the forces applied by the digits are highly coordinated to achieve precise positioning and to guarantee that the needed net force is applied. A challenge for performing a coordination task is to determine the distribution of the applied forces. This has been identified as a motor redundancy problem (also known as Bernstein’s problem) since the resultant grasping force can be achieved with an infinite number of combinations of individual finger forces when the number of constraints of the grasping task is smaller than the number of the employed fingers [29].
In attempting to discover the strategies behind multifinger force production, recent experimental results suggest that forces exerted by all digit pairs tend to be synchronized under a variety of task constraints [30, 31]. Studies have also found that the strength and consistency of finger force covariation may vary with task [32, 33], age [34–36], magnitude of the total force [37] and disease, such as stroke [38], carpal tunnel syndrome [39], olivopontocerebellar atrophy [40] and hemiplegic cerebral palsy [41].
A force plate, one of the most popular tools for quantitative assessment of postural stability, is a rectangular plate with force transducers that measure the applied force [42]. In similarity to finger forces that need to be coordinated with one another to sum to the total needed force, the measurements from the force transducers, which are typically located at the four corner of the force plate, also need to sum to the total GRF. This similarity can be demonstrated by comparing the task of holding a cup and the event of standing on a force plate. Just as the weight of the cup must be shared by the four fingers that are holding it, the weight of the human body must be balanced by the four GRF components measured by the force plate. A fundamental difference between these two weight sharing tasks is that finger forces are directly generated and controlled by human hands whereas the GRF components are indirectly produced and controlled by human body motion. Despite such a difference, since both tasks require the coordination of multiple humangenerated forces, it is valuable to investigate whether the results of extensive studies on the multidigit force coordination problems can help us understand the coordination mechanism of the GRF components during postural balance. Considering the sensitivity of the interfinger force coordination patterns to aging and other operating and physiological conditions, it seems reasonable to assume that there are associations between the GRF component coordination patterns and postural stability.
Compared to previous GRF based postural balance studies that extract stability features from the GRF resultant, this work looks into the components of the GRF. In particular, this work proposes features to characterize the strength, the changing speed and the stabilizing effect of the coordination patterns of the GRF components. To the best of our knowledge, this is the first attempt to use the GRF coordination properties to assess postural balance. The hypothesis is that the proposed features can detect the agingrelated changes in postural balance. The promise of this work is that the proposed features can help uncover ways in which the postural control system is compromised with the aging process and thereby provide valuable information in identifying persons at risk of falling.
Methods
Subjects
The tested subjects consisted of two adult age groups: an older age group (68.7 ± 2.96 years, range 65–73 years, BMI 23.9 ± 4.13 kg/m^{2}) and a younger age group (20.1 ± 1.29 years, range 18–24 years, BMI 22.5 ± 3.21 kg/m^{2}). Each group consisted of ten male and ten female healthy adults, none of whom had pathological conditions that would compromise their postural performance. All experimental procedures were approved by the Institutional Review Board of the Kaohsiung Medical University ChungHo Memorial Hospital, Kaohsiung, Taiwan.
Measurements
For a conventional force plate, shown in Fig. 1, load cells or piezoelectric force transducers are positioned at the four corners. In Fig. 1, the x, y and z axes correspond to the anterior–posterior (AP), mediallateral (ML) and vertical directions, respectively. The origin of the coordinate system is located at the center of the force plate. With F _{ i } signifying the force vector at the ith corner of the force plate, the GRF vector (denoted by f) can be expressed as
where i, j and k are unit vectors along the x, y and z axes and f _{ x }, f _{ y } and f _{ z } represent the GRF components in the AP, ML and vertical directions, respectively. Denoting the vertical component of F _{ i }, as F _{ iz },
Hereafter, we refer to f _{ z } as VGRF and F _{1z }, F _{2z }, F _{3z } and F _{4z } as VGRF components at the respective corners. Next, by defining F _{ ijα } as the sum of F _{ iα } and F _{ jα } (i, j = 1, 2, 3 or 4 and α = x, y or z), the following equations of the COP coordinates can be derived by using Newtonian mechanics:
As shown in Fig. 1, d _{ x } and d _{ y } represent the distances along the xaxis and yaxis, respectively, from the coordinate axes to the force sensors, and d _{ z } represents the distance along the zaxis from the coordinate system to the support surface of the force plate.
In this work, the measurement system consists of a force plate (Kistler 9286AA) connected to a PCbased signal processing system. This force plate system can provide the following signals: two COP coordinates (COP _{ x } and COP _{ y }) and eight GRF components (F _{1z }, F _{2z }, F _{3z }, F _{4z }, F _{12x }, F _{34x }, F _{14y } and F _{23y }), all of which were recorded at a sampling rate of 512 Hz and filtered by a sixthorder Butterworth filter with a cutoff frequency of 5 Hz. The data acquisition and display software was a custom program written in LabVIEW (National Instruments).
Each subject was tested on 2 days with two experimental sessions performed each day. Each session included three 80 s eyesopenclosed trials, during which the participants were barefoot and instructed to look straight ahead at a visual reference and stand quietly (with arms at the side) in a comfortable stance near the center of the force plate for the first 40 s and then maintaining the same posture the subjects closed their eyes in the remaining 40 s of the trial. The data collected from 5 to 35 s of the trials were used for this study. Trials and sessions were separated by approximately one and 5 min of rest, respectively. In summary, 240 (20 subjects × 3 trials × 4 sessions) measures were collected for both age groups. This work focuses on the eyesopen data.
Proposed postural stability features
By analyzing the interaction among GRF components measured by the four force sensors, this section proposes three set of features to characterize postural stability. Instead of studying the interdependency between F _{1}, F _{2}, F _{3} and F _{4}, this work focuses on the coupling between the VGRF components, namely, F _{1z }, F _{2z }, F _{3z } and F _{4z }. This decision is based on the following reasons. First, threedimensional force sensors are much more expensive than onedimensional force sensors. Therefore, developing postural stability features that can be derived directly from the VGRF components can significantly reduce the cost of the measurement system. Second, correlation measures between scalar variables are much simpler to develop than correlation measures between vectors. Third, during quiet standing, the VGRF is much larger than the GRF components in the AP and ML directions. As a result, the horizontal components of the GRF can be considered negligible in computing the COP for quiet standing [43, 44]. This is also the reason why that, even lacking the capability of measuring the horizontal components of the GRF, Nintendo Wii balance board has still been demonstrated to be a valid and reliable tool for assessing the performance of standing balance [45–47].
To quantitatively demonstrate the dominance of the VGRF in quiet standing, by denoting the energy of signal F _{14z }–F _{23z } as E _{ v } and the energy of f _{ x } as E _{ x }, the ratio E _{ v }/(E _{ v } + E _{ x }) can be used to quantify the contribution of the VGRF on the COP _{ x }. For the younger age group the sample mean of this ratio is about 0.991, whereas for the older age group the sample mean of this ratio is about 0.996. Similar results were obtained in the y direction. These data support the assumption that the position of the COP during quiet standing depends largely on the VGRF. Consequently, during the process of quiet standing, the coordinates of the COP can be approximated by
In studying finger force production, finger force correlation has been used to characterize the strength of interdigit coupling [48, 49]. Similarly, to characterize the degree of coupling among the VGRF components, this work uses the Pearson productmoment correlation coefficient to quantify the strength of coupling. Considering the possibility that such coupling may be timevarying, this work divides the time response of these VGRF components into a number of time subintervals. For the kth time subinterval, the following correlation coefficient was computed for every pair of F _{ iz } and F _{ jz } for i = 1, …, 4 and j = i + 1, …, 4.
where V(X, Y) represents the sample covariance between random variables X and Y, S(X) denotes the sample standard deviation of random variable X, and F _{ iz }[k] signifies the sampled data vector associated with the kth time subinterval of F _{ iz }. With a sampling rate of 512 Hz and 0.125 s subinterval length, F _{ iz }[k] is a 64dimensional vector. Note that the selection of 0.125 s was determined by a trialanderror process that tested many different subinterval lengths. The subinterval 0.125 s is chosen since it yields the best classification results for our experimental study.
To quantify the coupling strength between the VGRF components, the first feature set proposed in this work consists of the mean of the absolute value of ρ _{ ij } and can therefore be expressed as
where K is the number of time subintervals. In addition, the average value of these correlation strength features C _{ ij } is used to represent the overall coupling among the four VGRF components. This overall correlation strength feature, M _{ c }, can be written as
To help understand the neurophysiological and biomechanical bases of the proposed correlation strength features, we consider the analogy between the postural balance problem and the finger force coordination problem. As noted in the “Background” section, force sharing among multiple fingers is a typical motor redundancy problem since many different combinations of finger forces can give the same total force. Discovering the coordination strategies that can eliminate or utilize the redundant degrees of freedom (DOFs) has thus become a central issue of motor control [50]. A possible approach for reducing such redundancy is through constraining the biomechanical DOFs. For example, rather than controlling each DOF independently, the DOFs can be neurally coupled. Consequently, the occurrence of significant correlation between finger forces has been hypothesized to reflect control of the fingers more as a single unit. Based on this concept, it has been suggested that finger force correlation should be high when a reduction in the controlled DOFs is permitted by the task conditions. In contrast, finger force correlation should be lower when the task demands larger finger decoupling. Hence, finger force correlation can mirror results for the dynamical DOFs in finger forces [51, 52]. Based on this very idea, the proposed correlation strength features can be used to characterize the degree of interdependence among the VGRF components.
As described in the previous section, in order to maintain postural balance, the human balance control system needs to adjust the COP in order to regulate the COM. The moving direction of the COP therefore changes with time. Consequently, the GRF correlation patterns may also be time varying. For example, from Eq. (5), a possible coordination strategy to move the COP along the positive xdirection is to increase F _{1z } and F _{4z } and decrease F _{2z } and F _{3z }, simultaneously. This results in ρ _{12} < 0, ρ _{13} < 0, ρ _{14} > 0, ρ _{23} > 0, ρ _{24} < 0, ρ _{34} < 0. Similarly, from Eq. (6), the COP can be driven to move in the positive ydirection by increasing F _{1z } and F _{2z } and decreasing F _{3z } and F _{4z } at the same time. This leads to ρ _{12} > 0, ρ _{13} < 0, ρ _{14} < 0, ρ _{23} < 0, ρ _{24} < 0, ρ _{34} > 0. Note that in these two cases, the correlation coefficients ρ _{12}, ρ _{14}, ρ _{23} and ρ _{34} have different signs.
These correlation coefficients may vary with the movement of the COP. Therefore, the second feature set proposed in this work characterizes the changing speed of the correlation patterns. Specifically, by defining the zerocrossing points (ZCPs) as the time instant that the correlation coefficient ρ _{ ij } changes sign, the changing speed feature for ρ _{ ij } is defined as the number of ZCPs of ρ _{ ij }. By denoting the number of ZCPs of ρ _{ ij } as n _{ ij }, the overall correlation pattern changing speed feature is defined as
In the remaining part of the manuscript, n _{ ij } and N _{ C } will be referred to as features of the correlation pattern changing speed.
As addressed in the “Background” section, with COM as the output and COP as the control variables, the postural balance system can be modelled as a feedback control system. One of the key performance measures of a feedback control system is its speed of response which depends on the system bandwidth, the saturation limits of the control variables as well as the maximum changing rate of the control variables, etc. Since different COP movements require different GRF coordination patterns, therefore a rapidly varying COP motion requires a fast GRF coordination pattern changing rate. Consequently, the proposed correlation pattern changing speed features have the potential to assess the response speed of the postural balance system.
During multifinger tasks, previous studies have identified two types of finger interaction modes: enslaving and compensation [53, 54]. In the enslaving mode there is a positive correlation between the individual finger forces; in the compensation mode there is a negative correlation. The distinction between these two interaction modes can be illustrated by considering a simple force coordination problem where a person uses two fingers to push a cart. To stabilize the cart, the time responses of the two finger forces (f _{1} and f _{2}) should be negatively correlated and thus united into a forcestabilizing synergy [55, 56]. In contrast, if the goal is to efficiently change the cart’s acceleration by preventing these two forces from canceling each other’s effect, the time responses of these of two forces should be positively correlated. To generalize this concept to the postural balance problem, by setting G _{ ijz } = −F _{ ijz }, Eqs. (5) and (6) are rewritten as:
As shown by Eq. (11), by interacting in the compensation mode, negatively correlated F _{14z } and G _{23z } have the tendency to stabilize the COP _{ x }. In contrast, when interacting in the enslaving mode, positively correlated F _{14z } and G _{23z } tend to be more capable of producing rapid COP _{ x } movement than negatively correlated F _{14z } and G _{23z }.
To investigate the relative amount of time that the balance control system operates in these two interaction modes, the first member of the third proposed feature set is the ratio of time forces F _{14z } and G _{23z } interacted in the enslaving mode and is signified as SMR _{ x }. Similarly, the second member of the third feature set, SMR _{ y }, represents the ratio of time forces F _{12z } and G _{34z } is coordinated in the enslaving mode. Hereafter, these two features will be referred to as the enslaving mode ratio features.
For a multifinger interaction task, enslaving reflects mechanical and neural connections among fingers, while compensation results from synergic control of fingers to stabilize their net output [55]. Therefore, the proposed features can shed some light on the strategy employed by the postural balance control system in coordinating the VGRF components. This topic will be explored in more detail in the “Discussion” section.
Statistical and classification methods
Validity (accuracy) and reliability (repeatability) are two fundamental properties that need to assured for any measurement method. In this study, the test–retest reliability of the postural stability features was evaluated using the intraclass correlation coefficients (ICC) that compares withinsubject variability with betweensubject variability. In particular, among different versions of ICC, the ICC (2, 1) [57] was employed since it is the most frequently adopted ICC model for COP features [22].
To compare the accuracy of the tested features in detecting the aging effects, each of the tested features was independently used to classify the two age groups. In solving the binary classification problems, we chose the older age group as the positive class when the sample mean of the tested feature was larger in the older group than in the younger group. Otherwise, the younger group was selected as the positive class. In classifying any given sample data point, if its feature value was larger than a specified threshold, then the classifier assigned this data point to the positive class. Otherwise, this data point was classified to the negative class.
If a positive class sample data point was correctly classified as positive, it was counted as a true positive. Otherwise, it was counted as a false negative. Similarly, if a negative class sample data point was classified as negative, it was counted as a true negative. Otherwise, this classification result was counted as a false positive. The true and false positive rates of the overall classification results vary with the value of the threshold. By sweeping the threshold from the most positive to the most negative values, a receiver operating characteristic (ROC) curve can be generated by plotting the true positive rate as a function of the false positive rate. Because the area under the curve (AUC) of the ROC represents an estimate of the probability that the classifier ranks a randomly chosen positive example higher than a negative example, AUC has been considered to be a good performance index for classification methods [58]. In addition to the AUC, this study calculated the best accuracy of the ROC curve as the second efficacy measure for the tested features.
As noted previously, there were twelve measurements from twelve experimental trials (3 measurements/session × 4 sessions) for each subject. Since the variability of a sample mean can be reduced by increasing the sample size, the reliability of the tested features can be improved by replacing the original single measurement result with the average of several measurement results. In addition to improving reliability, this multipletrial average method can improve the performance of the classifiers since the generalization capability of a classifier depends strongly on the training set size. As an example, by using the 3trial averaging method, 220 sample data points can be generated for each subject (the number of possible combinations for selecting 3 objects from 12 objects is 220). Therefore, the 3trial average method increases the training set size from 240 (12 measures/person × 20 persons) to 44,000 (220 measures/person × 20 persons) data points for both classes. Another advantage of the multipletrial average method is that by comparing results obtained from different number of trials the influences of intertrial variability on classification results can be investigated. Unless otherwise indicated, the experimental results presented in this manuscript were obtained using 3trial average method. The influence of the number of trials used by the multipletrial average method will be independently investigated in the final part of the experimental study.
Results
Tables 1 and 2 summarize the mean and standard deviations of the correlation strength features C _{ ij } and correlation pattern changing speed n _{ ij }, respectively. As shown in Table 1, the mean of the correlation strength features of the older group are all larger than those of the younger group. The results of Table 2 show that the mean of n _{ ij } of the younger group are all larger than those of the older group. Next, the proposed features are compared to one of the most popular and effective COP postural stability features.
The conventional COP features can be broadly divided into three groups, namely, displacementrelated, velocityrelated and frequencyrelated measures. This work only compares the proposed features to a COP velocity feature for the following reasons. First, the COP velocity is one of the most commonly used postural stability measures [22]. Second, numerous studies have demonstrated that the COP velocity is one of the most reliable and informative COP measures for postural steadiness [19, 21, 23, 59–62]. Third, COP velocity was found to be significantly correlated with the risk of falling among elderly [1, 4]. Traditionally, the COP velocity can be represented by the mean velocity (MV), mean velocity in the ML direction (MV _{ ML }) and mean velocity in the AP direction (MV _{ AP }). This study chose the MV _{ AP } as the benchmark reference since the MV _{ AP } is considered the most sensitive measure for postural control assessment related regulatory balance activity [13, 60, 63].
Considering the equally important roles of COP and COM in postural balance, this study also tests two recently proposed COM features which can clearly distinguish the difference between age groups [64]. By denoting COM _{ acc } as the ratio of the GRF in the AP direction to the body mass of the participant, the first tested COM feature is the standard deviation of COM _{ acc } divided by the body height. The second tested COM feature is the peak frequency domain power of the COM _{ acc } divided by squared body height. Hereafter, these two COM features will be referred as COM _{ acc }SD and COM _{ acc }PP, respectively.
To compare to these three COP and COM features, this work chose the overall correlation feature M _{ c }, the overall correlation changing speed feature N _{ c } and the enslaving mode ratio feature SMR _{ x } as the representative features for the first, second and third feature sets, respectively. In the remaining part of this manuscript, M _{ c }, N _{ c } and SMR _{ x } will be referred to as the core correlation features.
Table 3 summarizes the mean and standard deviations of MV _{ AP }, two tested COM _{ acc } features, the core correlation features and the enslaving mode ratio SMR _{ y }. In agreement with previous studies [13], the experimental results here indicate that aging increases the COP velocity. As shown in Table 3, the older age group has larger MV _{ AP }, COM _{ acc }SP, COM _{ acc }PP, M _{ c }, SMR _{ x } and SMR _{ y } but smaller N _{ C } than the younger group. The results of Table 4 demonstrate that all three core correlation features have higher accuracy and larger AUC than the three tested COP and COM features.
To investigate how many trials must be averaged to obtain reliable and accurate measurement results, the ICC values and the results of classification accuracy obtained by the multipletrial average method were plotted in Figs. 2 and 3, respectively. These two figures show a general trend that the reliability and classification accuracy improve as the number of trials increases. As shown in Fig. 2, the resulting ICC curves of MV _{ AP }, N _{ c } and SMR _{ x } are quite similar. In comparison, the reliability of M _{ c } is not as good as the other three features. The results of Fig. 3 show that, regardless of the number of trials, all three core correlation features were able to achieve better classification accuracy than MV _{ AP }.
Discussion
The results of Tables 1, 2 and 3 show that aging intensifies the coupling strength among the VGRF components, reduces the changing speed of the correlation patterns and increases the ratio of time of enslaving mode. To analyze the results associated with the correlation strength features, we have considered the analogy between the postural balance problem and the finger force coordination problem. As noted in the “Background” section, under the constraints of the task, higher finger force correlation is associated with the reduction in the controlled DOFs of fingers and lower finger force correlation is related to larger finger decoupling. Hence, finger force correlation can mirror results for the dynamical DOFs in finger forces [51, 52].
For the postural balance problem, the first two task constraints are the COP coordinate Eqs. (5) and (6). The third constraint comes from the force balance in the vertical direction,
where M is the mass of the human body, g is gravity and a _{ z } is the vertical acceleration of the mass center of the human body. With Eqs. (5), (6) and (13) as the task constraints, the coordination of the four VGRF components can therefore be regarded as a redundant DOF problem whose goal is to identify how a person copes with redundancy during quiet standing. As shown by the comparative results of the correlation strength features summarized in Table 1, aging increases the level of coupling among VGRF components. The results imply that aging leads to restricted dynamical DOFs and thus constrains the older age group’s capability in manipulating COP which in turn results in poor postural stability.
The results of Table 2 indicate that the features of correlation pattern changing speed of the younger group is larger than those of the older group. This may be associated with the agingrelated decline of the speed of response of the balance control system. Such a decline may be attributed to the deterioration of the speed and organization of muscle force [65], decline in nervous system speed [66] and progressive decrease of sensory input [67]. Since all these ageinduced factors can degrade the ability to respond quickly to a loss of balance and increase the risk of falling, the possible associations between these agingrelated impairments and the proposed features of correlation pattern changing speed warrant further study.
Regarding the behavior of the enslaving mode ratio features SMR _{ x } and SMR _{ y }, the results of Table 3 show that both age groups spent more time in the enslaving mode than in the compensation mode. The results of Table 3 also show that the older age group spent more time in the enslaving mode than the younger group. Such results can be interpreted by considering the connection between the COP and VGRF components. As addressed in the “Methods” section , any difference between the vertical projection of the COM and the COP will result in a moment that can destabilize the balance of the body. Therefore, the balance control system needs to regulate the COM by moving the COP. To keep the COM within a desired region between the two feet, the COP oscillates on either side of the moving COM at higher frequencies. Since the enslaving mode represents a more efficient choice for moving the COP than the compensation mode, it is reasonable that most of the time (about 75 % for the younger group and 88 % for the older group) forces F _{14z } and G _{23z } interacted in the enslaving mode. However, if the AP component of the COM coincides with COP _{ x } or the difference is too small to be detected by the human sensory system, then F _{14z } and G _{23z } interact in the compensation mode since negative covariation of these two forces produces a stabilizing effect on COP _{ x }. This stabilizing effect leads directly to small MV _{ AP }. A similar argument applies to ML direction. Since it is well known that a smaller COP velocity represents better balance stability, this explains why the younger age group spent more time in the compensation mode than the older group.
The results of Table 3 also show that SMR _{ y } is considerably less than SMR _{ x } (about 52 % for the younger group and 60 % for the older group). This indicates that, compared to the forces F _{14z } and G _{23z } of the AP direction, forces F _{12z } and G _{34z } of the ML direction spent more time in the compensation mode. Such a difference between SMR _{ x } and SMR _{ y } can be explained based on several previous findings. First, it is well known that the posture balance control strategies for the AP and ML directions are fundamentally different. By simultaneously loading one limb and unloading the opposite limb to shift the COP, the control policy of the ML direction has been identified as a hip “load/unload” strategy while the AP control was independent and controlled at the ankle [17]. Second, as the feet were positioned well apart, the stance in the ML direction was considered to be particularly stable [68, 69]. This observation is supported by the experimental results that the COP _{ x } typically has larger range, path length and velocity than those of the COP _{ y } [13, 70, 71]. It was also found that COP features in the AP direction were generally more variable than those in the ML direction [72]. The results on SMR _{ x } and SMR _{ y } are considered in agreement with these findings since by spending more time in the compensation mode, forces F _{12z } and G _{34z } provide a more stabilizing effect on the COP _{ y } than forces F _{14z } and G _{23z } can provide for the COP _{ x }.
Since human balance performance is not perfectly repeatable and thus may vary from trial to trial, the variability can degrade the performance of the postural stability features in assessing balance ability. Therefore, an important property that also needs to be examined is the test–retest reliability of the postural stability features. In this work, we used ICC (2, 1) as the reliability metric to quantify the degree of agreement between measurement results across trials. Theoretically, the value of ICC varies from 0 (no reliability) to 1 (perfect reliability). To determine the number of trials required by the multipletrial average method to achieve a satisfactory ICC, one can incrementally increase the number of trials until the value of ICC plateaus. Since adding extra trials increases the time and effort required for the experimental tests, an efficient feature should be able to satisfy the reliability requirement by using the averages of a small number of trials.
The value of ICC was interpreted using the following criteria: 0.00–0.39 poor, 0.40–0.59 fair, 0.60–0.74 good and 0.75–1.00 excellent [73]. With these criteria, the results of Fig. 2 show that the COP feature MV _{ AP } and the proposed features N _{ c } and SMR _{ x } can achieve excellent reliability with singletrial measurements. The proposed features N _{ c } and SMR _{ x } provide a level of reliability comparable to that of MV _{ AP }, which was considered as one of the most reliable COP features. In addition, for the four features shown in Fig. 2, their ICC values are all larger than 0.9 when the number of trials increases to 3 or higher.
Despite the fact that MV _{ AP }, N _{ c } and SMR _{ x } provide comparable ICC test results, comparing the classification performances of these features is still considered to be indispensable since reliability is only a necessary but not a sufficient condition for validity. The results of Fig. 3 show that, compared to MV _{ AP }, all three core correlation features can achieve higher classification accuracy with a smaller number of trials. In particular, as shown by Table 4, with singletrial measurements, the classification accuracy obtained by MV _{ AP }, M _{ c }, N _{ c } and SMR _{ x } are 0.733, 0.779, 0.801 and 0.864, respectively. In comparison, by increasing the number of trials to six, the multipletrial average method improved the classification accuracy of MV _{ AP } to 0.746. This is still inferior to the classification accuracy obtained by the core correlation features that used only single trial averages.
Figure 3 shows that in applying the multipleaverage method the proposed features N _{ c } and SMR _{ x } plateau when the number of trials increases to two. In contrast, to approach the plateau level, the number of trials required by M _{ c } and MV _{ AP } must be four or larger. Regarding reliability, with 2trial averages, the values of ICC(2, 1) for these two features are close to 0.9. These results indicate that in terms of the capability in distinguishing the aging effects on postural balance N _{ c } and SMR _{ x } are more robust to intertrial variability than M _{ c } and MV _{ AP }. This advantage of robustness is a very attractive property since one can reduce the time and cost of the measurement process by requiring fewer experimental trials.
One possible direction of future work is to perform indepth studies of the proposed features. These include, but are not limited to, investigating the associations between the proposed features and conventional postural stability features, studying the performance of the proposed features in characterizing the postural stability for different stance conditions and examining the connections between the proposed features and standard clinical balance measures such as BBS. Another direction for possible future work is to explore the potential applications of the proposed features. These may include predicting the risk of falling for the aged population, characterizing the severity of balance disorders, progressive monitoring of the outcomes of rehabilitation designed to improve the functional capabilities of balance control and studying the impacts of other pathologies on balance ability.
Conclusion
By extracting information from the correlation patterns of the VGRF components, this study proposes three sets of features to assess human postural stability during quiet standing. The first set of features quantifies the coupling strength among the VGRF components. The second set of features characterizes the changing speed of the correlation patterns. By comparing the amount of time the VGRF components interacted in enslaving (positive correlation pattern) and compensation modes (negative correlation pattern), the third set of features evaluates the COP stabilizing effect produced by the VGRF components. In addition to experimentally demonstrating the reliability of the proposed features, the efficacy of the proposed features has also been tested by using them to classify two age groups (18–24 and 65–73 years) in quiet standing. Experimental results show that aging intensifies the coupling strength among the VGRF components, reduces the changing speed of the correlation patterns and increases the ratio of time of enslaving mode. In addition, the proposed features are not only robust to intertrial variability but also more accurate than one of the most effective COP features and two recently proposed COM features in classifying the older and younger age groups. An additional advantage of the proposed approach is that it reduces the force sensing requirement from 3D to 1D, substantially reducing the cost of the force plate measurement system.
Abbreviations
 AP:

anterior–posterior
 AUC:

area under the curve
 BBS:

Berg balance scale
 COM:

center of mass
 COP:

center of pressure
 DOF:

degree of freedom
 GRF:

ground reaction force
 ICC:

intraclass correlation coefficients
 ML:

medial–lateral
 MV:

mean velocity
 MV_{AP} :

mean velocity in the AP direction
 MV_{ML} :

mean velocity in the ML direction
 ROC:

receiver operating characteristic
 VGRF:

vertical ground reaction force
 ZCP:

zerocrossing point
References
 1.
Fernie GR, Gryfe CI, Holliday PJ, Llewellyn A. The relationship of postural sway in standing to the incidence of falls in geriatric subjects. Age Ageing. 1982;11:11–6.
 2.
Thapa PB, Gideon P, Brockman KG, Fought RL, Ray WA. Clinical and biomechanical measures of balance as fall predictors in ambulatory nursing home residents. J Gerontol A Biol Sci Med Sci. 1996;51:239–46.
 3.
Piirtola M, Era P. Force platform measurements as predictors of falls among older peoplea review. Gerontology. 2006;52:1–16.
 4.
Pajala S, Era P, Koskenvuo M, Kaprio J, Törmäkangas T, Rantanen T. Force platform balance measures as predictors of indoor and outdoor falls in communitydwelling women aged 63–76 years. J Gerontol A Biol Sci Med Sci. 2008;63:171–8.
 5.
Kruz L, Oddsson I, Melzer I. Characteristics of balance control in older persons who fall with injury–a prospective study. J Electromyogr Kinesiol. 2013;23:814–9.
 6.
Centers for Disease Control and Prevention. (CDC). Selfreported falls and fallrelated injuries among persons aged ≥ 65 years—United States. Morb Mortal Wkly Rep. 2008;57:225–9.
 7.
Panel on Prevention of Falls in Older Persons, American Geriatrics Society, British Geriatrics Society. Summary of the Updated American Geriatrics Society/British Geriatrics Society clinical practice guideline for prevention of falls in older persons. J Am Geriatr Soc. 2011;59:148–57.
 8.
Johansson R, Magnusson M. Human postural dynamics. Crit Rev Biomed Eng. 1991;18:413–37.
 9.
Browne JE, O’Hare NJ. Review of the different methods for assessing standing balance. Physiotherapy. 2001;87:489–95.
 10.
Jancová J. Measuring the balance control system–review. Acta Medica (Hradec Kralove). 2008;51:129–37.
 11.
Panjan A, Sarabon N. Review of methods for the evaluation of human body balance. Sport Sci Rev. 2010;19:131–63.
 12.
Prieto TE, Myklebust JB, Myklebust BM. Characterization and modeling of postural steadiness in the elderly: a review. IEEE Trans Rehab Eng. 1993;1:26–34.
 13.
Prieto TE, Myklebust JB, Hoffmann RG, Lovett EG, Myklebust BM. Measures of postural steadiness: differences between healthy young and elderly adults. IEEE Trans Biomed Eng. 1996;43:956–66.
 14.
de Haart M, Geurts AC, Huidekoper SC, Fasotti L, van Limbeek J. Recovery of standing balance in postacute stroke patients: a rehabilitation cohort study. Arch Phys Med Rehabil. 2004;85:886–95.
 15.
Jonsson E, Seiger A, Hirschfeld H. Oneleg stance in healthy young and elderly adults: a measure of postural steadiness. Clin Biomech. 2004;19:688–94.
 16.
Soangra R, Lockhart TE. Comparison of intra individual physiological sway complexity from force plate and inertial measurement unit—biomed 2013. Biomed Sci Instrum. 2013;49:180–6.
 17.
Winter DA. Human balance and posture control during standing and walking. Gait Posture. 1995;3:193–214.
 18.
Maurer C, Peterka RJ. A new interpretation of spontaneous sway measures based on a simple model of human postural control. J. Neurophysiology. 2005;93:189–200.
 19.
Raymakers JA, Samson MM, Verhaar HJ. The assessment of body sway and the choice of the stability parameter(s). Gait Posture. 2005;21:48–58.
 20.
Rougier PR. What insights can be gained when analysing the resultant centre of pressure trajectory? Neurophysiol Clin. 2008;38:363–73.
 21.
Salavati M, Hadian MR, Mazaheri M, Negahban H, Ebrahimi I, Talebian S, Jafari AH, Sanjari MA, Sohani SM, Parnianpour M. Testretest reliability of center of pressure measures of postural stability during quiet standing in a group with musculoskeletal disorders consisting of low back pain, anterior cruciate ligament injury and functional ankle instability. Gait Posture. 2009;29:460–4.
 22.
Ruhe A, Fejer R, Walker B. The testretest reliability of centre of pressure measures in bipedal static task conditions–a systematic review of the literature. Gait Posture. 2010;32:436–45.
 23.
Moghadam M, Ashayeri H, Salavati M, Sarafzadeh J, Taghipoor KD, Saeedi A, Salehi R. Reliability of center of pressure measures of postural stability in healthy older adults: effects of postural task difficulty and cognitive load. Gait Posture. 2011;33:651–5.
 24.
Gasq D, Labrunée M, Amarantini D, Dupui P, Montoya R, Marque P. Betweenday reliability of centre of pressure measures for balance assessment in hemiplegic stroke patients. J Neuroeng Rehabil. 2014;11:39.
 25.
Karlsson A, Frykberg G. Correlations between force plate measures for assessment of balance. Clin Biomech. 2000;15:365–9.
 26.
Frykberg GE, Lindmark B, Lanshammar H, Borg J. Correlation between clinical assessment and force plate measurement of postural control after stroke. J Rehabil Med. 2007;39:448–53.
 27.
Minamisawa T, Sawahata H, Takakura K, Yamaguchi T. Characteristics of temporal fluctuation of the vertical ground reaction force during quiet stance in Parkinson’s disease. Gait Posture. 2012;35:308–11.
 28.
Zadpoor AA, Nikooyan AA. The relationship between lowerextremity stress fractures and the ground reaction force: a systematic review. Clin Biomech. 2011;26:23–8.
 29.
Sporns O, Edelman GM. Solving Bernstein’s problem: a proposal for the development of coordinated movement by selection. Child Dev. 1993;64:960–81.
 30.
Santello M, Soechting JF. Force synergies for multifingered grasping. Exp Brain Res. 2000;133:457–67.
 31.
Rearick MP, Santello M. Force synergies for multifingered grasping: effect of predictability in object center of mass and handedness. Exp Brain Res. 2002;144:38–49.
 32.
Rearick MP, Casares A, Santello M. Taskdependent modulation of multidigit force coordination patterns. J Neurophysiol. 2003;89:1317–26.
 33.
Shim JK, Olafsdottir H, Zatsiorsky VM, Latash ML. The emergence and disappearance of multidigit synergies during forceproduction tasks. Exp Brain Res. 2005;164:260–70.
 34.
Shim JK, Lay BS, Zatsiorsky VM, Latash ML. Agerelated changes in finger coordination in static prehension tasks. J Appl Physiol. 2004;97:213–24.
 35.
Diermayr G, McIsaac TL, Gordon AM. Finger force coordination underlying object manipulation in the elderlya minireview. Gerontology. 2011;57:217–27.
 36.
Park J, Sun Y, Zatsiorsky VM, Latash ML. Agerelated changes in optimality and motor variability: an example of multifinger redundant tasks. Exp Brain Res. 2011;212:1–18.
 37.
Olafsdottir H, Yoshida N, Zatsiorsky VM, Latash ML. Anticipatory covariation of finger forces during selfpaced and reaction time force production. Neurosci Lett. 2005;381:92–6.
 38.
Li S, Latash ML, Yue GH, Siemionow V, Sahgal V. The effects of stroke and age on finger interaction in multifinger force production tasks. Clin Neurophysiol. 2003;114:1646–55.
 39.
Zhang W, Johnston JA, Ross MA, Coakley BJ, Gleason EA, Dueck AC, Santello M. Effects of carpal tunnel syndrome on adaptation of multidigit forces to object mass distribution for wholehand manipulation. J Neuroeng Rehabil. 2012;9:83–95.
 40.
Park J, Lewis MM, Huang X, Latash ML. Effects of olivopontocerebellar atrophy (OPCA) on finger interaction and coordination. Clin Neurophysiol. 2013;124:991–8.
 41.
Bleyenheuft Y, Gordon AM. Precision grip control, sensory impairments and their interactions in children with hemiplegic cerebral palsy: a systematic review. Res Dev Disabil. 2013;34:3014–28.
 42.
Cross R. Standing, walking, running, and jumping on a force plate. Am J Phys. 1999;67:304–9.
 43.
Baratto L, Morasso PG, Re C, Spada G. A new look at posturographic analysis in the clinical context: swaydensity versus other parameterization techniques. Mot Control. 2002;6:246–70.
 44.
Duarte M, Freitas SM. Revision of posturography based on force plate for balance evaluation. Rev Bras Fisioter. 2010;14:183–92.
 45.
Clark RA, Bryant AL, Pua Y, McCrory P, Bennell K, Hunt M. Validity and reliability of the Nintendo Wii Balance Board for assessment of standing balance. Gait Posture. 2010;31:307–10.
 46.
Clark RA, McGough R, Paterson K. Reliability of an inexpensive and portable dynamic weight bearing asymmetry assessment system incorporating dual Nintendo Wii Balance Boards. Gait Posture. 2011;34:288–91.
 47.
Park DS, Lee G. Validity and reliability of balance assessment software using the Nintendo Wii balance board: usability and validation. J Neuroeng Rehabil. 2014;11:99.
 48.
Keogh J, Morrison S, Barrett R. Agerelated differences in interdigit coupling during finger pinching. Eur J Appl Physiol. 2006;97:76–88.
 49.
Li K, Marquardt TL, Li ZM. Removal of visual feedback lowers structural variability of interdigit force coordination during sustained precision pinch. Neurosci Lett. 2013;545:1–5.
 50.
Bernstein N. The coordination and regulation of movements. Oxford: Pergamon; 1967.
 51.
James EG. Dynamical degrees of freedom and correlations in isometric finger force production. Exp Brain Res. 2012;223:533–9.
 52.
Li K, Wei N, Yue S, Thewlis D, Fraysse F, Immink M, Eston R. Coordination of digit force variability during dominant and nondominant sustained precision pinch. Exp Brain Res. 2015;233:2053–60.
 53.
Zatsiorsky VM, Latash ML. Prehension synergies. Exerc Sport Sci Rev. 2004;32:75–80.
 54.
Martin JR, Latash ML, Zatsiorsky VM. Interaction of finger enslaving and error compensation in multiple finger force production. Exp Brain Res. 2009;192:293–8.
 55.
Martin JR, Zatsiorsky VM, Latash ML. Multifinger interaction during involuntary and voluntary single finger force changes. Exp Brain Res. 2011;208:423–35.
 56.
Wu YH, Pazin N, Zatsiorsky VM, Latash ML. Improving finger coordination in young and elderly persons. Exp Brain Res. 2013;226:273–83.
 57.
Shrout PE, Fleiss JL. Intraclass correlation: uses in assessing rater reliability. Psychol Bull. 1979;86:420–8.
 58.
Fawcett T. An introduction to ROC analysis. Pattern Recogn Lett. 2006;27:861–74.
 59.
Lafond D, Corriveau H, Hébert R, Prince F. Intrasession reliability of center of pressure measures of postural steadiness in healthy elderly people. Arch Phys Med Rehabil. 2004;85:896–901.
 60.
CornilleauPérès V, Shabana N, Droulez J, Goh JC, Lee GS, Chew PT. Measurement of the visual contribution to postural steadiness from the COP movement: methodology and reliabilit y. Gait Posture. 2005;22:96–106.
 61.
Lin D, Seol H, Nussbaum MA, Madigan ML. Reliability of COPbased postural sway measures and agerelated differences. Gait Posture. 2008;28:337–42.
 62.
Qiu H, Xiong S. Centerofpressure based postural sway measuresreliability and ability to distinguish between age, fear of falling and fall history. Int J Indus Ergon. 2015;47:37–44.
 63.
Kouzaki M, Masani K. Postural sway during quiet standing is related to physiological tremor and muscle volume in young and elderly adults. Gait Posture. 2012;35:11–7.
 64.
Oba N, Sasagawa S, Yamamoto A, Nakazawa K. Difference in postural control during quiet standing between young children and adults: assessment with center of mass acceleration. PLoS One. 2015;10:e0140235.
 65.
Sturnieks DL, St George R, Lord SR. Balance disorders in the elderly. Neurophysiol Clin. 2008;38:467–78.
 66.
Laughton CA, Slavin M, Katdare K, Nolan L, Bean JF, Kerrigan DC, Phillips E, Lipsitz LA, Collins JJ. Aging, muscle activity, and balance control: physiologic changes associated with balance impairment. Gait Posture. 2003;18:101–8.
 67.
Berger L, Chuzel M, Buisson G, Rougier P. Undisturbed upright stance control in the elderly: part 1. Agerelated changes in undisturbed upright stance control. J Mot Behav. 2005;37:348–58.
 68.
MoeNilssen R, Helbostad JL. Trunk accelerometry as a measure of balance control during quiet standing. Gait Posture. 2002;16:60–8.
 69.
Jeter PE, Wang J, Gu J, Barry MP, Roach C, Corson M, Yang L, Dagnelie G. Intrasession testretest reliability of magnitude and structure of center of pressure from the Nintendo Wii Balance Board™ for a visually impaired and normally sighted population. Gait Posture. 2015;41:482–7.
 70.
Winter DA, Patla AE, Prince F, Ishac M, GieloPerczak K. Stiffness control of balance in quiet standing. J Neurophysiol. 1998;80:1211–21.
 71.
Mayagoitia RE, Lötters JC, Veltink PH, Hermens H. Standing balance evaluation using a triaxial accelerometer. Gait Posture. 2002;16:55–9.
 72.
Baltich J, von Tscharner V, Zandiyeh P, Nigg BM. Quantification and reliability of center of pressure movement during balance tasks of varying difficulty. Gait Posture. 2014;40:327–32.
 73.
Fleiss JL. The design and analysis of clinical experiments. New York: John Wiley & Sons; 1986.
Authors’ contributions
All authors participated in the research. CYH, MLN CWY developed the proposed approach. LYG and RS designed the algorithm framework, CYH and JLS wrote the computer programs and conducted the experiments. MLN and CWY prepared and finished the manuscript. All authors read and approved the final manuscript.
Acknowledgements
We thank Dr. BoLin Su for his assistance in the experimental study. We also thank the Department of Sports Medicine of Kaohsiung Medical University for providing the laboratory space.
Competing interests
The authors declare that they have no competing interests.
Availability of data and supporting materials
The dataset of this manuscript can be obtained by contacting the corresponding author.
Consent of publication
Not applicable.
Ethics approval and consent to participate
All experimental procedures were approved by the Institutional Review Board of the Kaohsiung Medical University ChungHo Memorial Hospital, Kaohsiung, Taiwan [register number KMUHIRB20120807(I)].
Funding
This work was founded by National Sun YatSen University and Kaohsiung Medical University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Author information
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Received
Accepted
Published
DOI
Keywords
 Ground reaction force
 Postural balance
 Quiet standing
 Force plate