• Journal of the Korean Data and Information Science Society
• /
• v.26 no.3
• /
• pp.569-580
• /
• 2015

Log-rank is widely used for testing equality of two survival functions, and this method is efficient only under the proportional hazard assumption. However, crossing survival functions are common in practice. Therefore, many approaches have been suggested to test equality of them. This study considered several methods; Renyi type test, modified Kolmogorov-Smirnov and Cramer-von Mises test, and weighted Log-rank test, which can be applied when the survival functions cross, and simulated power of those methods. Based on the simulation results, we provide the useful information to choose a suitable approach in a given situation.

• The Korean Journal of Applied Statistics
• /
• v.22 no.6
• /
• pp.1345-1353
• /
• 2009

We introduce a simple methodology, so-called generalized gamma-polynomial approximation, based on moment-matching technique to approximate survival and hazard functions in the context of parametric survival analysis. We use the generalized gamma-polynomial approximation to approximate the density and distribution functions of convolutions and finite mixtures of random variables, from which the approximated survival and hazard functions are obtained. This technique provides very accurate approximation to the target functions, in addition to their being computationally efficient and easy to implement. In addition, the generalized gamma-polynomial approximations are very stable in middle range of the target distributions, whereas saddlepoint approximations are often unstable in a neighborhood of the mean.

• Communications for Statistical Applications and Methods
• /
• v.24 no.2
• /
• pp.115-127
• /
• 2017

This paper proposes a Bayesian test for the equivalence of survival functions in multiple groups. Proposed Bayesian test use the model of Cox's regression with time-varying coefficients. B-spline expansions are used for the time-varying coefficients, and the proposed test use only the partial likelihood, which provides easier computations. Various simulations of the proposed test and typical tests such as log-rank and Fleming and Harrington tests were conducted. This result shows that the proposed test is consistent as data size increase. Specifically, the power of the proposed test is high despite the existence of crossing hazards. The proposed test is based on a Bayesian approach, which is more flexible when used in multiple tests. The proposed test can therefore perform various tests simultaneously. Real data analysis of Larynx Cancer Data was conducted to assess applicability.

• Journal of the Korean Data and Information Science Society
• /
• v.13 no.2
• /
• pp.31-38
• /
• 2002

In the case of multiple survival times which might be censored at each covariate vector, we study the regression quantile estimations in this paper. The estimations are based on the empirical distribution functions of the censored times and the sample quantiles of the observed survival times at each covariate vector and the weighted least square method is applied for the estimation of the regression quantile. The estimators are shown to be asymptotically normally distributed under some regularity conditions.

• Journal of Korean Society of Forest Science
• /
• v.88 no.3
• /
• pp.320-326
• /
• 1999

In this study, 15 survival functions in integral and difference forms for forest plantation were derived based on assumptions for the number of surviving trees and the differential forms of the mortality rate model. Then, performance of the models was evaluated by fitting to remeasurement data of unthinned white pine(Pinus strobes) forest plantation. As a result, three equations associated with a power function of age, $t^{\beta}$, are somewhat more suitable for describing the effect of self-thinning over time. On the other hand, a general survival function for Japanese larch(Larix leptolepts) forest plantation was derived in order to exam the effect of site quality on self-thinning procedures. The results indicate that the $N_{min}$ is negatively correlated with site index and, even though the same initial stand density was assumed, the survival function curves differ in shapes associated with site index values.

• Genomics & Informatics
• /
• v.17 no.4
• /
• pp.41.1-41.12
• /
• 2019

Survival analysis mainly deals with the time to event, including death, onset of disease, and bankruptcy. The common characteristic of survival analysis is that it contains "censored" data, in which the time to event cannot be completely observed, but instead represents the lower bound of the time to event. Only the occurrence of either time to event or censoring time is observed. Many traditional statistical methods have been effectively used for analyzing survival data with censored observations. However, with the development of high-throughput technologies for producing "omics" data, more advanced statistical methods, such as regularization, should be required to construct the predictive survival model with high-dimensional genomic data. Furthermore, machine learning approaches have been adapted for survival analysis, to fit nonlinear and complex interaction effects between predictors, and achieve more accurate prediction of individual survival probability. Presently, since most clinicians and medical researchers can easily assess statistical programs for analyzing survival data, a review article is helpful for understanding statistical methods used in survival analysis. We review traditional survival methods and regularization methods, with various penalty functions, for the analysis of high-dimensional genomics, and describe machine learning techniques that have been adapted to survival analysis.

• Journal of the Korean Statistical Society
• /
• v.23 no.1
• /
• pp.179-185
• /
• 1994

To test the equality of survival distributions in the presence of arbitrary right censorship, the choice of weights which are functions of the number of individuals at risk at the time of each death is very important in increasing the power of the test. In this paper a weight by probit scale is derived and the efficiencies relative to the other weight's are also investigated.

• Journal of Applied Reliability
• /
• v.16 no.1
• /
• pp.56-63
• /
• 2016

Purpose: The purpose of this study is to point out that the Kaplan-Meier method is not valid to calculate the survival probability or failure probability (risk) in the presence of competing risks and to introduce more valid method of cumulative incidence function. Methods: Survival analysis methods have been widely used in biostatistics division. However the same methods have not been utilized in reliability division. Especially competing risks cases, where several causes of failure occur and the occurrence of one event precludes the occurrence of the other events, are scattered in reliability field. But they are not noticed in the realm of reliability expertism or they are analysed in the wrong way. Specifically Kaplan-Meier method which assumes that the censoring times and failure times are independent is used to calculate the probability of failure in the presence of competing risks, thereby overestimating the real probability of failure. Hence, cumulative incidence function is introduced and sample competing risks data are analysed using cumulative incidence function and some graphs. Finally comparison of cumulative incidence functions and regression type analysis are mentioned briefly. Results: Cumulative incidence function is used to calculate the survival probability or failure probability (risk) in the presence of competing risks and some useful graphs depicting the failure trend over the lifetime are introduced. Conclusion: This paper shows that Kaplan-Meier method is not appropriate for the evaluation of survival or failure over the course of lifetime. In stead, cumulative incidence function is shown to be useful. Some graphs using the cumulative incidence functions are also shown to be informative.

• Communications for Statistical Applications and Methods
• /
• v.28 no.5
• /
• pp.425-445
• /
• 2021

A generalization of the log-logistic (LL) distribution called exponentiated log-logistic (ELL) distribution on lines of exponentiated Weibull distribution is considered. In this paper, based on progressive type-II censored samples, we have derived the maximum likelihood estimators and Bayes estimators for three parameters, the survival function and hazard function of the ELL distribution. Then, under the balanced squared error loss (BSEL) and the balanced linex loss (BLEL) functions, their corresponding Bayes estimators are obtained using Lindley's approximation (see Jung and Chung, 2018; Lindley, 1980), Tierney-Kadane approximation (see Tierney and Kadane, 1986) and Markov Chain Monte Carlo methods (see Hastings, 1970; Gelfand and Smith, 1990). Here, to check the convergence of MCMC chains, the Gelman and Rubin diagnostic (see Gelman and Rubin, 1992; Brooks and Gelman, 1997) was used. On the basis of their risks, the performances of their Bayes estimators are compared with maximum likelihood estimators in the simulation studies. In this paper, research supports the conclusion that ELL distribution is an efficient distribution to modeling data in the analysis of survival data. On top of that, Bayes estimators under various loss functions are useful for many estimation problems.

• Communications for Statistical Applications and Methods
• /
• v.20 no.3
• /
• pp.169-174
• /
• 2013

The ROC curve is drawn with two conditional cumulative distribution functions (or survival functions) of the univariate random variable. In this work, we consider joint cumulative distribution functions of k random variables, and suggest a ROC curve for multivariate random variables. With regard to the values on the line, which passes through two mean vectors of dichotomous states, a joint cumulative distribution function can be regarded as a function of the univariate variable. After this function is modified to satisfy the properties of the cumulative distribution function, a ROC curve might be derived; moreover, some illustrative examples are demonstrated.