Risk factors for the development of nosocomial pneumonia and mortality on intensive care units: application of competing risks models
Critical Care volume 12, Article number: R44 (2008)
Pneumonia is a very common nosocomial infection in intensive care units (ICUs). Many studies have investigated risk factors for the development of infection and its consequences. However, the evaluation in most of theses studies disregards the fact that there are additional competing events, such as discharge or death.
A prospective cohort study was conducted over 18 months in five intensive care units at one university hospital. All patients that were admitted for at least 2 days were included, and surveillance of nosocomial pneumonia was conducted. Various potential risk factors (baseline- and time-dependent) were evaluated in two competing risks models: the acquisition of nosocomial pneumonia and discharge (dead or alive; model 1) and for the risk of death in the ICU and discharge alive (model 2).
Patients from 1,876 admissions were included. A total of 158 patients developed nosocomial pneumonia. The main risk factors for nosocomial pneumonia in the multivariate analysis in model 1 were: elective surgery (cause-specific hazard ratio = 1.95; 95% CI 1.33 to 2.85) or emergency surgery (1.59; 95% CI 1.10 to 2.28) prior to ICU admission, usage of a nasogastric tube (3.04; 95% CI 1.25 to 7.37) and mechanical ventilation (5.90; 95% CI 2.47 to 14.09). Nosocomial pneumonia prolonged the length of ICU stay but was not directly associated with a fatal outcome (p = 0.55).
More studies using competing risk models, which provide more accurate data compared to naive survival curves or logistic models, should be carried out to verify the impact of risk factors and patient characteristics for the acquisition of nosocomial infections and infection-associated mortality.
Nosocomial pneumonia (NP) is the most commonly reported infection in intensive care units (ICUs), especially in mechanically ventilated patients with an incidence of about 15 infections per 1,000 ventilation days . This infection is associated with a significantly increased length of hospital stay and may have a considerable impact on morbidity and mortality .
Endpoints, possible risk factors for the acquisition of NP and the clinical outcome after the infection has occurred have been addressed in numerous studies. However, many of these studies did not take into account the fact that there are other possible endpoints competing with the event of interest [3, 4]. For example 'death' or 'discharge' are competing events for the onset of infection. A competing risks methodology allows for a better understanding of why NP increases mortality. Unlike logistic regression, it allows modelling of the time-dependency of certain procedures (for example intubation), thereby avoiding biased results. For this, multi-state models are a more accurate approach in order to consider competing events [5, 6]. We present here the results of a competing risks analysis to address two major objectives: (1) to identify potential risk factors for NP in ICUs, considering discharge (dead or alive without prior NP) as the competing event, and (2) to investigate several risk factors, including blood stream infection, NP and other lower respiratory tract infections as time-dependent risks, for mortality in ICU patients with discharge (alive) as the competing endpoint.
Materials and methods
Patients and infections
The presenr study was conducted in five ICUs (one medical, one surgical, one neurosurgical and two interdisciplinary) at one German university hospital from February 2000 to July 2001 (a total study period of 18 months). All patients with a duration of ICU stay of at least 2 days were enrolled. Prospective surveillance of nosocomial infections was performed by trained staff of the German Nosocomial Infection Surveillance System (KISS)  using the standardized US Centers for Disease Control and Prevention (CDC) definitions for NP . The method of surveillance remained unchanged over the study period. As all investigations represented routine diagnostic procedures, the Institutional Board on the Ethics of Clinical Studies waived the need for informed consent. Further details on the setting of the study are described elsewhere [9, 10].
Analysis of risk factors for the acquisition of NP (model 1)
In model 1, we studied risk factors for NP acquisition as well as the competing risk 'discharge (dead or alive without prior NP)' (Figure 1). After admission to the ICU (event 0) the patient may (event 1) or may not (event 2) acquire NP. The impact of the following baseline risk factors were investigated: age, gender, simplified acute physiology score (SAPS) II, intubation at ICU admission, infection present already at the time point of ICU admission (pneumonia, urinary tract infection and other infections), hospitalization prior to ICU admission, elective or emergency surgery before ICU admission (for example, head trauma, multiple trauma, vascular surgery and neurosurgery), underlying diseases (cardial/pulmonal, gastrointestinal, neurological, and metabolic/renal) and other underlying diseases (including sepsis, malignancies or alcoholism). The impact of the following time-dependent risk factors were investigated as time-dependent covariates (which start with value = 0 and may increase to 1): ventilation, chest drainage, colostomy, enterostomy, jejunostomy, nasogastric tube and urinary catheter. Age and SAPS II score were included in the model as continuous variables; all other factors were binary variables only.
Analysis of risk factors for mortality (model 2)
In model 2 we studied competing risks for mortality and discharge (Figure 1). After admission to the ICU (event 0) the patient may either die during their ICU stay (event 1) or be discharged from the ICU (event 2). Here, we are mainly interested in NP as a time-dependent risk factor for death in the ICU. The same baseline and time-dependent risk factors as described for model 1 were also applied in model 2. We also checked for lower respiratory tract (LRT) infections other than pneumonia on admission as baseline, and for nosocomial LRT and nosocomial blood stream infection as time-dependent variables.
For both models 1 and 2 a competing risk analysis was performed using cause-specific hazards [11, 12]. This analysis follows separate Cox models for each event assuming proportional hazards. In such competing risks analyses with two endpoints, it is possible to interpret both cause-specific hazard ratios (CSHRs) simultaneously for each risk factor. Cumulative incidence functions have been displayed for each endpoint. The proportional hazard assumptions were assessed by study of the graphs of the Schoenfeld's residuals; this technique is especially suitable for time-dependent covariates . The correlation matrices of each Cox model were considered in order to check whether there are correlations among the risk factors, respectively. Risk factors with a p value ≤ 0.157 for at least one of the CSHRs from the univariate analysis were included in a consecutive multivariate analysis. This benchmark corresponds to the well established Akaike information criterion for model selection . A p value ≤ 0.05 was considered statistically significant. For all analyses the R 2.4.1 software was used (R Foundation, Vienna, Austria), especially the R functions coxph, cuminc and cox.zph, from the survival and cmprsk libraries.
Additional data file 1 contains information on the required data format and SAS and R calculations for cause-specific hazard ratios in a competing risks analysis with time-dependent covariates represented.
Patients and infections
A total of 7,269 patients were admitted to the ICUs (35,817 patient days) during the study period; of those, 1,876 admissions (28,498 patient days) required treatment of ≥ 48 h. Only those patients were included in this study. In all, 158 (8.4%) of the included patients developed NP; 132 of these (83.5% of all NP) were ventilator-associated NP (incidence of 8.5 per 1,000 ventilator days) and 33 of these (20.9% of all NP cases) died in the ICU. Overall, in 214 of the 1,876 admissions (11.4%) the patient died in the ICU. More details of risk factors and outcomes are shown in Table 1.
Analysis of risk factors for the acquisition of nosocomial pneumonia (model 1)
Detailed information on the CSHRs of baseline and time-dependent risk factors of model 1 are shown in Table 2. According to this model, significant risk factors for the acquisition of NP in our patient population were (1) pneumonia at admission (CSHR = 0.02), whereas this risk factor also had a reducing effect on the competing event discharge (CSHR = 0.66), (2) undergoing elective surgery prior to ICU admission (CSHR = 1.95), and this effect was accentuated since the CSHR was reduced for discharge (CSHR = 0.54), (3) undergoing emergency surgery prior to ICU admission (CSHR = 1.59), with no significant effect on discharge (CSHR = 1.08), (4) use of a nasogastric tube (CSHR = 3.04), without effect on discharge (CSHR = 0.89), and (5) mechanical ventilation of the patient (CSHR = 5.90), which also significantly reduced the CSHR for discharge from the ICU (CSHR = 0.53; 95% CI 0.45 to 0.62).
In addition to the analysis of model 1, we considered a model with three competing events: nosocomial pneumonia, discharge (alive) and death in the ICU. The CSHRs for pneumonia are the same as in model 1 with the combined competing event. However, the following risk factors had an opposite influence on discharge (alive) and death in the ICU: SAPS II, other infections on admission, surgical patients, metabolic/renal underlying disease and other underlying diseases. This is in line with the results for model 2.
Cumulative incidence functions (CIF)(model 1)
In addition to CSHR, cumulative incidence functions are suitable to illustrate the results of a competing risk analysis. This was exemplarily performed for the risk factors of elective surgery and pneumonia on admission. The CIF of pneumonia starts to increase at an earlier time point for patients with elective surgery, but later for the competing endpoint death/discharge (Figure 2a).
There is only a very low cause-specific risk to acquire nosocomial pneumonia if the patient already had pneumonia on admission (Figure 2b). Regarding discharge (dead or alive) as the endpoint, the cumulative incidence function of the patient group with pneumonia on admission is below the function of the group without until about 40 days in the ICU, but above afterwards.
Analysis of risk factors for mortality (model 2)
Detailed information on the CSHRs of baseline and time-dependent risk factors of model 2 are shown in Table 3. The baseline variables of age, SAPS II and other underlying diseases significantly increased the CSHR for a fatal outcome. No nosocomial infection was significantly associated with the CSHR for death. However, patients with nosocomial pneumonia stay significantly longer in the ICU (CSHR = 0.59); a similar effect was seen for patients with nosocomial LRT (CSHR = 0.56). The CSHRs with regard to death in the ICU were not significant for these nosocomial infections.
Cumulative incidence functions (model 2)
Although patients with an elective surgery had a lower cause-specific risk of death (CSHR = 0.43), they tended to stay longer in the ICU compared to those patients without an elective surgery (CSHR = 0.56). This effect can also be seen in Figure 3a: the cumulative incidences of both endpoints start at a later time point for patients with elective surgery.
Patients with pneumonia on admission stay longer in the ICU (CSHR = 0.61); the CSHR for death was not significant. However, that also means that patients with pneumonia on admission die more frequently. This effect can be viewed in Figure 3b: the cause-specific risk of death decreased for patients with pneumonia on admission at the beginning of their ICU stay, but increased if they stay longer; the curves intersect.
Correlations among risk factors
The following time-dependent risk factors were highly correlated among each other: colostomy, enterostomy and jejunostomy (absolute values range between 0.6 to 0.9). There was a low correlation of the baseline risk factor 'intubated on admission' and the SAPS II score (0.5). All other correlation coefficients ranged between -0.4 and 0.4.
Many patient characteristics and significant risk factors for ventilator-associated pneumonia have been published. These include age, male gender, hospitalization prior to ICU admission, length of ICU stay, treatment in large hospitals, a low Glasgow Coma Scale (GCS), a poor Acute Physiology and Chronic Health Evaluation (APACHE) II or SAPS II score, respiratory failure, congestive heart failure, acute renal failure and dialysis, bronchoscopy, tracheotomy, re-intubation, duration of mechanical ventilation, detection of certain multi drug resistant pathogens, use of central vein catheters, bacteraemia, enteral feeding, and application of sucralfat or corticosteroids, [4, 15–24].
However, in most of these studies the time-dependent issue of nosocomial infections was ignored, that is, the time-dependent exposure was analysed as being known at time origin. This results in time-dependent bias . In addition, competing events such as discharge or death were not explicitly modelled. Recently, Resche-Rigon and co-authors point out that ICU discharge should be considered a competing event, when estimating the mortality of ICU patients . In this context, Schoenfeld argued that one should better focus on whether patients die rather then when they die, and therefore mortality should be analysed as a binary variable (30-day mortality) using a logistic regression . But that means that the time-dependent nature of nosocomial infections is ignored and it is impossible to consider time-dependent risk factors as for example, ventilation. In the present paper we applied multi-state models in order to accurately take these two important issues (that is, time-dependent risk-factors and competing events) into account.
The competing risks situation at hand, however, requires careful interpretation of the results: for example, in model 2 we find that pneumonia on admission has a (non-significant reducing) effect on the cause-specific hazard ratio of death, and an even more reducing (and significant) effect on the CSHR of discharge. This suggests that pneumonia on admission prolongs ICU stay; however, as the death hazard is not reduced as much as the discharge hazard is, there will eventually be more patients who are deceased . Thus, the competing risks model explains how pneumonia on admission contributes to mortality: pneumonia on admission prolongs ICU stay; each day, such a patient is again exposed to the (not significantly altered) risk of dying. As a consequence, there will be more patients with pneumonia on admission, who stay longer and die in the ICU.
In this study, we could show that elective surgery increases the CSHR for nosocomial pneumonia (model 1). Although nosocomial pneumonia is a risk factor for death, patients with elective surgery have a lower cause-specific risk of dying (model 2). However, these patients stay longer in the ICU. There are two possible explanations for this: firstly, there is an effect independent of whether they acquire NP during their ICU stay, and secondly via a nosocomial pneumonia which extends their ICU stay as well.
Our data from a competing risk model 1 confirmed mechanical ventilation as the key risk factor for the development of NP, with an increase in the CSHR of 5.90 (Table 2); this effect is accentuated by the parallel competing risks analysis of CSHR for direct discharge, which is significantly reduced by mechanical ventilation. Additional significant factors in our study were some form of surgery prior to ICU stay and the use of a nasogastric tube, though as a limitation it should be remembered that we did not consider all of the above-mentioned factors from previous works. Patients with diagnosed pneumonia on admission were much less likely to develop NP (CSHR = 0.02; Table 2). Our interpretation of this is that very few patients resolve from the initial pneumonia, thus they cannot acquire an additional NP afterwards.
There is little doubt that the acquisition of NP increases the length of ICU stay and the overall health care costs [18, 28]. However it is controversial whether NP also influences ICU mortality. Some studies found an increase in mortality due to NP, while other did not or found an increase for certain pathogens only . When comparing and evaluating these findings the possibility of publication bias should be kept in mind. It is less likely that studies without a significant increase in mortality will get published. None of the studies carried out previously have ever used a model of time-dependent variables to address the question of the mortality attributable to NP. Our competing risk model 2 did not show an increase of the CSHR for a fatal outcome after NP (CSHR = 0.87; p = 0.55; Table 3). However, as stated above, patients with NP require longer treatment in the ICU on average. This was confirmed by our findings (CSHR for discharge = 0.59; p < 0.01; Table 3). As a consequence patients with NP are exposed to the (not significantly altered) risk of dying in the ICU for a longer time period compared to patients without NP. As a result of this, more patients will die after NP. This is a typical competing risks phenomenon, which is discussed in detail by Beyersmann et al. .
More studies using competing risk models should be carried out to re-evaluate the impact of risk factors (especially time-dependent variables) on the occurrence of nosocomial infections and patient outcomes thereafter.
Nosocomial infections are time-dependent risk factors and should be analysed as such.
Ignoring the time-dependency of nosocomial infections leads to biased conclusions.
If the time to acquisition of a nosocomial infection is of interest, discharge/death is a competing event.
Whenever the length of ICU stay is of interest, death in the ICU is a competing event.
Only appropriate time-to-event analysis methods such as multi-state models can take the time-dependency of risk factors and competing events into account.
= Centers for Disease Control and Prevention
= cause-specific hazard ratio
= intensive care unit
= German Nosocomial Infection Surveillance System
= lower respiratory tract
= nosocomial pneumonia
= simplified acute physiology score.
Cook DJ, Walter SD, Cook RJ, Griffith LE, Guyatt GH, Leasa D, Jaeschke RZ, Brun-Buisson C: Incidence of and risk factors for ventilator-associated pneumonia in critically ill patients. Ann Intern Med 1998, 129: 433-440.
Safdar N, Dezfulian C, Collard HR, Saint S: Clinical and economic consequences of ventilator-associated pneumonia: a systematic review. Crit Care Med 2005, 33: 2184-2193. 10.1097/01.CCM.0000181731.53912.D9
Meric M, Willke A, Caglayan C, Toker K: Intensive care unit-acquired infections: incidence, risk factors and associated mortality in a Turkish university hospital. Jpn J Infect Dis 2005, 58: 297-302.
Apostolopoulou E, Bakakos P, Katostaras T, Gregorakos L: Incidence and risk factors for ventilator-associated pneumonia in 4 multidisciplinary intensive care units in Athens, Greece. Respir Care 2003, 48: 681-688.
Putter H, Fiocco M, Geskus RB: Tutorial in biostatistics: competing risks and multi-state models. Stat Med 2007, 26: 2389-2430. 10.1002/sim.2712
Andersen PK, Abildstrom SZ, Rosthoj S: Competing risks as a multi-state model. Stat Methods Med Res 2002, 11: 203-215. 10.1191/0962280202sm281ra
Gastmeier P, Geffers C, Sohr D, Dettenkofer M, Daschner F, Ruden H: Five years working with the German nosocomial infection surveillance system (Krankenhaus Infektions Surveillance System). Am J Infect Control 2003, 31: 316-321. 10.1067/mic.2003.66
Garner JS, Jarvis WR, Emori TG, Horan TC, Hughes JM: CDC definitions for nosocomial infections, 1988. Am J Infect Control 1988, 16: 128-140. 10.1016/0196-6553(88)90053-3
Grundmann H, Barwolff S, Tami A, Behnke M, Schwab F, Geffers C, Halle E, Gobel UB, Schiller R, Jonas D, Klare I, Weist K, Witte W, Beck-Beilecke K, Schumacher M, Rüden H, Gastmeier P: How many infections are caused by patient-to-patient transmission in intensive care units? Crit Care Med 2005, 33: 946-951. 10.1097/01.CCM.0000163223.26234.56
Beyersmann J, Gastmeier P, Grundmann H, Barwolff S, Geffers C, Behnke M, Ruden H, Schumacher M: Use of multistate models to assess prolongation of intensive care unit stay due to nosocomial infection. Infect Control Hosp Epidemiol 2006, 27: 493-499. 10.1086/503375
Andersen P, Borgan A, Gill D, Keiding N: Statistical models based on counting processes. New York: Springer; 1993.
Tai BC, Machin D, White I, Gebski V, On behlf of the EOI (The European Osteosarcoma Intergroup): Competing risks analysis of patients with osteosarcoma: a comparison of four different approaches. Stat Med 2001, 20: 661-684. 10.1002/sim.711
Therneau T, Grambsch P: Modeling survival data: extending the Cox model (statistics for biology and health). New York: Springer; 2000.
Teräsvirta T, Mellin I: Model selection criteria and model selection tests in regerssion models. Scandinavian J Statistics 1986, 13: 159-171.
Sofianou DC, Constandinidis TC, Yannacou M, Anastasiou H, Sofianos E: Analysis of risk factors for ventilator-associated pneumonia in a multidisciplinary intensive care unit. Eur J Clin Microbiol Infect Dis 2000, 19: 460-463. 10.1007/s100960000236
Georges H, Leroy O, Guery B, Alfandari S, Beaucaire G: Predisposing factors for nosocomial pneumonia in patients receiving mechanical ventilation and requiring tracheotomy. Chest 2000, 118: 767-774. 10.1378/chest.118.3.767
Ibrahim EH, Tracy L, Hill C, Fraser VJ, Kollef MH: The occurrence of ventilator-associated pneumonia in a community hospital: risk factors and clinical outcomes. Chest 2001, 120: 555-561. 10.1378/chest.120.2.555
Rello J, Ollendorf DA, Oster G, Vera-Llonch M, Bellm L, Redman R, Kollef MH, VAP Outcomes Scientific Advisory Group: Epidemiology and outcomes of ventilator-associated pneumonia in a large US database. Chest 2002, 122: 2115-2121. 10.1378/chest.122.6.2115
Pawar M, Mehta Y, Khurana P, Chaudhary A, Kulkarni V, Trehan N: Ventilator-associated pneumonia: incidence, risk factors, outcome, and microbiology. J Cardiothorac Vasc Anesth 2003, 17: 22-28. 10.1053/jcan.2003.4
Erbay RH, Yalcin AN, Zencir M, Serin S, Atalay H: Costs and risk factors for ventilator-associated pneumonia in a Turkish university hospital's intensive care unit: a case-control study. BMC Pulm Med 2004, 4: 3. 10.1186/1471-2466-4-3
Boots RJ, Lipman J, Bellomo R, Stephens D, Heller RF: Disease risk and mortality prediction in intensive care patients with pneumonia. Australian and New Zealand practice in intensive care (ANZPIC II). Anaesth Intensive Care 2005, 33: 101-111.
Myny D, Depuydt P, Colardyn F, Blot S: Ventilator-associated pneumonia in a tertiary care ICU: analysis of risk factors for acquisition and mortality. Acta Clin Belg 2005, 60: 114-121.
Tejerina E, Frutos-Vivar F, Restrepo MI, Anzueto A, Abroug F, Palizas F, Gonzalez M, D'Empaire G, Apezteguia C, Esteban A, Internacional Mechanical Ventilation Study Group: Incidence, risk factors, and outcome of ventilator-associated pneumonia. J Crit Care 2006, 21: 56-65. 10.1016/j.jcrc.2005.08.005
Gastmeier P, Sohr D, Geffers C, Behnke M, Ruden H: Risk factors for death due to nosocomial infection in intensive care unit patients: findings from the krankenhaus infektions surveillance system. Infect Control Hosp Epidemiol 2007, 28: 466-472. 10.1086/510810
van Walraven C, Davis D, Forster AJ, Wells GA: Time-dependent bias was common in survival analyses published in leading clinical journals. J Clin Epidemiol 2004, 57: 672-682. 10.1016/j.jclinepi.2003.12.008
Resche-Rigon M, Azoulay E, Chevret S: Evaluating mortality in intensive care units: contribution of competing risks analyses. Crit Care 2006, 10: R5. 10.1186/cc3921
Schoenfeld D: Survival methods, including those using competing risk analysis, are not appropriate for intensive care unit outcome studies. Crit Care 2006, 10: 103. 10.1186/cc3949
Schulgen G, Kropec A, Kappstein I, Daschner F, Schumacher M: Estimation of extra hospital stay attributable to nosocomial infections: heterogeneity and timing of events. J Clin Epidemiol 2000, 53: 409-417. 10.1016/S0895-4356(99)00182-1
Beyersmann J, Dettenkofer M, Bertz H, Schumacher M: A competing risks analysis of bloodstream infection after stem-cell transplantation using subdistribution hazards and cause-specific hazards. Stat Med 2007, 26: 5360-5369. 10.1002/sim.3006
We would like to thank all people that were involved in the German SIR-3 study.
The authors declare that they have no competing interests.
HG and PG initiated the SIR-3 study. MB created the database and online platform for the KISS system. SB and CG participated in collecting of the data. MW, JB and MS participated in the statistical analysis of the data. RPV, PG and HR participated in interpreting the data and drafting of the manuscript. All authors read and approved the final manuscript.
About this article
Cite this article
Wolkewitz, M., Vonberg, R.P., Grundmann, H. et al. Risk factors for the development of nosocomial pneumonia and mortality on intensive care units: application of competing risks models. Crit Care 12, R44 (2008). https://doi.org/10.1186/cc6852