Research Article  Open Access
Richard E. Hughes, "Using a Bayesian Network to Predict L5/S1 Spinal Compression Force from Posture, Hand Load, Anthropometry, and Disc Injury Status", Applied Bionics and Biomechanics, vol. 2017, Article ID 2014961, 7 pages, 2017. https://doi.org/10.1155/2017/2014961
Using a Bayesian Network to Predict L5/S1 Spinal Compression Force from Posture, Hand Load, Anthropometry, and Disc Injury Status
Abstract
Stochastic biomechanical modeling has become a useful tool most commonly implemented using Monte Carlo simulation, advanced mean value theorem, or Markov chain modeling. Bayesian networks are a novel method for probabilistic modeling in artificial intelligence, risk modeling, and machine learning. The purpose of this study was to evaluate the suitability of Bayesian networks for biomechanical modeling using a static biomechanical model of spinal forces during lifting. A 20node Bayesian network model was used to implement a wellestablished static twodimensional biomechanical model for predicting L5/S1 compression and shear forces. The model was also implemented as a Monte Carlo simulation in MATLAB. Mean L5/S1 spinal compression force estimates differed by 0.8%, and shear force estimates were the same. The model was extended to incorporate evidence about disc injury, which can modify the prior probability estimates to provide posterior probability estimates of spinal compression force. An example showed that changing disc injury status from false to true increased the estimate of mean L5/S1 compression force by 14.7%. This work shows that Bayesian networks can be used to implement a wholebody biomechanical model used in occupational biomechanics and incorporate disc injury.
1. Introduction
Stochastic modeling has become a useful analytical method in biomechanics over the last twenty years. Applications have included lumbar region muscle forces [1], musculoskeletal injury risk prediction [2–4], upper extremity joint mechanics [5–8] muscle modeling [9], population morphological modeling [10], probabilistic sensitivity analyses [11], tissue repair [12], and orthopaedic implant design [13, 14]. Past stochastic biomechanical models use Monte Carlo simulation, advanced mean value theorem, or Markov chain models. Unfortunately, these methods are all “forward” simulation methods that produce probability estimates of outcomes based on model inputs represented as random variables.
Bayesian networks are well established in artificial intelligence [15, 16], risk analysis [17], and reliability engineering [18–21] and have the potential to enhance biomechanical modeling. A Bayesian network is a graphical probabilistic model containing nodes and directed edges that can be used to compute probabilities when evidence has been entered at a variety of nodes in the graph, some of which are descendants of the primary variable of interest [22]. Consider the problem of predicting lumbar spine compression force during lifting, which is a common task in occupational biomechanics. Existing biomechanical models require taskrelated inputs like mass in the hands, anthropometry, and joint angles to predict compression force [23, 24]. Physiologically driven models also require electromyographic measurements [25–28]. The results of these models are predictions of spinal compression and shear forces. In theory, a Bayesian network could be used to also incorporate intervertebral disc injury status, which existing models cannot.
Therefore, the purpose of this project was to evaluate the feasibility of implementing a wellestablished static twodimensional biomechanical model of lifting as a Bayesian network and extend it to include disc prolapse as a model input.
2. Methods
The design of this study consisted of five steps: (1) implement a wellestablished deterministic static twodimensional biomechanical model for predicting L5/S1 shear and compression force as a Bayesian network, (2) extend the Bayesian network to include uncertainty in model inputs, (3) verify the implementation, (4) augment model to include disc injury status, and (5) evaluate effect of including disc injury status on L5/S1 compression force estimates.
A wellestablished deterministic model for predicting L5/S1 compression force during lifting was selected as the basis for this model [29] (pp. 130131). This model computes intersegmental reaction forces and moments at the elbow, shoulder, and L5/S1 disc. Body segment lengths and masses are scaled from status and total body mass, respectively. The location of the center of mass for each segment is scaled to segment length. All necessary body segment anthropometric parameters are described in Chaffin et al. [29] (pp. 41–47). The model allows for a downward directed hand force representing the holding of an object in the hands. The model is bilaterally symmetric. A 50thpercentile female (161.8 cm stature and 65.6 kg body mass) was used for all simulations. The erector spinae moment arm was 5.3 cm according to Chaffin et al.’s recommendation for models not incorporating intraabdominal pressure (p. 134).
The deterministic and stochastic versions of the lifting model were also coded in MATLAB R2014a (The Mathworks, Natick, MA) for purposes of validation. Monte Carlo simulation was used for the stochastic simulation, and each simulation contained 10^{7} iterations. Hand calculations were also performed for selected deterministic cases for code verification.
The next step was to implement the same model as a Bayesian network [30] using AgenaRisk software (Agena, Cambridge, UK). A Bayesian network is a graphical probabilistic modeling framework developed initially in artificial intelligence [15, 16]. A Bayesian network is a directed acyclic graph in which the directed edges represent conditional independence assumptions and node probability tables are associated with each node. Each node represents a random variable. The probability table associated with a node consists of a full statement of conditional probabilities for all nodes that are parents of the node. It is also possible to code deterministic relationships between nodes. Probability distributions can be defined for nodes, as well as observed values of the random variables. Junction tree methods are used to propagate probability distributions through the network [31].
In the Bayesian network implementation of the lifting model, nodes for input data were defined: mass in hands, elbow angle, shoulder angle, torso angle, knee angle, and ankle angle (Figure 1). Note that in this paper, variables are indicated in italics. Intermediate nodes were also defined for purposes of computation. Nodes were introduced for the reaction forces and moments at the elbow (elbow reaction force and elbow reaction moment), shoulder (shoulder reaction force and shoulder reaction moment), and L5/S1 level (L5/S1 reaction force and L5/S1 reaction moment). A node was created for the erector spinae moment arm, and another for erector spine force. The effect if intraabdominal pressure was excluded because of the nonlinearity of the relationship between hip moment and pressure (p. 132). A node for L5/S1 disc angle was included, as were three nodes (β, T, and K) to implement a regression model for predicting deviation of disc angle from 40 degrees [32], which is how Chaffin et al. compute disc angle (p. 132). Two output nodes were defined: (1) L5/S1 compression force and (2) L5/S1 shear force. Edges were introduced to link the nodes in a way that represented mechanics. In AgenaRisk, each arc was given the deterministic relationships linking nodes. This defined the deterministic version of the lifting model.
The model was extended to a stochastic model by modifying the input nodes to be normally distributed random variables rather than deterministic values. Liu et al. [33] showed that assessment of body segment angles for a lifting model using common ergonomics tools has measurement error and these error estimates were used to define the standard deviation for each joint angle. Since Liu et al. [33] reported mean absolute errors, normality was assumed and the result of Geary [34] was used to compute standard deviation for the elbow, shoulder, and torso angles. The standard deviation for the error across all joints was also computed, and this was used for the knee and ankle.
The method was evaluated on a lifting task obtained from a study of heavy work that has been previously reported [35]. The task involved lifting chunks of carbon that had broken off large anodes in an aluminum smelter. Posture was collected from video of the job taken from a single viewpoint. The postural measurements were subject to error as shown by Liu et al. [33], and hand load was random because the carbon pieces broke off from anodes by a process that created a wide distribution of masses. Therefore, both posture and hand load were modeled as normally distributed random variables. Table 1 provides the mean (and standard deviation) joint angles. The mass of the object in the hands was modeled as a normally distributed random variable having mean and standard deviation of 53.1 kg and 12.8 kg, respectively.

A modified model was developed that incorporates intervertebral disc prolapse status (a socalled “injuryaugmented” model). It is the full DAG shown in Figure 1. It was developed by adding a probabilistic inference model of structural failure in which failure occurs when the applied load exceeds the compressive strength. Both applied load and strength were modeled as random variables [36]. Therefore, a node for intervertebral disc compressive strength was added. The distribution of compressive strength was assumed to be normally distributed. The mean (5448 N) and standard deviation (2366 N) values were taken from reports of cadaver testing [37]. A disc prolapse Boolean variable node was added that took on a true value if and only if the compression force exceeded the disc compressive strength. For the numerical example analyzed to illustrate the method, separate simulations were performed for three injury conditions: (1) no disc injury evidence entered, (2) disc injury Boolean variable set to false, and (3) disc injury Boolean variable set to true. The model also included modeling postures and hand loads as random variables as described above.
3. Results
The deterministic model implemented as a Bayesian network produced L5/S1 compression force estimates that were within 1.4 N (0.04%) of the MATLAB implementation predictions (Table 2). Similarly, both deterministic model versions produced the same L5/S1 shear force estimate. This demonstrates verification of the deterministic implementation. The Bayesian network basic model in which the joint angles were modeled as random variables produced greater average L5/S1 compression force than the MATLAB implementation. However, the difference in means was only 28.5 N (0.8%). When hand load was simulated as a random variable, the difference in means was only 3.8 N (0.1%). For the random posture case, the 25th percentile, median, and 75th percentile of the compression force distribution were within 290.2 N (4.5%), 1.6 N (0.05%), and 261.9 N (7.0%) of the MATLAB code Monte Carlo implementation results.

The injuryaugmented model’s predicted L5/S1 compression force changed as evidence was entered at the injury node. When evidence of no injury (disc injury = false) was entered at the injury node, the mean compression force was 3190.2 N (SD 789.3). However, when evidence of injury (disc injury = true) was entered, the predicted mean compression force increased to 3660.9 N (SD 820.0). This was an increase of 470.7 N (14.7%). Moreover, the injury evidence shifted the entire distribution. The 25th percentile value went from 2664.1 N to 3095.9 N when disc injury was changed from false to true; the 75th percentile similarly increased from 3690.1 N to 4180.1 N going from false to true. Predicted mean L5/S1 shear force also changed from 673.1 N (SD 108.9) to 715.7 N (SD 111.0) when disc injury was changed from false to true.
4. Discussion
This project demonstrated the ability of implementing a biomechanical force prediction as a Bayesian network. The error between mean L5/S1 compression force estimates was small (<1%), and it was lower when modeling hand load as a random variable than when modeling joint angles as random variables. While Bayesian networks have had limited use in biomechanical data analysis, they have not been used to implement wholebody biomechanical models.
A Bayesian network is a directed acyclic graph (DAG) whose nodes encode conditional probability statements and edges reflect conditional independence assumptions. Directed edges can be thought of as indicating causal relationships [38]. Directed paths in the DAG can also represent causal structures. In the first model developed (Figure 1), a DAG is constructed so that directed paths go from input variables (mass in hands, elbow angle, etc.) and internal structural variables (erector spinae moment arm, L5/S1 disc angle, etc.) to low back kinetic variables (L5/S1 shear force and L5/S1 compression force). This “basic model,” which merely implements the standard model of Chaffin, is the DAG shown outside of the dashed rectangle in Figure 1. Calculations on this network are straightforward and do not require Bayes’ theorem. However, the full power of a Bayesian network is apparent when Bayes’ theorem is used to calculate probabilities backwards along directed paths. The injuryaugmented model (entire network in Figure 1) does this. Moreover, the injuryaugmented model uses both the input data required for the basic model (Figure 1) and disc prolapse status (disc injury node in Figure 1 denoting the Boolean variable describing whether the person being analyzed has a prolapsed disc or not). The difference was 14.7%, which demonstrates the impact incorporating disc injury status can have on spinal compression force estimates.
Another interesting feature of calculations performed on a Bayesian network is how the disc injury variable affects L5/S1 shear force. Note that the disc injury node is a descendant of the L5/S1 shear force node, yet changing disc injury from false to true increases the mean L5/S1 shear force from 673.1 N to 715.7 N. If it is a descendant, how can it affect L5/S1 shear force? It is because there is—in terminology from graphical causal modeling—a “trek” from disc injury to L5/S1 shear force that does not contain a collider node [39]. Changing the value of disc injury from false to true affects L5/S1 reaction force through L5/S1 compression force because of Bayes’ theorem. In turn, L5/S1 reaction force affects L5/S1 shear force. A similar effect could occur via the L5/S1 disc angle node.
Existing stochastic simulation methods are forward methods that propagate probability from distributions of input parameters to distributions of outputs. Previously, Monte Carlo simulation has been the most common tool for introducing probabilistic elements into deterministic biomechanical models [1–4, 7, 8, 13, 14, 40]. While intuitively easy to implement, Monte Carlo techniques are computationally intensive. Laz et al. [14] introduced the use of the advanced mean value theorem approach to probabilistic simulation in biomechanics. Markov chain models have also been used for stochastic modeling of softtissue failure [12]. In Bayesian network parlance, all of these methods place the input variables on nodes that are ancestors of the variable to be predicted. Variables that are graphtheoretic descendants of the predicted variable cannot influence the predictions. This makes it impossible to include injury status into the calculation of spinal compression and shear forces using existing probabilistic methods.
The most significant limitation of the Bayesian network implementation of the model was the difficulty in predicting the tails of the spinal compression force distribution. Specifically, there was a difference of 7% for the 75th percentile. An additional limitation was that the Bayesian network produced a better estimate of the mean compression force when the hand load was modeled as a random variable than when the joint angles were random. This is because the effect of joint angles on spinal compression force is nonlinear.
Comparisons of these results to other Bayesian network models in biomechanics are difficult because so few applications of Bayesian networks exist in biomechanics. Van Gestel et al. [41] and Lo et al. [42] used a Bayesian network for gait classification, but neither used the Bayesian network to implement a biomechanical model. Ma et al. [43] used a Bayesian network to analyze spinal kinematics in patients having low bone mineral density. Again, no biomechanical model was used. Bayesian networks have also been used in the analysis of clinical data related to degenerative disease of the lumbar spine. Takenaka and Aono [44] used a Bayesian network to develop a classifier for lower extremity muscle strength for patients presenting with drop foot following lumbar spinal decompression surgery.
A different comparison to existing literature can be made by noting the similarity to work done in the structural reliability area. The model developed here considers disc prolapse to occur when the applied load exceeds the compressive strength of the intervertebral disc to withstand load [36]. In this context, the current model can be viewed as using a stochastic biomechanical model to predict the load and cadaver testing to provide strength data. Thus, there is a conceptual connection to work implementing interference models in structural reliability models using Bayesian networks [45]. Donnell et al. [12] have also used the interference model of structural failure in biomechanics, although not implemented as a Bayesian network.
This model takes a simplistic view of spinal injury mechanics to highlight the Bayesian network modeling technique; future work should include subtle factors. The disc compression limit came from a study [37] that tested spinal motion segments in mild flexion, matching the lifting posture simulated. However, it used a single load to failure test paradigm, where other authors argue that cyclic loading is a very important factor in intervertebral disc failure [46, 47]. The current model also only considered compression in disc injury; the effects of shear should also be investigated [48].
In conclusion, this is a novel contribution to low back modeling because heretofore, biomechanical models predicted spinal compression force from posture and hand load without consideration of disc injury status. It is Bayes’ theorem that allows for disc injury, which is a graphtheoretic descendant of L5/S1 compression force on the causal directed path, to influence the estimate of L5/S1 compression force. This innovation may be useful in areas such as litigation, for example, where the analysis is being done to estimate whether compression force exceeds a recommended threshold for a specific person who has experienced a prolapsed intervertebral disc.
Conflicts of Interest
There are no competing interests.
Acknowledgments
This work was supported by the Department of Orthopaedic Surgery of the University of Michigan.
References
 G. A. Mirka and W. S. Marras, “A stochastic model of trunk muscle coactivation during trunk bending,” Spine, vol. 18, no. 11, pp. 1396–1409, 1993. View at: Google Scholar
 C. F. Lin, M. Gross, C. Ji et al., “A stochastic biomechanical model for risk and risk factors of noncontact anterior cruciate ligament injuries,” Journal of Biomechanics, vol. 42, no. 4, pp. 418–423, 2009. View at: Publisher Site  Google Scholar
 S. G. McLean, X. Huang, and A. J. van den Bogert, “Investigating isolated neuromuscular control contributions to noncontact anterior cruciate ligament injury risk via computer simulation methods,” Clinical Biomechanics, vol. 23, no. 7, pp. 926–936, 2008. View at: Publisher Site  Google Scholar
 P. L. Davidson, D. J. Chalmers, and B. D. Wilson, “Stochasticrheological simulation of freefall arm impact in children: application to playground injuries,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 7, no. 2, pp. 63–71, 2004. View at: Publisher Site  Google Scholar
 N. G. Flieg, C. J. Gatti, L. C. Doro, J. E. Langenderfer, J. E. Carpenter, and R. E. Hughes, “A stochastic analysis of glenoid inclination angle and superior migration of the humeral head,” Clinical Biomechanics, vol. 23, no. 5, pp. 554–561, 2008. View at: Publisher Site  Google Scholar
 J. E. Langenderfer, R. E. Hughes, and J. E. Carpenter, “A stochastic model of elbow flexion strength for subjects with and without long head biceps tear,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 8, no. 5, pp. 315–322, 2005. View at: Publisher Site  Google Scholar
 J. E. Langenderfer, J. E. Carpenter, M. E. Johnson, K. N. An, and R. E. Hughes, “A probabilistic model of glenohumeral external rotation strength for healthy normals and rotator cuff tear cases,” Annals of Biomedical Engineering, vol. 34, no. 3, pp. 465–476, 2006. View at: Publisher Site  Google Scholar
 F. J. ValeroCuevas, M. E. Johanson, and J. D. Towles, “Towards a realistic biomechanical model of the thumb: the choice of kinematic description may be more critical than the solution method or the variability/uncertainty of musculoskeletal parameters,” Journal of Biomechanics, vol. 36, no. 7, pp. 1019–1030, 2003. View at: Publisher Site  Google Scholar
 J. N. ChoppHurley, J. E. Langenderfer, and C. R. Dickerson, “Probabilistic evaluation of predicted force sensitivity to muscle attachment and glenohumeral stability uncertainty,” Annals of Biomedical Engineering, vol. 42, no. 9, pp. 1867–1879, 2014. View at: Publisher Site  Google Scholar
 J. N. ChoppHurley, J. E. Langenderfer, and C. R. Dickerson, “A probabilistic orthopaedic population model to predict fatiguerelated subacromial geometric variability,” Journal of Biomechanics, vol. 49, 2016. View at: Publisher Site  Google Scholar
 J. E. Langenderfer, P. J. Rullkoetter, A. G. Mell, and P. J. Laz, “A multisubject evaluation of uncertainty in anatomical landmark location on shoulder kinematic description,” Computer Methods in Biomechanics and Biomedical Engineering, vol. 12, no. 2, pp. 211–216, 2009. View at: Publisher Site  Google Scholar
 D. M. S. Donnell, J. L. Seidelman, C. L. Mendias, B. S. Miller, J. E. Carpenter, and R. E. Hughes, “A stochastic structural reliability model explains rotator cuff repair retears,” International Biomechanics, vol. 1, no. 1, pp. 29–35, 2014. View at: Publisher Site  Google Scholar
 S. K. Easley, S. Pal, P. R. Tomaszewski, A. J. Petrella, P. J. Rullkoetter, and P. J. Laz, “Finite elementbased probabilistic analysis tool for orthopaedic applications,” Computer Methods and Programs in Biomedicine, vol. 85, no. 1, pp. 32–40, 2007. View at: Publisher Site  Google Scholar
 P. J. Laz, S. Pal, J. P. Halloran, A. J. Petrella, and P. J. Rullkoetter, “Probabilistic finite element prediction of knee wear simulator mechanics,” Journal of Biomechanics, vol. 39, no. 12, pp. 2303–2310, 2006. View at: Publisher Site  Google Scholar
 J. Pearl, “Fusion, propagation, and structuring in belief networks,” Artificial Intelligence, vol. 29, pp. 241–288, 1986. View at: Publisher Site  Google Scholar
 J. Pearl, Probabilistic Reasoning in Intelligent Systems: Networks of Plausible Inference, Morgan Kaufmann, San Francisco, CA, USA, 1988.
 M. Neil, N. Fenton, and M. Tailor, “Using Bayesian networks to model expected and unexpected operational losses,” Risk Analysis, vol. 25, no. 4, pp. 963–972, 2005. View at: Publisher Site  Google Scholar
 M. Neil, N. Fenton, S. Forey, N. Fenton, and R. Harris, “Using Bayesian belief networks to predict the reliability of military vehicles,” Computing and Control Engineering Journal, vol. 12, no. 1, pp. 11–20, 2001. View at: Publisher Site  Google Scholar
 A. Bobbio, L. Portinale, M. Minichino, and E. Ciancamerla, “Improving the analysis of dependable systems by mapping fault trees into Bayesian networks,” Reliability Engineering and System Theory, vol. 71, pp. 249–260, 2001. View at: Publisher Site  Google Scholar
 H. Langseth and L. Portinale, “Applications of Bayesian networks in reliability analysis,” in Bayesian Network Technologies: Applications and Graphical Models, A. Mittal and A. Kassim, Eds., pp. 84–102, IGI Publishing, New York, NY, USA, 2007. View at: Publisher Site  Google Scholar
 H. Langseth and L. Portinale, “Bayesian networks in reliability,” Reliability Engineering & System Safety, vol. 92, no. 1, pp. 92–108, 2007. View at: Publisher Site  Google Scholar
 D. Koller and N. Friedman, Probabilistic Graphical Models: Principles and Techniques, MIT Press, Cambridge, MA, USA, 2009.
 D. B. Chaffin, “A computerized biomechanical model—development of and use in studying gross body actions,” Journal of Biomechanics, vol. 2, pp. 429–441, 1969. View at: Publisher Site  Google Scholar
 D. B. Chaffin and M. Erig, “Threedimensional biomechanical static strength prediction model sensitivity to postural and anthropometric inaccuracies,” IIE Transactions, vol. 23, no. 3, pp. 215–227, 1991. View at: Publisher Site  Google Scholar
 K. P. Granata and W. S. Marras, “An EMGassisted model of trunk loading during freedynamic lifting,” Journal of Biomechanics, vol. 28, no. 11, pp. 1309–1317, 1995. View at: Publisher Site  Google Scholar
 M. A. Nussbaum and D. B. Chaffin, “Lumbar muscle force estimation using a subjectinvariant 5parameter EMGbased model,” Journal of Biomechanics, vol. 31, pp. 667–672, 1998. View at: Publisher Site  Google Scholar
 J. Cholewicki, S. M. McGill, and R. W. Norman, “Comparison of muscle forces and joint load from an optimization and EMG assisted lumbar spine model: towards development of a hybrid approach,” Journal of Biomechanics, vol. 28, no. 3, pp. 321–331, 1995. View at: Publisher Site  Google Scholar
 S. M. McGill and R. W. Norman, “1986 Volvo award in biomechanics: partitioning of the L4  L5 dynamic moment into disc, ligamentous, and muscular components during lifting,” Spine, vol. 11, no. 7, pp. 666–678, 1986. View at: Publisher Site  Google Scholar
 D. B. Chaffin, G. B. J. Andersson, and B. J. Martin, Occupational Biomechanics, WileyInterscience, Hoboken, NJ, USA, 2006.
 N. Fenton and M. Neil, Risk Assessment and Decision Analysis with Bayesian Networks, CRC Press, Boca Raton, FL, USA, 2013.
 S. L. Lauritzen and D. J. Spiegelhalter, “Local computations with probabilities on graphical structures and their application to expert systems,” Journal of the Royal Statistical Society, Series B, vol. 50, no. 2, pp. 157–224, 1988. View at: Google Scholar
 C. K. Anderson, D. B. Chaffin, and G. D. Herrin, “A study of lumbosacral orientation under varied static loads,” Spine, vol. 11, no. 5, pp. 456–462, 1986. View at: Publisher Site  Google Scholar
 Y. Liu, X. Zhang, and D. Chaffin, “Preception and visualization of human posture information for computeraided ergonomic analysis,” Ergonomics, vol. 40, no. 8, pp. 818–833, 1997. View at: Publisher Site  Google Scholar
 R. C. Geary, “The ratio of the mean deviation to the standard deviation as a test of normality,” Biometrika, vol. 27, no. 3/4, pp. 310–332, 1935. View at: Publisher Site  Google Scholar
 R. E. Hughes, “Choice of optimization models for predicting spinal forces in a threedimensional analysis of heavy work,” Ergonomics, vol. 38, no. 12, pp. 2476–2484, 1995. View at: Publisher Site  Google Scholar
 P. ThoftChristensen and M. J. Baker, Structural Reliability Theory and its Applications, Springer, New York, NY, USA, 1982. View at: Publisher Site
 M. A. Adams and W. C. Hutton, “Prolapsed intervertebral disc. A hyperflexion injury 1981 Volvo award in basic science,” Spine, vol. 7, no. 3, pp. 184–191, 1982. View at: Google Scholar
 J. Pearl, Causality: Models, Reasoning, and Inference, Cambridge University Press, New York, NY, USA, 2009.
 P. Spirtes, C. Glymour, and R. Scheines, Causation, Prediction, and Search, MIT Press, London, 2000.
 P. J. Laz, J. Q. Stowe, M. A. Baldwin, A. J. Petrella, and P. J. Rullkoetter, “Incorporating uncertainty in mechanical properties for finite elementbased evaluation of bone mechanics,” Journal of Biomechanics, vol. 40, no. 13, pp. 2831–2836, 2007. View at: Publisher Site  Google Scholar
 L. Van Gestel, T. De Laet, E. Di Lello et al., “Probabilistic gait classification in children with cerebral palsy: a Bayesian approach,” Research in Developmental Disabilities, vol. 32, no. 6, pp. 2542–2552, 2011. View at: Publisher Site  Google Scholar
 B. Lo, J. Pansiot, and G.Z. Yang, “Bayesian analysis of subplantar ground reaction force with BSN,” in Sixth International Workshop on Wearable and Implantable Body Sensor Networks, 2009. BSN 2009, Berkeley, CA, USA, 2009, The Institute of Electrical and Electronics Engineers. View at: Publisher Site  Google Scholar
 H. T. Ma, J. F. Griffith, Z. Yang, A. W. Kwok, P. C. Leung, and R. Y. Lee, “Kinematics of the lumbar spine in elderly subjects with decreased bone mineral density,” Medical & Biological Engineering & Computing, vol. 47, no. 7, pp. 783–789, 2009. View at: Publisher Site  Google Scholar
 S. Takenaka and H. Aono, “Prediction of postoperative clinical recovery of drop foot attributable to lumbar degenerative diseases, via a Bayesian network,” Clinical Orthopaedics and Related Research, vol. 475, no. 3, pp. 872–880, 2017. View at: Publisher Site  Google Scholar
 D. Straub and A. Der Kiureghian, “Bayesian network enhanced with structural reliability methods: methodology,” Journal of Engineering Mechanics, vol. 136, no. 10, pp. 1248–1258, 2010. View at: Publisher Site  Google Scholar
 S. Gallagher, W. S. Marras, A. S. Litsky, and D. Burr, “An exploratory study of loading and morphometric factors associated with specific failure modes in fatigue testing of lumbar motion segments,” Clinical Biomechanics, vol. 21, no. 3, pp. 228–234, 2006. View at: Publisher Site  Google Scholar
 R. J. Parkinson and J. P. Callaghan, “The role of dynamic flexion in spine injury is altered by increasing dynamic load magnitude,” Clinical Biomechanics, vol. 24, no. 2, pp. 148–154, 2009. View at: Publisher Site  Google Scholar
 V. R. Yingling and S. M. McGill, “Mechanical properties and failure mechanics of the spine under posterior shear load: observations from a porcine model,” Journal of Spinal Disorders, vol. 12, no. 6, pp. 501–508, 1999. View at: Google Scholar
Copyright
Copyright © 2017 Richard E. Hughes. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.