1
|
Debnath M, Chang J, Bhandari K, Nagy DJ, Insperger T, Milton JG, Ngu AHH. Pole balancing on the fingertip: model-motivated machine learning forecasting of falls. Front Physiol 2024; 15:1334396. [PMID: 38638278 PMCID: PMC11024436 DOI: 10.3389/fphys.2024.1334396] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 03/12/2024] [Indexed: 04/20/2024] Open
Abstract
Introduction: There is increasing interest in developing mathematical and computational models to forecast adverse events in physiological systems. Examples include falls, the onset of fatal cardiac arrhythmias, and adverse surgical outcomes. However, the dynamics of physiological systems are known to be exceedingly complex and perhaps even chaotic. Since no model can be perfect, it becomes important to understand how forecasting can be improved, especially when training data is limited. An adverse event that can be readily studied in the laboratory is the occurrence of stick falls when humans attempt to balance a stick on their fingertips. Over the last 20 years, this task has been extensively investigated experimentally, and presently detailed mathematical models are available. Methods: Here we use a long short-term memory (LTSM) deep learning network to forecast stick falls. We train this model to forecast stick falls in three ways: 1) using only data generated by the mathematical model (synthetic data), 2) using only stick balancing recordings of stick falls measured using high-speed motion capture measurements (human data), and 3) using transfer learning which combines a model trained using synthetic data plus a small amount of human balancing data. Results: We observe that the LTSM model is much more successful in forecasting a fall using synthetic data than it is in forecasting falls for models trained with limited available human data. However, with transfer learning, i.e., the LTSM model pre-trained with synthetic data and re-trained with a small amount of real human balancing data, the ability to forecast impending falls in human data is vastly improved. Indeed, it becomes possible to correctly forecast 60%-70% of real human stick falls up to 2.35 s in advance. Conclusion: These observations support the use of model-generated data and transfer learning techniques to improve the ability of computational models to forecast adverse physiological events.
Collapse
Affiliation(s)
- Minakshi Debnath
- Department of Computer Science, Texas State University, San Marcos, TX, United States
| | - Joshua Chang
- Department of Neurology, Dell Medical School, The University of Texas at Austin, Austin, TX, United States
- Oden Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, TX, United States
| | - Keshav Bhandari
- Department of Computer Science, Texas State University, San Marcos, TX, United States
| | - Dalma J. Nagy
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
- HUN-REN–BME Dynamics of Machines Research Group, Budapest, Hungary
| | - John G. Milton
- Department of Neurology, Dell Medical School, The University of Texas at Austin, Austin, TX, United States
| | - Anne H. H. Ngu
- Department of Computer Science, Texas State University, San Marcos, TX, United States
| |
Collapse
|
2
|
Nagy DJ, Milton JG, Insperger T. Controlling stick balancing on a linear track: Delayed state feedback or delay-compensating predictor feedback? BIOLOGICAL CYBERNETICS 2023; 117:113-127. [PMID: 36943486 PMCID: PMC10160210 DOI: 10.1007/s00422-023-00957-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/26/2022] [Accepted: 02/18/2023] [Indexed: 05/06/2023]
Abstract
A planar stick balancing task was investigated using stabilometry parameters (SP); a concept initially developed to assess the stability of human postural sway. Two subject groups were investigated: 6 subjects (MD) with many days of balancing a 90 cm stick on a linear track and 25 subjects (OD) with only one day of balancing experience. The underlying mechanical model is a pendulum-cart system. Two control force models were investigated by means of numerical simulations: (1) delayed state feedback (DSF); and (2) delay-compensating predictor feedback (PF). Both models require an internal model and are subject to certainty thresholds with delayed switching. Measured and simulated time histories were compared quantitatively using a cost function in terms of some essential SPs for all subjects. Minimization of the cost function showed that the control strategy of both OD and MD subjects can better be described by DSF. The control mechanism for the MD subjects was superior in two aspects: (1) they devoted less energy to controlling the cart's position; and (2) their perception threshold for the stick's angular velocity was found to be smaller. Findings support the concept that when sufficient sensory information is readily available, a delay-compensating PF strategy is not necessary.
Collapse
Affiliation(s)
- Dalma J Nagy
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
| | - John G Milton
- W. M. Keck Science Center, Claremont Colleges, Claremont, CA, 91711, USA
| | - Tamas Insperger
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary.
- ELKH-BME Dynamics of Machines Research Group, Budapest, Hungary.
| |
Collapse
|
3
|
Critical delay as a measure for the difficulty of frontal plane balancing on rolling balance board. J Biomech 2022; 138:111117. [DOI: 10.1016/j.jbiomech.2022.111117] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Revised: 03/29/2022] [Accepted: 04/28/2022] [Indexed: 11/19/2022]
|
4
|
Nagy DJ, Insperger T. Predictor feedback models for stick balancing with delay mismatch and sensory dead zones. CHAOS (WOODBURY, N.Y.) 2022; 32:053108. [PMID: 35649988 DOI: 10.1063/5.0087019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 04/18/2022] [Indexed: 06/15/2023]
Abstract
Human stick balancing is investigated in terms of reaction time delay and sensory dead zones for position and velocity perception using a special combination of delayed state feedback and mismatched predictor feedback as a control model. The corresponding mathematical model is a delay-differential equation with event-driven switching in the control action. Due to the sensory dead zones, initial conditions of the actual state cannot always be provided for an internal-model-based prediction, which indicates that (1) perfect prediction is not possible and (2) the delay in the switching condition cannot be compensated. The imperfection of the predictor is described by the delay mismatch, which is treated as a lumped parameter that creates a transition between perfect predictor feedback (zero delay mismatch) and delayed state feedback (mismatch equal to switching delay). The maximum admissible switching delay (critical delay) is determined numerically based on a practical stabilizability concept. This critical delay is compared to a realistic reference value of 230 ms in order to assess the possible regions of the threshold values for position and velocity perception. The ratio of the angular position and angular velocity for 44 successful balancing trials by 8 human subjects was used to validate the numerical results. Comparison of actual human stick balancing data and numerical simulations based on the mismatched predictor feedback model provided a plausible range of parameters: position detection threshold 1°, velocity detection threshold between 4.24 and 9.35°/s, and delay mismatch around 100-150 ms.
Collapse
Affiliation(s)
- Dalma J Nagy
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Tamás Insperger
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
5
|
Kovacs BA, Insperger T. Virtual stick balancing: skill development in Newtonian and Aristotelian dynamics. J R Soc Interface 2022; 19:20210854. [PMID: 35232278 PMCID: PMC8889188 DOI: 10.1098/rsif.2021.0854] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Human reaction delay significantly limits manual control of unstable systems. It is more difficult to balance a short stick on a fingertip than a long one, because a shorter stick falls faster and therefore requires faster reactions. In this study, a virtual stick balancing environment was developed where the reaction delay can be artificially modulated and the law of motion can be changed between second-order (Newtonian) and first-order (Aristotelian) dynamics. Twenty-four subjects were separated into two groups and asked to perform virtual stick balancing programmed according to either Newtonian or Aristotelian dynamics. The shortest stick length (critical length, Lc) was determined for different added delays in six sessions of balancing trials performed on different days. The observed relation between Lc and the overall reaction delay τ reflected the feature of the underlying mathematical models: (i) for the Newtonian dynamics Lc is proportional to τ2; (ii) for the Aristotelian dynamics Lc is proportional to τ. Deviation of the measured Lc(τ) function from the theoretical one was larger for the Newtonian dynamics for all sessions, which suggests that, at least in virtually controlled tasks, it is more difficult to adopt second-order dynamics than first-order dynamics.
Collapse
Affiliation(s)
- Balazs A. Kovacs
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| |
Collapse
|
6
|
Zelei A, Milton J, Stepan G, Insperger T. Response to perturbation during quiet standing resembles delayed state feedback optimized for performance and robustness. Sci Rep 2021; 11:11392. [PMID: 34059718 PMCID: PMC8167093 DOI: 10.1038/s41598-021-90305-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2021] [Accepted: 05/05/2021] [Indexed: 12/03/2022] Open
Abstract
Postural sway is a result of a complex action–reaction feedback mechanism generated by the interplay between the environment, the sensory perception, the neural system and the musculation. Postural oscillations are complex, possibly even chaotic. Therefore fitting deterministic models on measured time signals is ambiguous. Here we analyse the response to large enough perturbations during quiet standing such that the resulting responses can clearly be distinguished from the local postural sway. Measurements show that typical responses very closely resemble those of a critically damped oscillator. The recovery dynamics are modelled by an inverted pendulum subject to delayed state feedback and is described in the space of the control parameters. We hypothesize that the control gains are tuned such that (H1) the response is at the border of oscillatory and nonoscillatory motion similarly to the critically damped oscillator; (H2) the response is the fastest possible; (H3) the response is a result of a combined optimization of fast response and robustness to sensory perturbations. Parameter fitting shows that H1 and H3 are accepted while H2 is rejected. Thus, the responses of human postural balance to “large” perturbations matches a delayed feedback mechanism that is optimized for a combination of performance and robustness.
Collapse
Affiliation(s)
- Ambrus Zelei
- MTA-BME Research Group on Dynamics of Machines and Vehicles, Budapest, 1111, Hungary.,MTA-BME Lendület Human Balancing Research Group, Budapest, 1111, Hungary
| | - John Milton
- The Claremont Colleges, W. M. Keck Science Center, Claremont, CA, 91711, USA
| | - Gabor Stepan
- MTA-BME Research Group on Dynamics of Machines and Vehicles, Budapest, 1111, Hungary.,Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, 1111, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, 1111, Hungary. .,MTA-BME Lendület Human Balancing Research Group, Budapest, 1111, Hungary.
| |
Collapse
|
7
|
Molnar CA, Zelei A, Insperger T. Rolling balance board of adjustable geometry as a tool to assess balancing skill and to estimate reaction time delay. J R Soc Interface 2021; 18:20200956. [PMID: 33784884 DOI: 10.1098/rsif.2020.0956] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
The relation between balancing performance and reaction time is investigated for human subjects balancing on rolling balance board of adjustable physical parameters: adjustable rolling radius R and adjustable board elevation h. A well-defined measure of balancing performance is whether a subject can or cannot balance on balance board with a given geometry (R, h). The balancing ability is linked to the stabilizability of the underlying two-degree-of-freedom mechanical model subject to a delayed proportional-derivative feedback control. Although different sensory perceptions involve different reaction times at different hierarchical feedback loops, their effect is modelled as a single lumped reaction time delay. Stabilizability is investigated in terms of the time delay in the mechanical model: if the delay is larger than a critical value (critical delay), then no stabilizing feedback control exists. Series of balancing trials by 15 human subjects show that it is more difficult to balance on balance board configuration associated with smaller critical delay, than on balance boards associated with larger critical delay. Experiments verify the feature of the mechanical model that a change in the rolling radius R results in larger change in the difficulty of the task than the same change in the board elevation h does. The rolling balance board characterized by the two well-defined parameters R and h can therefore be a useful device to assess human balancing skill and to estimate the corresponding lumped reaction time delay.
Collapse
Affiliation(s)
- Csenge A Molnar
- Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, Hungary.,MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| | - Ambrus Zelei
- MTA-BME Research Group on Dynamics of Machines and Vehicles, Budapest, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest, Hungary.,MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| |
Collapse
|
8
|
Gyebrószki G, Csernák G, Milton JG, Insperger T. The effects of sensory quantization and control torque saturation on human balance control. CHAOS (WOODBURY, N.Y.) 2021; 31:033145. [PMID: 33810721 DOI: 10.1063/5.0028197] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2020] [Accepted: 02/25/2021] [Indexed: 06/12/2023]
Abstract
The effect of reaction delay, temporal sampling, sensory quantization, and control torque saturation is investigated numerically for a single-degree-of-freedom model of postural sway with respect to stability, stabilizability, and control effort. It is known that reaction delay has a destabilizing effect on the balancing process: the later one reacts to a perturbation, the larger the possibility of falling. If the delay is larger than a critical value, then stabilization is not even possible. In contrast, numerical analysis showed that quantization and control torque saturation have a stabilizing effect: the region of stabilizing control gains is greater than that of the linear model. Control torque saturation allows the application of larger control gains without overcontrol while sensory quantization plays a role of a kind of filter when sensory noise is present. These beneficial effects are reflected in the energy demand of the control process. On the other hand, neither control torque saturation nor sensory quantization improves stabilizability properties. In particular, the critical delay cannot be increased by adding saturation and/or sensory quantization.
Collapse
Affiliation(s)
- Gergely Gyebrószki
- Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest 1111, Hungary
| | - Gábor Csernák
- Department of Applied Mechanics, Budapest University of Technology and Economics, Budapest 1111, Hungary
| | - John G Milton
- The Claremont Colleges, W. M. Keck Science Center, Claremont, California 91711, USA
| | - Tamás Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, Budapest 1111, Hungary
| |
Collapse
|
9
|
Suzuki Y, Nakamura A, Milosevic M, Nomura K, Tanahashi T, Endo T, Sakoda S, Morasso P, Nomura T. Postural instability via a loss of intermittent control in elderly and patients with Parkinson's disease: A model-based and data-driven approach. CHAOS (WOODBURY, N.Y.) 2020; 30:113140. [PMID: 33261318 DOI: 10.1063/5.0022319] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/21/2020] [Accepted: 10/28/2020] [Indexed: 06/12/2023]
Abstract
Postural instability is one of the major symptoms of Parkinson's disease. Here, we assimilated a model of intermittent delay feedback control during quiet standing into postural sway data from healthy young and elderly individuals as well as patients with Parkinson's disease to elucidate the possible mechanisms of instability. Specifically, we estimated the joint probability distribution of a set of parameters in the model using the Bayesian parameter inference such that the model with the inferred parameters can best-fit sway data for each individual. It was expected that the parameter values for three populations would distribute differently in the parameter space depending on their balance capability. Because the intermittent control model is parameterized by a parameter associated with the degree of intermittency in the control, it can represent not only the intermittent model but also the traditional continuous control model with no intermittency. We showed that the inferred parameter values for the three groups of individuals are classified into two major groups in the parameter space: one represents the intermittent control mostly for healthy people and patients with mild postural symptoms and the other the continuous control mostly for some elderly and patients with severe postural symptoms. The results of this study may be interpreted by postulating that increased postural instability in most Parkinson's patients and some elderly persons might be characterized as a dynamical disease.
Collapse
Affiliation(s)
- Yasuyuki Suzuki
- Graduate School of Engineering Science, Osaka University, Osaka 5608531, Japan
| | - Akihiro Nakamura
- Graduate School of Engineering Science, Osaka University, Osaka 5608531, Japan
| | - Matija Milosevic
- Graduate School of Engineering Science, Osaka University, Osaka 5608531, Japan
| | - Kunihiko Nomura
- Department of Information Technology and Social Sciences, Osaka University of Economics, Osaka 5338533, Japan
| | - Takao Tanahashi
- Department of Neurology, Osaka Rosai Hospital, Osaka 5918025, Japan
| | - Takuyuki Endo
- Department of Neurology, Osaka Toneyama Medical Center, Osaka 5608552, Japan
| | - Saburo Sakoda
- Department of Neurology, Osaka Toneyama Medical Center, Osaka 5608552, Japan
| | - Pietro Morasso
- Center for Human Technologies, Istituto Italiano di Tecnologia, Genoa 16163, Italy
| | - Taishin Nomura
- Graduate School of Engineering Science, Osaka University, Osaka 5608531, Japan
| |
Collapse
|
10
|
State-space intermittent feedback stabilization of a dual balancing task. Sci Rep 2020; 10:8470. [PMID: 32439947 PMCID: PMC7242428 DOI: 10.1038/s41598-020-64911-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2019] [Accepted: 03/13/2020] [Indexed: 12/25/2022] Open
Abstract
Balancing the body in upright standing and balancing a stick on the fingertip are two examples of unstable tasks that, in spite of strong motor and sensory differences, appear to share a similar motor control paradigm, namely a state-space intermittent feedback stabilization mechanism. In this study subjects were required to perform the two tasks simultaneously, with the purpose of highlighting both the coordination between the two skills and the underlying interaction between the corresponding controllers. The experimental results reveal, in particular, that upright standing (the less critical task) is modified in an adaptive way, in order to facilitate the more critical task (stick balancing), but keeping the overall spatio-temporal signature well known in regular upright standing. We were then faced with the following question: to which extent the physical/biomechanical interaction between the two independent intermittent controllers is capable to explain the dual task coordination patterns, without the need to introduce an additional, supervisory layer/module? By comparing the experimental data with the output of a simulation study we support the former hypothesis, suggesting that it is made possible by the intrinsic robustness of both state-space intermittent feedback stabilization mechanisms.
Collapse
|
11
|
Buza G, Milton J, Bencsik L, Insperger T. Establishing metrics and control laws for the learning process: ball and beam balancing. BIOLOGICAL CYBERNETICS 2020; 114:83-93. [PMID: 31955261 PMCID: PMC7062859 DOI: 10.1007/s00422-020-00815-z] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/27/2019] [Accepted: 01/04/2020] [Indexed: 06/02/2023]
Abstract
Understanding how dexterity improves with practice is a fundamental challenge of motor control and neurorehabilitation. Here we investigate a ball and beam implementation of a dexterity puzzle in which subjects stabilize a ball at the mid-point of a beam by manipulating the angular position of the beam. Stabilizability analysis of different biomechanical models for the ball and beam task with time-delayed proportional-derivative feedback identified the angular position of the beam as the manipulated variable. Consequently, we monitored the changes in the dynamics with learning by measuring changes in the control parameters. Two types of stable motion are possible: node type (nonoscillatory) and spiral type (oscillatory). Both types of motion are observed experimentally and correspond to well-defined regions in the parameter space of the control gains. With practice the control gains for each subject move close to or on the portion of the boundary which separates the node-type and spiral-type solutions and which is associated with the rightmost characteristic exponent of smallest real part. These observations suggest that with learning the control gains for ball and beam balancing change in such a way that minimizes overshoot and the settling time. This study provides an example of how mathematical analysis together with careful experimental observations can shed light onto the early stages of skill acquisition. Since the difficulty of this task depends on the length of the beam, ball and beam balancing tasks may be useful for the rehabilitation of children with dyspraxia and those recovering from a stroke.
Collapse
Affiliation(s)
- Gergely Buza
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
- MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| | - John Milton
- W. M. Keck Science Center, Claremont Colleges, Claremont, CA 91711 USA
| | - Laszlo Bencsik
- MTA-BME Research Group on Dynamics of Machines and Vehicles, Budapest, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Faculty of Mechanical Engineering, Budapest University of Technology and Economics, Budapest, Hungary
- MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| |
Collapse
|
12
|
Fu C, Suzuki Y, Morasso P, Nomura T. Phase resetting and intermittent control at the edge of stability in a simple biped model generates 1/f-like gait cycle variability. BIOLOGICAL CYBERNETICS 2020; 114:95-111. [PMID: 31960137 DOI: 10.1007/s00422-020-00816-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/09/2019] [Accepted: 01/04/2020] [Indexed: 06/10/2023]
Abstract
The 1/f-like gait cycle variability, characterized by temporal changes in stride-time intervals during steady-state human walking, is a well-documented gait characteristic. Such gait fractality is apparent in healthy young adults, but tends to disappear in the elderly and patients with neurological diseases. However, mechanisms that give rise to gait fractality have yet to be fully clarified. We aimed to provide novel insights into neuro-mechanical mechanisms of gait fractality, based on a numerical simulation model of biped walking. A previously developed heel-toe footed, seven-rigid-link biped model with human-like body parameters in the sagittal plane was implemented and expanded. It has been shown that the gait model, stabilized rigidly by means of impedance control with large values of proportional (P) and derivative (D) gains for a linear feedback controller, is destabilized only in a low-dimensional eigenspace, as P and D decrease below and even far below critical values. Such low-dimensional linear instability can be compensated by impulsive, phase-dependent actions of nonlinear controllers (phase resetting and intermittent controllers), leading to the flexible walking with joint impedance in the model being as small as that in humans. Here, we added white noise to the model to examine P-value-dependent stochastic dynamics of the model for small D-values. The simulation results demonstrated that introduction of the nonlinear controllers in the model determined the fractal features of gait for a wide range of the P-values, provided that the model operates near the edge of stability. In other words, neither the model stabilized only by pure impedance control even at the edge of linear stability, nor the model stabilized by specific nonlinear controllers, but with P-values far inside the stability region, could induce gait fractality. Although only limited types of controllers were examined, we suggest that the impulsive nonlinear controllers and criticality could be major mechanisms for the genesis of gait fractality.
Collapse
Affiliation(s)
- Chunjiang Fu
- Graduate School of Engineering Science, Osaka University, Osaka, 5608531, Japan
- Honda R&D Innovative Research Excellence, Wako, Japan
| | - Yasuyuki Suzuki
- Graduate School of Engineering Science, Osaka University, Osaka, 5608531, Japan
| | - Pietro Morasso
- Center for Human Technologies, Istituto Italiano di Tecnologia, 16152, Genoa, Italy
| | - Taishin Nomura
- Graduate School of Engineering Science, Osaka University, Osaka, 5608531, Japan.
| |
Collapse
|
13
|
Kovacs BA, Milton J, Insperger T. Virtual stick balancing: sensorimotor uncertainties related to angular displacement and velocity. ROYAL SOCIETY OPEN SCIENCE 2019; 6:191006. [PMID: 31827841 PMCID: PMC6894588 DOI: 10.1098/rsos.191006] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/04/2019] [Accepted: 11/01/2019] [Indexed: 05/22/2023]
Abstract
Sensory uncertainties and imperfections in motor control play important roles in neural control and Bayesian approaches to neural encoding. However, it is difficult to estimate these uncertainties experimentally. Here, we show that magnitude of the uncertainties during the generation of motor control force can be measured for a virtual stick balancing task by varying the feedback delay, τ. It is shown that the shortest stick length that human subjects are able to balance is proportional to τ 2. The proportionality constant can be related to a combined effect of the sensory uncertainties and the error in the realization of the control force, based on a delayed proportional-derivative (PD) feedback model of the balancing task. The neural reaction delay of the human subjects was measured by standard reaction time tests and by visual blank-out tests. Experimental observations provide an estimate for the upper boundary of the average sensorimotor uncertainty associated either with angular position or with angular velocity. Comparison of balancing trials with 27 human subjects to the delayed PD model suggests that the average uncertainty in the control force associated purely with the angular position is at most 14% while that associated purely with the angular velocity is at most 40%. In the general case when both uncertainties are present, the calculations suggest that the allowed uncertainty in angular velocity will always be greater than that in angular position.
Collapse
Affiliation(s)
- Balazs A. Kovacs
- Department of Applied Mechanics, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| | - John Milton
- W. M. Keck Science Department, The Claremont Colleges, Claremont, CA 91711, USA
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, Budapest, Hungary
| |
Collapse
|
14
|
McKee KL, Neale MC. Direct estimation of the parameters of a delayed, intermittent activation feedback model of postural sway during quiet standing. PLoS One 2019; 14:e0222664. [PMID: 31527893 PMCID: PMC6748412 DOI: 10.1371/journal.pone.0222664] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2019] [Accepted: 09/04/2019] [Indexed: 01/12/2023] Open
Abstract
Human postural sway during quiet standing has been characterized as a proportional-integral-derivative controller with intermittent activation. In the model, patterns of sway result from both instantaneous, passive, mechanical resistance and delayed, intermittent resistance signaled by the central nervous system. A Kalman-Filter framework was designed to directly estimate from experimental data the parameters of the model’s stochastic delay differential equations with discrete dynamic switching conditions. Simulations showed that all parameters could be estimated over a variety of possible data-generating configurations with varying degrees of bias and variance depending on their empirical identification. Applications to experimental data reveal distributions of each parameter that correspond well to previous findings, suggesting that many useful, physiological measures may be extracted from sway data. Individuals varied in degree and type of deviation from theoretical expectations, ranging from harmonic oscillation to non-equilibrium Langevin dynamics.
Collapse
Affiliation(s)
- Kevin L. McKee
- Virginia Commonwealth University, Virginia Institute of Psychiatric and Behavioral Genetics, Richmond, Virginia, United States of America
- * E-mail:
| | - Michael C. Neale
- Virginia Commonwealth University, Virginia Institute of Psychiatric and Behavioral Genetics, Richmond, Virginia, United States of America
| |
Collapse
|
15
|
Milton J, Insperger T. Acting together, destabilizing influences can stabilize human balance. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2019; 377:20180126. [PMID: 31329069 PMCID: PMC6661324 DOI: 10.1098/rsta.2018.0126] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 05/30/2019] [Indexed: 05/20/2023]
Abstract
The causes of falling in the elderly are multi-factorial. Three factors that influence balance stability are the time delay, a sensory dead zone and the maximum ankle torque that can be generated by muscular contraction. Here, the effects of these contributions are evaluated in the context of a model of an inverted pendulum stabilized by time-delayed proportional-derivative (PD) feedback. The effect of the sensory dead zone is to produce a hybrid type of control in which the PD feedback is switched ON or OFF depending on whether or not the controlled variable is larger or smaller than the detection threshold, Π. It is shown that, as Π increases, the region in the plane of control parameters where the balance time (BT) is greater than 60 s is increased slightly. However, when maximum ankle torque is also limited, there is a dramatic increase in the parameter region associated with BTs greater than 60 s. This increase is due to the effects of a torque limitation on over-control associated with bang-bang type switching controllers. These observations show that acting together influences, which are typically thought to destabilize balance, can actually stabilize balance. This article is part of the theme issue 'Nonlinear dynamics of delay systems'.
Collapse
Affiliation(s)
- John Milton
- W. M. Keck Science Center, The Claremont Colleges, Claremont, CA 91711, USA
- e-mail:
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology, and MTA-BME Lendület Human Balancing Research Group, 1111 Budapest, Hungary
| |
Collapse
|
16
|
Franklin DW, Cesonis J, Franklin S, Leib R. A Technique for Measuring Visuomotor Feedback Contributions to the Control of an Inverted Pendulum. ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY. ANNUAL INTERNATIONAL CONFERENCE 2019; 2019:1513-1516. [PMID: 31946181 DOI: 10.1109/embc.2019.8857119] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
We developed a new technique to measure the contributions of rapid visuomotor feedback responses to the stabilization of a simulated inverted pendulum. Human participants balanced an inverted pendulum simulated on a robotic manipulandum. At a random time during the balancing task, the visual representation of the tip of the pendulum was shifted by a small displacement to the left or right while the motor response was measured. This response was either the exerted force against a fixation position, or the motion to re-stabilize the pendulum in the free condition. Our results demonstrate that rapid involuntary visuomotor feedback responses contribute to the stabilization of the pendulum.
Collapse
|
17
|
Zhang L, Stepan G, Insperger T. Saturation limits the contribution of acceleration feedback to balancing against reaction delay. J R Soc Interface 2019; 15:rsif.2017.0771. [PMID: 29386400 DOI: 10.1098/rsif.2017.0771] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2017] [Accepted: 01/09/2018] [Indexed: 11/12/2022] Open
Abstract
A nonlinear model for human balancing subjected to a saturated delayed proportional-derivative-acceleration (PDA) feedback is analysed. Compared to the proportional-derivative (PD) controller, it is confirmed that the PDA controller improves local stability even for large feedback delays. However, it is shown that the saturated PDA controller typically introduces subcritical Hopf bifurcation into the system, which can also lead to falling for large enough perturbations. The subcriticality becomes stronger as the acceleration feedback gain increases or the saturation torque limit decreases. These explain some features of human balancing failure related to the increased reaction delay of inactive elderly people.
Collapse
Affiliation(s)
- Li Zhang
- State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, People's Republic of China
| | - Gabor Stepan
- Department of Applied Mechanics, Budapest University of Technology and Economics, 1521 Budapest, Hungary
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics, 1521 Budapest, Hungary.,Economics and MTA-BME Lendület Human Balancing Research Group, 1521 Budapest, Hungary
| |
Collapse
|
18
|
Morasso P, Nomura T, Suzuki Y, Zenzeri J. Stabilization of a Cart Inverted Pendulum: Improving the Intermittent Feedback Strategy to Match the Limits of Human Performance. Front Comput Neurosci 2019; 13:16. [PMID: 31024281 PMCID: PMC6461063 DOI: 10.3389/fncom.2019.00016] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Accepted: 03/11/2019] [Indexed: 11/13/2022] Open
Abstract
Stabilization of the CIP (Cart Inverted Pendulum) is an analogy to stick balancing on a finger and is an example of unstable tasks that humans face in everyday life. The difficulty of the task grows exponentially with the decrease of the length of the stick and a stick length of 32 cm is considered as a human limit even for well-trained subjects. Moreover, there is a cybernetic limit related to the delay of the multimodal sensory feedback (about 230 ms) that supports a feedback stabilization strategy. We previously demonstrated that an intermittent-feedback control paradigm, originally developed for modeling the stabilization of upright standing, can be applied with success also to the CIP system, but with values of the critical parameters far from the limiting ones (stick length 50 cm and feedback delay 100 ms). The intermittent control paradigm is based on the alternation of on-phases, driven by a proportional/derivative delayed feedback controller, and off-phases, where the feedback is switched off and the motion evolves according to the intrinsic dynamics of the CIP. In its standard formulation, the switching mechanism consists of a simple threshold operator: the feedback control is switched off if the current (delayed) state vector is closer to the stable than to the unstable manifold of the off-phase and is switched on in the opposite case. Although this simple formulation is effective for explaining upright standing as well as CIP balancing, it fails in the most challenging configuration of the CIP. In this work we propose a modification of the standard intermittent control policy that focuses on the explicit selection of switching times and is based on the phase reset of the estimated state vector at each switching time and on the simulation of an approximated internal model of CIP dynamics. We demonstrate, by simulating the modified intermittent control policy, that it can match the limits of human performance, while operating near the edge of instability.
Collapse
Affiliation(s)
- Pietro Morasso
- Robotics, Brain and Cognitive Sciences Department, Center for Human Technologies, Italian Institute of Technology, Genoa, Italy
| | - Taishin Nomura
- Mechanical Science and Bioengineering Department, Graduate School of Engineering Science, Osaka University, Toyonaka, Japan
| | - Yasuyuki Suzuki
- Mechanical Science and Bioengineering Department, Graduate School of Engineering Science, Osaka University, Toyonaka, Japan
| | - Jacopo Zenzeri
- Robotics, Brain and Cognitive Sciences Department, Center for Human Technologies, Italian Institute of Technology, Genoa, Italy
| |
Collapse
|
19
|
Morasso P, Cherif A, Zenzeri J. Quiet standing: The Single Inverted Pendulum model is not so bad after all. PLoS One 2019; 14:e0213870. [PMID: 30897124 PMCID: PMC6428281 DOI: 10.1371/journal.pone.0213870] [Citation(s) in RCA: 47] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Accepted: 03/01/2019] [Indexed: 11/19/2022] Open
Abstract
In the study of balance and postural control the (Single) Inverted Pendulum model (SIP) has been taken for a long time as an acceptable paradigm, with the implicit assumption that only ankle rotations are relevant for describing and explaining sway movements. However, more recent kinematic analysis of quiet standing revealed that hip motion cannot be neglected at all and that ankle-hip oscillatory patterns are characterized by complex in-phase and anti-phase interactions, suggesting that the SIP model should be substituted by a DIP (Double Inverted Pendulum) model. It was also suggested that DIP control could be characterized as a kind of optimal bi-axial active controller whose goal is minimizing the acceleration of the global CoM (Center of Mass). We propose here an alternative where active feedback control is applied in an intermittent manner only to the ankle joint, whereas the hip joint is stabilized by a passive stiffness mechanism. The active control impulses are delivered to the ankle joint as a function of the delayed state vector (tilt rotation angle + tilt rotational speed) of a Virtual Inverted Pendulum (VIP), namely a pendulum that links the ankle to the CoM, embedded in the real DIP. Simulations of such DIP/VIP model, with the hybrid control mechanism, show that it can reproduce the in-phase/anti-phase interaction patterns of the two joints described by several experimental studies. Moreover, the simulations demonstrate that the DIP/VIP model can also reproduce the measured minimization of the CoM acceleration, as an indirect biomechanical consequence of the dynamic interaction between the active control of the ankle joint and the passive control of the hip joint. We suggest that although the SIP model is literally false, because it ignores the ankle-hip coordination, it is functionally correct and practically acceptable for experimental studies that focus on the postural oscillations of the CoM.
Collapse
Affiliation(s)
- Pietro Morasso
- Department of Robotics, Brain and Cognitive Sciences, Istituto Italiano di Tecnologia, Genoa, Italy
| | - Amel Cherif
- Department of Robotics, Brain and Cognitive Sciences, Istituto Italiano di Tecnologia, Genoa, Italy
- Department of Informatics, Bioengineering, Robotics, and System Engineering, University of Genoa, Genoa, Italy
| | - Jacopo Zenzeri
- Department of Robotics, Brain and Cognitive Sciences, Istituto Italiano di Tecnologia, Genoa, Italy
- * E-mail:
| |
Collapse
|
20
|
Milton JG, Insperger T, Cook W, Harris DM, Stepan G. Microchaos in human postural balance: Sensory dead zones and sampled time-delayed feedback. Phys Rev E 2018; 98:022223. [PMID: 30253531 DOI: 10.1103/physreve.98.022223] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2017] [Indexed: 06/08/2023]
Abstract
Models for the stabilization of an inverted pendulum figure prominently in studies of human balance control. Surprisingly, fluctuations in measures related to the vertical displacement angle for quietly standing adults with eyes closed exhibit chaos. Here we show that small-amplitude chaotic fluctuations ("microchaos") can be generated by the interplay between three essential components of human neural balance control, namely time-delayed feedback, a sensory dead zone, and frequency-dependent encoding of force. When the sampling frequency of the force encoding is decreased, the sensitivity of the balance control to changes in the initial conditions increases. The sampled, time-delayed nature of the balance control may provide insights into why falls are more common in the very young and the elderly.
Collapse
Affiliation(s)
- John G Milton
- W. M. Keck Science Center, The Claremont Colleges, Claremont, California 91711, USA
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, 1111 Budapest, Hungary
| | - Walter Cook
- W. M. Keck Science Center, The Claremont Colleges, Claremont, California 91711, USA
| | - David Money Harris
- Department of Engineering, Harvey Mudd College, Claremont, California 91711, USA
| | - Gabor Stepan
- Department of Applied Mechanics, Budapest University of Technology and Economics, 1111 Budapest, Hungary
| |
Collapse
|
21
|
Stepan G, Milton JG, Insperger T. Quantization improves stabilization of dynamical systems with delayed feedback. CHAOS (WOODBURY, N.Y.) 2017; 27:114306. [PMID: 29195339 DOI: 10.1063/1.5006777] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We show that an unstable scalar dynamical system with time-delayed feedback can be stabilized by quantizing the feedback. The discrete time model corresponds to a previously unrecognized case of the microchaotic map in which the fixed point is both locally and globally repelling. In the continuous-time model, stabilization by quantization is possible when the fixed point in the absence of feedback is an unstable node, and in the presence of feedback, it is an unstable focus (spiral). The results are illustrated with numerical simulation of the unstable Hayes equation. The solutions of the quantized Hayes equation take the form of oscillations in which the amplitude is a function of the size of the quantization step. If the quantization step is sufficiently small, the amplitude of the oscillations can be small enough to practically approximate the dynamics around a stable fixed point.
Collapse
Affiliation(s)
- Gabor Stepan
- Department of Applied Mechanics, Budapest University of Technology and Economics, 1111 Budapest, Hungary
| | - John G Milton
- W. M. Keck Science Center, The Claremont Colleges, Claremont, California 91711, USA
| | - Tamas Insperger
- Department of Applied Mechanics, Budapest University of Technology and Economics and MTA-BME Lendület Human Balancing Research Group, 1111 Budapest, Hungary
| |
Collapse
|
22
|
Wang J, Kuske R. The influence of parametric and external noise in act-and-wait control with delayed feedback. CHAOS (WOODBURY, N.Y.) 2017; 27:114319. [PMID: 29195303 DOI: 10.1063/1.5006776] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
We apply several novel semi-analytic approaches for characterizing and calculating the effects of noise in a system with act-and-wait control. For concrete illustration, we apply these to a canonical balance model for an inverted pendulum to study the combined effect of delay and noise within the act-and-wait setting. While the act-and-wait control facilitates strong stabilization through deadbeat control, a comparison of different models with continuous vs. discrete updating of the control strategy in the active period illustrates how delays combined with the imprecise application of the control can seriously degrade the performance. We give several novel analyses of a generalized act-and-wait control strategy, allowing flexibility in the updating of the control strategy, in order to understand the sensitivities to delays and random fluctuations. In both the deterministic and stochastic settings, we give analytical and semi-analytical results that characterize and quantify the dynamics of the system. These results include the size and shape of stability regions, densities for the critical eigenvalues that capture the rate of reaching the desired stable equilibrium, and amplification factors for sustained fluctuations in the context of external noise. They also provide the dependence of these quantities on the length of the delay and the active period. In particular, we see that the combined influence of delay, parametric error, or external noise and on-off control can qualitatively change the dynamics, thus reducing the robustness of the control strategy. We also capture the dependence on how frequently the control is updated, allowing an interpolation between continuous and frequent updating. In addition to providing insights for these specific models, the methods we propose are generalizable to other settings with noise, delay, and on-off control, where analytical techniques are otherwise severely scarce.
Collapse
Affiliation(s)
- Jiaxing Wang
- Department of Mathematics, University of British Columbia, Vancouver, BC V6T1Z2 Canada
| | - Rachel Kuske
- Department of Mathematics, University of British Columbia, Vancouver, BC V6T1Z2 Canada
| |
Collapse
|
23
|
Tanabe H, Fujii K, Kouzaki M. Intermittent muscle activity in the feedback loop of postural control system during natural quiet standing. Sci Rep 2017; 7:10631. [PMID: 28878227 PMCID: PMC5587544 DOI: 10.1038/s41598-017-10015-8] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2017] [Accepted: 08/02/2017] [Indexed: 12/02/2022] Open
Abstract
The origin of continual body oscillation during quiet standing is a neural-muscular-skeletal closed feedback loop system that includes insufficient joint stiffness and a time delay. Thus, muscle activity and joint oscillations are nonlinear during quiet standing, making it difficult to demonstrate the muscular-skeletal relationship experimentally. Here we experimentally revealed this relationship using intermittent control theory, in which non-actuation works to stabilize the skeletal system towards equilibrium. We found that leg muscles were activated/inactivated when the state point was located in the opposite/same direction as the direction of anatomical action, which was associated with joint torque actuating the body towards equilibrium. The derivative values of stability index defined in the phase space approximately 200 ms before muscle inactivation were also larger than those before activation for some muscles. These results indicate that bipedal standing might be achieved by monitoring the rate of change of stability/instability components and generating joint torque to stabilize the body. In conclusion, muscles are likely to activate in an event-driven manner during quiet standing and a possible metric for on/off switching is SI dot, and our methodology of EMG processing could allows us to extract such event-driven intermittent muscle activities.
Collapse
Affiliation(s)
- Hiroko Tanabe
- Graduate School of Arts and Sciences, The University of Tokyo, 3-8-1 Komaba, Meguro, Tokyo, 153-8902, Japan.
| | - Keisuke Fujii
- Center for Advanced Intelligence Project, Institute of Physical and Chemical Research, 6-2-3 Furuedai, Suita, Osaka, 565-0874, Japan
| | - Motoki Kouzaki
- Graduate School of Human and Environmental Studies, Kyoto University, Yoshida-nihonmatsu, Sakyo-ku, Kyoto, 606-8501, Japan
| |
Collapse
|
24
|
Abstract
When we knock on a door, we perceive the impact as a collection of simultaneous events, combining sound, sight, and tactile sensation. In reality, information from different modalities but from a single source is flowing inside the brain along different pathways, reaching processing centers at different times. Therefore, interpreting different sensory modalities which seem to occur simultaneously requires information processing that accounts for these different delays. As in a computer-based robotic system, does the brain use some explicit estimation of the time delay, to realign the sensory flows? Or does it compensate for temporal delays by representing them as changes in the body/environment mechanics? Using delayed-state or an approximation for delayed-state manipulations between visual and proprioceptive feedback during a tracking task, we show that tracking errors, grip forces, and learning curves are consistent with predictions of a representation that is based on approximation for delay, refuting an explicit delayed-state representation. Delayed-state representations are based on estimating the time elapsed between the movement commands and their observed consequences. In contrast, an approximation for delay representations result from estimating the instantaneous relation between the expected and observed motion variables, without explicit reference to time.
Collapse
|
25
|
A negative group delay model for feedback-delayed manual tracking performance. J Comput Neurosci 2016; 41:295-304. [DOI: 10.1007/s10827-016-0618-4] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2016] [Revised: 07/27/2016] [Accepted: 08/03/2016] [Indexed: 11/25/2022]
|