Nonparametric Estimation of the Multivariate Distribution Function in a Censored Regression Model with Applications

In a regression model with univariate censored responses, a new estimator of the joint distribution function of the covariates and response is proposed, under the assumption that the response and the censoring variable are independent conditionally to the covariates. This estimator is based on the conditional Kaplan–Meier estimator of Beran (1981 Beran , R. ( 1981 ). Nonparametric regression with randomly censored survival data. Technical Report, University of California, Berkeley, California . [Google Scholar]), and happens to be an extension of the multivariate empirical distribution function used in the uncensored case. We derive asymptotic i.i.d. representations for the integrals with respect to the measure defined by this estimated distribution function. These representations hold even in the case where the covariates are multidimensional under some additional assumption on the censoring. Applications to censored regression and to density estimation are considered.     

An Estimator for the Survival Function on the Basis of Observable Censoring Times

Abstract

Chiu [Chiu, S. N. (1999 Chiu, S. N. 1999. An unbiased estimator for the survival function of censored data. Commun. Statist. - Theory Meth., 28(9): 2249–2260. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). An unbiased estimator for the survival function of censored data. Commun. Statist. - Theory Meth. 28(9):2249–2260.] proposed a nonparametric estimator for the survival function which is based on observable censoring times in the general censoring model. His estimator is less efficient than the Product-Limit estimator. Considering an informative censoring model this drawback can partially be overcome. This is shown by a nonparametric, uniformly consistent estimator based on observable censoring times within the simple Koziol–Green model. Some asymptotic properties of the new estimator are investigated and it is compared with the well-known ACL-estimator.     

Maximum Likelihood Estimation in Compound Geometric Failure Model with Changing Parameters From Type-I Two-Stage Progressively Censored and Group Censored Samples

Boardman and Kendell (1970 Boardman , T. J. , Kendell , P. J. ( 1970 ). Estimation in compound failure models . Technometrics 12 : 891 – 908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered the problem of estimation with respect to Type-I censoring when an item is subjected to only one of the two causes of failure assuming exponential model. Patel and Gajjar (1992 Patel , M. N. , Gajjar , A. V. ( 1992 ). Maximum likelihood estimation in compound exponential failure model with changing failure rates from Type-I progressively censored and group censored samples . Commun. Statist. Theor. Meth. 21 ( 10 ): 2899 – 2908 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) considered extension of the Boardman and Kendell's results in case of two-stage progressive censoring. Here we have considered geometric competing risk failure model with two independent causes of failures. Maximum likelihood estimation of the parameters is carried out using Type-I two-stage progressively censored and group censored samples. Asymptotic standard errors of the estimators are obtained for both the cases. Two illustrative examples are cited for ungroup and group competing risk models.     

On goodness-of-fit tests for Aalen's additive model based on left-truncated right-censored or doubly censored data

ABSTRACT

Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) proposed goodness-of-fit tests for Aalen's additive risk model. In this article, we demonstrate that the approach of Gandy and Jensen (2005 Gandy, A., Jensen, U. (2005). On goodness-of-fit tests for Aalen's additive risk model. Scan. J. Stat. 32:425–445.[Crossref], [Web of Science ®] , [Google Scholar]) can be applied to left-truncated right-censored (LTRC) data and doubly censored data. A simulation study is conducted to investigate the performance of the proposed tests. The proposed tests are illustrated using heart transplant data.     

Estimation for exponential distribution based on competing risk middle censored data

Abstract

For a general censoring scheme called “middle censoring” scheme which was proposed by Jammalamadaka and Mangalam (2003 Jammalamadaka, S.R., Mangalam, V. (2003). Nonparametric estimation for middle censored data. J. Nonparametric Statist. 15:253–265.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in nonparametric set up. In this article, point and interval estimation problems are considered for the exponential distribution when the failure data is middle censored with two independent competing failure risks. Different methods are introduced to estimate the unknown model parameters such as maximum likelihood estimation, midpoint approximation, equivalent quantities estimation. The Bayesian estimation is also considered with gamma priors. Two numerical examples are analyzed to show the performance of the proposed methods.     

A proportional hazards model for the analysis of doubly censored competing risks data

ABSTRACT

Competing risks data are common in medical research in which lifetime of individuals can be classified in terms of causes of failure. In survival or reliability studies, it is common that the patients (objects) are subjected to both left censoring and right censoring, which is refereed as double censoring. The analysis of doubly censored competing risks data in presence of covariates is the objective of this study. We propose a proportional hazards model for the analysis of doubly censored competing risks data, using the hazard rate functions of Gray (1988 Gray, R.J. (1988). A class of k-sample tests for comparing the cumulative incidence of a competing risk. Ann. Statist. 16:1141–1154.[Crossref], [Web of Science ®] , [Google Scholar]), while focusing upon one major cause of failure. We derive estimators for regression parameter vector and cumulative baseline cause specific hazard rate function. Asymptotic properties of the estimators are discussed. A simulation study is conducted to assess the finite sample behavior of the proposed estimators. We illustrate the method using a real life doubly censored competing risks data.     

Estimation of Survivorship Function Based on Doubly Censored Data with Application to HIV/AIDS

Left censoring concept has been defined in different ways in statistical applications. Turnbull (1974 Turnbull , B. W. ( 1974 ). Nonparametric estimation of a survivorship function with doubly censored data . J. Amer. Statist. Assoc. 69 : 169 – 173 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) defines it in a particular way. Whereas in recent literature, especially in epidemiological studies, it has been defined in another way. This difference between the two approaches is the main reason that despite simplicity, Turnbull method cannot be applicable in all cases of doubly censored data. In this article we present a modified Turnbull method for analysis of doubly censored data adequate with recent definition. Comparison has been done with other statistical methods, including imputation estimator, full likelihood-based and conditional likelihood-based approach using Iranian HIV data.     

Chi-Squared Type Test for the AFT-Generalized Inverse Weibull Distribution

The generalized inverse Weibull distribution is a newlife time probability distribution which can be used to model a variety of failure characteristics. It has several desirable properties and nice physical interpretations which enable them to be used frequently. In this article, we present a chi-squared goodness-of-fit test for an accelerated failure time (AFT) model with generalized inverse Weibull distribution (GIW) as the baseline distribution, in both of complete and censored data. This test is based on a modification of the NRR (Nikulin-Rao-Robson) statistic Y², proposed by Bagdonavicius and Nikulin (2011 Bagdonavicius, V., Nikulin, M. (2011). Chi-squared tests for general composite hypotheses from censored samples. Comptes Rendus Mathematique, Ser. I, 349(3–4): 219–223.[Crossref], [Web of Science ®] , [Google Scholar]), for censored data. Two applications of real data are given to illustrate the potentiality of the proposed test.     

Analysis of discrete lifetime data under middle-censoring and in the presence of covariates

‘Middle censoring’ is a very general censoring scheme where the actual value of an observation in the data becomes unobservable if it falls inside a random interval (L, R) and includes both left and right censoring. In this paper, we consider discrete lifetime data that follow a geometric distribution that is subject to middle censoring. Two major innovations in this paper, compared to the earlier work of Davarzani and Parsian [3 N. Davarzani and A. Parsian, Statistical inference for discrete middle-censored data, J. Statist. Plan. Inference 141 (2011), pp. 1455–1462. doi: 10.1016/j.jspi.2010.10.012[Crossref], [Web of Science ®] , [Google Scholar]], include (i) an extension and generalization to the case where covariates are present along with the data and (ii) an alternate approach and proofs which exploit the simple relationship between the geometric and the exponential distributions, so that the theory is more in line with the work of Iyer et al. [6 S.K. Iyer, S.R. Jammalamadaka, and D. Kundu, Analysis of middle censored data with exponential lifetime distributions, J. Statist. Plan. Inference 138 (2008), pp. 3550–3560. doi: 10.1016/j.jspi.2007.03.062[Crossref], [Web of Science ®] , [Google Scholar]]. It is also demonstrated that this kind of discretization of life times gives results that are close to the original data involving exponential life times. Maximum likelihood estimation of the parameters is studied for this middle-censoring scheme with covariates and their large sample distributions discussed. Simulation results indicate how well the proposed estimation methods work and an illustrative example using time-to-pregnancy data from Baird and Wilcox [1 D.D. Baird and A.J. Wilcox, Cigarette smoking associated with delayed conception, J, Am. Med. Assoc. 253 (1985), pp. 2979–2983. doi: 10.1001/jama.1985.03350440057031[Crossref], [Web of Science ®] , [Google Scholar]] is included.     

A modified multiple comparisons with a control for exponential location parameters based on doubly censored sample under heteroscedasticity

ABSTRACT

In this paper, a modified one-stage multiple comparison procedures with a control for exponential location parameters based on the doubly censored sample under heteroscedasticity is proposed. A simulation study is done and the results show that the proposed procedures have shorter confidence length with coverage probabilities closer to the nominal ones compared with the one proposed in Wu (2017 Wu, S. F. 2017. Multiple comparisons of exponential location parameters with a control based on doubly censored sample under heteroscedasticity. Communications in Statistics: Simulation and Computation 46 (3):1858–1870. doi: 10.1080/03610918.2015.1017582.[Taylor & Francis Online] , [Google Scholar]). At last, an example of comparing the duration of remission for four drugs as the treatment of leukemia is given to demonstrate the proposed procedures.     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

Closeness of k-records to Progressive Type-II Censored Order Statistics for Location-scale Families

In this article, the Pitman closeness of upper and lower k-records to progressive Type-II censored order statistics for location-scale families is investigated. In each case, the special properties of the probability of Pitman closeness are obtained and the corresponding monotonicity properties are discussed. Moreover, the closest k-record to a specific progressive Type-II censored data is obtained. Finally, for the standard exponential and standard uniform distributions, explicit expressions for the probability of Pitman closeness are derived. For various censoring schemes, the results of the numerical computations are displayed in tables. Most of the results in Ahmadi and Balakrishnan (2013) Ahmadi, J., Balakrishnan, N. (2013). On the nearness of record values to order statistics from Pitman measure of closeness. Metrika 76:521–541.[Crossref], [Web of Science ®] , [Google Scholar] can be achieved as special cases.     

Semiparametric analysis of transformation models with doubly censored data

Double censoring arises when T represents an outcome variable that can only be accurately measured within a certain range, [L, U], where L and U are the left- and right-censoring variables, respectively. In this note, using Martingale arguments of Chen et al. [3 Chen, K., Jin, Z. and Ying, Z. 2002. Semiparametric analysis of transformation models with censored data. Biometrika, 89: 659–668. [Crossref], [Web of Science ®] [Google Scholar]], we propose an estimator (denoted by ?β) for estimating regression coefficients of transformation model when L is always observed. Under Cox proportional hazards model, the proposed estimator is equivalent to the partial likelihood estimator for left-truncated and right-censored data if the left-censoring variables L were regarded as left-truncated variables. In this case, the estimator ?β can be obtained by the standard software. A simulation study is conducted to investigate the performance of ?β. For the purpose of comparison, the simulation study also includes the estimator proposed by Cai and Cheng [2 Cai, T. and Cheng, S. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] [Google Scholar]] for the case when L and U are always observed.     

Inferences on partial linear models with right censored data by empirical likelihood

ABSTRACT

This article considers inference for partial linear models with right censored data. We use empirical likelihood based on the Buckley and James (1979 Buckley, J., James, I. (1979). Linear regression with censored data. Biometrika 66:429–436.[Crossref], [Web of Science ®] , [Google Scholar]) estimating equation to derive the confidence region for the regression parameter. We introduce an adjusted empirical likelihood ratio statistic for the parameter of interest and show that its limiting distribution is standard chi-square. A simulation is carried out to compare our method with the synthetic data approach in Wang and Li (2002 Wang, Q.-H., Li, G. (2002). Empirical Likelihood Semiparametric Regression Analysis under Random Censorship. J. Multivariate Anal. 83:469–486.[Crossref], [Web of Science ®] , [Google Scholar]).     

Locally Maximin Test for the Mixture Proportion with Type-I Censoring

Mixtures of distributions arise frequently in the context of life testing experiments, where one also encounters the problem of censored data. In this article, we derive the locally most powerful (LMP) test for testing the mixing proportion in a general mixture model based on type-I censored data. We also prove the additional properties of unbiasedness and locally maximin using a novel approach. To this end, we prove an extension of a standard lemma in the testing literature (Lehmann, 1986 Lehmann , E. L. ( 1986 ). Testing Statistical Hypotheses . New York : Wiley .[Crossref] , [Google Scholar]) relating to families with monotone likelihood ratio.     

Comparison of cumulative incidence functions of current status competing risks data with discrete observation times

Abstract

Competing risks data with current status censoring arise frequently from transversal studies in demography, epidemiology and reliability theory; where the only information about lifetime is whether the event of interest has occurred or not before a monitoring time. In practice, the monitoring times are discrete, but most of the studies consider them as continuous in nature. In the present paper, we propose a non parametric test for comparing cumulative incidence functions of current status competing risks data while the observation (monitoring) times are discrete. Asymptotic distribution of the test statistic is also derived. A simulation study is conducted to assess the finite sample behavior of the test statistic. The practical utility of the procedure is well demonstrated using a real-life data set on menopausal history of 2423 women given in Jewell, van der Laan, and Henneman (2003 Jewell, N. P., M. van der Laan, and T. Henneman. 2003. Nonparametric estimation from current status data with competing risks. Biometrika 90 (1):183–197. doi: 10.1093/biomet/90.1.183.[Crossref], [Web of Science ®] , [Google Scholar]).     

Inferences for Mixtures of Distributions for Centrally Censored Data with Partial Identification

In this article, several methods to make inferences about the parameters of a finite mixture of distributions in the context of centrally censored data with partial identification are revised. These methods are an adaptation of the work in Contreras-Cristán, Gutiérrez-Peña, and O'Reilly (2003 Contreras-Cristán , A. , Gutiérrez-Peña , E. , O'Reilly , F. ( 2003 ). Inferences using latent variables for mixtures of distributions for censored data with partial identification . Comm. Stat. Theor. Meth. 32 ( 4 ): 749 – 774 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in the case of right censoring. The first method focuses on an asymptotic approximation to a suitably simplified likelihood using some latent quantities; the second method is based on the expectation-maximization (EM) algorithm. Both methods make explicit use of latent variables and provide computationally efficient procedures compared to non-Bayesian methods that deal directly with the full likelihood of the mixture appealing to its asymptotic approximation. The third method, from a Bayesian perspective, uses data augmentation to work with an uncensored sample. This last method is related to a recently proposed Bayesian method in Baker, Mengersen, and Davis (2005 Baker , P. , Mengersen , K. , Davis , G. ( 2005 ). A Bayesian solution to reconstructing centrally censored distributions . J. Agr. Biol. Environ. Stat. 1 : 61 – 84 . [Google Scholar]). Our proposal of the three adapted methods is shown to provide similar inferential answers, thus offering alternative analyses.     

Simultaneous Confidence Intervals for the Parameters of Pareto Distribution under Progressive Censoring

In this article, we consider the progressive Type II right censored sample from Pareto distribution. We introduce a new approach for constructing the simultaneous confidence interval of the unknown parameters of this distribution under progressive censoring. A Monte Carlo study is also presented for illustration. It is shown that this confidence region has a smaller area than that introduced by Ku? and Kaya (2007 Ku? , C. , Kaya , M. F. ( 2007 ). Estimation for the parameters of the Pareto distribution under progressive censoring . Commun. Statist. Theor. Meth. 36 : 1359 – 1365 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

Survival Weibull regression model for mismeasured outcomes

In some survival studies, the exact time of the event of interest is unknown, but the event is known to have occurred during a particular period of time (interval-censored data). If the diagnostic tool used to detect the event of interest is not perfectly sensitive and specific, outcomes may be mismeasured; a healthy subject may be diagnosed as sick and a sick one may be diagnosed as healthy. In such cases, traditional survival analysis methods produce biased estimates for the time-to-failure distribution parameters (Paggiaro and Torelli 2004 Paggiaro, A., and N. Torelli. 2004. The effect of classification errors in survival data analysis. Statistical Methods and Applications 13:213–25.[Crossref] , [Google Scholar]). In this context, we developed a parametric model that incorporates sensitivity and specificity into a grouped survival data analysis (a case of interval-censored data in which all subjects are tested at the same predetermined time points). Inferential aspects and properties of the methodology, such as the likelihood function and identifiability, are discussed in this article. Assuming known and non differential misclassification, Monte Carlo simulations showed that the proposed model performed well in the case of mismeasured outcomes; the estimates of the relative bias of the model were lower than those provided by the naive method that assumes perfect sensitivity and specificity. The proposed methodology is illustrated by a study related to mango tree lifetimes.