Disentangling the contribution of hospitals and municipalities for understanding patient level differences in one-year mortality risk after hip-fracture: A cross-classified multilevel analysis in Sweden

Abstract

Background

One-year mortality after hip-fracture is a widely used outcome measure when comparing hospital care performance. However, traditional analyses do not explicitly consider the referral of patients to municipality care after just a few days of hospitalization. Furthermore, traditional analyses investigates hospital (or municipality) variation in patient outcomes in isolation rather than as a component of the underlying patient variation. We therefore aimed to extend the traditional approach to simultaneously estimate both case-mix adjusted hospital and municipality comparisons in order to disentangle the amount of the total patient variation in clinical outcomes that was attributable to the hospital and municipality level, respectively.

Methods

We determined 1-year mortality risk in patients aged 65 or above with hip fractures registered in Sweden between 2011 and 2014. We performed cross-classified multilevel analysis with 54,999 patients nested within 54 hospitals and 290 municipalities. We adjusted for individual demographic, socioeconomic and clinical characteristics. To quantify the size of the hospital and municipality variation we calculated the variance partition coefficient (VPC) and the area under the receiver operator characteristic curve (AUC).

Results

The overall 1-year mortality rate was 25.1%. The case-mix adjusted rates varied from 21.7% to 26.5% for the 54 hospitals, and from 18.9% to 29.5% for the 290 municipalities. The VPC was just 0.2% for the hospital and just 0.1% for the municipality level. Patient sociodemographic and clinical characteristics were strong predictors of 1-year mortality (AUC = 0.716), but adding the hospital and municipality levels in the cross-classified model had a minor influence (AUC = 0.718).

Conclusions

Overall in Sweden, one-year mortality after hip-fracture is rather high. However, only a minor part of the patient variation is explained by the hospital and municipality levels. Therefore, a possible intervention should be nation-wide rather than directed to specific hospitals or municipalities.

About the Authors:

Pia Kjær Kristensen

Contributed equally to this work with: Pia Kjær Kristensen, Juan Merlo

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Validation, Visualization, Writing – original draft, Writing – review & editing

* E-mail: [email protected]

Affiliations Department of Clinical Epidemiology, Aarhus University Hospital, Aarhus, Denmark, Department of Orthopedic Surgery, Regional Hospital Horsens, Horsens, Denmark

ORCID logo http://orcid.org/0000-0001-5473-9386

Raquel Perez-Vicente

Roles Data curation, Formal analysis, Software, Validation, Writing – review & editing

Affiliation: Research Unit of Social Epidemiology, Clinical Research Centre, Faculty of Medicine, Lund University, Malmö, Sweden

George Leckie

Roles Conceptualization, Methodology, Software, Validation, Writing – review & editing

Affiliation: Centre for Multilevel Modelling, School of Education, University of Bristol, United Kingdom

ORCID logo http://orcid.org/0000-0003-1486-745X

Søren Paaske Johnsen

Roles Supervision, Validation, Writing – review & editing

Affiliation: Danish Center for Clinical Health Services Research, Department of Clinical Medicine, Aalborg University, Aalborg, Denmark

Juan Merlo

Contributed equally to this work with: Pia Kjær Kristensen, Juan Merlo

Roles Data curation, Funding acquisition, Methodology, Project administration, Resources, Software, Supervision, Validation, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Research Unit of Social Epidemiology, Clinical Research Centre, Faculty of Medicine, Lund University, Malmö, Sweden

ORCID logo http://orcid.org/0000-0001-8379-9708

Introduction

Increasing demands for healthcare fueled by aging populations, increasing prevalence of multimorbid patients and technological advances combined with restricted financial resources constitute a major challenge for health care systems globally [1]. To address this challenging task many countries have started to streamline their healthcare systems by improving care and reducing length of stay in hospitals, which includes earlier discharge to care in community settings or at home with professional support (i.e., homecare) [2, 3]. Hip fracture care is an example of this development where surgery is performed at the hospital while postoperative rehabilitation and further care of potential sequelae are done in the primary health care sector. Therefore, evaluating the quality of care after hip fracture needs to simultaneously consider both the hospital and the municipality settings. When evaluating quality of care after hip fracture, one-year mortality is frequently used as an outcome quality indicator in many healthcare systems [4]. From this perspective, identifying unwarranted variability between not only hospitals but also municipalities has considerable health policy relevance [5–7].

It is well-established that multilevel modelling is a suitable methodology when comparing different types of cluster (e.g., hospital, municipalities) for both statistical and conceptual reasons [8–18]. However, existing multilevel model studies evaluating healthcare performance after hip fracture have primarily focused on hospital level variation [12, 19–22]. Less work has been carried out analysing the potential variability among municipalities [23]. Cross-classified multilevel analysis allows us to decompose the total individual level variation into its hospital and its municipality level components [24, 25]. To our knowledge no previous studies have applied measures of components of variance and of discriminatory accuracy to examine individual variation in one-year mortality with respect to both hospitals and municipalities. Therefore, we use cross-classified multilevel analysis to evaluate one-year mortality in all patients suffering from hip fracture in Sweden.

Population and methods

Databases and study population

This historical follow-up study was based on prospectively collected data available from medical registries in Sweden (about 9.3 million inhabitants by 2010) with free access to medical care. At birth or upon immigration, all citizens in Sweden are assigned a unique registration number through which all contact with the healthcare system is recorded. This allows unambiguous record linkage between registries.

We obtained information on diagnoses, medication use and mortality from the Swedish Patient Register, the Swedish Prescribed Drug Register [26], and from the Cause of Death Register [27] administered by the National Board of Health and Welfare. We obtained demographic and socioeconomic information from the Population Register [28] and the Longitudinal Integration Database for Health Insurance and Labour Market Studies (LISA) register (http://www.scb.se/lisa), which are administrated by Statistics Sweden. To ensure the anonymity of the subjects, the Swedish authorities transformed the personal identification numbers of the individuals into arbitrary numbers before delivering the research databases to us, and we linked the databases using the anonymized identification number.

We defined hip fracture as the presence of any discharge diagnosis with a fracture of the femur coded according the International Classification of Diseases, 10th revision (ICD codes) as fracture of the femoral neck (S72.0), pertrochanteric (S72.1) or subtrochanteric (S72.2). From the Swedish Patient Registry, we identified all 56,161 patients being 65 years or older and residing in Sweden by 31th December 2010 who were discharged from the Swedish hospitals between 2011 and 2014 with a diagnosis of hip (i.e., femoral) fracture and surgery code. Next, we excluded 1,162 patients because they were residing less than five years in Sweden (n = 173), they have missing information on education (n = 976) or erroneous information on death date (n = 2). We also excluded 11 patients treated at facilities with less than 10 hip fracture patients during the study period (Fig 1). Patients which previous hip fracture compose of 11.95% (6,575/48,424) of our population.

[Figure omitted. See PDF.]

Fig 1. Flowchart patient inclusion.

https://doi.org/10.1371/journal.pone.0234041.g001

Ethics and data accessibility

The Regional Ethics Review Board in southern Sweden (Dnr: 2014/856) as well as the data safety committees from the National Board of Health and Welfare and from Statistics Sweden approved the construction of the database used in this study.

The database we analysed is not publicly available for ethical and data safety reasons. However, the same dataset can be constructed by request to the Swedish National Board of Health and Welfare and Statistics Sweden after approval of the research project by an Ethical Committee and by the data safety committees of the Swedish Authorities. The study also needs to be performed in collaboration with Swedish researchers. [29]

Assessment of variables

Mortality.

For each patient we ascertained all-cause mortality within one year from the admission date to hospital.

Sociodemographic and clinical characteristics of the patients.

Part of the initial differences between hospitals or municipalities may relate to differences in case mix. To make the observational measurement of hospital and municipality effects as valid as possible we therefore adjusted for potential differences in patient sociodemographic and clinical characteristics as well as use of medication (Table 1).

[Figure omitted. See PDF.]

Table 1. Characteristic of the hip fracture population.

Values are percentages (number of patients) if not otherwise indicated.

https://doi.org/10.1371/journal.pone.0234041.t001

Sociodemographic characteristics.

We use the age of the individuals in years as a continuous variable and included it in the analysis as a quadratic function. We categorized the patients as immigrants if they were born in another country and as natives if they were born in Sweden.

We calculated household individualized disposable income by dividing the total disposable income of the family of the patient by the number of family members, considering the different consumption weights of adults and children, according to Statistics Sweden (http://www.scb.se/lisa). We did so for the complete Swedish population in three occasions 2010, 2005 and 2000, and summarized the three occasion by computing the cumulative income. Finally, using the tertile values of the cumulative income distribution we divided the study population into high, medium and low income. In this form, we categorized each patient according his/her income category in the total population rather than the population of patients.

We classified education achievement as compulsory schooling of nine years or less, secondary education three years or less, and higher education. We defined low education achievement as compulsory school or less and use the value of the member of the family with the highest educational achievement for all member in the same family.

For defining those patients living alone, we first identified all patients who were cohabiting. That is, married couple, in a registered partnership or in an unregistered partnership with a common child as well as those living in a household with at least another adult. We group all other patients into the living alone category.

We combined sociodemographic characteristics into a single sociodemographic risk score for all-cause mortality. We did this via a logistic regression analysis modelling all-cause mortality as a function of age, gender, income, cohabiting, born in Sweden and education (see S1 Table). The individual predicted probabilities were then categorized into four groups by quartiles as low, medium, high, and very high. The low-risk score group was then used as the reference in the comparisons in the subsequent multilevel analyses.

Clinical characteristics.

For the purpose of patient-mix adjustment and using previous knowledge on risk factors for mortality in patients with hip fractures [30, 31] as well as considering the variables included in the Charlson comorbidity index (CCI), we selected a number of diseases (ICD-10 codes) identified during the five years before hospital admission (see S2 Table). We then applied logistic regression to obtain a clinical risk score (i.e., individual predicted probability) of one-year mortality based on the following variables: Chronic kidney disease (N18), Acute myocardial infarction (I21), Heart failure (I50), Other peripheral vascular diseases (I73), Cerebrovascular diseases (I60-I69), Atherosclerosis- Aortic aneurysm and dissection (I70-I71), Dementia (F01-F03), Chronic lower respiratory diseases (J40-J47), Other disorders of the skin and subcutaneous tissue (L80-L99), Peptic ulcer (K27), Diseases of liver (K70-K77), Diabetes mellitus (E08-E13), Hemiplegia (G81), Neoplasms (C00-D49), Human Immunodeficiency Virus (HIV) (B20), osteoporosis (M80-M81), previous Hip fracture (S70-S72). Finally, we created four categories of clinical risk scores using the quartile values of the risk score distribution and considering the group with the lowest risk score as reference in the comparisons.

Use of medication.

Using the Anatomical Therapeutic Chemical (ATC) classification from the Swedish Prescribed Drug Register we also obtained information on previous use of Analgesics (ATC cod: N02), Psycholeptics (N05), Psychoanaleptics (N06) and Bisphosphonates (M05BA) as these medications may influence the vital prognosis of the patients [32, 33].

Statistical analysis

We applied single-level (conventaional) and multilevel logistic regression analyses of discriminatory accuracy as decribed elsewhere [18]. We developed three consecutive logistic regression analyses to model one-year mortality.

The first model (model 1) is a single-level logistic regression including the socioeconomic risk score in four groups. This model aimed to evaluate the influence of patients’ demographical characteristic on one-year mortality.

The second (model 2) is also a single-level logistic regression that expands model 1 to include the risk score for clinical factors as well as use of medication.

In the final, third model (model 3), we expanded model 2 by adding two random effects, one for the hospital level and the other for municipality level. In this way the model was converted into a two-way cross-classified multilevel model with the patients nested within the 54 hospitals and the 290 municipalities. As described elsewhere [18], we calculated the hospital and municipality general contextual effects as expressed by the variance partition coefficient (VPC) which informs on the share of the total individual variance in the propensity of one-year mortality that is at the hospital (VPCH) and at the municipality level (VPCM).

For all models, we calculated the Area Under the receiving operator characteristics Curve (AUC) or C-statistics [34] as a measure of discriminatory accuracy. Hosmer and Lemeshow [35] suggest that an AUC of 0.70 to 0.80 could be considered as 'acceptable', 0.80 to 0.90 as 'excellent' and 0.90 or above as 'outstanding', while an area under the ROC curve of 0.50 suggests no discrimination between the outcome groups (i.e., similar as tossing a coin to decide group membership). In model 3 we separately calculated the AUC when adding the random effect for the hospital (AUCH) and the random effect for the municipality (AUCH). Then, we calculated the increment in the AUC [18] when going from model 1 to model 2 (AUCΔ2−1) and from model 2 to model 3 (AUCΔ3−2). The increment in AUC measures the improvement in the ability of the model to correctly classify individuals with or without the outcome (i.e., one-year mortality) when considering the hospital or the municipality of the patients. The AUC provides analogous information as the VPC [24] as both measures inform on the general contextual effect of the hospital and of the municipality levels in relation to the patients’ one-year mortality. As it has been explained elsewhere [12, 13, 18, 36–41] the general contextual effect, expresses the relevance of the hospital/municipality context for understanding patients’ differences in one-year mortality.

In model 3 we use the predicted hospital and municipality random effects (shrunken residuals) to obtain case-mix adjusted and reliability-weighted average mortality rates and their 95% credible intervals and created league tables illustrating the ranking of the hospitals and of the municipalities.

In the supplementary material S1 File we provided an extended explanation of the methodology.

Results

Characteristics of the hip fracture population

The Swedish cohort of hip fracture patients consists of 54,999 patients. The unadjusted population one-year mortality was 25.1%. The hip fracture patients were mainly women, and most of the patients were living alone. Table 1 describes the additional characteristics of the population.

Patient effects

High sociodemographic scores and high clinical scores were strong predictors for one-year mortality (Table 2). The AUC1 in model 1, which informs on the discriminatory accuracy of the sociodemographic information, had a value of 0.667 (95% CI: 0.662–0.672). Including the clinical score of the patients and use of medicines increased the AUC to 0.716 (95% CI: 0.711–0.720) (model 2). The use of analgesics, psycholeptics or psychoanaleptics was associated with an increased risk of one-year mortality (Table 2). However, use of bisphosphonates was not associated with one-year mortality.

[Figure omitted. See PDF.]

Table 2. Analysis of 1-year mortality after hip fracture in the Swedish hospitals.

https://doi.org/10.1371/journal.pone.0234041.t002

Hospital effects

In the period from 2011 to 2014, 54 hospitals treating patients with hip fracture were included in our study. These hospitals treated between 158 and 3,724 patients. The hospital unadjusted one-year mortality expressed as percentage varied between 19.6% and 29.8%. The league table obtained from the multilevel analysis shows that after adjustment for patient case-mix and municipality effects the hospital mortality rates extended from 21.7% to 26.5% (Fig 2). The general contextual effect of the hospital on the patients' risk for one-year mortality was, however, low as only 0.2% of the adjusted individual variation in one-year mortality laid at the hospital level. Also, the increase in the discriminative accuracy when adding the hospital level to the model including only patient level variables was only 0.002 units (model 2).

[Figure omitted. See PDF.]

Fig 2. League table ranking the 54 hospitals according to their adjusted absolute risk of 1-year mortality with 95% confidence intervals obtained from the cross-classified multilevel model.

https://doi.org/10.1371/journal.pone.0234041.g002

Municipality effects

After admission the hip fracture patients were discharged to 290 different municipalities in Sweden. The number of hip fracture patients in each municipality varied from 15 to 4,687 hip fracture patients in the three years. The unadjusted mortality rate expressed as percentage varies between 7.1% and 37.9%. After adjustment for the difference in case-mix, and the hospital effects in the multilevel analysis the municipality differences in one-year mortality extended between 18.9% and 29.5% (Fig 3). The general contextual effect of the municipalities on the patients' risk for one-year mortality was also low as only 0.1% of the total individual variance in the propensity of dying was at the municipality level. Also, the increase in the discriminative accuracy when adding the municipality level to the model including only patient level variables was only 0.002 units (model 2).

[Figure omitted. See PDF.]

Fig 3. League table ranking the 290 municipalities according to their adjusted absolute risk of 1-year mortality with 95% confidence intervals obtained from the cross-classified multilevel model.

https://doi.org/10.1371/journal.pone.0234041.g003

Discussion

In this nationwide population-based study of Swedish hip fracture patients, the crude one-year mortality rate was 25.1% during the period 2011–2014. Patient sociodemographic and clinical characteristics were predictors of one-year mortality (AUC = 0.718), whereas the cross-classified multilevel analysis revealed that the AUC of the regressions model including hospitals and municipality as random effects was only marginally higher compared to the single-level model. Accordingly, hospital and municipality level variation corresponded to less than 1% of the overall individual variation in the underlying propensity of death within one year.

For evaluating the performance of the hospital and municipality health care performance we considered two measures, the overall unadjusted one-year mortality rate after hip fracture and the size of the general contextual effects. In our study the overall unadjusted one-year mortality rate of 25.08% is comparable to a previous Danish register study, which in the period 2011 to 2014 reported a one-year mortality of 26% [42]. However, our mortality rate was higher than that found in a previous study on seven European countries, which found a mortality rate of 22.3% in Sweden in 2007 [43]. In this European study the lowest mortality, 19%, was observed in Italy. Those differences might be explained by differences in study design and changing comorbidity patterns across time as well as country differences in healthcare provision [42, 44].

After adjustment for patient case-mix, the size of the general context effect is assumed to inform on the influence of the hospital and municipality health care levels on the vital prognosis of the patients [24]. The VPCs we found in this study were small, being 0.1% for municipalities and 0.2% for hospitals. Meaning that once we have adjusted for patient characteristics, patients treated at the same hospital or at the same municipality have very little in common according to their propensity of death within one year. That is, there were very small differences between hospitals and between municipalities (se elsewhere for an extended explanation on this concept [24])

To our knowledge no previous studies have applied cross-classified multilevel analyses to hip fracture data and so a direct comparison is therefore not possible. However, our results are in line with a previous Danish multilevel study examining hospital variation in 30-day mortality among hip fracture patient (but ignoring municipality variation), which found a VPC for hospitals of 0.9%.

The AUC of the cross-classified multilevel analysis was slightly higher compared to the single-level model, which still indicated that the cross-classified multilevel analysis was a better model than the single-level model even though the VPC was low for both the hospital and the municipality levels. However, the AUCs of the models were overall moderate, ranging from 0.667 to 0.718.

In our study bisphosphonates were not associated with reduced mortality. In randomized controlled trials, use of bisphosphonates have been shown to be effective for preventing both osteoporotic fractures [45–47] and mortality after hip fracture [48]. The protective effect of bisphosphonates on mortality after hip fracture has also been observed in numerous observational studies. However, recent research has suggested that the reduction in mortality occurs within the first weeks after treatment and, therefore, could express confounding [33].

Our study was strengthened by the cross-classified multilevel study design, as well as by the nationwide population-based coverage with complete follow-up on one-year mortality due to linkage to Swedish registries. However, our study only analysed the Swedish hospitals and municipalities and, therefore, cannot be generalized to other country contexts.

The treatment, care and rehabilitation of hip fracture patients at both hospitals and municipalities are fundamentally important. We found that hospital and mortality differences in average mortality risks were only a minor share of the total individual variance in the propensity to die or that, in other words, adding the hospital and municipality level does not increase the discriminatory accuracy obtained by patient level information only. This means that the hospital and municipality performances are homogenous overall in the country. Therefore, our results indicate that we will not lower mortality after hip fracture in Sweden by focusing on specific hospitals or specific municipalities with high average mortality. Rather, special efforts to reduce mortality after hip fracture should be focused on vulnerable patient groups of hip fracture patients wherever they are. From this reason, however, hospitals and municipalities with a higher number of vulnerable patients than other will need more intense interventions. Future multilevel studies among hip fracture patients need to investigate differences in provided care and intermediate outcomes such as complications [49].

In conclusion, using cross-classified multilevel regression analysis, we observed that overall in Sweden, one-year mortality after hip-fracture at 25.1% was rather high during the period 2011–2014. We also observed that the average one-year adjusted mortality varied between 21.7% and 26.5% for hospitals and between 18.9% and 29.5% for municipalities. However, while the patient socioeconomic and clinical characteristics appears relevant for predicting mortality, only a minor part of the patient variation was explained by the hospital and municipality levels. Therefore, future interventions should focus on identifying high risk patient groups and be nationwide rather than directed at specific hospital or municipalities.

Supporting information

[Figure omitted. See PDF.]

S1 Table. Sociodemographic risk score.

https://doi.org/10.1371/journal.pone.0234041.s001

(DOCX)

S2 Table. Clinical risk score.

https://doi.org/10.1371/journal.pone.0234041.s002

(DOCX)

S1 File. Statistical appendix.

https://doi.org/10.1371/journal.pone.0234041.s003

(DOCX)

Acknowledgments

We thank the staff of the hospital departments caring for patients with hip fracture for their continuous effort and contribution to acquisition of the data in the registries.

Citation: Kristensen PK, Perez-Vicente R, Leckie G, Johnsen SP, Merlo J (2020) Disentangling the contribution of hospitals and municipalities for understanding patient level differences in one-year mortality risk after hip-fracture: A cross-classified multilevel analysis in Sweden. PLoS ONE 15(6): e0234041. https://doi.org/10.1371/journal.pone.0234041

References

1. Bungener M, Pierret J. De l'influence du chômage sur l'état de santé. In: Erès E, editor. Trajectoires sociales et inégalités1994.

2. Clarke A, Rosen R. Length of stay. How short should hospital care be? European journal of public health. 2001;11(2):166–70. pmid:11420803

3. Nordstrom P, Gustafson Y, Michaelsson K, Nordstrom A. Length of hospital stay after hip fracture and short term risk of death after discharge: a total cohort study in Sweden. BMJ (Clinical research ed). 2015;350:h696. pmid:25700551

4. Saez-Lopez P, Branas F, Sanchez-Hernandez N, Alonso-Garcia N, Gonzalez-Montalvo JI. Hip fracture registries: utility, description, and comparison. Osteoporosis international: a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 2017;28(4):1157–66. Epub 2016/11/23. pmid:27872956.

5. Corallo AN, Croxford R, Goodman DC, Bryan EL, Srivastava D, Stukel TA. A systematic review of medical practice variation in OECD countries. Health Policy. 2014;114(1):5–14. Epub 2013/09/24. pmid:24054709.

6. Fisher ES, Wennberg DE, Stukel TA, Gottlieb DJ, Lucas FL, Pinder EL. The implications of regional variations in Medicare spending. Part 1: the content, quality, and accessibility of care. Annals of Internal Medicine. 2003;138(4):273–87. 200302180–00006 [pii]. pmid:12585825

7. Fisher ES, Wennberg DE, Stukel TA, Gottlieb DJ, Lucas FL, Pinder EL. The implications of regional variations in Medicare spending. Part 2: health outcomes and satisfaction with care. Annals of Internal Medicine. 2003;138(4):288–98. 200302180–00007 [pii]. pmid:12585826

8. Merlo J. Changing analytical approaches in European epidemiology—a short comment on a recent article. European Journal of Epidemiology. 2005;20(8):737o–. WOS:000231732800014. pmid:16151890

9. Goldstein H, Spiegelhalter D. League Tables and Their Limitations: Statistical Issues in Comparisons of Institutional Performance. Journal of the Royal Statistical Society, SocA. 1996;159:385–443.

10. Leyland AH, Boddy FA. League tables and acute myocardial infarction. Lancet. 1998;351(9102):555–8. pmid:9492774.

11. Merlo J, Ostergren P, Broms K, Bjorck-Linne A, Liedholm H. Survival after initial hospitalisation for heart failure: a multilevel analysis of patients in Swedish acute care hospitals. Journal of Epidemiology and Community Health. 2001;55(5):323–9. WOS:000168138600010. pmid:11297650

12. Kristensen PK, Merlo J, Ghith N, Leckie G, Johnsen SP. Hospital differences in mortality rates after hip fracture surgery in Denmark. Clinical Epidemiology. 2019;11:605. pmid:31410068

13. Ghith N, Wagner P, Frolich A, Merlo J. Short Term Survival after Admission for Heart Failure in Sweden: Applying Multilevel Analyses of Discriminatory Accuracy to Evaluate Institutional Performance. PLoS One. 2016;11(2):e0148187. Epub 2016/02/04. pmid:26840122; PubMed Central PMCID: PMC4739586.

14. Ash AS, Fienberg SE, Louis TA, Normand SLT, Stukel TA, Utts J. STATISTICAL ISSUES IN ASSESSING HOSPITAL PERFORMANCE. Commissioned by the Committee of Presidents of Statistical Societies. The COPSS-CMS White Paper Committee:. 2011.

15. Adab P, Rouse AM, Mohammed MA, Marshall T. Performance league tables: the NHS deserves better. Bmj. 2002;324(7329):95–8. Epub 2002/01/12. pmid:11786455; PubMed Central PMCID: PMC64507.

16. Abel G, Lyratzopoulos G. Ranking hospitals on avoidable death rates derived from retrospective case record review: methodological observations and limitations. BMJ Qual Saf. 2015;24(9):554–7. Epub 2015/07/05. pmid:26141503; PubMed Central PMCID: PMC4552920.

17. Abel G, Elliott MN. Identifying and quantifying variation between healthcare organisations and geographical regions: using mixed-effects models. BMJ Qual Saf. 2019. Epub 2019/09/20. pmid:31533954.

18. Merlo J, Wagner P, Ghith N, Leckie G. An Original Stepwise Multilevel Logistic Regression Analysis of Discriminatory Accuracy: The Case of Neighbourhoods and Health. PloS one. 2016;11(4):e0153778. pmid:27120054

19. Hakkinen U, Rosenqvist G, Iversen T, Rehnberg C, Seppala TT, Euro Hsg. Outcome, Use of Resources and Their Relationship in the Treatment of AMI, Stroke and Hip Fracture at European Hospitals. Health Economics. 2015;24 Suppl 2:116–39. pmid:26633872

20. Paula FdL Cunha GMd, Leite IdC Pinheiro RS, Valente JG. Readmission of older patients after hospital discharge for hip fracture: a multilevel approach. Revista de Saúde Pública. 2016;50. pmid:27143616

21. Hekkert K, Kool RB, Rake E, Cihangir S, Borghans I, Atsma F, et al. To what degree can variations in readmission rates be explained on the level of the hospital? a multilevel study using a large Dutch database. BMC Health Serv Res. 2018;18(1):999. Epub 2018/12/29. pmid:30591058; PubMed Central PMCID: PMC6307249.

22. van Dishoeck AM, Lingsma HF, Mackenbach JP, Steyerberg EW. Random variation and rankability of hospitals using outcome indicators. BMJ Qual Saf. 2011;20(10):869–74. Epub 2011/06/07. pmid:21642443.

23. Tomioka S, Fujino Y, Nakaya T, Ohtani M, Kubo T, Matsuda S. Equality of Treatment for Hip Fracture Irrespective of Regional Differences in Socioeconomic Status: Analysis of Nationwide Hospital Claims Data in Japan. The Tohoku Journal of Experimental Medicine. 2019;247(3):161–71. pmid:30867369

24. Merlo J, Wagner P, Leckie G. A simple multilevel approach for analysing geographical inequalities in public health reports: The case of municipality differences in obesity. Health Place. 2019;58:102145. Epub 2019/06/14. pmid:31195211.

25. Leckie G. Cross-Classified Multilevel Models—Concepts. LEMMA VLE Module 12, 1–60. 2013 [cited 2020 04.02.2020]. Available from: http://www.bristol.ac.uk/cmm/learning/online-course/.

26. Wettermark B, Hammar N, Fored CM, Leimanis A, Otterblad Olausson P, Bergman U, et al. The new Swedish Prescribed Drug Register—opportunities for pharmacoepidemiological research and experience from the first six months. Pharmacoepidemiol Drug Saf. 2007;16(7):726–35. pmid:16897791.

27. Brooke HL, Talback M, Hornblad J, Johansson LA, Ludvigsson JF, Druid H, et al. The Swedish cause of death register. Eur J Epidemiol. 2017;32(9):765–73. Epub 2017/10/07. pmid:28983736; PubMed Central PMCID: PMC5662659.

28. Ludvigsson JF, Almqvist C, Bonamy AK, Ljung R, Michaelsson K, Neovius M, et al. Registers of the Swedish total population and their use in medical research. Eur J Epidemiol. 2016;31(2):125–36. Epub 2016/01/16. pmid:26769609.

29. Public Access to Information Secrecy Act. In: Justice Mo, editor. Stockholm2009.

30. Hu F, Jiang C, Shen J, Tang P, Wang Y. Preoperative predictors for mortality following hip fracture surgery: a systematic review and meta-analysis. Injury. 2012;43(6):676–85. pmid:21683355

31. Schousboe JT, Paudel ML, Taylor BC, Kats AM, Virnig BA, Dowd BE, et al. Pre-fracture individual characteristics associated with high total health care costs after hip fracture. Osteoporosis international: a journal established as result of cooperation between the European Foundation for Osteoporosis and the National Osteoporosis Foundation of the USA. 2016. pmid:27743069

32. Jantzen C, Madsen CM, Abrahamsen B, Van Der Mark S, Duus BR, Howland J, et al. Pre-fracture medication use as a predictor of 30-day mortality in hip fracture patients: an analysis of 141,201 patients. Hip Int. 2020;30(1):101–6. Epub 2019/03/02. pmid:30821177.

33. Bergman J, Nordström A, Hommel A, Kivipelto M, Nordström P. Bisphosphonates and mortality: confounding in observational studies? Osteoporosis International. 2019;30(10):1973–82. pmid:31367949

34. Pepe MS, Janes H, Longton G, Leisenring W, Newcomb P. Limitations of the odds ratio in gauging the performance of a diagnostic, prognostic, or screening marker. American Journal of Epidemiology. 2004;159(9):882–90. pmid:15105181

35. Hosmer DW, Lemeshow S, Sturdivant RX. Applied logistic regression. Third edition / ed. Hoboken, New Jersey: Wiley; 2013. xvi, 500 pages p.

36. Merlo J, Ostergren PO, Broms K, Bjorck-Linne A, Liedholm H. Survival after initial hospitalisation for heart failure: a multilevel analysis of patients in Swedish acute care hospitals. J Epidemiol Community Health. 2001;55(5):323–9. Epub 2001/04/12. pmid:11297650; PubMed Central PMCID: PMC1731888.

37. Ghith N, Merlo J, Frolich A. Albuminuria measurement in diabetic care: a multilevel analysis measuring the influence of accreditation on institutional performance. BMJ Open Qual. 2019;8(1):e000449. Epub 2019/02/08. pmid:30729192; PubMed Central PMCID: PMC6340563.

38. Ghith N, Frolich A, Merlo J. The role of the clinical departments for understanding patient heterogeneity in one-year mortality after a diagnosis of heart failure: A multilevel analysis of individual heterogeneity for profiling provider outcomes. PLoS One. 2017;12(12):e0189050. pmid:29211785; PubMed Central PMCID: PMC5718563.

39. Hjerpe P, Ohlsson H, Lindblad U, Bostrom KB, Merlo J. Understanding adherence to therapeutic guidelines: a multilevel analysis of statin prescription in the Skaraborg Primary Care Database. Eur J Clin Pharmacol. 2010. Epub 2010/12/31. pmid:21190018.

40. Ohlsson H, Librero J, Sundquist J, Sundquist K, Merlo J. Performance evaluations and league tables: do they capture variation between organizational units? An analysis of 5 Swedish pharmacological performance indicators. Med Care. 2011;49(3):327–31. Epub 2011/01/26. pmid:21263360.

41. Merlo J, Leckie G, Wagner P. A simple multilevel approach for evaluating geographical differences in public health reports: the case of obesity in the municipalities of Skåne, Sweden. Health & Place (in press). 2019.

42. Pedersen AB, Ehrenstein V, Szepligeti SK, Lunde A, Lagerros YT, Westerlund A, et al. Thirty-five-year Trends in First-time Hospitalization for Hip Fracture, 1-year Mortality, and the Prognostic Impact of Comorbidity: A Danish Nationwide Cohort Study, 1980–2014. Epidemiology. 2017;28(6):898–905. Epub 2017/08/03. pmid:28767515.

43. Medin E, Goude F, Melberg HO, Tediosi F, Belicza E, Peltola M. European Regional Differences in All-Cause Mortality and Length of Stay for Patients with Hip Fracture. Health Econ. 2015;24 Suppl 2:53–64. Epub 2015/12/04. pmid:26633868.

44. Turesson E, Ivarsson K, Thorngren KG, Hommel A. The impact of care process development and comorbidity on time to surgery, mortality rate and functional outcome for hip fracture patients: a retrospective analysis over 19 years with data from the Swedish National Registry for hip fracture patients, RIKSHOFT. BMC Musculoskelet Disord. 2019;20(1):616. Epub 2019/12/28. pmid:31878903; PubMed Central PMCID: PMC6933741.

45. Black DM, Cummings SR, Karpf DB, Cauley JA, Thompson DE, Nevitt MC, et al. Randomised trial of effect of alendronate on risk of fracture in women with existing vertebral fractures. Fracture Intervention Trial Research Group. Lancet (London, England). 1996;348(9041):1535–41. Epub 1996/12/07. pmid:8950879.

46. McClung MR, Geusens P, Miller PD, Zippel H, Bensen WG, Roux C, et al. Effect of risedronate on the risk of hip fracture in elderly women. Hip Intervention Program Study Group. N Engl J Med. 2001;344(5):333–40. Epub 2001/02/15. pmid:11172164.

47. Lyles KW, Colon-Emeric CS, Magaziner JS, Adachi JD, Pieper CF, Mautalen C, et al. Zoledronic acid and clinical fractures and mortality after hip fracture. The New England journal of medicine. 2007;357(18):1799–809. pmid:17878149

48. Peng J, Liu Y, Chen L, Peng K, Xu Z, Zhang D, et al. Bisphosphonates can prevent recurrent hip fracture and reduce the mortality in osteoporotic patient with hip fracture: A meta-analysis. Pakistan journal of medical sciences. 2016;32(2):499–504. Epub 2016/05/18. pmid:27182270; PubMed Central PMCID: PMC4859053.

49. Rogmark C. Further refinement of surgery will not necessarily improve outcome after hip fracture. Acta Orthop. 2020:1–3. Epub 2020/01/09. pmid:31913741.