- Open Access
Changes in hospital mortality for United States intensive care unit admissions from 1988 to 2012
Critical Carevolume 17, Article number: R81 (2013)
A decrease in disease-specific mortality over the last twenty years has been reported for patients admitted to United States (US) hospitals, but data for intensive care patients are lacking. The aim of this study was to describe changes in hospital mortality and case-mix using clinical data for patients admitted to multiple US ICUs over the last 24 years.
We carried out a retrospective time series analysis of hospital mortality using clinical data collected from 1988 to 2012. We also examined the impact of ICU admission diagnosis and other clinical characteristics on mortality over time. The potential impact of hospital discharge destination on mortality was also assessed using data from 2001 to 2012.
For 482,601 ICU admissions there was a 35% relative decrease in mortality from 1988 to 2012 despite an increase in age and severity of illness. This decrease varied greatly by diagnosis. Mortality fell by > 60% for patients with chronic obstructive pulmonary disease, seizures and surgery for aortic dissection and subarachnoid hemorrhage. Mortality fell by 51% to 59% for six diagnoses, 41% to 50% for seven diagnoses, and 10% to 40% for seven diagnoses. The decrease in mortality from 2001 to 2012 was accompanied by an increase in discharge to post-acute care facilities and a decrease in discharge to home.
Hospital mortality for patients admitted to US ICUs has decreased significantly over the past two decades despite an increase in the severity of illness. Decreases in mortality were diagnosis specific and appear attributable to improvements in the quality of care, but changes in discharge destination and other confounders may also be responsible.
Over the last two decades there has been a decrease in mortality for patients admitted to US hospitals. Among Medicare patients who may or may not have been admitted to an ICU, risk adjusted mortality fell by 18% to 46% between 1994 and 2007 [1, 2]. Disease specific decreases in the 30-day mortality rate over time have also been reported for patients hospitalized for congestive heart failure (50%, 1993 to 2008) , community acquired pneumonia (28%, 1987 to 2005)  and surgery for subarachnoid hemorrhage (50%, 1980 to 2005), coronary artery bypass graft (21%, 1999 to 2008) and other high-risk procedures . Changes in hospital mortality have also been reported for US patients with disorders commonly managed in ICUs including sepsis (9.9%, 1979 to 2000) , acute lung injury (9%, 1996 to 2005)  and surgery for aortic dissection (17%, 1979 to 2003) .
To our knowledge temporal changes in hospital mortality have not been reported for patients admitted to US ICUs. In contrast, multi-institutional data have shown mortality reductions over time in Australia and New Zealand (4%, 1993 to 2003)  and England (13.4%, 1998 to 2006) . Studies focused on ICU prognostic systems have also reported decreases in hospital mortality over the last two decades. These studies have repeatedly demonstrated 'model fade': the over prediction of mortality when risk adjustment models are applied to more recent data . As a result, the Simplified Acute Physiology Score (SAPS) , Mortality Probability Model (MPM) , and Acute Physiology and Chronic Health Evaluation (APACHE)  have required repeated updating. The developers of these systems have attributed mortality over prediction to increased treatment effectiveness, improved care before ICU admission, and more frequent discharge to post-acute care facilities [13–15].
Patients who are admitted to ICUs and survive hospitalization have a high mortality rate in the six months after discharge . Many of these post-discharge deaths are among patients transferred to other acute-care hospitals  or long-term acute care facilities [18, 19]. These findings have raised concerns that hospital discharge practices may partially account for reductions in hospital mortality and adversely impact performance assessments that use predicted hospital mortality [17, 20].
This study reports a time series analysis of hospital mortality among patients admitted to multiple US ICUs over a 24-year period. The objectives of the study are: to describe changes in mortality and case-mix among ICU admissions from 1988 to 2012; to describe mortality changes among ICU patients with specific diagnoses over the same period; and to present data about the potential impact of patient discharge destination on hospital mortality. Data comparisons used a Chi-square test for categorical measures and an analysis-of-variance test for continuous measures.
Materials and methods
We performed a retrospective analysis using data obtained from ICUs that had installed an APACHE system. Data for patients admitted before 2001 were obtained from published studies of APACHE. This emanated from two time periods: patient data for 1988 to 1989 were obtained from institutions randomly selected to represent US hospitals  and for 1993 to 1996 from a nonrandomized cohort of US institutions . For admissions from 2001 onwards data were obtained for individual patients in the electronic APACHE database and grouped into three-year intervals. Multiple Institutional Review Boards (George Washington University, University of Maryland, Baystate Medical Center) waived the need for informed consent, deeming the use of these databases to not be human subject research requiring approval [15, 21–24].
Patient data were generated as a result of patient care and collected for consecutive unselected ICU admissions. Demographic and clinical data were recorded on ICU day one, and vital status at ICU and hospital discharge. We did not collect data for patients with burns, ICU stay <4 hours, age <16 years and excluded second or subsequent ICU admissions. Further details about data collection, variables and the reliability of APACHE data have been described elsewhere [15, 21, 22].
Changes in mortality over time
We examined changes in hospital mortality for all eligible first ICU admissions from May 1988 to June 2012. Concurrently we tracked changes in mean age and the day 1 acute physiology score (APS) of APACHE III. We repeated these analyses for patients who were mechanically ventilated during ICU day 1. We did not compare standardized mortality ratios for each time interval due to the numerous updates made to the APACHE mortality prediction over this 24- year interval [15, 22].
To account for differences in diagnostic coding during the study, we examined only diagnoses with a consistent definition across all time intervals. We excluded diagnoses with <10% mortality during 1988 to 1989 to eliminate diagnoses that might have a small absolute change in mortality, which appears as a large relative mortality difference. These criteria produced 22 diagnoses. Using mortality rates in 1988 to 1989 as the baseline, we calculated the percent change in mortality in 2010 to 2012 for the 22 consistently defined diagnoses, and grouped them according to the percent decline in mortality: >60%, 51% to 59%, 41% to 50%, and 10% to 40%. Age and severity of illness (APS) data were available in aggregate but not for individual patients in the 1988 to 1989 and 1993 to 1996 cohorts were not available, thus we could not track diagnosis-specific changes in these variables.
Changes in discharge destination over time
For each ICU admission the discharge destination was recorded using a computerized pick list; this information was only available for admissions from 2001 to 2012. Discharge destinations included death, home (with or without home health care), another hospital (including long term acute care facilities), nursing home (including assisted living facilities), skilled nursing facility (including inpatient rehabilitation facilities) and other destinations. Discharge to a post-acute care (PAC) facility was defined by discharge to skilled nursing, inpatient rehabilitation or assisted living facilities [24, 25].
To assess the potential impact of changes in discharge destination on hospital mortality we examined the percentage of patients discharged home for each 3-year time interval from 2001 to 2012. We then repeated this analysis for each of the 22 diagnoses described above. To eliminate large changes in discharge destination due to small sample size we excluded diagnoses with <2,500 admissions in 2010 to 2012, the interval with the smallest number of patients.
We did not perform statistical analysis of differences between time periods. Analysis was not performed because significance was assured by the large size of each cohort.
The number of hospitals, ICUs and patients included in each data collection interval from May 1988 to June 2012 are shown in Table 1. The total number of ICU admissions was 482,601 and each time interval included ≥40 hospitals and ≥40 ICUs. With the exception of 1988 to 1989, all intervals contained at least 35,000 first admissions. The characteristics of the hospitals and ICUs are shown in Table 2. In general, hospital bed-size, teaching status, and geographic region were well dispersed, except for a smaller proportion of hospitals in the Northeast. The modal type of ICU was mixed medical-surgical, but most unit types were represented across the study.
The characteristics of patients admitted from 1988 to 2012 are shown in Table 3. For each cohort, patient age and severity of illness (APS) increased over time. There was a declining trend in the percentage of patients who were admitted post-operatively, from floors or by direct admission. Conversely, an increasing percentage of patients were admitted from another hospital or a step down unit. During the overall study period ICU length of stay decreased by 23% and hospital length of stay decreased by 38%.
Changes in patient mortality
Temporal trends from the1988/1989 to the 2010/2012 time interval for aggregate mean hospital mortality, age, and APS for the 482,601 admissions are shown in Figure 1. There was a relative decrease in mean hospital mortality of 35% from 1988/1989 to 2010/2012 (P <0.001). Age increased slightly over this time period, but severity of illness (APS) increased by 7.4% (P <0.001). When stratifying the mortality trends by whether or not a patient was ventilated on ICU day 1, the changes in mortality were parallel between the two strata. These changes are shown as supplementary information [see Additional File 1, Figure S1].
Table 4 shows the 22 ICU admission diagnoses that were consistently defined throughout the study and had > 10% mortality during 1988/1989. There were 14 medical and 8 postoperative diagnoses. For 19 (90.5%) of the 22 diagnoses there was a >20% mortality decline from 1988/1989 to 2010/2012. The remaining three diagnoses had a 0 to 20% decline. Compared to the 1988/1989 baseline there was a > 60% reduction in mortality for four diagnostic groups by 2010/2012 (Figure 2); they included surgery for aortic dissection (84%), surgery for subarachnoid hemorrhage (77%), seizures (72%), and chronic obstructive pulmonary disease (66%). There was a 51% to 59% mortality reduction for six diagnostic groups (Figure 3); they included surgery for intracerebral hemorrhage (58%), sepsis, urinary tract (57%), acute myocardial infarction (57%), surgery for gastrointestinal (GI) malignancy (57%), congestive heart failure (55%) and sepsis, non-urinary tract (51%).
There was a 41% to 50% mortality reduction for seven diagnostic groups from 1988/1989 to 2010/2012 (Figure 4); they included intracerebral hemorrhage (48%), stroke (47%), GI bleeding, upper-non-variceal (47%), viral pneumonia (46%), surgery for multiple-trauma, including the head (41%), surgery for GI perforation (41%) and GI bleeding from varices (41%). There was a 10% to 40% reduction for five diagnostic groups (Figure 5); they were surgery for GI obstruction (32%), acute respiratory distress syndrome (31%), surgery for GI bleeding (20%), cardiac arrest (19%) and head trauma with chest, abdomen, pelvis or spine injury (10%).
Changes in hospital discharge destination
Hospital death rates and the percentage of hospital discharges by destination for 427,533 ICU admissions from 2001 to 2012 are shown in Figure 6. Hospital death rates and the proportion of home discharges decreased over time. Conversely, the proportion of discharges to other hospitals and PAC facilities increased.
There were six diagnoses with sufficient sample size for analysis (>2,500 admissions in 2010 to 2012): cardiac arrest, congestive heart failure, GI bleeding (upper, non-variceal) stroke, acute myocardial infarction and sepsis (urinary tract). Changes in hospital discharge destination from 2001 to 2012 for patients with four of these diagnoses included an increased frequency of discharge to PAC facilities accompanied by a decreased hospital mortality rate and discharge to home. The exceptions were acute myocardial infarction (increase in discharge to home) and congestive heart failure (no change in discharge to home). These diagnosis specific results are shown graphically as supplementary information [see Additional File 1, Figures S2 to S7].
In a clinical database from multiple US ICUs, we found that mortality at hospital discharge had a relative decrease of 35% (17.3% to 12.4%) from 1988/1989 to 2010/2012. This mortality reduction was accompanied by an increase in severity of illness and number of chronic health conditions, factors that typically pose an increase in mortality risk. There was also a >50% mortality reduction across ten diagnostic categories over the same 24-year interval. Most of these dramatic relative decreases in hospital mortality rate occurred between 1988 to 1989 and 1993 to 1996. To the best of our knowledge, this is the first longitudinal analysis of hospital mortality rates in a clinical (non-administrative) database for a large number of admissions to multiple US ICUs.
The decrease in hospital mortality was not uniform across diagnostic categories. This may be because medical advances are typically related to the development of new therapies for specific diseases. The dramatic reductions in our study are consistent with prior reports of improved mortality for patients with stroke [1, 2], acute myocardial infarction [1, 2], congestive heart failure [1, 3], chronic obstructive pulmonary disease  and sepsis . Dramatic mortality reductions have also been reported after surgery for aortic dissection  and subarachnoid hemorrhage . It is not our intention, however, to compare these disease specific mortality reductions with prior reports or relate them to new therapies.
Compared to reports from ICUs in other countries the changes in demographics and severity of illness were similar but reductions in hospital mortality were larger for these US ICU patients [10,11). Unfortunately, documented changes in the accuracy of APACHE mortality predictions over time [15, 22] precluded evaluation of the impact of case-mix on changes in observed versus predicted mortality rates between 1988 and 2012. Our study's retrospective design, differences in patient selection across time periods, lack of information about ICU structure, process and use of disease-specific treatments, also make it impossible to assess the reasons for and timing of these mortality decreases. We are, therefore, unable to comment on the impact of differing ICU types, staffing by critical care specialists, ventilator and sepsis bundles, and changes in disease-specific therapy on mortality.
We are unable to account for the impact of Medicare's prospective payment system in reducing mortality at hospital discharge. This legislation created a substantial shift in the time US patients spent in acute care hospitals to PAC facilities between 1984 and 2002 [25, 29]. These changes in discharge destination have previously been reported to shift hospital deaths to PACs [17, 20], long-term acute care facilities [17, 18] or other acute care hospitals [17, 19]. Because we lacked data about discharge destination before 2001 we are unable to assess the impact of discharge destination on the dramatic fall in mortality between 1988 to 1989 and 1993 to 1996. Reductions in mortality might also be related to the 1998 Joint Commission performance measurement requirements for patients with acute myocardial infarction, heart failure, and pneumonia and subsequent public reporting in 2004 .
The modest reduction in hospital mortality (10.5%) from 2001 to 2012 may have been influenced, at least in part, by the increasing frequency of discharge to PACs and other hospitals (Figure 6). It is also possible, however, that severely ill patients discharged home in 2001 to 2003 might be going to PACs in 2010 to 2012 and, thus, not confounding the reduction in mortality. This emphasizes the need to obtain long-term survival data to properly assess temporal changes in mortality for patients surviving to hospital discharge [16–18].
The reduction in hospital mortality over this 24-year time interval explains why APACHE III (developed using 1988 to 1989 data) and APACHE IV (developed using 2002 to 2003 data) and other prognostic systems over predict mortality for contemporary ICU admissions. The over prediction of mortality is even greater for patients with specific diagnoses that had the largest mortality reductions between 1988/1989 and 2010/2012. The magnitude of these disease-specific changes may explain why recalibrating models that do not include specific diagnostic information often fails to improve calibration when tested in new populations [31–33].
Our study is subject to several limitations. First, our data are not necessarily representative of all US ICUs because they were obtained only from units with an APACHE system. However, there was great diversity in hospital size, teaching status, region, as well as ICU type. Second, our data are subject to a selection bias because data from 1988 to 1989 and 1993 to 1996 were population based, whereas data from 2001 onwards was from a self-selected sample of ICUs. In addition, the data for this study were not serially collected from the same hospitals over the 24-year interval, and we cannot exclude the possibility that changes in outcome were related to ICU based differences in quality of care. Third, we lacked individual patient data for 1988 to 1989 and 1993 to 1996 admissions. For specific diagnostic groups, this hampered our ability to evaluate the impact of changes in severity of illness and other confounders on mortality. Finally, we did not have information about mortality after hospital discharge. For this reason we were unable to evaluate accurately the impact of changes in discharge destination on hospital mortality.
Hospital mortality for patients admitted to intensive care units has decreased significantly over the past quarter century despite an increase in the severity of patient illness. Decreases in mortality are diagnosis-specific and might be, in part, attributable to improvements in quality of care, but changing patterns in discharge destination and other confounders may also be responsible.
Hospital mortality for patients admitted to US intensive care units had a relative decrease of 35% between 1988 and 2012.
Reductions in hospital mortality were greater and more varied for patients within specific diagnostic groups.
Changes in hospital mortality over time explain why prognostic scoring systems over predict mortality when applied to more contemporary patient data.
Long-term outcome data are needed to assess the impact of discharge to post-acute care facilities on hospital mortality rates.
Acute Physiology and Chronic Health Evaluation
Acute Physiology Score
Mortality Probability Model
Simplified Acute Physiology Score.
Andrews RM, Russo CA, Pancholi M: Trends in hospital risk-adjusted mortality for select diagnoses and procedures, 1994-2004. In HCUP Statistical Brief #38. Agency for Healthcare Research and Quality; Rockville, MD; 2007:1-7.
Hines A, Stranges E, Andrews R: Trends in hospital risk-adjusted mortality for select diagnoses by patient subgroups, 2000-2007. In HCUP Statistical Brief #98. Agency for Healthcare Research and Quality, Rockville, MD; 2010:1-11.
Bueno H, Ross JS, Wang Y, Chen J, Vidan MT, Normand SLT, Curtis JP, Drye EE, Lichtman JH, Keenan PS, Kosiborod M, Krumholz HM: Trends in length of stay and short-term outcomes among Medicare patients hospitalized for heart failure: 1993-2008. JAMA 2010, 303: 2141-2147. 10.1001/jama.2010.748
Ruhnke GW, Coca-Perraillon M, Kitch BT, Cutler DM: Marked reduction in 30-day mortality among elderly patients with community-acquired pneumonia. Am J Med 2011, 124: 171-178. 10.1016/j.amjmed.2010.08.019
Lovelock CE, Rinkel GJE, Rothwell PM: Time trends in outcome of subarachnoid hemorrhage. Neurology 2010, 74: 1494-1501. 10.1212/WNL.0b013e3181dd42b3
Finks JF, Osborne NH, Birkmeyer JD: Trends in hospital volume and operative mortality for high-risk surgery. N Engl J Med 2011, 364: 2128-2837. 10.1056/NEJMsa1010705
Martin GS, Mannino DM, Eaton S, Moss M: The epidemiology of sepsis in the United States from 1979 through 2000. N Engl J Med 2003, 348: 1546-1554. 10.1056/NEJMoa022139
Erickson SE, Martin GS, Davis JL, Matthay MA, Eisner MD: Recent trends in acute lung injury mortality: 1996-2005. Crit Care Med 2009, 37: 1574-1579. 10.1097/CCM.0b013e31819fefdf
Stevens LM, Madsen JC, Isselbacher EM, Khairy P, MacGillivray TE, Hilgenberg AD, Agnihotri AK: Surgical management and long-term outcomes for acute ascending aortic dissection. J Thorac Cardiovasc Surg 2009, 138: 1349-1357. 10.1016/j.jtcvs.2009.01.030
Moran JL, Bristow P, Solomon PJ, George C, Hart GK: Mortality and length-of-stay outcomes, 1993-2003, in the binational Australian and New Zealand intensive care adult patient database. Crit Care Med 2008, 36: 46-61. 10.1097/01.CCM.0000295313.08084.58
Hutchings A, Durand MA, Grieve R, Harrison D, Rowan K, Green J, Carins J, Black N: Evaluation of modernization of adult critical care services in England: time series and cost effectiveness analysis. BMJ 2009, 339: b4353. 10.1136/bmj.b4353
Kramer AA: Predictive mortality models are not like fine wine. Crit Care 2005, 9: 636-637. 10.1186/cc3899
Moreno RP, Metnitz PGH, Almeida E, Jordan B, Bauer P, Campos RA, Iapichino G, Edbrooke D, Capuzzo M, Le Gall JR: SAPS 3 - from evaluation of the patient to evaluation of the intensive care unit. Part 2: development of a prognostic model for hospital mortality at ICU admission. Intensive Care Med 2005, 31: 1345-1355. 10.1007/s00134-005-2763-5
Higgins TL, Teres D, Copes WS, Nathanson BH, Stark M, Kramer AA: Assessing contemporary intensive care unit outcome: an updated Mortality Probability Admission Model (MPM O -III). Crit Care Med 2007, 35: 827-835. 10.1097/01.CCM.0000257337.63529.9F
Zimmerman JE, Kramer AA, McNair DS, Malila FM: Acute Physiology and Chronic Health Evaluation (APACHE) IV: hospital mortality assessment for today's critically ill patients. Crit Care Med 2006, 34: 1297-1310. 10.1097/01.CCM.0000215112.84523.F0
Wunsch H, Guerra C, Barnato AE, Angus DC, Li G, Linde-Zwirble WT: Three-year outcomes for Medicare beneficiaries who survive intensive care. JAMA 2010, 303: 849-856. 10.1001/jama.2010.216
Vasilevskis EE, Kuzniewicz MW, Dean ML, Clay T, Vittinghoff E, Rennie DJ, Dudley RA: Relationship between discharge practices and intensive care unit in-hospital mortality performance: evidence of a discharge bias. Med Care 2009, 47: 803-812. 10.1097/MLR.0b013e3181a39454
Hall WB, Willis LE, Medvedev S, Carson SS: The implications of long-term acute care hospital transfer practices for measures of in-hospital mortality and length of stay. Am J Respir Crit Care Med 2012, 185: 53-57. 10.1164/rccm.201106-1084OC
Kelly KB, Koeppel ML, Como JJ, Carter JW, McCoy AM, Claridge JA: Continued rationale of why hospital mortality is not an appropriate measure of trauma outcomes. Am J Surg 2012, 203: 366-369. 10.1016/j.amjsurg.2011.10.004
Kahn JM, Kramer AK, Rubenfeld GD: Transferring critically ill patients out of hospital improves the standardized mortality ratio. Chest 2007, 131: 68-75. 10.1378/chest.06-0741
Knaus WA, Wagner DP, Draper EA, Zimmerman JE, Bergner M, Bastos PG, Sirio CA, Murphy DJ, Lotring T, Damiano A, Harrell FE: The APACHE III prognostic system: risk prediction of hospital mortality for critically ill hospitalized adults. Chest 1991, 100: 1619-1636. 10.1378/chest.100.6.1619
Zimmerman JE, Wagner DP, Draper EA, Wright L, Alzola C, Knaus WA: Evaluation of Acute Physiology and Chronic Health Evaluation III predictions of hospital mortality in an independent database. Crit Care Med 1998, 26: 1317-1326. 10.1097/00003246-199808000-00012
Kramer AA, Higgins TL, Zimmerman JE: Intensive care unit readmissions in U.S. hospitals: patient characteristics, risk factors, and outcomes. Crit Care Med 2012, 40: 3-10. 10.1097/CCM.0b013e31822d751e
Kramer AA, Zimmerman JE: Institutional variations in frequency of discharge of elderly intensive care survivors to post-acute care facilities. Crit Care Med 2010, 38: 2319-2328. 10.1097/CCM.0b013e3181fa02e4
Buntin MB: Access to postacute rehabilitation. Arch Phys Med Rehabil 2007, 88: 1488-1493. 10.1016/j.apmr.2007.07.023
Ford ES, Mannino DM, Zhao G, Li C, Croft JB: Changes in mortality among US adults with COPD in two national cohorts recruited from 1971-1975 and 1988-1994. Chest 2012, 141: 101-110. 10.1378/chest.11-0472
Kumar G, Kumar N, Taneja A, Kaleekal T, Tarima S, McGinley E, Jimenez E, Mohan A, Kahn RA, Whittle J, Jacobs E, Nanchal R: Nationwide trends of severe sepsis in the 21stcentury (2000-2007). Chest 2011, 140: 1223-1231. 10.1378/chest.11-0352
Qureshi AI, Suri MF, Nasar A, Kirmani JF, Divani AA, He W, Hopkins LN: Trends in hospitalization and mortality for subarachnoid hemorrhage and unruptured aneurysms in the United States. Neurosurgery 2005, 57: 1-8.
Chan L, Koepsell TD, Devo RA, Esselman PC, Haselkorn JK, Lowery JK, Stolov WC: The effect of Medicare's payment system for rehabilitation hospitals on length of stay, charges, and total payments. N Engl J Med 1997, 337: 978-985. 10.1056/NEJM199710023371406
Chassin MR, Loeb JM, Schmaltz SP, Wachter RM: Accountability measures - using measurement to promote quality improvement. N Engl J Med 2010, 363: 683-688. 10.1056/NEJMsb1002320
Keegan MT, Gajic O, Afessa B: Severity of illness scoring systems in the intensive care unit. Crit Care Med 2011, 39: 163-169. 10.1097/CCM.0b013e3181f96f81
Breslow MJ, Badawi O: Severity scoring in the critically ill: part 1 - interpretation and accuracy of outcome prediction scoring systems. Chest 2012, 141: 245-252. 10.1378/chest.11-0330
Kuzniewicz MW, Vasilevskis EE, Lane R, Dean ML, Trivedi NG, Rennie DJ, Clay T, Kotler PL, Dudley RA: Variation in ICU risk-adjusted mortality: impact of methods of assessment and potential cofounders. Chest 2008, 133: 1319-1327. 10.1378/chest.07-3061
JEZ is a paid consultant and receives research support from Cerner Corporation. AAK is a full-time employee and stockholder in Cerner Corporation. WAK declares that he has no competing interests. Cerner Corporation, Kansas City, MO, supported the study. Cerner markets the APACHE Outcomes Clinical Information System and owns databases used for this study.
All authors took part in the planning, design, interpretation and drafting of the study. AAK acquired the data and performed the analyses. All authors have read and approved the manuscript for publication.