• The Korean Journal of Applied Statistics
• /
• v.29 no.7
• /
• pp.1231-1246
• /
• 2016

Cause-specific hazard model (Prentice et al., 1978) and subdistribution hazard model (Fine and Gray, 1999) are mostly used for the right censored survival data with competing risks. Some other models for survival data with competing risks have been subsequently introduced; however, those models have not been popularly used because the models cannot provide reliable statistical estimation methods or those are overly difficult to compute. We introduce simple and reliable competing risk regression models which have been recently proposed as well as compare their methodologies. We show how to use SAS and R for the data with competing risks. In addition, we analyze survival data with two competing risks using five different models.

• Communications for Statistical Applications and Methods
• /
• v.25 no.4
• /
• pp.385-396
• /
• 2018

Competing risks are commonly encountered in biomedical research. Regression models for competing risks data can be developed based on data routinely collected in hospitals or general practices. However, these data sets usually contain the covariate missing values. To overcome this problem, multiple imputation is often used to fit regression models under a MAR assumption. Here, we introduce a multivariate imputation in a chained equations algorithm to deal with competing risks survival data. Using pseudo-observations, we make use of the available outcome information by accommodating the competing risk structure. Lastly, we illustrate the practical advantages of our approach using simulations and two data examples from a coronary artery disease data and hepatocellular carcinoma data.

• The Korean Journal of Applied Statistics
• /
• v.28 no.6
• /
• pp.1209-1216
• /
• 2015

Competing-risks events are often observed in a clustered clinical study such as a multi-center clinical trial. We propose a joint modelling approach via a shared frailty term for competing risks survival data from a cluster. For the inference we use the hierarchical likelihood (or h-likelihood), which avoids an intractable integration. We derive the corresponding h-likelihood procedure. The proposed method is illustrated via the analysis of a practical data set.

• Journal of the Korean Data and Information Science Society
• /
• v.27 no.3
• /
• pp.711-724
• /
• 2016

Under variable complex operating environment, various factors can affect the lifetimes of systems. In this research, we study bivariate reliability models having multiple dynamic competing risks. As competing risks, in addition to the natural failure, we consider the increased stress caused by the failure of one component, external shocks, and the level of stress of the working environment at the same time. Considering two reliability models which take into account all of these competing risks, we derive bivariate life distributions. Furthermore, we compare these two models and also compare the distributions of maximum and minimum statistics in the two models.

• Asian Pacific Journal of Cancer Prevention
• /
• v.17 no.3
• /
• pp.1193-1196
• /
• 2016

Background: Colorectal cancer (CRC) is the commonest malignancy in the lower gastrointestinal tract in both men and women. It is the third leading cause of cancer-dependent death in the world. In Iran the incidence of colorectal cancer has increased during the last 25 years. Materials and Methods: In this article we analyzed the survival of 447 colorectal patients of Taleghani hospital in Tehran using parametric competing-risks models. The cancers of these patients were diagnosed during 1985 - 2012 and followed up to 2013. The purpose was to assess the association between survival of patients with colorectal cancer in the presence of competing-risks and prognostic factors using parametric models. The analysis was carried out using R software version 3.0.2. Results: The prognostic variables included in the model were age at diagnosis, tumour site, body mass index and sex. The effect of age at diagnosis and body mass index on survival time was statistically significant. The median survival for Iranian patients with colorectal cancer is about 20 years. Conclusions: Survival function based on Weibull model compared with Kaplan-Meier survival function is smooth. Iranian data suggest a younger age distribution compared to Western reports for CRC.

• The Korean Journal of Applied Statistics
• /
• v.33 no.4
• /
• pp.379-393
• /
• 2020

A terminal event such as death may censor an intermediate event such as relapse, but not vice versa in semi-competing risks data, which is often seen in medicine, public health, and epidemiology. We propose a Weibull regression model with a normal frailty to analyze semi-competing risks data when all three transition times of the illness-death model are possibly interval-censored. We construct the conditional likelihood separately depending on the types of subjects: still alive with or without the intermediate event, dead with or without the intermediate event, and dead with the intermediate event missing. Optimal parameter estimates are obtained from the iterative quasi-Newton algorithm after the marginalization of the full likelihood using the adaptive importance sampling. We illustrate the proposed method with extensive simulation studies and PAQUID (Personnes Agées Quid) data.

• The Korean Journal of Applied Statistics
• /
• v.29 no.7
• /
• pp.1311-1327
• /
• 2016

We propose a multi-state model for analyzing semi-competing risks data with interval-censored or missing intermediate events. This model is an extension of the 'illness-death model', which composes three states, such as 'healthy', 'diseased', and 'dead'. The state of 'diseased' can be considered as an intermediate event. Two more states are added into the illness-death model to describe missing events caused by a loss of follow-up before the end of the study. One of them is a state of 'LTF', representing a lost-to-follow-up, and the other is an unobservable state that represents the intermediate event experienced after LTF occurred. Given covariates, we employ the Cox proportional hazards model with a normal frailty and construct a full likelihood to estimate transition intensities between states in the multi-state model. Marginalization of the full likelihood is completed using the adaptive Gaussian quadrature, and the optimal solution of the regression parameters is achieved through the iterative Newton-Raphson algorithm. Simulation studies are carried out to investigate the finite-sample performance of the proposed estimation procedure in terms of the empirical coverage probability of the true regression parameter. Our proposed method is also illustrated with the dataset adapted from Helmer et al. (2001).

• Journal of Korean Society of Industrial and Systems Engineering
• /
• v.40 no.3
• /
• pp.92-98
• /
• 2017

The rapid growth of engineering technology and the emergence of systemized and large-scale engineering systems have resulted in complexity and uncertainty throughout the lifecycle activities of engineering systems. This complex and large-scale engineering system consists of numerous components, but system failure can be caused by failure of any one of a number of components. There is a real difficulty in managing such a complex and large-scale system as a part. In order to efficiently manage the system and have high reliability, it is necessary to structure a system with a complex structure as a sub-system. Also, in the case of a system in which cause of failures exist at the same time, it is required to identify the correlation of the components lifetime and utilize it for the design policy or maintenance activities of the system. Competitive risk theory has been used as a theory based on this concept. In this study, we apply the competitive risk theory to the models with combined structure of series and parallel which is the basic structure of most complex engineering systems. We construct a competing risks model and propose a mathematical model of net lifetime and crude lifetime for each cause of failure, assuming that the components consisting a parallel system are mutually dependent. In addition, based on the constructed model, the correlation of cause of failure is mathematically analyzed and the hazard function is derived by dividing into net lifetime and crude lifetime.

• The Korean Journal of Applied Statistics
• /
• v.31 no.4
• /
• pp.529-538
• /
• 2018

Tuberculosis causes high morbidity and mortality. However, Korea still has the highest tuberculosis (TB) incidence and mortality among OECD countries despite decreasing incidence and mortality due to the development of modern medicine. Korea has now implemented various policy projects to prevent and control tuberculosis. This study analyzes the effects of public-private mix (PPM) tuberculosis control program on treatment outcomes and identifies the factors that affecting the success of TB treatment. We analyzed 130,000 new tuberculosis patient cohort from 2012 to 2015 using data of tuberculosis patient reports managed by the Disease Control Headquarters. A cumulative incidence function (CIF) compared the cumulative treatment success rates for each factor. We compared the results of the analysis using two popular types of competition risk models (cause-specific Cox's proportional hazards model and subdistribution hazard model) that account for the main event of interest (treatment success) and competing events (death).

• The Korean Journal of Applied Statistics
• /
• v.30 no.4
• /
• pp.539-553
• /
• 2017

We propose a multi-state model to analyze semi-competing risks data with interval-censored or missing intermediate events. This model is an extension of the three states of the illness-death model: healthy, disease, and dead. The 'diseased' state can be considered as the intermediate event. Two more states are added into the illness-death model to incorporate the missing events, which are caused by a loss of follow-up before the end of a study. One of them is a state of the lost-to-follow-up (LTF), and the other is an unobservable state that represents an intermediate event experienced after the occurrence of LTF. Given covariates, we employ the Lin and Ying additive hazards model with log-normal frailty and construct a conditional likelihood to estimate transition intensities between states in the multi-state model. A marginalization of the full likelihood is completed using adaptive importance sampling, and the optimal solution of the regression parameters is achieved through an iterative quasi-Newton algorithm. Simulation studies are performed to investigate the finite-sample performance of the proposed estimation method in terms of empirical coverage probability of true regression parameters. Our proposed method is also illustrated with a dataset adapted from Helmer et al. (2001).