Modified FRAX Score for Prediction of Fall-induced Hip Fractures; The Added Value of Fall Energy, Number, and Social Vulnerability Index

Introduction

Hip fractures are a debilitating health concern worldwide that is estimated to rise from 1.26 million cases globally in 1990 to 4.5 million cases by 2050.^{1, 2} An estimated 37.3 million falls occur annually, contributing to over 95% of hip fractures.^1,3,4 According to the literature, a total of 1.75 million disability-adjusted life years (DALYs) are lost globally, accounting for 0.1% of the worldwide disease burden and 1.4% of the disease burden among women in established market economies.⁵ Therefore, predicting and implementing early intervention strategies is of crucial value to decrease the burden of hip fractures associated with falls.

To accomplish this goal, the Fracture Risk Assessment Tool (FRAX) generated by the World Health Organization (WHO) has been developed to predict the 10-year probability of an adult having a hip fracture.⁶ FRAX utilizes the most important risk factors known for major osteoporotic fractures, including age, sex, weight, height, alcohol use, glucocorticoid use, femoral neck bone mineral density (BMD), history of previous fracture in the patients and their parents, smoking status, and rheumatoid arthritis.⁷ With 6 million yearly calculations done in 173 countries, FRAX has been immensely helpful in predicting fracture risks. However, it is frequently criticized for lacking a detailed and standardized description of each risk factor, including the recency of fracture and glucocorticoid dose.⁸ Furthermore, since falls are one of the most critical factors leading to hip fractures in older adults, a common criticism of FRAX is that it does not factor in fall history.^{8, 9}

Since the official launch of FRAX in 2008, concerns have arisen regarding its inability to consider the impact of falls.^10,11 In 2010, the Task Force acknowledged this concern, stating that “Fracture probability may be underestimated in individuals with a history of frequent falls, but quantification of this risk is not currently possible.” The primary reason for this challenge was the lack of sufficient evidence available at that time. Since then, the body of literature on this topic has expanded, with controversial results suggesting a complex relationship between FRAX and falls and highlighting the need for a more careful approach when integrating the two.^12–14 Furthermore, while much of the research focuses on this aspect, other fall characteristics seem not well investigated. Among these, the mechanics of the falls seem to play a significant role in the risk of hip fractures.¹⁵ Hayes et al. pointed out that the energy or impact of the fall is found to be a crucial factor when assessing the risk factors for hip fractures.¹⁶ This suggests that not just the frequency of falls but also their nature and severity might be essential factors in assessing fall-induced hip fracture risk.

Given the multifactorial nature of hip fractures, efforts have been made to point out the relevant demographic and other non-medical risk factors in hip fractures.^17–20 Some factors, including age and sex, are already addressed in the FRAX, while others, like race, ethnicity, and socioeconomic factors, are not. The Social Vulnerability Index (SVI), developed by the Centers for Disease Control and Prevention (CDC), is intended to measure the vulnerability of communities to environmental hazards, disease outbreaks, and other public health emergencies.²¹ It encompasses various subscales, including Socioeconomic Status, Household Characteristics, Racial and Ethnic Minority Status, Housing Type and Transportation. When used in medicine, it enables capturing the pre-existing social inadequacies that make individuals or communities disproportionately vulnerable to a particular medical condition or disease.²² By understanding these correlations, healthcare providers and policymakers can develop more targeted interventions and resource allocations to reduce the risk of fall-induced injuries in the most vulnerable populations.

In this study, we hypothesized that fall frequency and energy, along with the SVI, are of predictive value in fall-induced hip fractures. Moreover, in order to have a more comprehensive and accurate assessment of hip fracture risk, we aim to develop a modified version of the FRAX that incorporates additional variables found to be significant risk factors for fall-induced hip fractures in our study.

Methods and Materials

I. Study Design and Population

We conducted a retrospective case-control study through a tertiary institution located in Boston, Massachusetts. After getting approval from the institutional review board (IRB # 2023P000741), data were queried through the institution’s data repository, the Research Patient Data Registry (RPDR). The relevant ICD-9, ICD-10, and CPT codes, along with string searches for clinician notes, were used for the query. The initial results were screened, and participants aged ≥ 40 years with a confirmed diagnosis of fall-induced hip fracture were included as cases and those who experienced falls without a resultant hip injury as controls. A total of 2,132 patients with a history of falls who were assessed for hip fractures were screened. Cases (patients with hip fracture) and controls (without hip fracture) were age-matched for the final analysis. The study excluded individuals who did not have a mention of the mechanism of falls on their charts or experienced falls resulting from violent encounters and animal attacks, significant external forces such as car or motor vehicle accidents, high-impact sports like skiing, and those with fractures caused by underlying pathological conditions to avoid the influence of confounding injuries.²³ We also excluded patients who did not have enough information on their charts for the calculation of the FRAX score or for the assessment of their falling injury (Figure 1).

• Download figure
• Open in new tab

Figure 1.

Receiver Operator Characteristics (ROC) Curve for the Original FRAX Model

• Download figure
• Open in new tab

Figure 2

Receiver Operator Characteristics (ROC) Curve for the Best Developed Model. Including the FRAX, Race, Fall Frequency, and Housing Domain of the SVI

Results

A total of 261 patients per group were included, and their demographic data are depicted in Table 1. Race, BMI, FRAX, Racial and Ethnic Minority Status, Housing type and Transportation showed significant differences between the hip fracture and no fracture groups.

The outcome of the logistic regression test on various factors that showed significant differences between the fracture and no-fracture groups, as well as fall characteristics, including energy and number of falls, are shown in Table 2. The white race also showed a higher rate of fall-induced fracture (OR=2.85; CI:1.56-5.20), although there was a discrepancy between the number of whites and non-whites (86% vs 13%, respectively). Housing Type and Transportation were shown to decrease the chance of experiencing the injury significantly (OR=0.30, CI:0.12-0.80).

The FRAX tool and the prediction model developed in this study (modified FRAX) were evaluated using the data of the cases and the control group. The performances of the FRAX tool, as well as the modified FRAX, are shown in Table 3. The final prediction model was developed using the coefficients of the best model and deployed on the following web address https://fixus.mgh.harvard.edu/fall/

Discussions

Almost half of the elderly population experience a fall, and up to 10% have reported a skeletal injury.³ A strong association between the propensity to fall and future fracture is mentioned in various studies. These studies showed a greater risk in individuals with a history of multiple falls.^{26, 27} However, the history and frequency of falls have not been included in current prediction models such as the FRAX tool. Falls and hip fractures share many similar risk factors, including increasing age, smoking, alcohol consumption, and frailty; thus, hypothesizing that FRAX and any other prediction model that could predict fractures could also be modified to predict the risk of fall-induced fractures would be reasonable.^12,28 Including patient-specific characteristics such as race, country of origin, medications, and medical conditions affecting bone quality was considered in several studies to improve the accuracy, validity, and reliability of the FRAX tool¹² In this study, we investigated the added predictive value of fall characteristics, basic demographics, SVI, and its domains to the FRAX score for the prediction of fall-induced hip fractures. Our results showed that the number of falls, race, Housing type and Transportation subscale of the SVI play a considerable role in predicting the risk of fall-induced hip fractures.

In the nuanced realm of assessing fracture risk, the relationship between fall frequency and the FRAX score emerges as a complex and evolving topic.¹² In 2015, Harvey et al. conducted a study in a 1,836 male population comparing the predictive value of FRAX and fall frequency for future falls.²⁷ They suggested a surprising twist by concluding that the FRAX score is capable of predicting future falls (Hazard Ratio=1.64; CI:1.36-1.97). This finding implied that since falls and fractures share a lot of risk factors, FRAX has an inherent mention of falls. However, this perspective was challenged by subsequent research from the same team in a large-scale study with 7,857 male patients.¹⁴ Here, they discovered that a patient’s fall frequency indeed has its independent predictive value for future osteoporotic fractures, regardless of their FRAX score. This key insight brought to light the value of including fall history in assessing fracture risk, suggesting that the FRAX score’s effectiveness could be enhanced by integrating this factor.

Our study builds upon this existing body of knowledge, indicating the improved performance of the FRAX score when prior falls are included in the model. While Harvey et al. focused explicitly on male subjects, in the present study, we did not observe any changes when the model was adjusted for sex.²⁹ In a study by Kanis et al. involving 21,116 male and female participants, a progressive increase in Hazard Ratios was observed with the rising number of falls.³⁰ When these ratios were adjusted for age, sex, and FARX variables, a slight decrease was noted. However, the present study focused specifically on patients with hip fractures caused by falls, providing distinct insights relevant to that population. Moreover, in our study, mainly involving elderly subjects, we did not find the energy of falls to be a significant predictor within the FRAX model.^{31, 32} Notably, we observed a low incidence of high-energy falls among both the cases and controls, aligning with existing literature that low-energy falls are more common in older adults.^{33, 34} Including fall energy in the FRAX model may not significantly affect fracture risk assessment in older adults, but its impact on younger populations warrants further investigation across various age groups.

In predicting fall-induced hip fractures, the significance of incorporating patient sociodemographic characteristics is increasingly recognized. Previous studies have underscored the impact of factors such as race, geographic location, age, economic conditions, and housing on the incidence and management of hip fractures. For instance, disparities in hip fracture rates and outcomes have been observed across different racial and ethnic groups,³⁵ while geographic variations highlight the influence of environmental and healthcare access differences.³⁶ This understanding has led to the adaptation of the FRAX in various countries, with each version tailored to reflect the specific characteristics and needs of the local population. Our study revealed a notable association between white race and an increased risk of fall-induced injuries, aligning with findings by Ellis et al., who reported that white individuals experienced approximately double the rate of fall-induced hip fractures compared to Hispanic and Black populations and triple the rate observed in Asian/Pacific Islanders.³⁷ It is imperative, however, to reflect on the demographic makeup of our study cohort, predominantly composed of white participants, which underscores the need for caution in extrapolating our findings to a broader populations.^38–40

Whether the broader sociodemographic factors captured by the SVI are of added value in fall-induced injuries is a topic that is not directly addressed in the literature. In this study, we found a significant inverse role for the SVI’s Housing Type and Transportation domain, meaning the less vulnerable individuals had a higher chance of sustaining a fall-induced hip fracture. This counterintuitive result warrants a deeper investigation into the dynamics of the Housing Type and Transportation domain, which includes the number of units in a building, the presence of mobile homes, the number of occupants per room, the presence of a vehicle in a household, and group living arrangements.⁴¹ One possible explanation could be that a higher number of occupants, typically associated with a higher vulnerability score, might paradoxically confer a protective effect against falls through heightened physical engagement.⁴² Bernhart et al. showed that older adults residing in households with two or more individuals were twice as likely to adhere to aerobic guidelines for physical activity compared to those living in single-adult households.⁴³ Similarly, Meghani et al concluded that residing with family plays a crucial role in assisting older adults to stay active and manage their personal responsibilities effectively.⁴⁴ While other socioeconomic factors seem to have a remarkable effect on hip fractures and falls separately, in this study, we did not find a role for other domains of the SVI in fall-induced hip fractures. In contrast, using the Housing-based Index of Socioeconomic Status (HOUSES), Ryu et al. concluded that a higher socioeconomic status reduces the likelihood of falls.⁴⁵ It is important to note, however, that their study included individuals across all age groups. Similarly, a study conducted in the UK highlighted the impact of material deprivation on the incidence of hip fractures among young adults, with a relative risk (RR=1.64; CI:1.57-1.72) in poorer electoral wards.⁴⁶ However, they reported no association in individuals aged 85 and over (RR=0.94, CI:0.87-1.01), suggesting the decreased effect of socioeconomic factors on the elderly, which is consistent with the outcomes of our study.

This study has a few limitations to be considered. While previous studies have shown different results with specific elderly age groups, in this study, because of a relatively small sample size, we could not conduct a group-wise analysis by age. Also, the energy of falls was measured based on a height criterion, and since oftentimes the height was not explicitly mentioned in clinicians’ notes, we relied on the interpretation of the clinical notes made by expert observers. However, if the expert observers were not confident about the interpretation of the fall energy, we would not include that factor for the patient in assessing our prediction models. Despite these limitations, this study is unique in its explicit inclusion of individuals with fall-induced hip fractures. Another strength is the direct extraction of prior fall history from patient records, which reduces the likelihood of recall bias that is often associated with fall questionnaires commonly used in similar studies. Furthermore, this study was the first of its kind, especially in the New England region of the United States, to consider different socioeconomic and demographic characteristics of patients and their relevance to fall-induced injuries.

Conclusion

This study highlighted the importance of the history and frequency of prior falls for the risk assessment of fall-induced hip fractures. Our results revealed the role of the Housing and Transportation subscale of the SVI, which also happens to be a modifiable risk factor. Adding these patient-specific factors to the widely used and accepted FRAX model has led to a modified FRAX that outperformed the FRAX tool in the prediction of fall-induced hip fractures. However, there is still potential for further improvement in the FRAX model. We recommend future research focus on age group stratification and validation of the impact of social determinants of health and socioeconomic characteristics of patients. This should include, but not be limited to, factors like the SVI and the Area Deprivation Index, as well as the characteristics of the falls themselves. These prediction models should be built using large and granular data and should be externally validated to achieve higher validity and generalizability.