An inception cohort analysis to predict nonunion in tibia and 17 other fracture locations

The epidemiology of fracture nonunion has been characterized so it is potentially possible to predict nonunion using patient-related risk factors. However, prediction models are currently too cumbersome to be useful. We test a hypothesis that nonunion can be predicted with ≤10 variables, retaining the predictive accuracy of a full model with 42 variables.

We sought to predict nonunion with prospectively-acquired inception cohort data for 18 different bones, using the smallest possible number of variables that did not substantially decrease prediction accuracy. An American nationwide claims database of ∼90.1 million participants was used, which included medical and drug expenses for 2011-2012. Continuous enrollment was required for 12 months after fracture, to allow sufficient time to capture a nonunion diagnosis. Health claims were evaluated for 309,330 fractures. A training dataset used a random subset of 2/3 of these fractures, while the remaining fractures formed a validation dataset. Multivariate logistic regression and stepwise logistic regression were used to identify variables predictive of nonunion. P value and the Akaike Information Criterion (AIC) were used to select variables for reduced models. Area-under-the-curve (AUC) was calculated to characterize the success of prediction.

Nonunion rate in 18 fracture locations averaged 4.93%. Algorithms to predict nonunion in 18 locations in the full-model validation set had average AUC=0.680 (±0.034). In the reduced models, average validation set AUC=0.680 (±0.033) and the average number of risk factors required for prediction was 7.6. There was agreement across training set, validation set, and reduced set; in tibia, reduced model validation AUC=0.703, while the full-model validation AUC=0.709. Certain risk factors were important for predicting nonunion in ≥10 bones, including open fracture, multiple fracture, osteoarthritis, surgical treatment, and use of certain medications, including anticoagulants, anticonvulsants, or analgesics.

Nonunion can be predicted in 18 fracture locations using parsimonious models with <10 patient demography-related risk factors. The model reduction approach used results in simplified models that have nearly the same AUC as the full model. Reduced algorithms can predict nonunion because risk factors important in the full models remain important in the reduced models. This prognostic inception cohort study provides Level I evidence.