Regression analysis of current status data with latent variables.

This paper discusses regression analysis of current status data with dependent censoring, a problem that often occurs in many areas such as cross-sectional studies, epidemiological investigations and tumorigenicity experiments. Copula model-based methods are commonly employed to tackle this issue. However, these methods often face challenges in terms of model and parameter identification. The primary aim of this paper is to propose a copula-based analysis for dependent current status data, where the association parameter is left unspecified. Our method is based on a general class of semiparametric linear transformation models and parametric copulas. We demonstrate that the proposed semiparametric model is identifiable under certain regularity conditions from the distribution of the observed data. For inference, we develop a sieve maximum likelihood estimation method, using Bernstein polynomials to approximate the nonparametric functions involved. The asymptotic consistency and normality of the proposed estimators are established. Finally, to demonstrate the effectiveness and practical applicability of our method, we conduct an extensive simulation study and apply the proposed method to a real data example.

Current status data arise when each study subject is observed only once and the survival time of interest is known only to be either less or greater than the observation time. Such data often occur in, for example, cross-sectional studies, demographical investigations and tumorigenicity experiments and several semi-parametric and non-parametric methods for their analysis have been proposed. However, most of these methods deal only with the situation where observation time is independent of the underlying survival time completely or given covariates. This paper discusses regression analysis of current status data when the observation time may be related to the underlying survival time and inference procedures are presented for estimation of regression parameters under the additive hazards regression model. The procedures can be easily implemented and are applied to two motivating examples.

ABSTRACTMany methods have been developed in the literature for regression analysis of current status data with noninformative censoring and also some approaches have been proposed for semiparametric regression analysis of current status data with informative censoring. However, the existing approaches for the latter situation are mainly on specific models such as the proportional hazards model and the additive hazard model. Corresponding to this, in this paper, we consider a general class of semiparametric linear transformation models and develop a sieve maximum likelihood estimation approach for the inference. In the method, the copula model is employed to describe the informative censoring or relationship between the failure time of interest and the censoring time, and Bernstein polynomials are used to approximate the nonparametric functions involved. The asymptotic consistency and normality of the proposed estimators are established, and an extensive simulation study is conducted and indicates that the proposed approach works well for practical situations. In addition, an illustrative example is provided.

This paper discusses regression analysis of current status data, a type of failure time data where each study subject is observed only once, in the presence of dependent censoring. Furthermore, there may exist a cured subgroup, meaning that a proportion of study subjects are not susceptible to the failure event of interest. For the problem, we develop a sieve maximum likelihood estimation approach with the use of latent variables and Bernstein polynomials. For the determination of the proposed estimators, an EM algorithm is developed and the asymptotic properties of the estimators are established. Extensive simulation studies are conducted and indicate that the proposed method works well for practical situations. A motivating application from a tumorigenicity experiment is also provided.

On efficient estimation in additive hazards regression with current status data

This paper considers joint analysis of current status and marker data using a threshold model based on first hitting times. A failure time is defined as the time at which a subject's latent health status process first decreases to zero. We extend the bivariate Wiener process model in Whitmore et al. (1998) to the case when only current status data are available. We develop maximum likelihood estimation procedures and provide simulation studies. We apply our methods to a motivating example involving liver tumors in mice.

ABSTRACTIn this paper, we introduce a partially linear single-index additive hazards model with current status data. Both the unknown link function of the single-index term and the cumulative baseline hazard function are approximated by B-splines under a monotonicity constraint on the latter. The sieve method is applied to estimate the nonparametric and parametric components simultaneously. We show that, when the nonparametric link function is an exact B-spline, the resultant estimator of regression parameter vector is asymptotically normal and achieves the semiparametric information bound and the rate of convergence of the estimator for the cumulative baseline hazard function is optimal. Simulation studies are presented to examine the finite sample performance of the proposed estimation method. For illustration, we apply the method to a clinical dataset with current status outcome.

SUMMARY A great deal of recent attention has focused on the estimation of survival distributions based on current status data, an extreme form of interval censored data. This particular data structure arises in a wide variety of applications where cross-sectional observation either naturally occurs or is preferred to more traditional forms of follow-up. Here we consider current status data in the context of competing risks. We briefly consider simple parametric models as a backdrop to nonparametric procedures. We make some brief comparisons and remarks regarding the nonparametric maximum likelihood estimator. The ideas are illustrated on the data of Krailo & Pike (1983) which considers estimation of the age distribution at both natural and operative menopause. We also consider the case where there is exact observation of failure times due to one of the competing risks when failure occurs prior to the monitoring time. Basic techniques of survival analysis focus on estimation of the survival distribution based on various forms of censored data. An extreme form of censoring arises where the only information on studied individuals is their survival status at a single monitoring time. This particular data structure is known as current status data. Nonparametric estimation of the survival function and semiparametric techniques for related regression models, based on current status data, have been much studied of late. Here, we consider estimation of sub-survival functions based on current status observations in the presence of competing risks. For simplicity, we describe only two competing risks, but the ideas are easily extended when there are more than two possible causes of failure. Let )j denote the cause-specific hazard function for cause j = 1, 2 (Kalbfleisch & Prentice, 2002, p. 252). If T and J are the random variables that measure time to failure and cause of failure, respectively, the two sub-distribution functions of interest are

Recent attention has focused on the length of the post partum period of infecundity in those societies where child-spacing has traditionally been achieved through post partum sexual abstinence or prolonged breast-feeding. The tabulation of proportions of women abstaining at interview by the age of the youngest child has been suggested as a technique free of problems of digital preference and memory lapse, while being otherwise equivalent to the distribution of the proportion of women abstaining by completed months since their last confinement, calculated from retrospectively-collected durations of abstinence. Both stochastic and deterministic models demonstrate that the techniques are comparable only when they are applied to non-recurrent events such as first marriage. Childbearing, and the relaxation of the proscription on post partum sexual intercourse, is a recurrent event, and the distribution of proportions of currently abstaining women according to the age of their youngest child is heavily biased. This finding invalidates the use of the post partum current status distribution.

Recent attention has focused on the length of the post partum period of infecundity in those societies where child-spacing has traditionally been achieved through post partum sexual abstinence or prolonged breast-feeding. The tabulation of proportions of women abstaining at interview by the age of the youngest child has been suggested as a technique free of problems of digital preference and memory lapse, while being otherwise equivalent to the distribution of the proportion of women abstaining by completed months since their last confinement, calculated from retrospectively-collected durations of abstinence. Both stochastic and deterministic models demonstrate that the techniques are comparable only when they are applied to non-recurrent events such as first marriage. Childbearing, and the relaxation of the proscription on post partum sexual intercourse, is a recurrent event, and the distribution of proportions of currently abstaining women according to the age of their youngest child is heavily biased. This finding invalidates the use of the post partum current status distribution.

We consider a varying-coefficient partially linear proportional odds model with current status data. This model enables one to examine the extent to which some covariates interact nonlinearly with an exposure variable, while other covariates present linear effects. B-spline approach and sieve maximum likelihood estimation method are used to get an integrated estimate for the linear coefficients, the varying-coefficient functions and the baseline function. The proposed parameter estimators are proved to be semiparametrically efficient and asymptotically normal, and the estimators for the nonparametric functions achieve the optimal rate of convergence. Simulation studies and a real data analysis are used for assessment and illustration.

The current status and panel count data frequently arise from cancer and tumorigenicity studies when events currently occur. A common and widely used class of two sample tests, for current status and panel count data, is the permutation class. We manipulate the double saddlepoint method to calculate the exact mid-p-values of the underlying permutation distributions of this class of tests. Permutation simulations are replaced by analytical saddlepoint computations which provide extremely accurate mid-p-values that are exact for most practical purposes and almost always more accurate than normal approximations. The method is illustrated using two real tumorigenicity panel count data. To compare the saddlepoint approximation with the normal asymptotic approximation, a simulation study is conducted. The speed and accuracy of the saddlepoint method facilitate the calculation of the confidence interval for the treatment effect. The inversion of the mid-p-values to calculate the confidence interval for the mean rate of development of the recurrent event is discussed.

This study develops a new joint modeling approach to simultaneously analyze longitudinal and time-to-event data with latent variables. The proposed model consists of three components. The first component is a hidden Markov model for investigating a longitudinal observation process and its underlying transition process as well as their potential risk factors and dynamic heterogeneity. The second component is a factor analysis model for characterizing latent risk factors through multiple observed variables. The third component is a proportional hazards model for examining the effects of observed and latent risk factors on the hazards of interest. A shared random effect is introduced to allow the longitudinal and time-to-event outcomes to be correlated. A Bayesian approach coupled with efficient Markov chain Monte Carlo methods is developed to conduct statistical inference. The performance of the proposed method is evaluated through simulation studies. An application of the proposed model to a general health survey study concerning cognitive impairment and mortality for Chinese elders is presented.

We investigate empirical likelihood for the additive hazards model with current status data. An empirical log-likelihood ratio for a vector or subvector of regression parameters is defined and its limiting distribution is shown to be a standard chi-squared distribution. The proposed inference procedure enables us to make empirical likelihood-based inference for the regression parameters. Finite sample performance of the proposed method is assessed in simulation studies to compare with that of a normal approximation method, it shows that the empirical likelihood method provides more accurate inference than the normal approximation method. A real data example is used for illustration.

In survival analysis, estimating the failure time distribution is an important and difficult task, since usually the data is subject to censoring. Specifically, in this paper we consider current status data, a type of data where the failure time cannot be directly observed. The format of the data is such that the failure time is restricted to knowledge of whether or not the failure time exceeds a random monitoring time. We propose a flexible kernel machine approach for estimation of the failure time expectation as a function of the covariates, with current status data. In order to obtain the kernel machine decision function, we minimize a regularized version of the empirical risk with respect to a new loss function. Using finite sample bounds and novel oracle inequalities, we prove that the obtained estimator converges to the true conditional expectation for a large family of probability measures. Finally, we present a simulation study and an analysis of real-world data that compares the performance of the proposed approach to existing methods. We show empirically that our approach is comparable to current state of the art, and in some cases is even better.