This study provides proof of concept that imminent risk of fracture can be assessed by evaluating recent fracture, age, sex, race, medically significant falls, and psychoactive medications.
Objectives: To develop and validate predictive models for imminent fracture risk in a Medicare population.
Study Design: This retrospective administrative claims (Humana Research Database) study assessed imminent risk in Humana’s Medicare Advantage and Prescription Drug plan members.
Methods: Individuals (aged 67-87 years on January 1, 2015 [index]) with 1 year or more of history were followed for 3 months to up to 2 years, with censoring at death/disenrollment. The cohort was split into training and validation samples (1:1). Cox regression models assessed demographics, fracture history, medically significant falls, osteoporosis-related factors, frailty markers, and selected medications and comorbidities for independent predictors (P <.001) of incident nontraumatic clinical fractures in 12 and 24 months. A 6-variable model of 12-month risk used a published method for the risk-scoring point system.
Results: Of 1,287,354 individuals (mean age, 74.3 years; 56% female; 84% white), 3.8% had at least 1 fragility fracture at 12-month follow-up; 6.6% experienced fracture at 24 months (women vs men: 12 months, 4.8% vs 2.5%; 24 months, 8.3% vs 4.4%; both P <.01). At 12 months, recent fracture conferred approximately 3-fold-higher fracture risk (vs no recent fracture). Older age, white race, female sex, osteoporosis-related screening/diagnosis/medication, antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications, history of falls, fracture history, and respiratory conditions also increased risk (all P <.0001). The simplified model (recent fracture, age, sex, race, falls, antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications) performed well (C statistic = 0.71).
Conclusions: Recent fracture, older age, female sex, white race, falls, and antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications predict imminent fracture risk in an older-adult Medicare Advantage population. Imminent fracture risk can be assessed using 6 easily quantified factors.
Am J Manag Care. 2020;26(3):e91-e97. https://doi.org/10.37765/ajmc.2020.42641
Up to 25% of men and 50% of women will experience a fracture in their remaining lifetime after age 50 years.1 Osteoporosis-related fractures confer significant direct medical costs, quality-of-life decrements and functional disability, and indirect costs (lost productivity and caregiver burden).2-7 By 2025, fracture-related direct treatment costs are expected to reach $25.3 billion in the United States.2,6 The impact of fractures to individuals may be long-lasting and result in increased use of healthcare services, a greater likelihood of hospitalization and skilled nursing care, reduced economic circumstances, and increased mortality risk for up to 10 years post fracture.3,4,7,8
Various antiresorptive and bone-forming therapies have received FDA approval based on their ability to reduce the risk of fracture in men and women with osteoporosis.9 Fracture risk-assessment tools help clinicians estimate the risk of osteoporosis-related fracture in their patients, and these estimates inform treatment recommendations. The most commonly used tools, such as FRAX,10 focus exclusively on long-term fracture risk (ie, 10-year time horizon). The literature suggests that the imminent risk (ie, in next 1-2 years) of fracture also has important implications for treatment decisions.10-15 Imminent risk is responsive to changes in bone mineral density (BMD), fracture, and fall history, as well as factors that influence physical or cognitive functioning in ways that are not captured in existing tools.10,13-17 Furthermore, patients with osteoporosis may not perceive fracture risk accurately, and long-term risk assessments, in conjunction with patients’ potentially incomplete understanding of the relationship between osteoporosis and fracture risk, may not provide sufficient motivation for treatment initiation.18-20
Elevated imminent risk may require a different treatment approach from fracture risk that accumulates or extends over a longer time frame.10,14 The objective of this study was to develop and validate predictive models for imminent risk of fracture in a population of older adults enrolled in Humana’s Medicare Advantage and Prescription Drug (MAPD) plans. This paper describes the primary full model, which includes an extensive set of predictors for the outcome of any fracture within a 12-month follow-up, and a parsimonious model, which was developed to provide proof of concept for a risk-assessment tool that could easily be used for hand-calculating risk during a routine physician office visit.
Data Source and Population
This retrospective database study used a cohort design and multivariate analysis to identify predictors of imminent fracture risk in older adults enrolled in Humana MAPD plans. The Humana Research Database (Humana; Louisville, Kentucky), which contains administrative data for Humana’s fully insured commercial and Medicare individuals, was used. Data on MAPD enrollment, medical claims, and outpatient pharmacy claims were linked for each individual. The study period—January 1, 2014, through December 31, 2016—was split into a 1-year baseline period (January 1, 2014, through December 31, 2014) and a 2-year follow-up period (January 1, 2015, through December 31, 2016).
All study individuals were aged 67 to 87 years on the index date (January 1, 2015) and had continuous enrollment in an MAPD plan throughout the baseline period and for at least the first 3 months of the follow-up period. Follow-up was censored at the date of death or plan disenrollment for individuals lost to follow-up before the end of study. Individuals with an International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) or International Classification of Diseases, Tenth Revision, Clinical Modification (ICD-10-CM) diagnosis of either nonmelanomatous cancer or Paget disease of the bone and individuals enrolled in administrative services—only health plans or in a Humana health plan which is contractually excluded from research studies were excluded.
The sample was split into a randomly generated training (development) sample and a validation sample, with a 1:1 ratio. The primary study outcome was defined as incident fragility fracture. These fractures were identified using a combination of ICD-9-CM and ICD-10-CM diagnosis codes. Clinical fragility (nontraumatic) fractures with an osteoporosis likelihood score of 8 or greater, as defined by Warriner et al,21 were identified during 12- and 24-month follow-up intervals. This published algorithm was developed to enhance the accuracy of identifying osteoporosis-related fractures in administrative claims data and was used in the current study to reduce the potential for misclassification of the primary study outcome. Fracture codes recorded during the follow-up period but within 60 days of a baseline fracture were considered as a continuation of care for the baseline fracture and, therefore, not considered in the analysis of fracture outcomes during follow-up.
Covariates considered as potential independent variables for the predictive models included patient demographics (age, sex, race, and geographic region), history of fracture, falls that result in medical encounters or that require treatment (ie, medically significant falls), BMD testing, comorbidities, medication use, and markers of frailty. Comorbidities were quantified by the Deyo-Charlson Comorbidity Index22 score and yes/no indicator variables for selected comorbidities, including those associated with an increased likelihood of fractures/medically significant falls. Osteoporosis treatments, as well as medications associated with an increased likelihood of fractures/medically significant falls (eg, antihypertensives, anticonvulsants, antidepressants, antipsychotics), were assessed. Markers of frailty included data on difficulty walking, use of durable medical equipment or ambulance/life support, paralysis, weakness, and podiatric care.
The training sample was used to construct risk-prediction models for any fragility fracture occurring in the 12- and 24-month follow-up intervals. The primary analysis focused on incident fractures that occurred in the 12-month follow-up, and data from the 24-month follow-up provided a sensitivity analysis. Here we present predictive models for any fragility fracture in these 2 follow-up periods, as well as results from a model specific to hip fracture risk during a 12-month follow-up.
To accommodate variable length of follow-up for individuals in the training and validation samples, Cox regression models23 were used to estimate adjusted hazard ratios (HRs) for potential predictor variables. Variables significantly associated with fracture risk (P <.001), after controlling for all other predictors in the model, were considered to be independent predictors. Using the validation sample, the predictor coefficients from the finalized predictive models of fracture risk for hip fracture were applied to the patient baseline characteristics to obtain the likelihood of fracture risk for each patient. C statistics and other measures (false-positive rate, false-negative rate, sensitivity, specificity, positive predictive value, negative predictive value, and efficiency [ie, percentage of accurately predicted patients in the study sample]) were computed to assess model performance when data from the validation sample were used. Models are typically considered effective when the C statistic is greater than 0.7 and strong when it exceeds 0.8.24,25 Therefore, a minimum threshold of 0.7 for the C statistic was considered when comparing the models. Additional performance metrics were assessed at fracture probabilities of 5% or greater, 10% or greater, 15% or greater, 20% or greater, and 25% or greater.
A parsimonious version of the primary predictive model was developed to adapt the results of the analysis for use as an imminent (12-month) fracture risk assessment for an individual patient in the routine care setting. The goal was to create a parsimonious model with approximately 5 variables and to incorporate a risk-scoring point system based on the Sullivan methodology26 adapted for use with 1-year incidence estimates. In the parsimonious model, increasing numbers of points were assigned to increasing levels of risk (eg, increasing age was associated with a greater number of points), and the totaled points corresponded to an individual’s predicted fracture probability during 12 months of follow-up.
From Humana’s MAPD population, 1,287,354 individuals (mean age, 74.3 years; 56% female; 84% white race) met all criteria for inclusion in this study (eAppendix Figure 1 [eAppendix available at ajmc.com]; Table 1). Only 6% of these individuals had a history of falls during the baseline period. Approximately 4% of the individuals identified had a claim for a recent fracture (ie, within the baseline year) and fewer than 1% had a historical fracture at any time in recorded medical history. Although 19% of individuals had a diagnosis of osteoporosis, prescription fills for osteoporosis medications during the baseline period were relatively uncommon (6%). By contrast, 72% of individuals had been given a diagnosis of hypertension, and 64% of the total population used antihypertensives. Approximately 15% to 20% of individuals had used other medications of interest.
In the 12-month follow-up period, 3.8% of individuals experienced at least 1 fragility fracture, and the fracture incidence increased to 6.6% over the 24-month follow-up period (Table 2). Fracture incidence was higher among women than men over both follow-up intervals (12 months, 4.8% vs 2.5%; 24 months, 8.3% vs 4.4%; P <.01 for both). The incidence of hip fractures in men and women combined was 0.7% in the 12-month follow-up period and 1.4% in the 24-month follow-up period.
In the primary analysis, recent fracture (ie, within the past year prior to the start of the follow-up period) had the strongest association with 1-year fracture risk (Table 3); this factor conferred nearly a 3-fold-higher risk among the affected individuals compared with the fracture risk in individuals who had no recent fracture. Risk also increased steadily with age and was highest among individuals aged 85 to 87 years. White race and female sex were also strongly associated with increased fracture risk, as were osteoporosis-related screening and diagnosis, the use of osteoporosis-related medications, use of antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications, falls, history of fracture, and respiratory conditions (all P <.0001). The primary model performed well, with a C statistic of 0.73, false-positive rates ranging from 74% to 90%, and false-negative rates ranging from 2% to 4% across fracture probability cut-offs ranging from 5% or greater to 25% or greater (eAppendix Table 1 [A]; eAppendix Figure 2 [A]). Across these same fracture probability categories, the model’s specificity was better than its sensitivity. Higher areas under the curve and C statistics were observed with the 12-month model compared with the 24-month model developed as a sensitivity analysis (eAppendix Table 1 [B]; eAppendix Figure 2 [B]). Risk factors for hip fracture within 12 months exhibited significant overlap with the predictors identified for all fragility fractures during that time interval (Table 3). In particular, older age and recent fracture exerted the greatest influence on risk, with 2- to 3-fold increases in risk associated with these factors. The C statistic for the 12-month hip fracture model was 0.77.
The parsimonious model for any fracture, which included recent fracture, age, sex, race, falls, and antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications as predictors, exhibited good performance with a C statistic of 0.71 (Table 4). The scored points and corresponding predicted probability of fracture were compared with actual fracture rates in 2015 for individuals with each point total (Table 5). These rates increased as the number of points increased.
This study identified predictors of the imminent risk of fracture (ie, within 12-24 months) among older adults enrolled in Humana’s MAPD plans. Of the 75 patient characteristics initially considered for the predictive models, recent fracture (ie, fracture occurring during the 1-year baseline period), advanced age, and white race were most strongly associated with increased imminent risk of fractures. Other key factors included female sex, indicators of osteoporosis (screening, diagnoses, and medications), any history of fracture, falls, and antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications.
The risk factors identified in this study are consistent with those identified in prior observational studies that assessed imminent fracture risk in different populations, including individuals with and without osteoporosis. In osteoporosis, generally, white/Caucasian women are at greater risk of fracture than nonwhite women or men of similar age.27 Low BMD, prior fracture, and falls are also well-documented risk factors for fragility fracture.28 The growing body of research on risk factors for imminent fractures indicates that these factors, along with certain comorbidities and medications that affect physical and cognitive functioning, also influence the risk of fracture in the near term.10-16 The model results for hip fracture illustrate this last point well, with previous research indicating that various mechanisms may influence the association between dementia and hip fractures, including common risk factors for both and the role of dementia and its treatment as an indirect risk factor that increases fracture risk through intermediate risk factors (eg, falls, vitamin D deficiency, treatment adverse effects).29
We found that antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications had a significant impact on the imminent risk of fracture. It is possible, as recent literature suggests, that these medications increase the risk of falls, thereby indirectly increasing the risk for fragility fractures.30-32 Although the 2 variables—antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications and falls—were not highly correlated (Pearson correlation coefficient = 0.13), this may be because falls are vastly underreported in medical claims data, which would make it difficult to show this association in the available data.
The accumulating literature on short-term or imminent risk highlights the fact that fracture risk is not linear but varies over time as a patient’s characteristics and exposures change, and investigators have begun identifying factors that contribute specifically to increased risk within the next 1 to 2 years.10,14,33-35 An important aspect of several risk factors identified in our full models of any fracture and in previous studies of imminent risk is that they are dynamic and their influence on fracture risk occurs quickly. Prior fracture illustrates this concept well. There is strong evidence that the risk of fragility fracture is approximately doubled after an initial fracture and that the greatest increase in risk occurs within the first year post fracture.17,28,36 Our results confirm this pattern of risk by demonstrating that although any history of fracture is a significant predictor of imminent fracture risk (adjusted HR, 1.26; 95% CI, 1.17-1.36), the impact of recent fracture is more than twice as great (adjusted HR, 2.81; 95% CI, 2.70-2.92).
Fracture risk may be particularly variable in older individuals because events that we have identified as predictors of imminent risk, such as the occurrence of a stroke or diagnosis of a neurologic or psychologic disorder (eg, Alzheimer disease, depression, anxiety), generally become more common with age. Similarly, changes in physical or cognitive functioning, which may be related to the onset of new comorbidities or the initiation of new medications, may either increase fracture risk directly or act on fracture risk by increasing a patient’s risk of falls as mentioned previously.10,14,30-32 This inherent malleability of fracture risk underscores the importance of being able to quantify imminent risk as a patient’s medical status changes and to help patients understand the need for timely intervention to reduce the risk of fracture and its related sequelae.
This dynamic nature of fracture risk suggests a need for tools that can be used to regularly assess the changing risk of imminent fracture in older adults. In this population, such risk assessments should be done regularly, and more frequent use of a tool such as the one developed in this study might be useful for patients whose health status is in transition or particularly modifiable. Quantifying imminent risk in an ongoing, real-time fashion requires a simple risk-assessment tool. With this in mind, we developed a parsimonious model based on our full model and adapted this model into a prototype for a points-based risk-assessment tool. This model uses only 6 variables (recent fracture, age, sex, race, falls, and antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications) and was adapted to a point-based system, with total points (risk score) easily computed during a routine physician office visit. Given the importance of falls as a predictor of long-term fracture risk and their association with imminent risk of fracture, additional parsimonious models were tested to assess a composite physical factors variable as a potential predictor of imminent risk. Because physical factors may influence fall risk, in one model, the physical factor variable was simply substituted for falls in the existing parsimonious model, whereas in a second model, falls and physical factors were combined and tested with age, sex, race, and recent fracture. The model performance was similar in both of these versions and resulted in the same number of points being attributed to this fifth variable.
This work supports the feasibility of a formal risk assessment that can quickly and easily be used during physician office visits to identify patients who are at high risk for an imminent fracture. To guide the use of this tool, input could be obtained from physician focus groups to help determine clinically meaningful point score thresholds for fracture prevention initiatives. Such initiatives would aid in the identification of patients in the greatest need of timely osteoporosis management ranging from BMD screening to interventions including prescription medications and initiatives to reduce the risk of falls.
This study is subject to limitations common to administrative claims data research. Identifying fragility fractures and distinguishing incident fractures from prevalent fractures can be challenging with the information available in claims. To address this issue, we excluded all fractures that were recorded as trauma related and examined the timing of fracture-related care to identify and exclude claims most likely to represent care for prevalent fractures that occurred in the baseline period. We also limited the analysis to closed fractures with the highest likelihood of being classified as fragility fractures based on a previously published methodology designed for use with administrative claims data.21 Specifically, we included only fractures with an osteoporosis likelihood score of 8 or greater to reduce bias related to potential misclassification of fracture outcomes. There was also potential for selection bias, given the requirement for at least 12 months of continuous insurance coverage in the baseline period. Although the multivariate regression modeling used in this study would have reduced the potential for bias among measured covariates, bias associated with unmeasured covariates would have remained. We also note that results may not generalize to the full US population, as the geographic mix in the sample was not balanced across all regions.
In summary, recent fracture, older age, female sex, white race, falls, and use of antidepressant/antipsychotic/sedative hypnotic/muscle relaxant medications were strong predictors of the imminent risk of fracture in Medicare-enrolled older adults. A patient’s level of imminent risk may change over time, and this work provides proof of concept that the imminent risk of fracture can be assessed by evaluating a small number of factors. Administering such a risk-assessment tool regularly may improve physicians’ opportunities to identify elderly patients at high risk of imminent fractures and enhance their ability to provide timely interventions (screening, medications, and fall-prevention strategies). Additional research is needed to set clinically meaningful risk-score thresholds for intervention, determine the applicability and accuracy of this risk-score model to specific patient subgroups, and evaluate the best assessment schedule and optimal ways to implement this tool within the existing physician workflow.
The authors would like to thank Sally Wade (Wade Outcomes Research and Consulting) for medical writing assistance.Author Affiliations: Humana Healthcare Research, Inc, Humana Inc (RLS, LS, MKP), Louisville, KY; Amgen Inc (RLB), Thousand Oaks, CA.
Source of Funding: Amgen Inc (Thousand Oaks, CA).
Author Disclosures: Mr Sheer is an employee of and stockholder in Humana Inc. Mr Barron was employed at Amgen until May 2018 and owned stock in Amgen as an employee. Ms Sudharshan was employed at Humana at the time of this study; Humana received consulting fees from Amgen for the conduct of this study. Dr Pasquale owns stock in Amgen and is an employee of Humana Health Research, which received funding from Amgen to conduct the study.
Authorship Information: Concept and design (RLS, RLB, LS, MKP); analysis and interpretation of data (RLS, RLB, LS, MKP); drafting of the manuscript (RLS, MKP); critical revision of the manuscript for important intellectual content (RLS, RLB, LS, MKP); statistical analysis (RLS, LS); provision of patients or study materials (LS); obtaining funding (RLB); administrative, technical, or logistic support (RLB, LS); and supervision (RLB, MKP).
Address Correspondence to: Margaret K. Pasquale, PhD, Humana Healthcare Research, Inc, Humana Inc, 500 W Main St, Louisville, KY 40202. Email: firstname.lastname@example.org.REFERENCES
1. What is osteoporosis and what causes it? National Osteoporosis Foundation website. nof.org/patients/what-is-osteoporosis. Accessed April 3, 2018.
2. Burge R, Dawson-Hughes B, Solomon DH, Wong JB, King A, Tosteson A. Incidence and economic burden of osteoporosis-related fractures in the United States, 2005-2025. J Bone Miner Res. 2007;22(3):465-475. doi: 10.1359/jbmr.061113.
3. Dyer SM, Crotty M, Fairhall N, et al; Fragility Fracture Network (FFN) Rehabilitation Research Special Interest Group. A critical review of the long-term disability outcomes following hip fracture. BMC Geriatr. 2016;16:158. doi: 10.1186/s12877-016-0332-0.
4. Fischer S, Kapinos KA, Mulcahy A, Pinto L, Hayden O, Barron R. Estimating the long-term functional burden of osteoporosis-related fractures. Osteoporos Int. 2017;28(10):2843-2851. doi: 10.1007/s00198-017-4110-4.
5. Nazrun AS, Tzar MN, Mokhtar SA, Mohamed IN. A systematic review of the outcomes of osteoporotic fracture patients after hospital discharge: morbidity, subsequent fractures, and mortality. Ther Clin Risk Manag. 2014;10:937-948. doi: 10.2147/TCRM.S72456.
6. Ström O, Borgström F, Kanis JA, et al. Osteoporosis: burden, health care provision and opportunities in the EU: a report prepared in collaboration with the International Osteoporosis Foundation (IOF) and the European Federation of Pharmaceutical Industry Associations (EFPIA). Arch Osteoporos. 2011;6:59-155. doi: 10.1007/s11657-011-0060-1.
7. Tajeu GS, Delzell E, Smith W, et al. Death, debility, and destitution following hip fracture. J Gerontol A Biol Sci Med Sci. 2014;69(3):346-353. doi: 10.1093/gerona/glt105.
8. Haentjens P, Magaziner J, Colón-Emeric CS, et al. Meta-analysis: excess mortality after hip fracture among older women and men. Ann Intern Med. 2010;152(6):380-390. doi: 10.7326/0003-4819-152-6-201003160-00008.
9. Medication and treatment adherence. National Osteoporosis Foundation website. nof.org/patients/treatment/medicationadherence. Accessed April 3, 2018.
10. Bonafede M, Shi N, Barron R, Li X, Crittenden DB, Chandler D. Predicting imminent risk for fracture in patients aged 50 or older with osteoporosis using US claims data. Arch Osteoporos. 2016;11(1):26. doi: 10.1007/s11657-016-0280-5.
11. Adams A, Johnson E, Zhou H, et al. Prediction of two-year risk of fracture among older US women. J Bone Miner Res. 2016;31(suppl 1):S94. Abstract FR0237. doi: 10.1002/jbmr.3107.
12. Balasubramanian A, Zhang J, Chen L, et al. High risk of second fracture within 1, 2, 5 years after prior fracture among women 65 years or older. J Bone Miner Res. 2016;31(suppl 1):S93. Abstract FR0233. doi: 10.1002/jbmr.3107.
13. Hannan M, Weycker D, McLean R, et al. Predictors of imminent risk of non-vertebral fracture in older women: the Framingham Osteoporosis Study. J Bone Miner Res. 2016;31(suppl 1):S365. Abstract MO0232. doi: 10.1002/jbmr.3107.
14. Weycker D, Edelsberg J, Barron R, et al. Predictors of near-term fracture in osteoporotic women aged ≥65 years, based on data from the study of osteoporotic fractures. Osteoporos Int. 2017;28(9):2565-2571. doi: 10.1007/s00198-017-4103-3.
15. Yusuf A, Hu Y, Chandler D, Crittenden B, Barron R. Characteristics of patients at high one-year fracture risk. J Bone Miner Res. 2016;31(suppl 1):S363. Abstract MO0223. doi: 10.1002/jbmr.3107.
16. Huntjens KM, Kosar S, van Geel TA, et al. Risk of subsequent fracture and mortality within 5 years after a non-vertebral fracture. Osteoporos Int. 2010;21(12):2075-2082. doi: 10.1007/s00198-010-1178-5.
17. Johnell O, Kanis JA, Odén A, et al. Fracture risk following an osteoporotic fracture. Osteoporos Int. 2004;15(3):175-179. doi: 10.1007/s00198-003-1514-0.
18. Deloumeau A, Moltó A, Roux C, Briot K. Determinants of short term fracture risk in patients with a recent history of low-trauma non-vertebral fracture. Bone. 2017;105:287-291. doi: 10.1016/j.bone.2017.08.018.
19. Sale JE, Gignac MA, Hawker G, et al. Patients do not have a consistent understanding of high risk for future fracture: a qualitative study of patients from a post-fracture secondary prevention program. Osteoporos Int. 2016;27(1):65-73. doi: 10.1007/s00198-015-3214-y.
20. Alami S, Hervouet L, Poiraudeau S, Briot K, Roux C. Barriers to effective postmenopausal osteoporosis treatment: a qualitative study of patients’ and practitioners’ views. PLoS One. 2016;11(6):e0158365. doi: 10.1371/journal.pone.0158365.
21. Warriner AH, Patkar NM, Curtis JR, et al. Which fractures are most attributable to osteoporosis? J Clin Epidemiol. 2011;64(1):46-53. doi: 10.1016/j.jclinepi.2010.07.007.
22. Deyo RA, Cherkin DC, Ciol MA. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J Clin Epidemiol. 1992;45(6):613-619. doi: 10.1016/0895-4356(92)90133-8.
23. Dunkler D, Plischke M, Leffondré K, Heinze G. Augmented backward elimination: a pragmatic and purposeful way to develop statistical models. PLoS One. 2014;9(11):e113677. doi: 10.1371/journal.pone.0113677.
24. Hosmer DW, Lemeshow S. Applied Logistic Regression. New York, NY: John Wiley & Sons; 1989.
25. Hosmer DW, Lemeshow S. Applied Logistic Regression. 2nd ed. New York, NY: John Wiley & Sons; 2000.
26. Sullivan LM, Massaro JM, D’Agostino RB Sr. Presentation of multivariate data for clinical use: the Framingham Study risk score functions. Stat Med. 2004;23(10):1631-1660. doi: 10.1002/sim.1742.
27. Taylor AJ, Gary LC, Arora T, et al. Clinical and demographic factors associated with fractures among older Americans [erratum in Osteoporos Int. 2011;22(4):1275-1276. doi: 10.1007/s00198-010-1347-6]. Osteoporos Int. 2011;22(4):1263-1274. doi: 10.1007/s00198-010-1300-8.
28. Klotzbuecher CM, Ross PD, Landsman PB, Abbott TA III, Berger M. Patients with prior fractures have an increased risk of future fractures: a summary of the literature and statistical synthesis. J Bone Miner Res. 2000;15(4):721-739. doi: 10.1359/jbmr.2000.15.4.721.
29. Friedman SM, Menzies IB, Bukata SV, Mendelson DA, Kates SL. Dementia and hip fractures: development of a pathogenic framework for understanding and studying risk. Geriatr Orthop Surg Rehabil. 2010;1(2):52-62. doi: 10.1177/2151458510389463.
30. Beunza-Sola M, Hidalgo-Ovejero ÁM, Martí-Ayerdi J, Sánchez-Hernández JG, Menéndez-García M, García-Mata S. Study of fall risk-increasing drugs in elderly patients before and after a bone fracture. Postgrad Med J. 2018;94(1108):76-80. doi: 10.1136/postgradmedj-2017-135129.
31. Coupland C, Hill T, Morriss R, Moore M, Arthur A, Hippisley-Cox J. Antidepressant use and risk of adverse outcomes in people aged 20-64 years: cohort study using a primary care database. BMC Med. 2018;16(1):36. doi: 10.1186/s12916-018-1022-x.
32. Thorell K, Ranstad K, Midlöv P, Borgquist L, Halling A. Is use of fall risk-increasing drugs in an elderly population associated with an increased risk of hip fracture, after adjustment for multimorbidity level: a cohort study. BMC Geriatr. 2014;14:131. doi: 10.1186/1471-2318-14-131.
33. Balasubramanian A, Zhang J, Chen L, et al. Risk of subsequent fracture after prior fracture among older women. Osteoporos Int. 2019;30(1):79-92. doi: 10.1007/s00198-018-4732-1.
34. Masud T, Binkley N, Boonen S, Hannan MT; FRAX Position Development Conference Members. Official positions for FRAX clinical regarding falls and frailty: can falls and frailty be used in FRAX? from Joint Official Positions Development Conference of the International Society for Clinical Densitometry and International Osteoporosis Foundation on FRAX. J Clin Densitom. 2011;14(3):194-204. doi: 10.1016/j.jocd.2011.05.010.
35. Roux C, Briot K. Imminent fracture risk. Osteoporos Int. 2017;28(6):1765-1769. doi: 10.1007/s00198-017-3976-5.
36. Johnell O, Oden A, Caulin F, Kanis JA. Acute and long-term increase in fracture risk after hospitalization for vertebral fracture. Osteoporos Int. 2001;12(3):207-214. doi: 10.1007/s001980170131.