Fracture risk assessment in the presence of competing risk of death

Purpose To identify the optimal statistical approach for predicting the risk of fragility fractures in the presence of competing event of death. Methods We used real-world data from the Dubbo Osteoporosis Epidemiology Study that has monitored 3035 elderly participants for bone health and mortality. Fragility fractures were ascertained radiologically. Mortality was confirmed by the State Registry. We considered four statistical models for predicting fracture risk: (i) conventional Cox’s proportional hazard model, (ii) cause-specific model, (iii) Fine-Gray sub-distribution model, and (iv) multistate model. These models were fitted and validated in the development (60% of the original sample) and validation (40%) subsets, respectively. The model performance was assessed by discrimination and calibration analyses. Results During a median follow-up of 11.3 years (IQR: 7.2, 16.2), 628 individuals (34.5%) in the development cohort fractured, and 630 (34.6%) died without a fracture. Neither the discrimination nor the 5-year prediction performance was significantly different among the models, though the conventional model tended to overestimate fracture risk (calibration-in-the-large index = − 0.24; 95% CI: − 0.43, − 0.06). For 10-year risk prediction, the multistate model (calibration-in-the-large index = − 0.05; 95% CI: − 0.20, 0.10) outperformed the cause-specific (− 0.23; − 0.30, − 0.08), Fine-Gray (− 0.31; − 0.46, − 0.16), and conventional model (− 0.54; − 0.70, − 0.39) which significantly overestimated fracture risk. Conclusion Adjustment for competing risk of death has minimum impact on the short-term prediction of fracture. However, the multistate model yields the most accurate prediction of long-term fracture risk and should be considered for predictive research in the elderly, who are also at high mortality risk. Summary Fracture risk assessment might be compromised by the competing event of death. This study, using real-world data found a multistate model was superior to the current competing risk methods in fracture risk assessment. A multistate model is considered an optimal statistical method for predictive research in the elderly.