- Research article
- Open Access
Non-steroidal anti-inflammatory drugs and the risk of a second hip fracture: a propensity-score matching study
BMC Musculoskeletal Disorders volume 17, Article number: 201 (2016)
Non-steroidal anti-inflammatory drugs (NSAIDs) are frequently prescribed for elderly patients, particularly after a hip fracture. However, we are not clear about the effect of NSAIDs on the risk of a second hip fracture because of confounding factors.
This was a Taiwan National Health Insurance Research Database-based study using propensity-score matching (PSM) to control for confounding. Enrollees were selected from patients with a hip fracture during 1996–2004 and followed longitudinally until December 2009. After PSM for comorbidities and bisphosphonate therapy, 94 patients with a second hip fracture were assigned to the Cases group and 461 without it to the Controls group. The target drugs are NSAIDs; paracetamol and dexamethasone are used for comparison.
The correlation between the mean daily-dose (MDD) ratios of NSAIDs and the probability values of the current statistical tests were highly negative (Pearson’s r = −0.920, P = 0.003), which indicated that the higher the MDD ratios, the greater the risks of a second hip fracture. A Kaplan-Meier survival analysis showed a time-dependent trend of increasing risk of a second hip fracture in patients taking NSAIDs (P < 0.001). Moreover, patients ≥60 years old had a higher risk of a second hip fracture than did those <60 and taking the NSAIDs diclofenac (P = 0.016) and celecoxib (P = 0.003) and the corticosteroid dexamethasone (P = 0.018), but not those taking analgesic paracetamol (P = 0.074).
We conclude that taking NSAIDs after a fragility hip fracture dose- and time-dependently significantly increases the risk of a second hip fracture, especially in elderly patients. To lower the risk of a second hip fracture, any underlying causes for excessively using NSAIDs should be treated and thus fewer NSAIDs prescribed after a first hip fracture.
The occurrence of a second hip fracture and subsequent mortality in patients with a fragility hip fracture is high [1–6]. Ryg et al., using Denmark’s National Hospital Discharge Register for the period 1997 to 2001, explored this sequela of osteoporosis . They found that patients with a fragility hip fracture had twice the risk of a second hip fracture, and that they had a 5-year mortality of approximately 60 %. Moreover, the estimated 1-year risk of a second hip fracture, whether in Western or Asian populations, should vary from 2 to 5 %, depending on the age of the patient . These findings highlight the importance of formulating and proposing a tertiary strategy for osteoporosis to prevent a subsequent hip fracture [3, 7].
We previously reported that, in addition to age, female gender, and comorbidities, the prolonged use of analgesics, e.g., paracetamol, and anti-inflammatory medications, e.g., dexamethasone and NSAIDs (non-steroidal anti-inflammatory drugs), is a significant risk factor for a second hip fracture after hip fracture surgery . Paracetamol, dexamethasone, and NSAIDs are commonly prescribed. Paracetamol is a mild analgesic and dexamethasone a potent anti-inflammatory corticosteroid. NSAIDs have analgesic effects and, in higher doses, anti-inflammatory effects because they inhibit the cyclooxygenase-2 (COX-2) activity and then reduce the synthesis of prostaglandins . Prostaglandins are potent, multifunctional regulators of bone metabolism, which can stimulate and inhibit bone resorption and formation [9–13]. In general, NSAIDs are supposed to reduce the rate of bone loss, and thus improve bone mineral density (BMD) and prevent a fragility hip fracture [14–16], but conflicting results have also been published [17, 18]. Additional analyses that use different approaches to control for confounding are necessary to clarify the degree of association and/or causal link between NSAIDs and the risk of a second hip fracture.
During the last decade, interest in determining the degree of association and/or causal link between medications and the risk of a fragility hip fracture has grown [17, 19]. However, establishing the degree of association and/or causal link is difficult because observational studies are notoriously vulnerable to the effect of different types of confounding . It is well recognized that the estimate of causality obtained by comparing cases with controls could be prejudiced because of problems such as selection bias or other systematic errors [20–22]. Rosenbaum and Rubin developed and popularized the propensity-score matching (PSM) method for observational (non-experimental) studies, in which only a small subset of controls comparable to cases must be selected ; PSM reduces the selection bias by balancing groups on the probability of being treated based on specific covariates . Many studies have used the PSM method to control for confounding.
In the current study, we used the PSM method to control for confounding and aimed to determine the strength of the association between NSAIDs and the risk for a second hip fracture in patients after hip fracture surgery. We also explored the age-specific risk of a second hip fracture in these patients.
The Taiwan National Health Insurance (NHI) program, a single-payer universal program that began in March 1995, now enrolls more than 99 % of the population. Its claims data for reimbursements provide one of the world’s largest datasets for health research. The National Health Insurance Research Database (NHIRD), derived from the payment system of the National Health Insurance Administration (NHIA) and maintained by the National Health Research Institute (NHRI) in Taiwan, provides vital information for research. The accuracy of the NHIRD diagnoses of major disorders, e.g., stroke and acute coronary syndrome, has been validated . The NHIRD includes patient demographics, disease diagnoses, medical care institution names, medical expenditure, and prescription claims data. For each medical expenditure reimbursement, the types of medical services and details of medical orders and costs are recorded. In the present study, data were obtained from the Longitudinal Health Insurance Database (LHID2000), a representative subset of 1 million patients from the NHIRD (http://nhird.nhri.org.tw/). All the identifiers of individual patients and medical care providers (medical professionals and institutions) are deleted by the NHIA before data are transferred to the NHIRD. Institutional review board (IRB) approval, an agreement to approve, monitor, and review biomedical and behavioral researches involving humans, is pre-approved by the NHRI for de-identified data.
Study design and participants
This study was a nationwide population-based observational study using PSM to control valuables that are measured at baseline for confounding. The study participants were selected from the NHIRD/LHID2000 registry covering the period from January 1996 to December 2004 and followed longitudinally until December 31, 2009. The study participants were identified from the database based on the following criteria: (1) a diagnosis code of hip fracture (ICD-9-CM diagnosis codes 820.0–820.9); (2) a procedure code of internal fixation or partial hip replacement (ICD-9-CM procedure codes 79.15, 79.35, and 81.52); and (3) patients 40 years old or older. The first admission date for a hip fracture was defined as the index date. After deleting the records of patients who died between 1996 and 2004 (n = 7), 1,545 patients with a hip fracture were included and followed longitudinally until the end of the study. During the follow-up period, patients with a second hip fracture were selected as Cases (n = 94). Four hundred sixty-one of 1,451 patients were selected as Controls using 1:5 matching on propensity scores for comorbidities and bisphosphonate therapy. By matching variables that would otherwise confound comparisons between groups, the PSM method effectively creates similar case and control sets from an existing dataset for an observational study . The Cases group consisted of 94 patients with a second hip fracture and the Controls group of 461 without. The flowchart of patient selection process is presented in Fig. 1.
Drug exposure and primary and secondary endpoints
NSAIDs, identified based on the anatomical therapeutic chemical (ATC) classification system , are known to work in two ways: (1) analgesic effects at lower doses and (2) anti-inflammatory effects at higher doses. To clarify the dose-response relationship between NSAIDs and the risk of a second hip fracture, we included paracetamol (a pure analgesic) and dexamethasone (a corticosteroid and potent anti-inflammatory drug) in the current study. The prescription history of each patient for every medication studied was summarized as either exposed or not exposed during the paired case and control periods. For each patient, to explore the possible dose response, we also estimated the mean daily dose (MDD) by calculating the cumulative dose of medications prescribed and then divided that number by the number of days in the period. The primary endpoint of this study was to determine the strength of association between drug exposure and the occurrence of a second hip fracture in patients after hip fracture surgery. For our secondary endpoint, information on the age-specific risk of a second hip fracture in these patients was collected for further analysis.
A χ 2 test was used to analyze categorical data. For numerical variables, Student’s t test was used for between-group comparisons. Univariate and multivariate analyses using a stepwise logistic regression model were done to detect significant predictive factors of a second hip fracture after hip fracture surgery. Between-group comparisons were done by estimating the odds ratio (OR) and the 95 % confidence interval (CI) in a logistic regression model. To examine the effects of MDD on the risk of a second hip fracture, the Pearson correlation test was used. The 15-year second hip fracture-free survival rate after hip fracture surgery was estimated using the Kaplan-Meier method and the log-rank test. The cumulative hazard function was also used and the proportional hazards assumption was met for each parameter in additional models. In Cox-regression proportional hazards survival analysis, potential predictors in the forced-entry model were used as covariates (e.g., age, gender, and comorbidities). Significance was set at p < 0.05 (two-sided). SAS 9.2 (SAS, Inc., Cary, NC, USA) was used for all analyses.
From 1996 through 2004, 555 patients ≥40 years who sustained a hip fracture and underwent hip fracture surgery were selected and enrolled in this study. There were no significant differences in comorbidities and bisphosphonate therapy between the Cases and Controls (Table 1). Patients with a second hip fracture after hip fracture surgery were older (mean age: 74.0 years vs. 69.5 years, P < 0.001). There were, however, no significant differences in gender distribution.
Patients in the Cases group had taken significantly more paracetamol and had a higher MDD than did those in the Controls group (35.8 % vs. 18.1 %, P = 0.026; 605.6 mg vs. 496.0 mg, P = 0.009, respectively); the same was true for diclofenac (41.7 % vs. 15.4 %, P < 0.001; 206.9 mg vs. 116.1 mg, P < 0.001), ibuprofen (28.1 % vs. 20.1 %, P < 0.001; 439.9 mg vs. 257.6 mg, P < 0.001), celecoxib (23.1 % vs. 9.5 %, P < 0.001; 305.3 mg vs. 198.2 mg, P < 0.001), and dexamethasone (39.4 % vs. 16.5 %, P < 0.001; 9.8 mg vs. 4.0 mg, P < 0.001). There were, however, no significant between-group differences for aspirin, naproxen, nabumetone, etodolac, or rofecoxib (all P ≥ 0.057) (Table 2). There was a highly negative correlation between the MDD ratios and the probability values of the current statistical tests (Pearson’s r = −0.920, P = 0.003 for NSAIDs only) (Fig. 2); thus, the higher the MDD ratios were, the greater risks of a second hip fracture the patients had.
Using the PSM method to control for confounding, age but not gender was identified as a significant predictor of a second hip fracture (Table 1). For every 1-year increase in age, there was at least a 2.4 % increase in the risk of a second hip fracture, and at least a 13.0 % increase for every 5-year increase in age (all P ≤ 0.044) (Table 3). Overall, the risk ratio (RR) was 14.8 % higher for female patients than for male patients (P < 0.001) (Table 4). Using 40–49-year-olds as the reference group, the RRs for a second hip fracture were significantly higher for those 60–69, 70–79, and ≥80 (all P ≤ 0.040), but not for those 50–59 (P = 0.627).
The 15-year survival analyses showed that male patients ≥60 were significantly more likely to have a second hip fracture after their hip fracture surgery than those <60 (P = 0.019). This was not true for the total group of patients or for the group of female patients (both P ≥ 0.092) (Fig. 3a–c). For the group as a whole, the RR for a second hip fracture was significantly (P < 0.001) and time-dependently higher for patients taking paracetamol, NSAIDs, and dexamethasone (Fig. 3d). Kaplan-Meier survival analyses also showed that the RR for a second hip fracture was significantly higher for patients ≥60 than for those <60 exposed to diclofenac (P = 0.016) (Fig. 3e), celecoxib (P = 0.003) (Fig. 3f), and dexamethasone (P = 0.018) (Fig. 3g), but not to paracetamol (P = 0.074) (Fig. 3h).
We confirmed that patients who take NSAIDs after hip fracture surgery have a significantly higher risk for a second hip fracture than do patients who do not take NSAIDs. Moreover, the positive association was dose- and time (years after first hip fracture surgery)-dependent for the group as a whole and dose- and age-dependent for patients ≥60. Furthermore, the higher the MDD of NSAIDs was, the greater the risk of a second hip fracture was. These findings highlight the importance of changing the treatment strategy of elderly patients after their first hip fracture.
To the best of our knowledge, this is the first report on confirming the association between taking NSAIDs and the risk of a second hip fracture in patients after hip fracture surgery. We previously reported that age, female gender, comorbidities, and the prolonged use of analgesics and anti-inflammatory medications are all significant risk factors for a second hip fracture after hip fracture surgery, but that bisphosphonate therapy was protective . However, in an observational study, there are many confounders that usually distort the relationship between an exposure and an outcome [20–22, 25]. The PSM method is useful in these circumstances because it provides a neutral weighting formula that yields unbiased estimates of the effects of treatment. In the current study, we thus used the propensity score model for controlling variables that are measured at baseline. After using the PSM for controlling comorbidity and bisphosphonate therapy, age but not gender remains to be a significant risk factor for the occurrence of a second hip fracture in patients after hip fracture surgery.
Although users of NSAIDs are reported to have higher BMD than do non-users [19, 26, 27], the effect of NSAIDs on the risk of a hip fracture is still not understood [14–18]. Two studies found that the hip fracture risk was lower in patients taking NSAIDs [14, 17]. In contrast, a report based on the Danish Osteoporosis Prevention Study (DOPS) showed that users of NSAIDs had more hip fractures than expected . Leaving aside the debate whether NSAIDs increase the risk of a first hip fracture, the authors highlight the issue of tertiary prevention for osteoporosis because of the healthcare resources constraint [2, 3]. Our data revealed that the higher the MDD ratios of NSAIDs, the greater the likelihood of a second hip fracture in patients after hip fracture surgery. One possible explanation is that, when prolonged, the anti-inflammatory effects of NSAIDs at higher doses elevate the risks of a second hip fracture. Another explanation is that NSAID users have an impaired gait because of the sequelae of their index fragility hip fracture. No matter what the mechanism is, we recommend that NSAIDs should be cautiously prescribed for patients who have just undergone hip fracture surgery, particularly for elderly patients. As with other drugs prescribed for elderly patients [17, 28, 29], the most judicious approach is to restrict NSAIDs to situations in which their benefits clearly overweigh their risks, and to use them only after any underlying causes for excessive drug use have been adequately treated and potentially safer alternatives have been tried. When treatment with NSAIDs is necessary, the lowest feasible dose should be used for the shortest duration required to achieve the desired effect [28, 29].
This study has some limitations. First, it is retrospective. Second, it might still have an uncontrolled selection and recall bias despite our using stratification and the PSM method to control for confounders. For example, the presence of a diagnosis of comorbidity is different from controlling for the severity of the disease, which may be expected to be different in older individuals. Third, the risk of a second hip fracture may be underestimated because we included the survival rate after the index hip fracture only and also excluded patients who died from any cause between 1996 and 2004. We could not account for patients who sustained a hip fracture 2 years before the index date and thus wrongly had their second hip fracture classified as an index hip fracture despite their having no hips at risk. Finally, we lacked information on what caused the higher exposure rate and higher MDD of NSAIDs in our patients after their hip fracture surgery.
In conclusion, taking NSAIDs after hip fracture surgery significantly increases the risk of a second hip fracture. The positive association is dose- and time-dependent for all patients and age-dependent for the elderly. Because of a global deterioration of health conditions in elderly patients after a fragility hip fracture , they are usually prescribed significantly more drugs, especially analgesics and anti-inflammatory agents. Physicians should limit prescribing NSAIDs to situations with a positive benefit-risk balance and use them only after the underlying cause for excessively using NSAIDs have been adequately treated and potentially safer alternatives have failed. The lowest feasible dose for the shortest duration required to achieve the desired effect should be considered. If more long-range treatment with NSAIDs is necessary and inevitable, more aggressive monitoring and prevention for another fragility fracture is warranted [30, 31], because long-term NSAIDs treatment might contribute to rises in the risk for a second hip fracture.
Ethic approval or consent to participate
Consent for publication
Availability of data and materials
The website of the Taiwan National Health Insurance Research Database is found at http://nhird.nhri.org.tw/.
anatomical therapeutic Chemical
bone mineral density
Danish Osteoporosis Prevention Study
9th revision of International Classification of Diseases, Clinical Modification
Institutional review board
Longitudinal Health Insurance Database
National Health Insurance
National Health Insurance Administration
National Health Insurance Resarch Database
Non-Steroidal Anti-Inflammatory Drugs
Kanis JA, Odén A, McCloskey EV, Johansson H, Wahl DA, Cooper C. A systematic review of hip fracture incidence and probability of fracture worldwide. Osteoporos Int. 2012;23(9):2239–56.
Ryg J, Rejnmark L, Overgaard S, Brixen K, Vestergaard P. Hip fracture patients at risk of second hip fracture: a nationwide population-based cohort study of 169,145 cases during 1977–2001. J Bone Miner Res. 2009;24(7):1299–307.
Shen SH, Huang KC, Tsai YH, et al. Risk analysis for second hip fracture in patients after hip fracture surgery: a nationwide population-based study. J Am Med Dir Assoc. 2014;15(10):725–31.
Lönnroos E, Kautiainen H, Karppi P, Hartikainen S, Kiviranta I, Sulkava R. Incidence of second hip fractures. A population-based study. Osteoporos Int. 2007;18(9):1279–85.
Khosla S, Amin S, Orwoll E. Osteoporosis in men. Endocr Rev. 2008;29(4):441–64.
Pietschmann P, Rauner M, Sipos W, Kerschan-Schindl K. Osteoporosis: an age-related and gender-specific disease - a mini-review. Gerontology. 2009;55(1):3–12.
Vestergaard P, Rejnmark L, Mosekilde L. Hip fracture prevention: Cost-effective strategies. Pharmacoeconomics. 2001;19(5 Pt 1):449–68.
Raisz LG. Potential impact of selective cyclooxygenase-2 inhibitors on bone metabolism in health and disease. Am J Med. 2001;110(Suppl 3A):43S–5S.
Kawaguchi H, Pilbeam CC, Harrison JR, Raisz LG. The role of prostaglandins in the regulation of bone metabolism. Clin Orthop Relat Res. 1995;313:36–46.
Raisz LG. Prostaglandins and bone: physiology and pathophysiology. Osteoarthritis Cartilage. 1999;7(4):1–3.
Kawaguchi H, Pilbeam CC, Vargas SJ, Morse EE, Lorenzo JA, Raisz LG. Ovariectomy enhances and estrogen replacement inhibits the activity of bone marrow factors that stimulate prostaglandin production in cultured mouse calvariae. J Clin Invest. 1995;96(1):539–48.
Jee WS, Ma YF. The in vivo anabolic actions of prostaglandins in bone. Bone. 1997;21(4):297–304.
Faye-Petersen OM, Johnson Jr WH, Carlo WA, et al. Prostaglandin E1-induced hyperostosis: clinicopathologic correlations and possible pathogenetic mechanisms. Pediatr Pathol Lab Med. 1996;16(3):489–507.
Rashiq S, Logan RF. Role of drugs in fractures of the femoral neck. BMJ. 1986;292(6524):861–3.
van Staa TP, Leufkens HG, Cooper C. Use of nonsteroidal anti-inflammatory drugs and risk of fractures. Bone. 2000;27(4):563–8.
Thompson DD, Rodan GA. Indomethacin inhibition of tenotomy-induced bone resorption in rats. J Bone Miner Res. 1988;3(4):409–14.
Rossini M, Viapiana O, Adami S, et al. Medication use before and after hip fractures: a population based cohort and case-control study. Drugs Aging. 2014;31(7):547–53.
Vestergaard P, Hermann P, Jensen JE, Eiken P, Mosekilde L. Effects of paracetamol, non-steroidal anti-inflammatory drugs, acetylsalicylic acid, and opioids on bone mineral density and risk of fracture: results of the Danish Osteoporosis Prevention Study (DOPS). Osteoporos Int. 2012;23(4):1255–65.
Carbone LD, Tylaysky FA, Cauley JA, et al. Association between bone mineral density and the use of nonsteroidal anti-inflammatory drugs and aspirin: impact of cyclooxygenase selectivity. J Bone Miner Res. 2003;18(10):1795–802.
Mann B, Wood E. Confounding in observational studies explained. Open Epidemiol J. 2012;5:18–20.
Rosenbaum PR, Rubin DB. The central role of the propensity score in observational studies for causal effects. Biometrika. 1983;70(1):41–55.
Dehejia RH, Wahba S. Propensity score-matching methods for nonexperimental causal studies. Rev Econ Statistics. 2002;84(1):151–61.
Cheng CL, Kao YH, Lin SJ, Lee CH, Lai ML. Validation of the national health insurance research database with ischemic stroke cases in Taiwan. Pharmacoepidemiol Drug Saf. 2011;20(3):236–42.
The Anatomical Therapeutic Chemical (ATC) classification system and the Defined Daily Dose (DDD). http://www.whocc.no/
Goodman SN. Toward evidence-based medical statistics. 1: The P value fallacy. Ann Intern Med. 1999;130(12):995–1004.
Bauer DC, Orwoll ES, Fox KM, et al. Aspirin and NSAID use in older women: effect on bone mineral density and fracture risk. Study of Osteoporotic Fractures Research Group. J Bone Miner Res. 1996;11(1):29–35.
Wilkin TJ, Devendra D. Bone densitometry is not a good predictor of hip fracture. BMJ. 2001;323(7316):795–7.
Rochon PA, Gurwitz JH. Optimising drug treatment for elderly people: the prescribing cascade. BMJ. 1997;315(7115):1096–9.
Ekman EF. The role of the orthopaedic surgeon in minimizing mortality and morbidity associated with fragility fractures. J Am Acad Orthop Surg. 2010;18(5):278–85.
Kannus P, Parkkari J, Niemi S, et al. Prevention of hip fracture in elderly people with use of a hip protector. N Engl J Med. 2000;343(21):1506–13.
McLellan AR, Gallacher SJ, Fraser M, McQuillian C. The fracture liaison service: success of a program for the evaluation and management of patients with osteoporotic fracture. Osteoporos Int. 2003;14(12):1028–34.
We thank Fei-Kai Syu; for her assistance with data collection and statistical analysis.
No resource of funding.
The authors declare that they have no competing interests.
PYC, SHS and KCH participated in he conception and design of the study. PYC, TYY and TWH performed the literature search and selection of studies. PYC, SHS and TYY extracted relevant data. PYC, SHS and KCH participated in the analysis and interpretation of data and in the preparation of the manuscript. All authors read and approved the final manuscript.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Chuang, PY., Shen, SH., Yang, TY. et al. Non-steroidal anti-inflammatory drugs and the risk of a second hip fracture: a propensity-score matching study. BMC Musculoskelet Disord 17, 201 (2016). https://doi.org/10.1186/s12891-016-1047-2
- Non-steroidal anti-inflammatory drugs
- Fragility hip fracture
- Second hip fracture
- Population-based study
- Propensity-score matching