‘Modeling and Analysis of Competing Risks Data’

there has been much research on analyzing various forms of competing risks data. Nevertheless, there are several occasions in survival studies, where the existing models and methodologies are inadequate for the analysis competing risks data. ldentifiabilty problem and various types of and censoring induce more complications in the analysis of competing risks data than in classical survival analysis. Parametric models are not adequate for the analysis of competing risks data since the assumptions about the underlying lifetime distributions may not hold well. Motivated by this, in the present study. we develop some new inference procedures, which are completely distribution free for the analysis of competing risks data.

Longitudinal outcomes of high-risk human papillomavirus (HPV) infections as competing-risks events following cervical HPV test at baseline visit in the *NIS-LAMS** cohort

Background: The complex natural history of human papillomavirus (HPV) infections following a single HPV test can be modeled as competing-risks events (i.e., no-, transient- or persistent infection) in a longitudinal setting. The covariates associated with these compet ng events have not been previously assessed using competing-risks regression models. Objectives: To gain further insights in the outcomes of cervical HPV infections, we used univariate- and multivariate competing-risks regression models to assess the covariaies associated with these competing events. Study Design and Methods: Covariates associated with three competing outcomes (no-, transient- or persistent HR-HPV infection) were analysed in a sub-cohort of 1,865 women prospectively followed-up in the NIS (n = 3,187) and LAMS Study (n = 12,114). Results: In multivariate competing-risks models (with two other outcomes as competing events), permanently HR-HPV negative outcome was significantly predicted only by the clearance of ASCUS+Pap during FU, while three independent covariates predicted transient HR-HPV infections: i) number of recent (< 12 months) sexual partners (risk increased), ii) previous Pap screening history (protective), and history of previous CIN (increased risk). The two most powerful predictors of persistent HR-HPV infections were persistent ASCUS+Pap (risk increased), and previous Pap screening history (protective). In pair-wise comparisons, number of recent sexual partners and previous CIN history increase the probability of transient HR-HPV infection against the HR-HPV negative competing event, while previous Pap screening history is protective. Persistent ASCUS+Pap during FU and no previous Pap screening history are significantly associated with the persistent HR-HPV outcome (compared both with i) always negative, and ii) transient events), whereas multiparity is protective. Conclusions: Different covariates are associated with the three main outcomes of cervical HPV infections. The most significant covariates of each competing events are probably distinct enough to enable constructing of a risk-profile for each main outcome.

Covariates of high-risk human papillomavirus (HPV) infections are distinct for incident CIN1, CIN2 and CIN3 as disclosed by competing-risks regression models

Background: In addition to the oncogenic human papillomavirus (HPV), several cofactors are needed in cervical carcinogenesis, but whether the HPV covariates associated with incident i) CIN1 are different from those of incident ii) CIN2 and iii) CIN3 needs further assessment. Objectives: To gain further insights into the true biological differences between CIN1, CIN2 and CIN3, we assessed HPV covariates associated with incident CIN1, CIN2, and CIN3. Study Design and Methods: HPV covariates associated with progression to CIN1, CIN2 and CIN3 were analysed in the combined cohort of the NIS (n = 3,187) and LAMS study (n = 12,114), using competing-risks regression models (in panel data) for baseline HR-HPV-positive women (n = 1,105), who represent a sub-cohort of all 1,865 women prospectively followed-up in these two studies. Results: Altogether, 90 (4.8%), 39 (2.1%) and 14 (1.4%) cases progressed to CIN1, CIN2, and CIN3, respectively. Among these baseline HR-HPV-positive women, the risk profiles of incident GIN I, CIN2 and CIN3 were unique in that completely different HPV covariates were associated with progression to CIN1, CIN2 and CIN3, irrespective which categories (non-progression, CIN1, CIN2, CIN3 or all) were used as competing-risks events in univariate and multivariate models. Conclusions: These data confirm our previous analysis based on multinomial regression models implicating that distinct covariates of HR-HPV are associated with progression to CIN1, CIN2 and CIN3. This emphasises true biological differences between the three grades of GIN, which revisits the concept of combining CIN2 with CIN3 or with CIN1 in histological classification or used as a common end-point, e.g., in HPV vaccine trials.

On modifications to the long-term survival mixture model in the presence of competing risks

A mixture model for long-term survivors has been adopted in various fields such as biostatistics and criminology where some individuals may never experience the type of failure under study. It is directly applicable in situations where the only information available from follow-up on individuals who will never experience this type of failure is in the form of censored observations. In this paper, we consider a modification to the model so that it still applies in the case where during the follow-up period it becomes known that an individual will never experience failure from the cause of interest. Unless a model allows for this additional information, a consistent survival analysis will not be obtained. A partial maximum likelihood (ML) approach is proposed that preserves the simplicity of the long-term survival mixture model and provides consistent estimators of the quantities of interest. Some simulation experiments are performed to assess the efficiency of the partial ML approach relative to the full ML approach for survival in the presence of competing risks.

An EM-based Semi-Parametric Mixture Model Approach to the Regression Analysis of Competing-Risks Data

We consider a mixture model approach to the regression analysis of competing-risks data. Attention is focused on inference concerning the effects of factors on both the probability of occurrence and the hazard rate conditional on each of the failure types. These two quantities are specified in the mixture model using the logistic model and the proportional hazards model, respectively. We propose a semi-parametric mixture method to estimate the logistic and regression coefficients jointly, whereby the component-baseline hazard functions are completely unspecified. Estimation is based on maximum likelihood on the basis of the full likelihood, implemented via an expectation-conditional maximization (ECM) algorithm. Simulation studies are performed to compare the performance of the proposed semi-parametric method with a fully parametric mixture approach. The results show that when the component-baseline hazard is monotonic increasing, the semi-parametric and fully parametric mixture approaches are comparable for mildly and moderately censored samples. When the component-baseline hazard is not monotonic increasing, the semi-parametric method consistently provides less biased estimates than a fully parametric approach and is comparable in efficiency in the estimation of the parameters for all levels of censoring. The methods are illustrated using a real data set of prostate cancer patients treated with different dosages of the drug diethylstilbestrol. Copyright (C) 2003 John Wiley Sons, Ltd.

Health investment complementarities under competing risks

Applying the competing--risks model to multi--cause mortality,this paper provides a theoretical and empirical investigation of the positive complementarities that occur between disease--specific policy interventions. We argue that since an individual cannot die twice, competing risks imply that individuals will not waste resources on causes that are not the most immediate, but will make health investments so as to equalize cause--specific mortality. However, equal mortality risk from a variety of diseases does not imply that disease--specific public health interventions are a waste. Rather, a cause--specific intervention produces spillovers to other disease risks, so that the overall reduction in mortality will generally be larger than the direct effect measured on the targeted disease. The assumption that mortality from non--targeted diseases remains the same after a cause--specific intervention under--estimates the true effect of such programs, since the background mortality is also altered as a result of intervention. Analyzing data from one of the most important public health programs ever introduced, the Expanded Program on Immunization (EPI) of the United Nations, we find evidence for the existence of such complementarities, involving causes that are not biomedically, but behaviorally, linked.

Competing risks to breast cancer mortality in Catalonia

Background: Breast cancer mortality has experienced important changes over the last century. Breast cancer occurs in the presence of other competing risks which can influence breast cancer incidence and mortality trends. The aim of the present work is: 1) to assess the impact of breast cancer deaths among mortality from all causes in Catalonia (Spain), by age and birth cohort and 2) to estimate the risk of death from other causes than breast cancer, one of the inputs needed to model breast cancer mortality reduction due to screening or therapeutic interventions. Methods: The multi-decrement life table methodology was used. First, all-cause mortality probabilities were obtained by age and cohort. Then mortality probability for breast cancer was subtracted from the all-cause mortality probabilities to obtain cohort life tables for causes other than breast cancer. These life tables, on one hand, provide an estimate of the risk of dying from competing risks, and on the other hand, permit to assess the impact of breast cancer deaths on all-cause mortality using the ratio of the probability of death for causes other than breast cancer by the all-cause probability of death. Results: There was an increasing impact of breast cancer on mortality in the first part of the 20th century, with a peak for cohorts born in 1945–54 in the 40–49 age groups (for which approximately 24% of mortality was due to breast cancer). Even though for cohorts born after 1955 there was only information for women under 50, it is also important to note that the impact of breast cancer on all-cause mortality decreased for those cohorts. Conclusion: We have quantified the effect of removing breast cancer mortality in different age groups and birth cohorts. Our results are consistent with US findings. We also have obtained an estimate of the risk of dying from competing-causes mortality, which will be used in the assessment of the effect of mammography screening on breast cancer mortality in Catalonia.

Nonparametric Estimation of Survivor Function in Bivariate Competing Risk Models’

So far, in the bivariate set up, the analysis of lifetime (failure time) data with multiple causes of failure is done by treating each cause of failure separately. with failures from other causes considered as independent censoring. This approach is unrealistic in many situations. For example, in the analysis of mortality data on married couples one would be interested to compare the hazards for the same cause of death as well as to check whether death due to one cause is more important for the partners’ risk of death from other causes. In reliability analysis. one often has systems with more than one component and many systems. subsystems and components have more than one cause of failure. Design of high-reliability systems generally requires that the individual system components have extremely high reliability even after long periods of time. Knowledge of the failure behaviour of a component can lead to savings in its cost of production and maintenance and. in some cases, to the preservation of human life. For the purpose of improving reliability. it is necessary to identify the cause of failure down to the component level. By treating each cause of failure separately with failures from other causes considered as independent censoring, the analysis of lifetime data would be incomplete. Motivated by this. we introduce a new approach for the analysis of bivariate competing risk data using the bivariate vector hazard rate of Johnson and Kotz (1975).

Competing regression models for longitudinal data

The choice of an appropriate family of linear models for the analysis of longitudinal data is often a matter of concern for practitioners. To attenuate such difficulties, we discuss some issues that emerge when analyzing this type of data via a practical example involving pretestposttest longitudinal data. In particular, we consider log-normal linear mixed models (LNLMM), generalized linear mixed models (GLMM), and models based on generalized estimating equations (GEE). We show how some special features of the data, like a nonconstant coefficient of variation, may be handled in the three approaches and evaluate their performance with respect to the magnitude of standard errors of interpretable and comparable parameters. We also show how different diagnostic tools may be employed to identify outliers and comment on available software. We conclude by noting that the results are similar, but that GEE-based models may be preferable when the goal is to compare the marginal expected responses.

Impact of age and comorbidities on long-term survival of patients with high-risk prostate cancer treated with radical prostatectomy: a multi-institutional competing-risks analysis

Survival after surgical treatment using competing-risk analysis has been previously examined in patients with prostate cancer (PCa). However, the combined effect of age and comorbidities has not been assessed in patients with high-risk PCa who might have heterogeneous rates of competing mortality despite the presence of aggressive disease.

Semiparametric Methods for Semi-competing Risks Problem with Censoring and Truncation

Studies of chronic life-threatening diseases often involve both mortality and morbidity. In observational studies, the data may also be subject to administrative left truncation and right censoring. Since mortality and morbidity may be correlated and mortality may censor morbidity, the Lynden-Bell estimator for left truncated and right censored data may be biased for estimating the marginal survival function of the non-terminal event. We propose a semiparametric estimator for this survival function based on a joint model for the two time-to-event variables, which utilizes the gamma frailty specification in the region of the observable data. Firstly, we develop a novel estimator for the gamma frailty parameter under left truncation. Using this estimator, we then derive a closed form estimator for the marginal distribution of the non-terminal event. The large sample properties of the estimators are established via asymptotic theory. The methodology performs well with moderate sample sizes, both in simulations and in an analysis of data from a diabetes registry.

The Impact of income support programs on transitions to employment

The income support programs are created with the purpose of fighting both, the poverty trap and the inactivity trap. The balance between both is fragile and hard to find. Thus, the goal of this work is to contribute to solve this issue by finding how income support programs, particularly the Portuguese RSI, affect transitions to employment. This is made through duration analysis, namely using Cox and Competing Risks models. A particular feature is introduced in this work as it incorporates the possibility of Defective Risks. The estimated hazard elasticity with respect to the amount of RSI received for individuals who move to employment is -0,41. More than a half of RSI receivers stays for more than a year and the probability of never leaving to employment is 44%. The results appear to indicate that RSI has affected negatively transitions to employment.

Adjuvant cisplatin-based combined chemotherapy for lymph node (LN)-positive urothelial carcinoma of the bladder (UCB) after radical cystectomy (RC): a retrospective international study of >1500 patients.

OBJECTIVE: To compare outcomes of patients with lymph node (LN)-positive urothelial carcinoma of the bladder (UCB) treated with or without cisplatin-based combined adjuvant chemotherapy (AC) after radical cystectomy (RC). PATIENTS AND METHODS: We retrospectively analysed 1523 patients with LN-positive UCB, who underwent RC with bilateral pelvic LN dissection. All patients had no evidence of disease after RC. AC was administered within 3 months. Competing-risks models were applied to compare UCB-related mortality. RESULTS: Of the 1523 patients, 874 (57.4%) received AC. The cumulative 1-, 2- and 5-year UCB-related mortality rates for all patients were 16%, 36% and 56%, respectively. Administration of AC was associated with an 18% relative reduction in the risk of UCB-related death (subhazard ratio 0.82, P = 0.005). The absolute reduction in mortality was 3.5% at 5 years. The positive effect of AC was detectable in patients aged ≤70 years, in women, in pT3-4 disease, and in those with a higher LN density and lymphovascular invasion. This study is limited by its retrospective and non-randomised design, selection bias, the absence of central pathological review and lack in standardisation of LN dissection and cisplatin-based protocols. CONCLUSION: AC seems to reduce UCB-related mortality in patients with LN-positive UCB after RC. Younger patients, women and those with high-risk features such as pT3-4 disease, a higher LN density and lymphovascular invasion appear to benefit most. Appropriately powered prospective randomised trials are necessary to confirm these findings.