This content is not included in your SAE MOBILUS subscription, or you are not logged in.
Statistical Simulations to Evaluate the Methods of the Construction of Injury Risk Curves
Published November 07, 2011 by The Stapp Association in United States
Annotation ability available
Several statistical methods are currently used to build injury risk curves in the biomechanical field. These methods include the certainty method (Mertz et al. 1996), Mertz/Weber method (Mertz and Weber 1982), logistic regression (Kuppa et al. 2003, Hosmer and Lemeshow 2000), survival analysis with Weibull distribution (Kent et al. 2004, Hosmer and Lemeshow 2000), and the consistent threshold estimate (CTE) (Nusholtz et al. 1999, Di Domenico and Nusholtz 2005). There is currently no consensus on the most accurate method to be used and no guidelines to help the user to choose the more appropriate one. Injury risk curves built for the WorldSID 50th side impact dummy with these different methods could vary significantly, depending on the sample considered (Petitjean et al. 2009). As a consequence, further investigations were needed to determine the fields of application of the different methods and to recommend the best statistical method depending on the biomechanical sample considered.
This study used statistical simulations on theoretical samples to address these questions. Two different theoretical distributions of injury thresholds were utilized to assess the five different methods of constructing injury risk curves. A normal distribution and a Weibull distribution, whose shape was not similar to a normal distribution, were selected. One hundred sets of “test subjects” were randomly chosen from each theoretical distribution, with sample sizes ranging from 10 to 50. A “stimulus value” was chosen for each “test subject.” The stimulus values were equally spaced, distributed tightly or loosely about the theoretical mean injury threshold, concentrated below the mean value, or concentrated above the mean value. An adaptive method was also used to assign stimulus values, based on the proportion of uninjured and injured in an early subset of the test subjects. The influences of 10%, 25%, and 50% exact data were compared to stimulus values that were either right or left censored. The test subject was considered to be uninjured if the stimulus was less than the subject's threshold or injured if the stimulus was equal to or greater than the subject's threshold. In all, 12,800 simulated data sets with both normal and Weibull distributions were used to construct injury risk curves by each of the five statistical methods. Cumulative errors of the constructed injury risk curves, compared to the theoretical curves, were calculated across the whole curve, as well as the portion of the theoretical curve up to 25% risk of injury. P-values were used to assess the significance of the differences in the errors.
The CTE and the survival analysis take into account the exact data whatever the theoretical distribution of injury threshold, while the logistic regression, the Mertz/Weber and the certainty methods do not. For left and right censored data, the logistic regression and/or the survival analysis lead to the lowest error. The survival analysis leads to the lowest error whatever the sample size, the level of censoring and the theoretical distribution evaluated. Increasing the sample size generally decreased the error. However, the benefit from increasing the sample size decreased when the sample size was already high. For the survival analysis, increasing the proportion of exact data decreased the error. The same way, the benefit from increasing the proportion of exact data decreased when the proportion of exact data was already high. Survival analysis may not converge for small sample size with left and right censored data. The number of simulations for which the survival analysis did not converge highly decreases with the increase of proportion of exact data and the increase of the sample size.
Therefore, it is recommended to use survival analysis with Weibull distribution to build risk curves compared to the four other statistical methods evaluated. The accuracy of survival analysis with other distributions (log-normal, log-logistic, etc) was not studied. There is no recommendation for the method to be used when survival analysis does not converge. The balance between maximal acceptable error and the need for an injury risk curve, even for a small dataset of poor quality, is not addressed.
|Technical Paper||How Few? Bayesian Statistics in Injury Biomechanics|
|Technical Paper||Three Applications of a Mathematical Model, PRAKIMOD, in Frontal Collisions|
|Technical Paper||The Application of Anthropometry to Automotive Design|
CitationPetitjean, A. and Trosseille, X., "Statistical Simulations to Evaluate the Methods of the Construction of Injury Risk Curves," SAE Technical Paper 2011-22-0015, 2011, https://doi.org/10.4271/2011-22-0015.
- Di Domenico, L., Nusholtz, G. (2003) Comparison of Parametric and Non-Parametric Methods for Determining Injury Risk. SAE Paper N°2003-01-1362, Society of Automotive Engineers, Warrendale, PA.
- Goldstein, S., Frankenburg, E., and Kuhn, J. (1993) Biomechanics of bone. In Accidental Injury: Biomechanics and Prevention, ed. Nahum, A.M. and Melvin, J.W., pp. 198-223. Springer-Verlag, New York.
- Di Domenico, L., Nusholtz, G. (2005) Risk Curve Boundaries. Traffic Injury Prevention, 6:86-94.
- Eppinger, R., Sun, E., Bandak, F., Haffner, M.,Khaewpong, N., Maltese, M., Kuppa, S., Nguyen, T., Takhounts, E., Tannous, R., Zhang, A., Saul, R., (1999) Development of Improved Injury Criteria for the Assessment of Advanced Automotive Restraint Systems - II, National Highway Traffic Safety administration, U.S. Department of Transportation, Docket No. NHTSA-99-6407, November 1999.
- Funk, J.R., Srinivasan, S.C.M., Crandall, J.R., Khaewpong, N., Eppinger, R.H., Jaffredo, A.S., Potier, P. (2002) The Effects of Axial Preload and Dorsiflexion on the Tolerance of the Ankle/Subtalar Joint to Dynamic Inversion and Eversion. Stapp Car Crash Journal 46: 245-265.
- Hosmer, D.W., Lemeshow, S. (2000) Applied Logistic Regression. Wiley, New York, 2nd ed. 2000.
- Kennedy, E.A., Hurst, W.J., Stitzel, J.D., Cormier, J.M., Hansen, G.A., Smith, E.P., Duma, S.M. (2004) Lateral and Posterior Dynamic Bending of the Mid-Shaft Femur: Fracture Risk Curves for the Adult Population. Stapp Car Crash Journal 48: 27-51.
- Kent, R.W., Funk, J.R. (2004) Data censoring and parametric distribution assignment in the development of injury risk functions from biomechanical data. SAE Paper N°2004-01-0317, Society of Automotive Engineers, Warrendale, PA.
- Kuppa, S., Eppinger, R., McKoy, F., Nguyen, T. (2003) Development of side impact thoracic injury criteria and their application to the modified ES-2 dummy with rib extensions (ES-2re). Stapp Car Crash Journal 47: 189-210.
- Laituri, T.R., Prasad, P., Kachnowski, B.P., Sullivan, K., Przybylo, P.A. (2003) Predictions of AIS3+ Thoracic Risks for Belted Occupants in Full-Engagement, Real-World Frontal Impacts: Sensitivity to Various Theoretical Risk Curves. SAE Paper N°2003-01-1355, Society of Automotive Engineers, Warrendale, PA.
- Mertz, H.J., Weber, D. A. (1982) Interpretations of the Impact Responses of a 3-year-old Child Dummy Relative to Child Injury Potential, SAE 826048, SP-736, Automatic Occupant Protection Systems, Feb 1988.
- Mertz, H. J., Prasad, P. and Nusholtz, G. (1996) Head Injury Risk Assessment for Forehead Impacts. SAE Paper N°960099, Society of Automotive Engineers, Warrendale, PA.
- Mertz, H.J., Prasad, P., Irwin, A.L. (1997) Injury risk curves for children and adults in frontal and rear collisions. Proc.41st Stapp Car Crash Conference, pp. 13-30. Society of Automotive Engineers, Warrendale, PA.
- Mertz, H.J., Prasad, P. (2000) Improved injury risk curves for tension and extension moment measurements of crash dummies. Proc.44th Stapp Car Crash Conference, pp. 59-75. Society of Automotive Engineers, Warrendale, PA.
- Mertz, H.J., Irwin, A.L., Prasad, P. (2003) Biomechanical and Scaling Bases for Frontal and Side Impact Injury Assessment Reference Values. Stapp Car Crash Journal 47:155-188.
- Nakahira, Y., Furrkawa, K., Niimi, H., Ishihara, T., Miki, K. (2000) A Combined Evaluation Method and A Modified Maximum Likelihood Method for Injury Risk Curves. Proceedings IRCOBI, International Conference of the Biomechanics of Impact, Montpellier, France, pp. 147-155.
- Nusholtz, G., Mosier, R. (1999) Consistent Threshold Estimate for Doubly Censored Biomechanical Data. SAE Paper N°1999-01-0714, Society of Automotive Engineers, Warrendale, PA.
- Petitjean, A., Trosseille, X., Petit, P., Irwin, A., Hassan, J., Praxl, N. (2009) Injury Risk Curves for the WorldSID 50th Male dummy. Stapp Car Crash Journal 53: 443-476.
- Therneau, T. (2011), R software, Package ‘survival’.
- The R development Core Team (2009) R: a language and Environment for Statistical Computing), Reference Index.
- Viano, D. (1989) Biomechanical Responses and Injuries in Blunt Lateral Impact. Proc.33rd Stapp Car Crash Conference, pp. 113-142. Society of Automotive Engineers, Warrendale, PA.
- Wang, L., Banglmaier, R., Prasad, P. (2003) Injury Risk Assessment of Several Crash Data Sets. SAE Paper N°2003-01-1214, Society of Automotive Engineers, Warrendale, PA.
- Zhou, M., Lee, L., Chen, K. (2010), R software, Package ‘dblcens’.