EJCTS Click here to go to Edwards website
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to Personal Folders
Right arrow Download to citation manager
Right arrow Author home page(s):
Eyal E. Porat
Anthony L. Estrera
Anders N. Vinnerkvist
Hazim J. Safi
Right arrow Permission Requests
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Miller, C. C.
Right arrow Articles by Safi, H. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Miller, C. C., III
Right arrow Articles by Safi, H. J.
Related Collections
Right arrow Education
Right arrow Great vessels

Eur J Cardiothorac Surg 2003;23:1023-1027
© 2003 Elsevier Science NL


Analysis of short-term multivariate competing risks data following thoracic and thoracoabdominal aortic repair

Charles C. Miller, III*, Eyal E. Porat, Anthony L. Estrera, Anders N. Vinnerkvist, Tam T.T. Huynh, Hazim J. Safi

Department of Cardiothoracic and Vascular Surgery, University of Texas Health Science Center at Houston, Medical School, Houston, TX 77030, USA

Received 12 September 2002; received in revised form 3 March 2003; accepted 6 March 2003.

* Corresponding author. Tel.: +1-713-500-5420 fax: +1-713-500-0656
e-mail: charles.c.miller{at}uth.tmc.edu


    Abstract
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
Objective: Estimating the overall successfulness of a treatment can be difficult when success is defined by freedom from multiple endpoints that are each subject to competing risks. We describe a method for modeling short-term competing outcomes. Methods: We used polytomous categorical variable modeling to describe the 30-day onset of renal failure, neurologic deficit, stroke or death (events) following repair of 841 thoracoabdominal aortic aneurysms. This was to determine whether common risk factors had a multivariate association with these outcomes, and whether predictor variables might be positively associated with some outcomes and negatively associated with others. The goal was to determine whether a single aggregate-endpoint logistic model could accurately predict the probability of good outcome 30 days following surgery. Results: When more than one event occurred in a single patient, the first (or most severe simultaneous) event was used for censoring. Five hundred and ninety-three out of 841 (70.5%) patients had no postoperative events. The most common event was renal failure. We detected five predictors that were significant for at least one of the four outcomes. These were age, poor preoperative renal function (RENAL), acute dissection, extent II aneurysm, and use of cerebrospinal fluid drainage and distal aortic perfusion (ADJUNCT). Only RENAL was significant for all outcomes. ADJUNCT was highly significant only for neurologic deficit in the polytomous analysis and dropped out of the aggregate-endpoint multiple logistic model. Conclusion: Polytomous-outcome multivariate categorical modeling can detect effects missed by aggregate models, and is a valuable and statistically powerful method for evaluating risk factor effects on multiple competing endpoints.

Key Words: Competing risks • Logistic regression • Multivariate model • Surgical outcome • Aortic aneurysm • Risk factors


    1. Introduction
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
Short-term risk analysis in the surgical literature has traditionally focused on the prediction of single events [1]. These events are usually limited to major postoperative complications such as renal failure, stroke or death, and their likelihood is described one outcome at a time, using separate analyses. What may be more important to patients than the likelihood of any specific undesirable event, however, is the overall probability of success – of going home in good condition after surgery.

Overall success is often cast as a competing risks problem. Traditional competing risks studies use failure-time statistics to look at long-term results, and draw a contrast between ‘actual’ and ‘actuarial’ analysis for describing the way outcomes are distributed over time [2,3]. Multivariate competing-risks analyses involving multiple predictor variables and multiple conditional failure-time distributions are a complex topic and are not the subject of this paper. Methods for modeling several outcomes simultaneously over the short term, such as the period of hospitalization, have received less attention in the literature. In this paper, we focus only on competing risks that are described by a short-term period prevalence endpoint.

One simple approach to estimating overall success might involve creation of a composite variable that is true when a good outcome occurs and false when any bad outcome occurs [4]. It would be termed a composite variable because several bad-outcome variables have to be false for ‘good outcome’ to be true. Accurate prediction of good outcome means predicting the absence of a range of undesirable events. However, these events may have entirely different risk factors, or may have common risk factors that nonetheless have dissimilar effects on outcome. This means that mathematical models attempting to predict a composite outcome with variables that are weighted differentially could suppress important predictor variables if their differential effects cancel or overpower one another. Consequently, we need a method to evaluate the possibility that a variety of risk factors imposing differential effects may form individual components or ‘layers’ of the composite outcome variable.

A composite endpoint is attractive for purposes of computing summary probability, but to evaluate global probability it is desirable to understand the multivariate effects of the risk factors on the component outcomes. We describe a simplified method for studying simultaneously the effects of predictor variables on multiple potential outcomes, following surgical repair of thoracoabdominal aortic aneurysm.


    2. Methods
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
The population for this study consisted of 854 patients who underwent repair of thoracic or thoracoabdominal aortic aneurysms between 1991 and, 2001 on our service at The Methodist Hospital or Memorial Hermann Hospital in Houston, Texas.

2.1. Surgical methods
The adjuncts cerebrospinal fluid (CSF) drainage and distal aortic perfusion play a critical role in our surgical technique, having been found to lower the incidence of paraplegia and paraparesis [5,6]. Briefly, the highlights of our techniques are as follows. For CSF drainage and for monitoring CSF pressure a catheter is placed in the third or fourth lumbar space. Pressure is maintained at less than 10 mmHg throughout the procedure and for 72 h postoperatively. Following the thoracoabdominal incision, the diaphragm is preserved, exposing only the aortic hiatus. The patient is anticoagulated with heparin. The left atrium is cannulated via the left pulmonary vein or the left atrial appendage. A BioMedicus (Minneapolis, MN) pump with an in-line heat exchanger is attached to this cannula and the arterial inflow is established through the left common femoral artery, or the descending thoracic aorta, if the femoral artery is not accessible. For extensive aneurysms, sequential clamping is used. Reattachment of patent lower intercostal arteries (T8–T12) is performed routinely. The celiac, superior mesenteric, and renal arteries are perfused with cold perfusate (blood at 4 °C). Renal temperature is maintained at approximately 15 °C, but core body temperature is kept between 32 and 33 °C to avoid hypothermia.

2.2. Statistical methods
Data were collected prospectively on a wide variety of demographic, physiologic and operative data by a trained Master's level research nurse and entered into our specially designed Microsoft Access database. Database queries were exported to SAS version 8.02 for analysis (SAS Institute Inc., Cary, NC, USA). Preoperative renal dysfunction was defined as a serum creatinine greater than 2 mg/dl, or a recent history of renal failure or hemodialysis. Postoperative renal dysfunction was defined by a rise in serum creatinine of 1 mg/dl per day for 2 consecutive days after surgery or a requirement for hemodialysis. Pulmonary complications were defined as any event leading to mechanical ventilation of greater than 72 h postoperatively. Postoperative neurologic deficit was defined as paraplegia or paraparesis, observed as the patient awoke from anesthesia, and regardless of severity or postoperative recovery. Postoperative stroke was defined as a focal neurologic deficit due to cerebral infarction, as identified by a neurologist following patient computed tomography scan or magnetic resonance imaging. Mortality rate was based on patient deaths occurring within 30 days of surgery.

Because we cast this as a competing risks problem, a censoring hierarchy was required. Ideally, the censoring event would always be identified as the first event to occur. But because our database is not designed to capture time-to-event for postoperative complications other than death, death was the only event subject to time censoring. Non-fatal events that occurred prior to death censored death in order of severity. The order of severity for non-fatal events, from least to most severe, was renal failure, paraplegia/paraparesis, and stroke. The most severe non-fatal event censored all other events that occurred in any individual and was the censoring event for the analysis – the event that made the complication variable true. Death within 30 days was the censoring event if not preceded by a non-fatal event. The study followed the historical cohort design, with a 30-day period prevalence endpoint.

We evaluated outcome in two ways. First, we used a polytomous outcome model. Each of the four undesirable outcomes were modeled on predictor variables individually but simultaneously, to highlight differences in the multivariate effects. Second, we used an aggregate-endpoint analysis that utilized a single composite outcome variable to model the overall probability of a good patient outcome. The outcomes considered were renal failure, stroke, paraplegia/paraparesis, or death within 30 days of surgery. Patients without any of these complications were considered to have had a good short-term outcome after surgery.

Although the ultimate goal of the aggregate analysis was to determine the probability of a good outcome, we actually defined the response variable as ‘poor outcome’ so that the regression coefficients for the aggregate-variable analysis would run in the same direction as the coefficients for the polytomous categorical model. The polytomous model described the relationship between risk factors and presence (rather than absence) of specific complications. The variable ‘poor outcome’ was true if any censoring variable was true. Variables considered as predictor variable candidates are shown in Table 2. Variables that were significant in at least one dimension of the polytomous model were retained in the multiple logistic analysis as well. For the aggregate-outcome analysis, we used multiple logistic regression, with a single outcome variable defined as described above. Continuous predictor variables such as age were kept in their native distribution for multivariable analysis. Indicator variables were created for dichotomous predictor variables such as preoperative renal dysfunction, extent II thoracoabdominal aneurysm, etc. Computations were performed using the SAS LOGISTIC procedure [7].


View this table:
[in this window]
[in a new window]
 
Table 2. Selected univariate predictors of the aggregate outcome

 
For the polytomous-outcome categorical model, an ordinal variable ‘censor’ was defined that had the value 1 if censoring was by renal failure, 2 by neurologic deficit, 3, by stroke, 4 by death. The variable had the value 5 if no event occurred. The model was estimated using maximum-likelihood computations with the predictor variables specified as direct. That is, the predictor variables were kept in their native distribution as in the logistic analysis rather than being converted to categorical variables, as is the default for this procedure. Computations were performed using the SAS CATMOD procedure. An advantage of the polytomous approach is that parallelism is not assumed as it would be with an ordinal logistic regression model. Instead, this type of model assumes that data follow a product multinomial distribution as though they were cells in a contingency table made up of samples of observed responses based on the predictor variable levels. The vector of the resulting cell proportions is transformed into a linear function using maximum-likelihood estimates [8,9].


    3. Results
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
Eight hundred and forty-one patients had complete data for all risk factors and outcomes. All results shown in the tables are based on the group of 841. Outcome frequencies are shown in Table 1. We attempted to include postoperative respiratory failure as an outcome as well, but it was so closely correlated with other outcome variables that, after censoring on severity, no significant risk factors for respiratory failure could be identified. We therefore considered respiratory failure to be a redundant component of the composite variable and dropped it from the models.


View this table:
[in this window]
[in a new window]
 
Table 1. Outcome frequency (30-day)

 
Predictor variables retained in the multivariable models were age, preoperative renal dysfunction, extent II thoracoabdominal aneurysm, acute aortic dissection, and the adjunct distal aortic perfusion+cerebrospinal fluid drainage. Univariate relationships between selected predictor variables and the aggregate outcome variable are shown in Table 2. Significant univariate risk factors were age, renal dysfunction, and extent II aneurysm. Chronic aortic dissection was inversely associated with the composite outcome.

In the polytomous model, where all of the complications were modeled individually but simultaneously, preoperative renal dysfunction was a significant predictor for every outcome (Table 3). Age and extent II aneurysm were significant for all outcomes but stroke. Acute dissection and adjunct were significant only for paraplegia/paraparesis. In the logistic model of the aggregate outcome, preoperative renal dysfunction, age and extent II aneurysm were highly significant predictors, but adjunct and acute dissection were not significant (Table 4).


View this table:
[in this window]
[in a new window]
 
Table 3. Polytomous model simultaneous multiple outcomes

 

View this table:
[in this window]
[in a new window]
 
Table 4. Multiple logistic regression model predictors of aggregate outcome

 
The polytomous analysis, which modeled all the outcomes simultaneously but with individual results, preserved relationships between two important predictor variables – acute aortic dissection and use of adjuncts – and the risk of paraplegia/paraparesis. These effects were lost in the aggregate endpoint analysis.


    4. Discussion
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
An understanding of the probability of a favorable outcome after surgery may be more important to patients than the risk of individual untoward events. Aggregate-endpoint analysis is one way to get at this basic idea. Unfortunately, aggregate-endpoint data analysis alone may not be sensitive enough to give an accurate portrait of risk factors that are highly specific to a small proportion of a composite outcome. The predictive power of risk factors for low-incidence outcomes may be diluted by summarization with the other potential outcomes. In this study, two of the risk factors for postoperative paraplegia/paraparesis were lost in the aggregate analysis. If we had not looked beyond the aggregate analysis, we would have been fooled into thinking that these risk factors were unimportant for overall outcome. But the conclusion that acute dissection and adjunct do not affect complication rate is false. Therefore, a more powerful multivariate analysis tool was required.

Disparity among outcome prevalences in the multiple components of a polytomous endpoint – which, as we have noted, causes major problems for composite endpoint analysis – is not completely without consequences for polytomous models. In general, one risk factor can be identified reliably without overdetermination for about every 10 events. In the polytomous case, this means ten events per outcome class. Therefore, model selection should consider the number of events in each outcome variable category, and this should be balanced against biological plausibility and what is generally known about the risk factors, in addition to statistical significance criteria.

The traditional approach of examining multiple endpoints with multiple separate analyses can be another informative way to investigate the relationships between risk factors and individual outcomes, if several conditions are met. It is possible to use individual single dependent-variable logistic regression equations in combination to arrive at combined outcome probabilities that accumulate to 100%. This was described by Blackstone and reported by Kirklin and Barratt-Boyes in their classic text on cardiac surgery [10]. Hosmer and Lemeshow demonstrate the computations [11]. These segmental ordinary logistic models can be made to provide identical results to those obtained by polytomous analysis under appropriate conditions. Interested readers are referred to the excellent chapters in Refs. [10] and [11] for a thorough discussion. Our general feeling is that this formulation is a bit more cumbersome than the polytomous approach, but as we have said, the two methods can be made to be equivalent. As a straightforward addition to composite endpoint analysis, in any event, polytomous multivariate analysis is a useful addition to the outcome analysis tool kit.

In our view, it makes sense to do things in several stages – univariate, multivariate, aggregate – so that small effects can be scrutinized and a clearer view of the story contained in the underlying data can be had.

The major finding of this study is that preoperative renal dysfunction is an important risk factor for all the components of a poor outcome following thoracoabdominal aortic aneurysm repair, and that this is best observed after adjustments for effects of multiple risk factors and multiple outcomes have been made. This finding has helped to formalize our clinical impressions, and provides the rationale for further study of this important risk factor and for future development of strategies for management of renal dysfunction.


    Footnotes
 
Presented at the 16th Annual Meeting of the European Association for Cardio-thoracic Surgery, Monte Carlo, Monaco, September 22–25, 2002.


    Appendix A. Conference discussion
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 
Dr J. Roquette (Lisbon, Portugal): What do you consider renal dysfunction in these groups of patients?

Dr Miller: Preoperatively it is a history of hemodialysis or a preoperative serum creatinine greater than 2.

Dr Roquette: Do you think this polytomous outcome should be used instead of another type of evaluation in this setting?

Dr Miller: I think polytomous analysis is another useful tool in the tool kit. The composite-endpoint analysis is useful, but I think that the polytomous analysis should be done in addition, because it allows us to see how the outcomes separate from one another and how they interact.

Dr F. Schmid (Regensburg, Germany): This is a huge number of patients you have analyzed, and I assume that the technical procedure evolved over the years. Is there an impact of the period of operation on that model?

Dr Miller: Practice has indeed changed over time, and period effect does exist for this reason. The polytomous categorical model is a generalization of the linear model, and as such, adjustments for period effect can be made if necessary.

Dr R. Mazhar (Doha, Qatar): You identified adjunct as a negative factor for poor outcome. Is this not the use of CSF tapping and distal cooling?

Dr Miller: It is beneficial, yes.


    References
 Top
 Abstract
 1. Introduction
 2. Methods
 3. Results
 4. Discussion
 Appendix A. Conference...
 References
 

  1. Miller C.C., III, Reardon M.J., Safi H.J. Risk stratification: a practical guide for clinicians. Cambridge: Cambridge University Press, 2001.
  2. Miller C.C., III, Safi H.J., Vinnerkvist A., Baldwin J.C. Actual versus actuarial analysis for cardiac valve complications: the problem of competing risks. Curr Opin Cardiol 1999;14:79-83.[Medline]
  3. Grunkemeier G.L., Wu Y. Interpretation of nonfatal events after cardiac surgery: actual versus actuarial reporting. J Thorac Cardiovasc Surg 2001;122:216-219.[Abstract/Free Full Text]
  4. LeMaire S.A., Miller C.C., 3rd, Conklin L.D., Schmittling Z.C., Koksoy C., Coselli J.S. A new predictive model for adverse outcomes after elective thoracoabdominal aortic aneurysm repair. Ann Thorac Surg 2001;71:1233-1238.[Abstract/Free Full Text]
  5. Safi H.J., Vinnerkvist A., Miller C.C., 3rd, Iliopoulos D.C., Reardon M.J., Espada R., Baldwin J.C. Effect of extended cross-clamp time during thoracoabdominal aortic aneurysm repair. Ann Thorac Surg 1998;66:1204-1209.[Abstract/Free Full Text]
  6. Safi H.J., Campbell M.P., Miller C.C., 3rd, Iliopoulos D.C., Khoynezhad A.K., Lestou G.V., Asimacopoulos P.J. Cerebral spinal fluid drainage and distal aortic perfusion decrease the incidence of neurologic deficit: the results of 343 descending and thoracoabdominal aortic aneurysm repairs. Eur J Vasc Endovasc Surg 1997;14:118-124.[CrossRef][Medline]
  7. SAS Institute. The logistic procedure. . SAS/STAT user's guide, 4th ed Cary, NC: SAS Institute, 1990.
  8. SAS Institute. The CATMOD procedure. . SAS/STAT user's guide, 4th ed Cary, NC: SAS Institute, 1990.
  9. Bishop Y.M., Fienberg S.E., Holland P.W. Discrete multivariate analysis: theory and practice. Cambridge, MA: MIT Press, 1977.
  10. Kirklin J.W., Barratt-Boyes B. Surgical concepts, research methods and data analysis and use [Chapter 6]. Cardiac surgery. New York: Churchill Livingstone, 1986.
  11. Hosmer D., Lemeshow S. Applied logistic regression. New York: Wiley, 1989.



This article has been cited by other articles:


Home page
Br J AnaesthHome page
E. C. Lases, M. A. Schepens, F. J. Haas, L. P. Aarts, H. T. ter Beek, E. P. van Dongen, H. P. Siegers, I. van der Tweel, and E. H. Boezeman
Clinical prospective study of biochemical markers and evoked potentials for identifying adverse neurological outcome after thoracic and thoracoabdominal aortic aneurysm surgery
Br. J. Anaesth., November 1, 2005; 95(5): 651 - 661.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Add to Personal Folders
Right arrow Download to citation manager
Right arrow Author home page(s):
Eyal E. Porat
Anthony L. Estrera
Anders N. Vinnerkvist
Hazim J. Safi
Right arrow Permission Requests
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Miller, C. C.
Right arrow Articles by Safi, H. J.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Miller, C. C., III
Right arrow Articles by Safi, H. J.
Related Collections
Right arrow Education
Right arrow Great vessels


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
ANN THORAC SURG ASIAN CARDIOVASC THORAC ANN EUR J CARDIOTHORAC SURG
J THORAC CARDIOVASC SURG ICVTS ALL CTSNet JOURNALS