Abstract
Abstract
Background:
Volume-infection relation studies have been published for high-risk surgical procedures, although the conclusions remain controversial. Inconsistent results may be caused by inconsistent categorization methods, the definitions of service volume, and different statistical approaches. The purpose of this study was to examine whether a relation exists between provider volume and coronary artery bypass graft (CABG) surgical site infection (SSI) using different categorization methods.
Methods:
A population-based cross-sectional multi-level study was conducted. A total of 10,405 patients who received CABG surgery between 2006 and 2008 in Taiwan were recruited. The outcome of interest was surgical site infection for CABG surgery. The associations among several patient, surgeon, and hospital characteristics was examined. The definition of surgeons' and hospitals' service volume was the cumulative CABG service volumes in the previous year for each CABG operation and categorized by three types of approaches: Continuous, quartile, and k-means clustering.
Results:
The results of multi-level mixed effects modeling showed that hospital volume had no association with SSI. Although the relation between surgeon volume and surgical site infection was negative, it was inconsistent among the different categorization methods.
Conclusions:
Categorization of service volume is an important issue in volume-infection study. The findings of the current study suggest that different categorization methods might influence the relation between volume and SSI. The selection of an optimal cutoff point should be taken into account for future research.
S
However, the findings of previous volume-infection studies are controversial. Some studies have found a negative relation between volume and infection [2–6]. Some found that the medium-volume group has better performance [7], whereas other studies have found no relation between provider's volume and infection [8,9]. Previous studies have shown the heterogeneity of defining what is high volume and what is low volume [10,11]. In addition, the neglect of cluster characterization of data that should use hierarchical models may result in biased estimation of the variation [12].
In recent years, most studies have adopted multi-level analysis to fix the data clustering issue, however, the heterogeneity of service volume definitions is still an issue that needs attention. Therefore, the purpose of this study was to compare the results of volume-infection studies under different service volume categorization methods, using coronary artery bypass graft (CABG) surgery as an example.
Patients and Methods
Study design
This retrospective study used a multi-level design to examine the relation between provider service volumes and SSI for CABG patients.
Database
We used the 2005–2009 Taiwan National Health Insurance Research Database. The database covers every episode of care provided to its 23 million Taiwanese enrollees (approximately 98% of the population). It is a de-identified secondary database containing patient-level demographic and administrative information. It is released for public access for research purposes [13].
Ethics statement
The protocol for this study was approved by the Institutional Review Board of the National Taiwan University Hospital (protocol #201001027R).
Study population
The study included all patients who underwent CABG surgery (ICD-9-OP code: 36.1x-36.2) surgery between 2006 and 2008 in all hospitals in Taiwan.
Exclusion criteria
Patients were excluded from analysis if they were: (1) Aged <20 y, (2) had post-operative surgical site due to prior operation, (3) death after surgery in hospital, and (4) patient's characteristics were incomplete.
Definition of Variables
Dependent variables
The SSI cases were identified by the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes. The SSI cases are divided into two categories: Index hospitalization events and post-discharge events (i.e., SSIs that occur within 1 y after discharge and require re-admission to a hospital and/or the use of ambulatory services). The ICD-9-CM codes for hospitalization events are 996.03, 996.61, 996.72, and 998.5. The ICD-9-CM codes for ambulatory services events are 038.0–038.4, 038.8, 038.9, 682.6, 682.9, 780.6, 790.7, 875.0, 875.1, 891.0, 891.1, 996.03, 996.61, 996.72, 998.3, and 998.5. Following the work by Wu et al. [6], this study used the secondary ICD-9-CM diagnosis codes for index hospitalization events, and the primary and secondary diagnosis codes for post-discharge events as criteria for SSI in order to avoid cases in which infection developed prior to hospitalization [6].
Independent variables
The definition of surgeons' and hospitals' service volume was the surgeons' and hospitals' cumulative CABG service volumes in the previous year for each CABG operation. We analyzed service volume as three different variables: (1) A continuous variable; (2) a subjective categorical variable based on first quartile and third quartile as cutoff points; and (3) a data-driven categorical variable based on k-means clustering approach. The subjective categorical variable was the most common method to categorize the service volume [5,14–16]. This study pre-determined the service volume be categorized into low, medium, and high groups. The cutoff points of surgeon's volume were 11 and 48 cases, and of hospital's volume were 73 and 176 cases. The data-driven categorical variable based on k-means clustering approach had been applied by Kuo et al. [12] to explore the relations among volume, mortality, and recurrence in breast cancer. This study also pre-determined the surgeon and hospital service volume into low-, medium-, and high-volume groups by k-means clustering method. The cutoff points determined by the statistical model were 25 and 64 cases for surgeon's volume, and 107 and181 cases for hospital's volume.
Control variables
In addition to service volume and SSI, the study collected patient-level data such as age, gender, Charlson/Romano comorbidity index, numbers of vessels obstructed (as a proxy indicator of duration of operation [17]), and underlying diseases, such as diabetes mellitus (ICD-9-CM code: 250), chronic obstructive pulmonary disease (ICD-9-CM code: 490-496), heart failure (ICD-9-CM code: 584-586), renal failure and renal insufficiency (ICD-9-CM code: 428), which are associated with surgical infections [6]. Regarding the surgeon-level data, the study collected surgeon's gender and age. For hospital level, this study also collected hospital ownership and accreditation level as control variables.
Statistical Analysis
The data were analyzed using multi-level logistic regression. Most previous studies have used conventional regression models, which ignore the possible correlations between outcomes for a given surgeon or hospital, resulting in the underestimation of a variable's standard error and produce incorrect statistical significance [13,18,19]. Therefore, multi-level logistic regression was applied to explore the relation between provider volume in CABG and SSI, adjusted for patient severity, comorbidities, and demographics, and characteristics of the surgeon/hospital. Patients (level 1) were considered to be nested within surgeons (level 2) and then within hospitals (level 3). A stepwise variable selection procedure was performed for avoiding multicollinearity. All statistical analyses were performed using SAS (version 9.2, SAS Inc., Cary, NC).
Results
The results showed that there were 10,405 CABG procedures performed by 214 surgeons in 54 hospitals during 2006–2008 in Taiwan. Most of the patients were male (>75%), with a mean age of 65 y. Regarding the status of the four diseases associated with surgical infections, more than 50% of the patients had diabetes mellitus, one-third had heart failure, approximately 15% had renal failure and renal insufficiency and more than 20% had chronic obstructive pulmonary disease. The mean number of vessels obstructed was approximately 1.6. The average comorbidity score was approximately 1.5. There were 403 patients (3.87%) identified as having SSIs during hospitalization and 1,079 patients (10.37%) were identified as index hospitalization and post-discharge event (Table 1).
Mean (standard deviation).
n (%).
Table 2 shows the characteristics of hospitals and surgeons that had ever performed CABG operations during 2006–2008. Among 54 hospitals, most were not-for-profit hospitals; 19 were medical centers (35.19%); 34 were regional hospitals (62.96%), and others were community hospitals (1.85%). The annual average service volume was from 149.84 cases in 2006 to 152.66 cases in 2008; the maximum volume was from 294 cases in 2006, to 392 cases in 2008. In terms of surgeons, more than 90% were male, with a mean age of approximately 42 y. The average service volume was from 52.53 cases in 2006, to 50.38 cases in 2008, the maximum volume was from 120 cases in 2006, to 169 cases in 2008.
n (%).
Mean (standard deviation).
As shown in Table 3, the incidence of SSI among female patients with four major comorbidities was significantly greater than that among males and those patients without four major comorbidities. The length of stay and the Charlson/Romano comorbidity index of SSI cases were also higher than non-SSI cases. Most patients underwent CABG surgery in medical centers and not-for-profit hospitals. Community hospitals had a lower incidence (5.88%) than medical centers (10.3%) and regional hospitals (10.53%). Not-for-profit hospitals had a greater incidence (10.84%) than other hospitals. However, the differences did not reach statistical significance. In terms of hospitals' and surgeons' service volume, the average of hospitals' service volume in the non-SSI group and the SSI group were 147.50 and 145.80 (p=0.717), the average of surgeons' service volume in non-SSI group and the SSI group were 50.62 and 44.14 (p<0.0001). The results also showed that the SSI incidence of high-volume hospitals and surgeons was the lowest among the groups, except for hospitals' volume categorized by the quartile method.
Mean (standard deviation).
n (%).
Hospital volume is defined as the number of coronary artery bypass graft operations performed at a particular hospital in the previous year before each operation.
Surgeon volume is defined as the number of coronary artery bypass graft operations performed by a particular surgeon in the previous year before each operation.
SSI=surgical site infection; COPD=chronic obstructive pulmonary disease.
All control variables were put into the model. After stepwise selection, diabetes mellitus, chronic obstructive pulmonary disease, heart failure, and renal failure and renal insufficiency were kept in the model. After that, this study put the independent variable with different definitions into the multi-level model. Table 4 summarizes our multi-level mixed effect model. Model 1 found that surgeon's volume had a negative relation with infection, however, hospital volume had no relation with infection. Model 2 found that patients treated by low-volume (odds ratio [OR]=1.458, p=0.017) and medium-volume surgeons (OR=1.308, p=0.059) had greater risk than by high-volume surgeons. In addition, the risk among high-, medium-, and low-volume hospitals was not statistically significant. Model 3 found that there was no relation between hospital volume and SSI infection, and also surgeon volume. In addition, four comorbidities were found to have a positive relation with SSI in all models. Also, model 3 had better performance in model fitting.
Adjusted by heart failure (reference: with heart failure), COPD (reference: with COPD), diabetes mellitus (reference: with diabetes), and renal insufficiency (reference: with renal insufficiency), all control variables were significant (p<0.001).
OR=Odds Ratio; CI=confidence interval; LL=lower limit; UL=upper limit; COPD=chronic obstructive pulmonary disease; SE=standard error. AIC=Akaike Information Criterion, one of the popular measures for mode fitting (small is better); BIC=Bayesian Information Criterion, one of the popular measures for model fitting (small is better).
Discussion
This nationwide population-based study was the first, to our knowledge, to use an uncategorized subjective method and a data-driven method to categorize service volume to explore the volume-infection relation. The surgeon and hospital annual average service volume were approximately 50 and 150, the service volume were less than that of the United States [20], but similar to previous experience in Taiwan [21]. In addition, unlike previous studies that calculated the accumulated service volume of hospitals or surgeons during the study period [3,5,6,8], this study considered the cumulative CABG service volumes in the previous year for each CABG operation, which was a similar design to that used by Solomon et al. [14] and Yasunaga et al. [22]. Therefore, this study could better identify whether “practice makes perfect” applies to the relation between CABG service volume and SSI.
There were three major findings of this study. First, the relation between hospital volume and SSI did not exist in these data; the relation between surgeon volume and SSI was partially consistent among different categorization methods. In addition, the data-driven categorization method had better model fitting than other methods. Second, the phenomenon of centralized service supply might exist. The hospital and surgeon annual average volume decreased slightly in 2006–2008, however, we also found that the standard deviation and maximum volume had increased. It might indicate that the phenomenon of CABG service centralization existed in Taiwan during these 3 y. Whereas we calculated the coefficient of variance, we found the level of deviation for surgeons was greater than hospitals (Table 2). The phenomenon of CABG service centralization was more obvious in surgeons than in hospitals. Third, patient's comorbidities might be more important than provider characteristics. In this study, patient, surgeon, and hospital characteristics were used as covariates, however, only four comorbidities remained after stepwise selection. These four comorbidities were also significant in all multi-level mixed effect models. It appeared the patient's comorbidities might be more important than surgeon and hospital characteristics.
In this study, we found the relation between provider's volume and SSI was ambiguous. In terms of surgeon's volume, the relation between surgeon volume and SSI was not consistent after adjusting other covariates. The result of the continuous model (model 1) showed that there was a negative relation between surgeon volume and SSI, whereas the result of quartile method model (model 2) demonstrated that the patients operated on by low-volume surgeons had higher risk than by high-volume surgeons. The risk of medium-volume surgeons was borderline significant higher than high-volume surgeons. However, the result of k-means model (model 3) revealed there was no relation between surgeon volume and SSI.
The controversial finding might be caused by the definition of service volume; it might also indicate that the cutoff point could be an important issue in volume-outcome studies [10,11,23]. Previous volume-outcome studies have categorized service volume in a subjective manner. Some of these studies have categorized service volume into three [6,8,24] or four groups [15,16,25] and others have used a specific case number [26] or percentage [2] for categorization. The categorization methods of previous studies and the current work were based on the linear assumption existed. Whereas a non-linear relation existed between volume and outcome, how to choose the optimal cutoff point should be dealt with carefully in future research.
In terms of hospital volume, there was no relation found in our study. Furthermore, our study contradicts that of Wu et al. [6] who found that the hospital's volume had a greater effect than surgeon's volume, and claimed that the findings may imply that a hospital's teamwork was more important than its individual surgeons. The data of our study were 5 y later and longer than that of Wu et al. [6]. During these 5 y, the severe acute respiratory syndrome (SARS) outbreak occurred in Taiwan in 2003. Subsequently, the hospital infection control system in Taiwan was reviewed and re-designed. These efforts might have also improved SSI control in hospitals and thus led to the inconsistent results with the study by Wu et al. [6].
Although a multi-level analysis was applied to manage the nesting issue in the data and three different methods were used to examine the relation between service volume and SSI, two common limitations of studies that use secondary data impacted this study. The first limitation was the accuracy of the outcome variable because data from health care claims data are not designed for infection surveillance. Previous studies indicated that the accuracy issue might exist in claim data because of an inaccurate code list or be distorted by payment scheme or other factors [27–29]. It was unavoidable in this study. The second limitation relates to unmeasured variables, such as length of stay before operation, infection condition, hair removal, duration of operation, the environment, surgical skills, use of post-operative drains, number of operations involved, and the surgical site and wound care [17]. Furthermore, the information regarding type (elective or urgent) and incision site of surgery was not available in Taiwan NHI claims data.
Conclusion
The findings of this study suggest that surgeon volume and patient comorbidity may play more important roles than hospital volume in SSI. However, because the relation between surgeon volume and infection was not consistent, how to choose the optimal cutoff point should be an important issue for future research.
Footnotes
Acknowledgments
This study was supported by the Taiwan Centers for Disease Control (DOH100-DC-1020).
T.H.Y developed study concept, analyzed the data, and drafted the manuscript; Y.C.T. was involved in collecting data and drafting the manuscript; and K.P.C. collected data, coordinated the study, and revised the manuscript. All of the authors have prepared, read, and approved the manuscript.
The authors would like to thank the Centers for Disease Control of Taiwan for financially supporting this project.
Author Disclosure Statement
The authors certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.
