Are subject-specific musculoskeletal models robust to the uncertainties in parameter identification?

• Published in: PloS one Vol. 9; no. 11; p. e112625
• Main Authors: Valente, Giordano, Pitto, Lorenzo, Testi, Debora, Seth, Ajay, Delp, Scott L, Stagni, Rita, Viceconti, Marco, Taddei, Fulvia
• Format: Journal Article
• Language: English
• Published: United States Public Library of Science 12.11.2014; Public Library of Science (PLoS)
• Subjects: Adult; Angles (geometry); Ankle; Bioengineering; Biology and Life Sciences; Biomechanical Phenomena - physiology; Biomechanics; Biomedical engineering; Computer and Information Sciences; Computer Simulation; Contact angle; Engineering and Technology; Gait; Gait - physiology; Geometry; Human mechanics; Humans; Inverse dynamics; Joint surgery; Joints - physiology; Knee Joint - physiology; Laboratories; Magnetic Resonance Imaging; Male; Mathematical models; Mechanical engineering; Medical research; Medical technology; Models, Biological; Monte Carlo simulation; Muscle, Skeletal - physiology; Muscles; Musculoskeletal diseases; Optimization; Parameter identification; Parameter sensitivity; Parameter uncertainty; Parameters; Posture; Probabilistic analysis; Sensitivity analysis; Statistical analysis; Tension; Walking; Walking - physiology
• ISSN: 1932-6203; 1932-6203
• DOI: 10.1371/journal.pone.0112625

Subject-specific musculoskeletal modeling can be applied to study musculoskeletal disorders, allowing inclusion of personalized anatomy and properties. Independent of the tools used for model creation, there are unavoidable uncertainties associated with parameter identification, whose effect on model predictions is still not fully understood. The aim of the present study was to analyze the sensitivity of subject-specific model predictions (i.e., joint angles, joint moments, muscle and joint contact forces) during walking to the uncertainties in the identification of body landmark positions, maximum muscle tension and musculotendon geometry. To this aim, we created an MRI-based musculoskeletal model of the lower limbs, defined as a 7-segment, 10-degree-of-freedom articulated linkage, actuated by 84 musculotendon units. We then performed a Monte-Carlo probabilistic analysis perturbing model parameters according to their uncertainty, and solving a typical inverse dynamics and static optimization problem using 500 models that included the different sets of perturbed variable values. Model creation and gait simulations were performed by using freely available software that we developed to standardize the process of model creation, integrate with OpenSim and create probabilistic simulations of movement. The uncertainties in input variables had a moderate effect on model predictions, as muscle and joint contact forces showed maximum standard deviation of 0.3 times body-weight and maximum range of 2.1 times body-weight. In addition, the output variables significantly correlated with few input variables (up to 7 out of 312) across the gait cycle, including the geometry definition of larger muscles and the maximum muscle tension in limited gait portions. Although we found subject-specific models not markedly sensitive to parameter identification, researchers should be aware of the model precision in relation to the intended application. In fact, force predictions could be affected by an uncertainty in the same order of magnitude of its value, although this condition has low probability to occur.