Reducing the Complexity of Musculoskeletal Models Using Gaussian Process Emulators

Musculoskeletal models (MSKMs) are used to estimate the muscle and joint forces involved in human locomotion, often associated with the onset of degenerative musculoskeletal pathologies (e.g., osteoarthritis). Subject-specific MSKMs offer more accurate predictions than their scaled-generic counterparts. This accuracy is achieved through time-consuming personalisation of models and manual tuning procedures that suffer from potential repeatability errors, hence limiting the wider application of this modelling approach. In this work we have developed a methodology relying on Sobol’s sensitivity analysis (SSA) for ranking muscles based on their importance to the determination of the joint contact forces (JCFs) in a cohort of older women. The thousands of data points required for SSA are generated using Gaussian Process emulators, a Bayesian technique to infer the input–output relationship between nonlinear models from a limited number of observations. Results show that there is a pool of muscles whose personalisation has little effects on the predictions of JCFs, allowing for a reduced but still accurate representation of the musculoskeletal system within shorter timeframes. Furthermore, joint forces in subject-specific and generic models are influenced by different sets of muscles, suggesting the existence of a model-specific component to the sensitivity analysis.