Flexible Modeling in the Koziol-Green Model by a Copula Function

In survival analysis, the classical Koziol-Green random censorship model is commonly used to describe informative censoring. Hereby, it is assumed that the distribution of the censoring time is a power of the distribution of the survival time. In this article, we extend this model by assuming a general function between these distributions. We determine this function from a relationship between the observable random variables which is described by a copula family that depends on an unknown parameter θ. For this setting, we develop a semi-parametric estimator for the distribution of the survival time in which we propose a pseudo-likelihood estimator for the copula parameter θ. As results, we show first the consistency and asymptotic normality of the estimator for θ. Afterwards, we prove the weak convergence of the process associated to the semi-parametric distribution estimator. Furthermore, we investigate the finite sample performance of these estimators through a simulation study and finally apply it to a practical data set on survival with malignant melanoma.     

Non-parametric Estimation of a Survival Function with Two-stage Design Studies

Abstract.  The two-stage design is popular in epidemiology studies and clinical trials due to its cost effectiveness. Typically, the first stage sample contains cheaper and possibly biased information, while the second stage validation sample consists of a subset of subjects with accurate and complete information. In this paper, we study estimation of a survival function with right-censored survival data from a two-stage design. A non-parametric estimator is derived by combining data from both stages. We also study its large sample properties and derive pointwise and simultaneous confidence intervals for the survival function. The proposed estimator effectively reduces the variance and finite-sample bias of the Kaplan–Meier estimator solely based on the second stage validation sample. Finally, we apply our method to a real data set from a medical device postmarketing surveillance study.     

A generalization of the compound rayleigh distribution: using a bayesian method on cancer survival times

In this paper the generalized compound Rayleigh model, exhibiting flexible hazard rate, is high¬lighted. This makes it attractive for modelling survival times of patients showing characteristics of a random hazard rate. The Bayes estimators are derived for the parameters of this model and some survival time parameters from a right censored sample. This is done with respect to conjugate and discrete priors on the parameters of this model, under the squared error loss function, Varian's asymmetric linear-exponential (linex) loss function and a weighted linex loss function. The future survival time of a patient is estimated under these loss functions. A Monte Carlo simu¬lation procedure is used where closed form expressions of the estimators cannot be obtained. An example illustrates the proposed estimators for this model.     

Joint Modeling of Survival and Longitudinal Ordered Data Using a Semiparametric Approach

Medical research frequently focuses on the relationship between quality of life (QoL) and survival time of subjects. QoL may be one of the most important factors that could be used to predict survival, making it worth identifying factors that jointly affect survival and QoL. We propose a semiparametric joint model that consists of item response and survival components, where these two components are linked through latent variables. Several popular ordinal models are considered and compared in the item response component, while the Cox proportional hazards model is used in the survival component. We estimate the baseline hazard function and model parameters simultaneously, through a profile likelihood approach. We illustrate the method using an example from a clinical study.     

Analysis of spatial frailty models by a weighted estimating equation

Spatially correlated survival data are frequently observed in ecological and epidemiological studies. An assumption in the clustered survival models is inter-cluster independence, which may not be adequate to model the dependence in spatial settings. For survival data, the likelihood function based on a spatial frailty may be complicated. In this paper, we develop a weighted estimating equation for spatially right-censored data. Some large sample properties for the estimate are developed. We also conduct simulations to compare estimation performance with other methods. A data set from a study of forest decline in Wisconsin is used to illustrate the proposed method.     

Estimation of the Average Survival Function Using a Censored Data Regression Model

In the presence of covariates information, assuming the linear relationship between a transformation of survival time and covariates, we propose a new estimator of survival function and show its consistency. In addition, a comparison of the proposed estimator with the product-limit estimator introduced by Kaplan and Meier (1958) is performed through Monte Carlo simulation studies. We illustrate the proposed estimator with the updated Stanford heart transplant data.     

An improved survival estimator for censored medical costs with a kernel approach

Cost assessment serves as an essential part in economic evaluation of medical interventions. In many studies, costs as well as survival data are frequently censored. Standard survival analysis techniques are often invalid for censored costs, due to the induced dependent censoring problem. Owing to high skewness in many cost data, it is desirable to estimate the median costs, which will be available with estimated survival function of costs. We propose a kernel-based survival estimator for costs, which is monotone, consistent, and more efficient than several existing estimators. We conduct numerical studies to examine the finite-sample performance of the proposed estimator.     

Maximum Penalized Likelihood Estimation in a Gamma-Frailty Model

The shared frailty models allow for unobserved heterogeneity or for statistical dependence between observed survival data. The most commonly used estimation procedure in frailty models is the EM algorithm, but this approach yields a discrete estimator of the distribution and consequently does not allow direct estimation of the hazard function. We show how maximum penalized likelihood estimation can be applied to nonparametric estimation of a continuous hazard function in a shared gamma-frailty model with right-censored and left-truncated data. We examine the problem of obtaining variance estimators for regression coefficients, the frailty parameter and baseline hazard functions. Some simulations for the proposed estimation procedure are presented. A prospective cohort (Paquid) with grouped survival data serves to illustrate the method which was used to analyze the relationship between environmental factors and the risk of dementia.     

Estimation of parameters under a random censorship model

The problem of estimation of parameters of a lifetime distribution is considered under the proportional hazards model of random censorship. Asymptotic variances of several estimators of survival function are compared in the eponential case.     

Nonparametric bayes estimation of the survival function from a record of failures and follow-ups

In some observational studies, we have random censoring model. However, the data available may be partially observable censored data consisting of the observed failure times and only those nonfailure times which are subject to follow-up. Suzuki (1985) discussed the problem of nonparametric estimation of the survival function from such partially observable censored data. In this article, we derive a nonparametric Bayes estimator of the survival function for such data of failures and follow-ups under a Dirichlet process prior and squared error loss. The limiting properties such as the mean square consistency, weak convergence and strong consistency of the Bayes estimator are studied. Finally, the procedures developed are illustrated by means of an example.     

Nonparametric change point estimation for survival distributions with a partially constant hazard rate

Lifetime Data Analysis - We present a new method for estimating a change point in the hazard function of a survival distribution assuming a constant hazard rate after the change point and a...     

A semi-parametric estimation of a survival function from incomplete and doubly censored data

The purpose of this paper is to present a semi-parametric estimation of a survival function when analyzing incomplete and doubly censored data. Under the assumption that the chance of censoring is not related to the individual's survivorship, we propose a consistent estimation of survival. The derived estimator treats the uncensored observations nonparametrically and uses parametric models for both right and left censored data. Some asymptotic properties and simulation studies are also presented in order to analyze the behavior of the proposed estimator.     

Detecting change in a hazard regression model with right-censoring

The hazard function plays an important role in survival analysis and reliability, since it quantifies the instantaneous failure rate of an individual at a given time point t, given that this individual has not failed before t. In some applications, abrupt changes in the hazard function are observed, and it is of interest to detect the location of such a change. In this paper, we consider testing of existence of a change in the parameters of an exponential regression model, based on a sample of right-censored survival times and the corresponding covariates. Likelihood ratio type tests are proposed and non-asymptotic bounds for the type II error probability are obtained. When the tests lead to acceptance of a change, estimators for the location of the change are proposed. Non-asymptotic upper bounds of the underestimation and overestimation probabilities are obtained. A short simulation study illustrates these results.     

Modeling Lifetimes by a Stochastic Process Hitting a Critical Point

First hitting times arise naturally in survival analysis where the underlying stochastic counting process represents the strength of the health of an individual. The patient experiences a clinical endpoint when this process reaches a critical point for the first time. We propose a very flexible and unified first hitting time density function in a stochastic carcinogenesis counting process, and its mathematical properties are investigated. The Poisson and negative binomial first hitting time models are addressed and two examples with real data are presented.     

Adjusted Hazard Rate Estimator Based on a Known Censoring Probability

In most reliability studies involving censoring, one assumes that censoring probabilities are unknown. We derive a nonparametric estimator for the survival function when information regarding censoring frequency is available. The estimator is constructed by adjusting the Nelson–Aalen estimator to incorporate censoring information. Our results indicate significant improvements can be achieved if available information regarding censoring is used. We compare this model to the Koziol–Green model, which is also based on a form of proportional hazards for the lifetime and censoring distributions. Two examples of survival data help to illustrate the differences in the estimation techniques.     

分数年龄假设的新方法：Kriging模型

在生存函数的计算中,生命表只提供了其在整数年龄上的值。当计算非整数年龄上的生存函数时就需要进行分数年龄假设。经典的分数年龄假设在数学上容易处理,但却容易导致死力函数不连续,更重要的是无法保证其在分数年龄上估计的精确性。分数年龄假设本质上是一种插值技术。本研究尝试将一种插值性能优越的Kriging模型引入到分数年龄假设中,对整数年龄上的生存函数进行插值,并基于良好拟合的生存函数进一步构建死力函数及平均余命函数。基于Kriging模型的分数年龄假设的有效性通过了Makeham法则下的生存函数的验证,其结果表明,Kriging模型的插值性能远胜过经典的分数年龄假设模型。     

Performance of nonparametric maximum likelihood estimator in a

When the probability of selecting an individual in a population is propor­tional to its lifelength, it is called length biased sampling. A nonparametric maximum likelihood estimator (NPMLE) of survival in a length biased sam­ple is given in Vardi (1982). In this study, we examine the performance of Vardi's NPMLE in estimating the true survival curve when observations are from a length biased sample. We also compute estimators based on a linear combination (LCE) of empirical distribution function (EDF) estimators and weighted estimators. In our simulations, we consider observations from a mix­ture of two different distributions, one from F and the other from G which is a length biased distribution of F. Through a series of simulations with vari­ous proportions of length biasing in a sample, we show that the NPMLE and the LCE closely approximate the true survival curve. Throughout the sur­vival curve, the EDF estimators overestimate the survival. We also consider a case where the observations are from three different weighted distributions, Again, both the NPMLE and the LCE closely approximate the true distribu­tion, indicating that the length biasedness is properly adjusted for. Finally, an efficiency study shows that Vardi's estimators are more efficient than the EDF estimators in the lower percentiles of the survival curves.     

A Bayesian nonparametric estimator of a multivariate survival function

Using reinforced processes related to beta-Stacy process and generalized Pólya urn scheme jointly with a structure assumption about dependence, a Bayesian nonparametric prior and a predictive estimator for a multivariate survival function are provided. This estimator can be computed through an easy implementation of a Gibbs sampler algorithm. Moreover consistency of the estimator is studied.     

A New Kernel Distribution Function Estimator Based on a Non-parametric Transformation of the Data

Abstract.  A new kernel distribution function (df) estimator based on a non-parametric transformation of the data is proposed. It is shown that the asymptotic bias and mean squared error of the estimator are considerably smaller than that of the standard kernel df estimator. For the practical implementation of the new estimator a data-based choice of the bandwidth is proposed. Two possible areas of application are the non-parametric smoothed bootstrap and survival analysis. In the latter case new estimators for the survival function and the mean residual life function are derived.     

Minimax estimation of a probability of success under LINEX loss

Minimax estimation of a binomial probability under LINEX loss function is considered. It is shown that no equalizer estimator is available in the statistical decision problem under consideration. It is pointed out that the problem can be solved by determining the Bayes estimator with respect to a least favorable distribution having finite support. In this situation, the optimal estimator and the least favorable distribution can be determined only by using numerical methods. Some properties of the minimax estimators and the corresponding least favorable prior distributions are provided depending on the parameters of the loss function. The properties presented are exploited in computing the minimax estimators and the least favorable distributions. The results obtained can be applied to determine minimax estimators of a cumulative distribution function and minimax estimators of a survival function.