- Research article
- Open Access
- Open Peer Review
Front-line intraperitoneal versus intravenous chemotherapy in stage III-IV epithelial ovarian, tubal, and peritoneal cancer with minimal residual disease: a competing risk analysis
BMC Cancervolume 16, Article number: 235 (2016)
In the analysis of survival data for cancer patients, the problem of competing risks is often ignored. Competing risks have been recognized as a special case of time-to-event analysis. The conventional techniques for time-to-event analysis applied in the presence of competing risks often give biased or uninterpretable results.
Using a prospectively collected administrative health care database in a single institution, we identified patients diagnosed with stage III or IV primary epithelial ovarian, tubal, and peritoneal cancers with minimal residual disease after primary cytoreductive surgery between 1995 and 2012. Here, we sought to evaluate whether intraperitoneal chemotherapy outperforms intravenous chemotherapy in the presence of competing risks. Unadjusted and multivariable subdistribution hazards models were applied to this database with two types of competing risks (cancer-specific mortality and other-cause mortality) coded to measure the relative effects of intraperitoneal chemotherapy.
A total of 1263 patients were recruited as the initial cohort. After propensity score matching, 381 patients in each arm entered into final competing risk analysis. Cumulative incidence estimates for cancer-specific mortality were statistically significantly lower (p = 0.017, Gray test) in patients receiving intraperitoneal chemotherapy (5-year estimates, 34.5 %; 95 % confidence interval [CI], 29.5–39.6 %, and 10-year estimates, 60.7 %; 95 % CI, 52.2–68.0 %) versus intravenous chemotherapy (5-year estimates, 41.3 %; 95 % CI, 36.2–46.3 %, and 10-year estimates, 67.5 %, 95 % CI, 61.6–72.7 %). In subdistribution hazards analysis, for cancer-specific mortality, intraperitoneal chemotherapy outperforms intravenous chemotherapy (Subdistribution hazard ratio, 0.82; 95 % CI, 0.70–0.96) after correcting other covariates.
In conclusion, results from this comparative effectiveness study provide supportive evidence for previous published randomized trials that intraperitoneal chemotherapy outperforms intravenous chemotherapy even eliminating the confounding of competing risks. We suggest that implementation of competing risk analysis should be highly considered for the investigation of cancer patients who have medium to long-term follow-up period.
Ovarian cancer is the sixth most common tumor in women. More than 200,000 new cases are diagnosed each year worldwide. Each year, it constitutes 4 % of all cancers diagnosed in women, and there are 6.6 new cases per 100,000 women per year . Epithelial ovarian cancer takes the predominant 90 % of cases among ovarian cancer patients.
Epithelial ovarian cancer typically spreads by intraperitoneal seeding or direct invasion or through the lymphatic and vascular circulation. Among the spreading routes, peritoneal spreading is the most common route of dissemination, and stage III disease is associated at best with a 5-year survival rate of approximately 32–47 % . Therefore, it is reasonable to consider the intraperitoneal chemotherapy for this disease entity . Intraperitoneal chemotherapy has been investigated for several decades, and there have been three large-scale randomized trials conducted in the US, all of which showed overall and/or progression-free survival benefit [4–6].
Most published studies for ovarian cancer use the time to some disease events as their primary outcome and hence, statistical methods developed for survival data are usually applied. Established methods for estimating and modelling these include the Kaplan–Meier estimator of the survival function and the Cox proportional hazards model for the hazard function [7, 8] An important assumption of these established survival analytical methods is that censoring is ‘independent’ . However, in some cases, several causes of failure are possible but the occurrence of one event precludes the occurrence of the other events (i.e., when failures are different causes of death, only the first one can be observed). This situation is known as competing risks. In a competing risk situation, standard techniques for survival analysis may lead to incorrect and biased results [10, 11]. In usual condition, ovarian cancer often presents a protracted disease course, and it is not uncommon to see a patient dies of other causes (e.g., heart failure and stroke), which precludes the occurrence of cancer-specific death.
In the current work, we conducted a competing risk analysis to investigate the therapeutic effects of intraperitoneal chemotherapy on stage III-IV epithelial ovarian, tubal, and peritoneal cancer with minimal residual disease using an administrative health care database constructed in a single tertiary care institution.
The study entailed a retrospective analysis of prospectively collected demographic, tumor profile, treatment, and comorbidity data, linking electronic data sources, including the cancer registry, administrative and clinical databases, and surgery records of consecutive patients with stage III or IV primary epithelial ovarian, tubal, and peritoneal cancers with minimal residual disease after primary cytoreductive surgery between January 1995 and December 2012. The standard patient informed consent for retrieval of personal information at the institution included an emphasized section describing the purpose of the study. This study was approved by the Institutional Review Board of Taipei Veterans General Hospital. The procedures used in this study were in accordance with the guidelines of the Helsinki Declaration on human experimentation.
Primary cytoreductive surgery and front-line chemotherapy
The standard cytoreductive surgery included total hysterectomy, bilateral salpingo-oophorectomy, infracolic omentectomy, and cytoreduction of all tumor nodules to a size of 1 cm or less in the greatest dimension.
Courses of front-line chemotherapy were repeated every 3 weeks for a total of six cycles, provided the serum creatinine concentration was less than or equal to 2.0 mg/dl, the white-cell count was higher than 3000/mm3, and the platelet count was higher than 80,000/mm3. Dosing schedules for intraperitoneal chemotherapy were either platinum-based or taxane-based. For platinum-based intraperitoneal chemotherapy, cisplatin at 100 mg/m2 or carboplatin at AUC 5 or 6 via Tenckhoff tubes was administered. For taxane-based intraperitoneal chemotherapy, protocol for GOG 172 was followed.6
For intravenous chemotherapy, both taxane (175 mg/m2) plus cisplatin (dosed at 50 mg/m2) or carboplatin (AUC = 5 or 6) were administered every 3 weeks for a total of six cycles.
In each cycle of chemotherapy, patients received physical examination, complete blood count, biochemical profiles, CA-125, and 24-h urine collection for measurement of clearance of creatinine. In the absence of clinical evidence of tumor progression, tumor evaluation by image studies including chest film, whole abdominal sonography, and CT scan (or MRI) was performed after six cycles of chemotherapy.
Continuous variables are presented as mean (± standard deviations) and were compared using Student’s unpaired t-test. Categorical variables are presented as counts and percentages and were compared with the χ2 test when appropriate (expected frequency > 5). Otherwise, Fisher’s exact test was used. Overall survival duration was calculated from disease diagnosis to event occurrence or last follow-up. Patients who succumbed to cancer-related death were classified as cancer-specific mortality (CSM), while patients who succumbed to other causes were classified as other-cause mortality (OCM).
Our statistical analyses consisted of two steps. In the first step, we attempted to adjust for the selection bias inherent in observational data by applying propensity score matching to balance the measured covariates between the intravenous and intraperitoneal chemotherapy group. The propensity score is a summary confounder score that is modeled using the exposure or treatment as the dependent variable . The propensity to receive intraperitoneal chemotherapy was calculated using a multivariable logistic regression model that adjusted to age at diagnosis, FIGO stage, tumor grade, histology, Charlson comorbidity index, and GOG performance score. We used the nearest neighbor matching with a caliper width of 0.2 of the standardized deviation of the logit to match cases. This optimizes the matching with minimal residual bias and highest precision . Covariate balance was evaluated using standardized differences of means (SDM), with SDM of < 0.1 (corresponding to <10 % difference between the arms) indicative of acceptable balance .
In the second step of analyses, competing risks analyses were conducted based on the propensity score-matched cohort. The crude cumulative incidence function were estimated for each type of competing risk using the method of Kalbfleisch and Prentice [15, 16] Cumulative incidence plots were used to graphically depict CSM and OCM rates. Statistical significance of differences in survival rates was assessed with the Gray test  Furthermore, we investigated differences in each cause of mortality (for CSM and OCM, respectively) using subdistribution hazard ratios estimated through Fine and Gray proportional hazards regression  Censoring time was set at December 31, 2012. The dataset supporting the conclusions of this article is listed in Additional file 1.
Finally, we performed a sensitivity analysis based on different sample size when extreme values of the propensity score were trimmed at different levels - a procedure shown to partly compensate for unobserved confounding .
All tests were performed two sided at the 5 % significance level. Statistical analyses were performed with IBM SPSS Statistics (version 20.0, IBM, Armonk, NY), and R software (version 2.15, R Foundation for Statistical Computing, Vienna, Austria) using the cmprsk, survival, and Matching packages .
A total of 1263 patients were recruited as the initial cohort. The patients’ baseline clinical characteristics of the initial cohort are summarized in the left part of Table 1. In terms of standardized difference of means, there were significant differences in age, stage, grade, histologic subtype, proportion of front-line regimen, GOG performance status, and Charlson comorbidity score. However, after propensity score-matching (summarized at right part of Table 1), the standardized difference of means in all the baseline clinical characteristics was less than 10 %, which indicates a high degree of similarity in the distribution of these covariates, with equal number of 381 patients receiving intravenous chemotherapy and intraperitoneal chemotherapy, respectively. The propensity score-matched cohort forms the basis for the following analyses. The median follow-up time for this cohort was 8.4 years.
Figure 1 shows the estimates and curves of cumulative incidence function stratified by CSM and OCM, and sub-stratified by type of chemotherapy (intravenous vs. intraperitoneal chemotherapy). Cumulative incidence estimates for CSM were statistically significantly lower (p = 0.017, Gray test) in patients receiving intraperitoneal chemotherapy (5-year estimates, 34.5 %; 95 % confidence interval [CI], 29.5–39.6 %, and 10-year estimates, 60.7 %; 95 % CI, 52.2–68.0 %) versus intravenous chemotherapy (5-year estimates, 41.3 %; 95 % CI, 36.2–46.3 %, and 10-year estimates, 67.5 %, 95 % CI, 61.6–72.7 %). In contrast, cumulative incidence estimates for OCM were comparable (p = 0.705, Gray test) between patients receiving intraperitoneal chemotherapy (5-year estimates, 6.8 %; 95 % CI, 4.4–9.8 %, and 10-year estimates, 21.8 %, 95 % CI, 16.5–27.5 %) and receiving intravenous chemotherapy (5-year estimate, 4.5 %; 95 % CI, 2.9–7.6 %, and 10-year estimate, 26.7 %; 95 % CI, 20.6–32.6 %).
Because ovarian cancer is now considered as a heterogeneous disease in which histologic phenotypes correlate with distinct genetic events, therefore we present cumulative incidence estimates for each histologic subtypes. For each histologic subtype, there remains a significant difference for CSM between intravenous chemotherapy and intraperitoneal chemotherapy. However, for OCM, there is no difference between intravenous chemotherapy and intraperitoneal chemotherapy (Fig. 2).
In conventional survival analysis, application of Cox proportional hazards regression uses multivariable hazard ratios can account for confounding. Nonetheless, Cox regression may lead to biased effect estimates in the presence of competing risks . Hence, we next conducted unadjusted and adjusted regression modelling proposed by Fine and Gray in the presence of competing risks . Unadjusted and adjusted subdistribution hazard ratios for CSM and OCM are given in Table 2. For CSM, several factors demonstrate statistically significant which include stage, grade, histology, and type of chemotherapy in the univariable analysis. However, after multivariable regression modelling, only some factors remain as an independent risk factor, including stage (IV vs. III, subdistribution hazard ratio [SDHR], 2.66; 95 % CI, 1.47–4.81]), grade (grade 3 vs. grade 1, SDHR, 1.89; 95 % CI, 1.04–3.43), histology (clear cell vs. serous, SDHR, 1.68; 95 % CI, 1.14–2.48, and others vs. serous, SDHR, 2.36; 95 % CI, 1.67–3.34), and type of chemotherapy (intraperitoneal vs. intravenous chemotherapy, SDHR, 0.82; 95 % CI, 0.70–0.96). Hence intraperitoneal chemotherapy remains a significant factor after multivariable regression modelling in competing risk analysis.
For OCM, two factors demonstrate statistically significant which include age and Charlson comorbidity index in the univariable analysis. After multivariable regression modelling, these two factors still remain as an independent factor: age (per 10 years increment, SDHR, 092; 95 % CI, 0.87–0.97) and Charlson comorbidity index (≥1 vs. 0, SDHR, 0.79; 95 % CI, 0.66–0.94).
Lastly we conducted a sensitivity analysis for CSM. Consistent results of subdistribution hazard ratio of intraperitoneal chemotherapy vs. intravenous chemotherapy were obtained when extreme values of the propensity score were trimmed at different levels − a procedure shown to partly compensate for unobserved confounding (Table 3).
This competing risks analyses based on a single-institution database of patients with stage III-IV epithelial ovarian, tubal, and peritoneal cancer with minimal residual disease demonstrate better survival outcome treated by intraperitoneal chemotherapy compared to intravenous chemotherapy with respect to CSM. Intraperitoneal chemotherapy has been investigated for a long time, and based on the positive results from a meta-analysis conducted by National Cancer Institute (NCI) and the Gynecologic Oncology Group, the NCI has released a clinical announcement encouraging the gynecological oncology community to consider intraperitoneal chemotherapy using cisplatin as the standard treatment for advanced ovarian cancer patients in whom the residual disease were debulked to 1 cm or less . The results of the current work can further consolidate the role of intraperitoneal chemotherapy in the treatment of advanced ovarian cancer.
Competing risks methods are common in biomedical research, particularly in cancer, where the need to deal with multiple potential outcomes is nearly ubiquitous. In fact, the concept of competing risks within clinical research was first introduced in the field of oncology . As treatment for cancer produced prolonged survival times, it became important to consider not only the effects of treatment on cancer-free survival, but also how competing risks, such as mortality from unrelated causes, might impact treatment decisions. For example, competing risks methods have been applied in the analyses of risk factors for breast cancer, either at the screening or treatment stage that may potentially impact the decision making [24–27].
Traditionally, when predict the unadjusted probability of a certain event of interest to occur, one can use the Kaplan–Meier (KM) method. However, in the presence of competing risks, using the KM method is problematic. This method can handle only one single event at one time, and all other events are treated as censored observations. Further, the complement of the KM estimate (1 − KM) is interpreted as the cumulative probability of the event of interest in a hypothetical world where no subject would experience the competing event. This kind of interpretation is not realistic in clinical practice .
While for the adjusted analysis of competing risk, Fine and Gray proposed a regression modeling applied directly on a cumulative incidence function for particular use in the competing risk analysis. For any event type, this approach focuses on the hazard associated with the cumulative incidence function.
Several limitations merit consideration in the current work. First, the composition of the initial cohort may be influenced by referral bias because of the monocentric design. The index hospital is a tertiary referral center with more patients harboring co-morbidities. Thus, the results of the present study may not be generalizable to the general population. Second, lack of laboratory data may influence the robustness of the results. For example, serum cancer antigen 125 (CA125) is widely used in ovarian cancer to monitor the effectiveness of therapy both in first line and recurrence [29, 30]. Lack of inclusion of CA125 into the calculation of propensity-score may `bias our results. Third, limitations in patient-level data collection were also present. For example, socioeconomic status have been found to be associated with cancer mortality  yet, our database did not capture this important factor which potentially affects the reliability of both propensity-score and competing risk modelling. Fourth, the protocol of follow-up was not strictly defined and censoring recordings may not be reflective of the true status. Fifth, the current analysis have not been validated with an external dataset which limit its clinical application.
Competing risks analyses allow disentangling the contribution of an intervention (e.g., intraperitoneal chemotherapy) on each event type separately. Intraperitoneal chemotherapy clearly shows therapeutic benefit in terms of CSM, but not OCM in the current work. In addition to intraperitoneal chemotherapy, stage IV (vs. III), grade 3 (vs. grade 1), and clear cell subtype (vs. serous subtype) also demonstrate as significant risk factors. Clear cell carcinoma accounts for 4 to 12 % of epithelial ovarian cancer in Western countries. Compared to serous adenocarcinoma, clear cell carcinoma is relatively resistant to conventional platinum, or taxane-based chemotherapy which is associated with its poor prognosis . The results of our work suggest that novel therapy should be developed for clear cell carcinoma in order to improve survival outcome.
From published reports, the mean or median age for ovarian cancer is between 55 and 60 years which implies that a significant proportion of patients are facing the problem of aging. During the course of treatment for cancer, the increased age and co-morbidities of older patients occasionally lead to death by the occurrence of competing events (e.g., stroke, infection, and accident). Clinical studies are often faced with the difficult problem of how to account for participants who die without experiencing the study outcome of interest. Conventional approaches to describe risk of disease like Kaplan-Meier survival analysis and Cox proportional hazards regression can overestimate risk of disease by failing to account for the competing risk of death .
In conclusion, results from this competing risk analysis provide supportive evidence for previous published randomized trials that intraperitoneal chemotherapy outperforms intravenous chemotherapy. We propose that implementation of competing risk analysis should be highly considered for the investigation of ovarian cancer patients who have medium to long-term follow-up period.
Jemal A, Siegel R, Xu J, Ward E. Cancer statistics, 2010. CA Cancer J Clin. 2010;60(5):277–300.
Ferlay J, Soerjomataram I, Dikshit R, Eser S, Mathers C, Rebelo M, et al. Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer. 2015;136(5):E359–86.
Fujiwara K, Armstrong D, Morgan M, Markman M. Principles and practice of intraperitoneal chemotherapy for ovarian cancer. Int J Gynecol Cancer. 2007;17(1):1–20.
Alberts DS, Liu PY, Hannigan EV, O'Toole R, Williams SD, Young JA, et al. Intraperitoneal cisplatin plus intravenous cyclophosphamide versus intravenous cisplatin plus intravenous cyclophosphamide for stage III ovarian cancer. N Engl J Med. 1996;335(26):1950–5.
Markman M, Bundy BN, Alberts DS, Fowler JM, Clark-Pearson DL, Carson LF, et al. Phase III trial of standard-dose intravenous cisplatin plus paclitaxel versus moderately high-dose carboplatin followed by intravenous paclitaxel and intraperitoneal cisplatin in small-volume stage III ovarian carcinoma: an intergroup study of the Gynecologic Oncology Group, Southwestern Oncology Group, and Eastern Cooperative Oncology Group. J Clin Oncol. 2001;19(4):1001–7.
Armstrong DK, Bundy B, Wenzel L, Huang HQ, Baergen R, Lele S, et al. Intraperitoneal cisplatin and paclitaxel in ovarian cancer. N Engl J Med. 2006;354(1):34–43.
Kalbfleisch JD, Prentice RL. The statistical analysis of failure time data. 2nd ed. Hoboken, New Jersey; John Wiley & Sons. Inc.; 2002.
Collett D. Modelling survival data in medical research: CRC press. 2003.
Lau B, Cole SR, Gange SJ. Competing risk regression models for epidemiologic data. Am J Epidemiol. 2009;170(2):244–56.
Helton JC, Oberkampf WL, Johnson JD. Competing failure risk analysis using evidence theory. Risk Anal. 2005;25(4):973-95
Southern DA, Faris PD, Brant R, Galbraith PD, Norris CM, Knudtson ML, et al. Kaplan-Meier methods yielded misleading results in competing risk scenarios. J Clin Epidemiol. 2006;59(10):1110–4.
Hullsiek KH, Louis TA. Propensity score modeling strategies for the causal analysis of observational data. Biostatistics. 2002;3(2):179–93.
Wang Y, Cai H, Li C, Jiang Z, Wang L, Song J, et al. Optimal caliper width for propensity score matching of three treatment groups: a Monte Carlo study. PLoS One. 2013;8(12), e81045.
Austin PC. Balance diagnostics for comparing the distribution of baseline covariates between treatment groups in propensity-score matched samples. Stat Med. 2009;28(25):3083–107.
Gooley TA, Leisenring W, Crowley J, Storer BE. Estimation of failure probabilities in the presence of competing risks: new representations of old estimators. Stat Med. 1999;18(6):695–706.
Putter H, Fiocco M, Geskus RB. Tutorial in biostatistics: competing risks and multi-state models. Stat Med. 2007;26(11):2389–430.
Gray RJ. A class of K-sample tests for comparing the cumulative incidence of a competing risk. Ann Stat. 1988;16(3):1141–54.
Fine JP, Gray RJ. A proportional hazards model for the subdistribution of a competing risk. J Am Stat Assoc. 1999;94(446):496–509.
Sturmer T, Rothman KJ, Avorn J, Glynn RJ. Treatment effects in the presence of unmeasured confounding: dealing with observations in the tails of the propensity score distribution--a simulation study. Am J Epidemiol. 2010;172(7):843–54.
Scrucca L, Santucci A, Aversa F. Regression modeling of competing risk using R: an in depth guide for clinicians. Bone Marrow Transplant. 2010;45(9):1388–95.
Satagopan JM, Ben-Porat L, Berwick M, Robson M, Kutler D, Auerbach AD. A note on competing risks in survival data analysis. Br J Cancer. 2004;91(7):1229–35.
NCI Clinical Announcement on Intraperitoneal Chemotherapy in Ovarian Cancer [http://ctep.cancer.gov/highlights/20060105_ovarian.htm] Accessed 2 December 2014.
Kalbfleisch JD, Prentice RL. Estimation of the average hazard ratio. Biometrika. 1981;68(1):105–12.
Harris HR, Bergkvist L, Wolk A. An estrogen-associated dietary pattern and breast cancer risk in the Swedish Mammography Cohort. International journal of cancer Journal international du cancer. 2015;137(9):2149–54.
Darabi H, Czene K, Zhao W, Liu J, Hall P, Humphreys K. Breast cancer risk prediction and individualised screening based on common genetic variation and breast density measurement. Breast Cancer Res. 2012;14(1):R25.
Tice JA, Miglioretti DL, Li CS, Vachon CM, Gard CC, Kerlikowske K. Breast Density and Benign Breast Disease: Risk Assessment to Identify Women at High Risk of Breast Cancer. J Clin Oncol. 2015;33(28):3137–43.
Vilaprinyo E, Gispert R, Martinez-Alonso M, Carles M, Pla R, Espinas JA, Rue M. Competing risks to breast cancer mortality in Catalonia. BMC Cancer. 2008;8:331.
Varadhan R, Weiss CO, Segal JB, Wu AW, Scharfstein D, Boyd C. Evaluating health outcomes in the presence of competing risks: a review of statistical methods and clinical applications. Med Care. 2010;48(6 Suppl):S96–105.
Pignata S, Cannella L, Leopardo D, Bruni GS, Facchini G, Pisano C. Follow-up with CA125 after primary therapy of advanced ovarian cancer: in favor of continuing to prescribe CA125 during follow-up. Ann Oncol. 2011;22 Suppl 8:viii40–4.
Crawford SM, Peace J. Does the nadir CA125 concentration predict a long-term outcome after chemotherapy for carcinoma of the ovary? Ann Oncol. 2005;16(1):47–50.
Albano JD, Ward E, Jemal A, Anderson R, Cokkinides VE, Murray T, et al. Cancer mortality in the United States by education level and race. J Natl Cancer Inst. 2007;99(18):1384–94.
Kajihara H, Yamada Y, Kanayama S, Furukawa N, Noguchi T, Haruta S, et al. Clear cell carcinoma of the ovary: potential pathogenic mechanisms (Review). Oncol Rep. 2010;23(5):1193–203.
Verduijn M, Grootendorst DC, Dekker FW, Jager KJ, le Cessie S. The analysis of competing events like cause-specific mortality--beware of the Kaplan-Meier method. Nephrol Dial Transplant. 2011;26(1):56–61.
We are indebted to the research panel: Lilly Wen, MD (Catholic Cardinal Tien Hospital, Taipei), Nan-Ni Chen, MD (Municipal Hospital for Women and Children, Taipei), and Kuo-Chang Wen, MD, Wei-Min Hu, MD, Hsiao-Wen Tsai, MD, Pi-Lin Sung, MD, Wei-Lun Hsu, MD, Chih-Yu Chen, MD, Jen-Yu Huang, MD, Chia-Ming Chang, MD, Peng-Hui Wang, MD, Nae-Fong Twu, MD, Hsiang-Tai Chao, MD, and Koun-Chung Chao, MD (Veterans General Hospital, Taipei).
Sources of support
This study was supported in part by National Science Council grant (NSC 104-NU-E-010-004-NU and MOST 103-2410-H-264-004).
The authors declare that they have no competing interests.
All authors contributed substantially to the conception and design of the study. YHC WPC and CMC contributed with acquisition and interpretation of data, reviewed the first and subsequent drafts of the article and revised it critically for intellectual content, and approved the final version before submission. WHL, YC, CWP, and CHC contributed with the collection of data, interpretation of the results, drafted the article, revised it and approved the final version before submission. All authors have read and approved the manuscript.
Dataset for computing competing risks. (CSV 18 kb)