Healthcare-associated infections (HAIs) are an important patient safety component of several key pieces of healthcare legislation in the United States. In 2008, the Centers for Medicare and Medicaid Services stopped reimbursing hospitals for several HAIs classified as “never events.”Reference Brown, Doloresco Iii and Mylotte 1 Additionally, as part of the Affordable Care Act (ACA), the hospital-acquired condition (HAC) reduction program imposes financial penalties for hospitals that perform poorly on measures of hospital-acquired adverse events, including HAIs. 2
The measurement of the cost savings associated with interventions to decrease the rates of HAIs is becoming an important tool for evaluating the cost-effectiveness of such policies. An important component of hospital costs attributable to HAI is the excess length of hospital stay (LOS).
Estimating the impact of HAI on LOS, however, is challenging because HAIs, by definition, are not present during a patient’s entire hospitalization.Reference Barnett, Beyersmann, Allignol, Rosenthal, Graves and Wolkewitz 3 , Reference Samore and Harbarth 4 Only those hospital days that accrue after the occurrence of the nosocomial infection are possibly secondary to infection. Van Walraven et alReference Beyersmann, Gastmeier, Wolkewitz and Schumacher 5 claimed, and Beyersmann et alReference van Walraven, Davis, Forster and Wells 6 provided mathematical proof, that treating a time-dependent exposure as time-fixed inflates estimates of the effect of that exposure on an outcome whether (1) the true effect prolongs the outcome, (2) the true effect accelerates the outcome, or (3) the exposure has no true effect on the outcome. For example, estimates of the excess LOS due to an HAI that do not take into account the fact that the HAI is not present during a patient’s entire hospital stay always overestimate the effect, regardless of what the true effect is. The concept of time-dependent bias is well appreciated, but less is known about the magnitude of its impact, particularly as it relates to outcomes attributable to HAIs.Reference De Angelis, Murthy, Beyersmann and Harbarth 7
The aim of this analysis was to summarize and quantify findings of studies of HAI-attributable LOS that compared methods that accounted for the timing of infection with those that do not account for the timing of infection. The purpose of this comparison is to inform the interpretation of literature on cost-saving impacts of preventing HAIs. An estimation of the magnitude of bias may be useful for decision analytic scientists to deflate the cost and LOS input parameters in cost-effectiveness analyses. The demonstration that this bias leads to nontrivial excess LOS estimates may also guide investigators when choosing between available epidemiologic study designs and statistical methods. Ultimately, policy makers may have more reasonable expectations regarding the effect of infection prevention efforts.
Time-Fixed Methods for Estimating Excess LOS Due to HAI
Estimates of excess LOS due to HAI are often generated by comparing the total LOS between patients who acquire HAIs during their hospitalization and those who do not.Reference Cosgrove, Kaye, Eliopoulous and Carmeli 8 – Reference Stewart and Hollenbeak 11 This comparison can be conducted as a simple comparison of means, as a multivariable regression model, or even using propensity score matching. However, because none of these methods accounts for the time-varying nature of the HAI, they are all subject to time-dependent bias. This bias occurs because in methods that treat the HAI as time fixed, the time spent in the hospital prior to the occurrence of the HAI is incorrectly attributed to the HAI. Therefore, the pre-HAI LOS will be added to the post-HAI LOS, thus making the attributable LOS appear to be larger than it actually is. In addition, patients are at greater risk for HAIs the longer they are in the hospital. This reverse causality can also lead to biased estimates of the true effect of HAI on LOS.
Methods for Overcoming Time-Dependent Bias
Two primary methods for overcoming time-dependent bias have been utilized in published studies: multistate model and matching on timing of infection. Both methods accomplish this task by taking into account the timing of infection. Multistate models treat an HAI as one of several mutually exclusive states through which a hospitalized patient can move. The method of matching on the timing of infection matches patients with HAIs to patients without HAIs with a requirement that the non-HAI patient LOS be at least as long as the time until the HAI in the infected case patient. In performing this matching, it is important that the selected uninfected patients are those who could still become infected at a later date. This method, called exposure density sampling, has been shown to reduce bias in the estimation of the effect of HAI because it is not conditional based on a future event.Reference Wolkewitz, Beyersmann, Gastmeier and Schumacher 12 By conditioning on the future event of having an HAI or not at a later date, matching methods that do not employ exposure density sampling are subject to similar biases that plague conventional estimation methods that do not take the timing of infection into account.
In the following section, we summarize the studies identified through a systematic literature search of PubMed for combinations of the following search terms: time-dependent bias, length of stay, attributable length of stay, excess length of stay, hospital-acquired infection, healthcare-associated infection, and nosocomial infection. From the list of published papers identified through these searches, we chose those that used one or both of these methods as well as methods in which an HAI was treated as a time-fixed variable using the same cohort of patients.
RESULTS
Search
Our literature search yielded 516 records; of these, 276 were duplicates. Of the 240 articles that remained after the removal of duplicates, 8 met the criteria described above and were included in our analysis and summary. To overcome time-dependent bias in the estimation of the excess LOS associated with HAIs, 4 studies used multistate models only, 2 studies matched HAI patients to control patients on the timing of infection, and 2 studies used both methods.
Multistate Models
(Table 1) In a landmark paper in 2000, SchulgenReference Schulgen, Kropec, Kappstein, Daschner and Schumacher 13 illustrated these concepts using 2 separate studies. The first study (Study 1) consisted of 765 patients admitted to the intensive care unit (ICU) after having undergone open heart, large-bowel, or biliary-tract surgery between November 1988 and September 1989. In Study 2, 756 patients were followed who were admitted to anesthesiological or medical ICUs in the same hospital between July 1991 and July 1993. The authors used 2 methods to estimate the impact of HAI on LOS in which the HAI was treated as a time-fixed variable. The first was a simple comparison of the unadjusted mean LOS for patients with and without an HAI. This second method involved a comparison of mean LOS after matching patients with and without an HAI on observable baseline confounders such as age, sex, diagnosis, and, for Study 1 only, the degree of contamination of the surgical wound. The multistate model used by SchulgenReference Schulgen, Kropec, Kappstein, Daschner and Schumacher 13 included 5 states: no nosocomial infection, nosocomial pneumonia, nosocomial sepsis, discharge, and death.
Published Estimates of the Magnitude of Time-Dependent Bias: Conventional Methods vs Multistate Models
NOTE. HAI, healthcare-associated infection; LOS, length of stay; CI, confidence interval; CLABSI, central line associated bloodstream infection; CAUTI, catheter-associated urinary tract infection; VAP, ventilator associated pneumonia.
a De Angelis (2011) and Macedo-Vinas (2011) stated that they identified HAIs due to methicillin-resistant Staphylococcus aureus. None of the other studies explicitly stated which organism was associated with the HAIs.
Roberts et alReference Roberts, Scott and Hota 14 examined outcomes for a cohort of 1,500 patients admitted to a large teaching hospital in the United States in the year 2000. Of these patients, 159 (10.6%) developed an HAI during their hospitalization. Treating HAI as a time-fixed variable, an ordinary least squares regression was used to estimate the impact of HAI on LOS. To overcome time-dependent bias, the investigators used multistate model with 3 states: no HAI, HAI, and discharge/death.
Using data from 9,545 patients (8.7% of whom had HAIs) admitted to the ICU in 1 of 11 hospitals in Buenos Aires, Argentina, Barnett et alReference Barnett, Beyersmann, Allignol, Rosenthal, Graves and Wolkewitz 3 estimated the impact of HAI on LOS assuming time-fixed HAI using a parametric γ model. This method yielded an excess LOS that was substantially larger than the estimate generated using a 3-state (no HAI, HAI, discharge/death) multistate model. The 700% relative difference between the 2 results was the largest relative difference of any among the studies in this analysis.
De Angelis et alReference De Angelis, Allignol and Murthy 15 used data from a prospective cohort studyReference Wolkewitz, Allignol, Harbarth, de Angelis, Schumacher and Beyersmann 16 of 976 patients admitted to the University of Geneva Hospital between July 2004 and May 2006, with 167 having a nosocomial MRSA infection. Wolkewitz et alReference Wolkewitz, Allignol, Harbarth, de Angelis, Schumacher and Beyersmann 16 used a fourth state in their multistate model that captured MRSA colonization as an intermediate state between admission and MRSA infection or admission and discharge/death. The estimates generated using methods that treated HAI as a time-fixed exposure are reported in the subsequent paper by Wolkewitz et al.Reference Wolkewitz, Allignol, Harbarth, de Angelis, Schumacher and Beyersmann 16
Using data from surgery patients admitted to the University of Geneva Hospital during 2009, Macedo-Vinas et alReference Macedo-Vinas, De Angelis and Rohner 17 estimated the excess LOS due to MRSA HAI using a multistate model with data from 167 infected patients and 25,766 uninfected patients. This estimate was slightly lower than the estimate obtained by taking the difference in mean LOS between the infected patients and a cohort of uninfected patients who were matched 2:1 (uninfected to infected) on age, sex, and diagnosis-related group.
Schumacher et alReference Schumacher, Allignol, Beyersmann, Binder and Wolkewitz 18 used data from a university hospital in Germany. Of the 1,876 patients admitted to the ICU during 2000–2001, 158 developed nosocomial pneumonia. These investigators contrast the excess LOS obtained using a Cox proportional hazards regression model with those from a multistate model with 4 states: no HAI, HAI, discharge, and death.
Matching on the Timing of Infection
In addition to a multistate model described in the earlier section, Schulgen also matched infected and non-infected patients on time to infection (Table 2).Reference Schulgen, Kropec, Kappstein, Daschner and Schumacher 13 The noninfected patients eligible to be matched were those who were discharged from the hospital without infections so exposure density sampling was not used. This method yielded estimates of the excess LOS that were larger than those generated from multistate models.
Published Estimates of the Magnitude of Time-Dependent Bias: Conventional Methods vs Matching on Timing of Infection

Note. HAI, healthcare-associated infection; LOS, length of stay; CI, confidence interval; UTI, urinary tract infection; BSI, bloodstream infection; SSI, surgical site infection; GI, gastrointestinal infection.
a None of the studies explicitly stated which organism was associated with the HAIs.
Vrijens et al,Reference Vrijens, Hulstaert, Van de Sande, Devriese, Morales and Parmentier 19 with a cohort of 1,839 patients admitted to 1 of 19 acute care hospitals in Belgium, examined the impact of bloodstream HAIs on LOS.Reference Vrijens, Hulstaert, Van de Sande, Devriese, Morales and Parmentier 19 Without taking the timing of infection into account, the investigators matched each infected patient to 1 noninfected patient regarding age, all-patient refined-diagnosis–related group (APR-DRG), principal diagnosis, and Charlson comorbidity index. The attributable LOS was calculated as the mean difference in LOS between the infected and matched uninfected patients. An additional analysis added the time of infection as a matching criterion. The authors do not state whether exposure density sampling was used.
In 2012, Vrijens et al Reference Vrijens, Hulstaert, Devriesse and van de Sande 20 extended their earlier work by including urinary tract, surgical site, and gastrointestinal infections along with bloodstream infections from 659 patients in 63 acute care hospitals in Belgium. These investigators demonstrated the importance of including the time prior to infection, which they refer to as exposure duration, as a matching variable when selecting noninfected control patients and matched infected and noninfected patients in a 3:1 ratio. The absolute difference in excess LOS estimates from matching on hospital, APR-DRG, and age compared to matching on these 3 variables plus exposure duration was 28.3 days, the highest absolute difference of any study in this analysis. Again, there is no indication as to whether exposure density sampling was used.
Finally, in addition to the multistate model described above, Schumacher et alReference Schumacher, Allignol, Beyersmann, Binder and Wolkewitz 18 also matched their infected patients to uninfected patients. This matching used exposure density sampling, and a hazard ratio resulted from this exercise that was very similar to the one generated from the multistate model but with a wider confidence interval (0.64; 95% CI, 0.55–0.73 vs 0.65 days; 95% CI, 0.51–0.83). However, for this method, the excess LOS was not reported.
Magnitude of Time-Dependent Bias
These summary estimates indicate that estimates of bias were greater when time-fixed methods were compared to multistate models than when they were compared with matched cohort analyses. This is true even within the same studies. For example, in Schumacher et al, the absolute (relative) differences in estimates of excess LOS were 15.7 (253.2%) when using a multistate model and 10.6 (93.8%) when matching on the timing of infection. Similarly, in the 2 studies described in Schulgen et al, the multistate model yielded an absolute (relative) difference of excess LOS of 9.0–10.0 days (91.8%–292.6%) and matching on the timing of infection resulted in an absolute (relative) difference of excess LOS of 5.2–7.4 days (62.7%–64.9%) compared with time-fixed approaches.
DISCUSSION
The results summarized here are useful for a number of purposes. First, with these estimates, policy and decision makers can make more realistic projections of the financial impact of infection prevention efforts. Second, when constructing models to perform cost-effectiveness analyses, analysts can use the estimates presented here to describe with greater accuracy the expected number of bed days saved from infection prevention programs. Not adjusting for sources of bias when estimating extra stay yields an overstatement of the cost savings from infection control programs and makes them appear to have greater value for money. Third, a demonstration that this bias leads to non-trivial difference in estimated LOS may further convince investigators of the importance of using alternative methods to the total LOS approach.
The magnitude of time-dependent bias in the estimation of excess LOS due to HAIs may depend on a number of factors, including the infection rate, the discharge rate, confounding, and model specification. In describing immortal time bias, SuissaReference Suissa 21 determined that the magnitude of this type of bias depends on the proportion of the exposed group’s time that is, in fact, unexposed (ie, immortal time). In other words, the immortal time bias is greater when more of the exposed patients’ time is unexposed. Suissa also showed that the size of this bias depends positively on the proportion of unexposed to exposed person-time.
Overestimation of the attributable LOS caused by treating a time-varying exposure as if it was present at baseline can occur in any of the many possible hospital-acquired adverse events.Reference Samore, Shen and Greene 22 Noninfectious disease-related events that are listed as part of the ACA HAC reduction program include pressure ulcers; iatrogenic pneumothorax; postoperative hip fractures, pulmonary embolism, and deep vein thrombosis; wound dehiscence; and accidental puncture and lacerations. 23 The reduction program will no doubt lead to efforts to reduce these adverse events, and economic analyses will be needed to evaluate the impacts of these efforts. However, these analyses will not be as useful as they could be because of inaccurate inputs resulting from time-dependent bias. Better estimates of the impact of non-HAI hospital-acquired adverse events are clearly needed.Reference Suissa 21 In the meantime, the estimates of the magnitude of time-dependent bias presented here can be used to deflate estimates for other conditions.
Estimates of the size of time-dependent bias may vary based on many factors. For instance, the bias may be reduced when the HAI occurs early during the hospitalization (eg, early-onset pneumonia after major trauma). The severity and type of infection as well as the infectious pathogen may also play a role, such as how soon after admission the HAI occurs. Finally, time-dependent bias may not be relevant for studies investigating the inpatient economic burden of HAIs that are identified after discharge. For example, surgical site infectionsReference McIntyre, Warner, Nester and Nathens 24 , Reference Sands, Vineyard and Platt 25 and hospital-acquired Clostridium difficile infectionsReference Chang, Krezolek, Johnson, Parada, Evans and Gerding 26 – Reference Palmore, Sohn, Malak, Eagan and Sepkowitz 29 often manifest after hospital discharge.
An important limitation to this study is publication bias. The estimates of time-dependent bias that we report come from published studies, many of which were focused on emphasizing the new methods for handling this bias. We were unable to find any published studies that reported zero or negative time-dependent bias. Furthermore, time-dependent bias could be less relevant in analyses of excess LOS after surgical site infection in patients with infection onset after discharge.
In conclusion, treating HAI as a time-fixed variable, as opposed to a time-varying variable, leads to substantial inflation of the effect of HAI on excess LOS. Future studies aiming to determine the excess costs and LOS of HAI should utilize one of several methods we discussed to reduce time-dependent bias and lead to more accurate estimates of HAI outcomes.
Acknowledgments
Financial support: The work in this paper was conducted through a contract with the Centers for Disease Control and Prevention. This material is the result, in part, of work supported with resources and the use of facilities of the George E. Wahlen Department of Veterans Affairs Medical Center, Salt Lake City, Utah.
Potential conflicts of interest: All authors report no conflicts of interest relevant to this article.

