New statistical software for the proportional hazards model with current status data

Currently, there are some statistical software packages such as intcox, survBayes, and BITE that are designed to analyze current status data (in which failure times are known to be either left- or right-censored) based on the proportional hazards model. They, however, either fail to directly provide standard errors for the estimated regression parameters or require frailty terms. As a result, practitioners often analyze their data using packages for right-censored data. By so doing, they mistreat left-censored observations as exact. This paper describes our newly developed statistical software for the proportional hazards model with current status data. The software is implemented in the R and C languages and consists of the following two simple steps: (a) find MLEs of the regression parameter and the cumulative hazard function; (b) compute the variance-covariance matrix of the regression parameter estimator by using the generalized missing information principle (GMIP) developed by Kim [Kim, J.S., 2003b. Efficient estimation for the proportional hazards model with left-truncated and Case 1 interval-censored data. Statista Sinica 13 (2), 519-537]. Our simulation study results show that our method works well in terms of bias, standard error, and power. By treating current status data as right-censored data, we also show the discrepancy in terms of bias, standard error, and power. Real examples are provided to illustrate the use of the software. This method can be extended to both general interval-censored data and truncated and interval-censored data.

[1]  T. Louis.  Finding the Observed Information Matrix When Using the EM Algorithm , 1982 .

[2]  B. Turnbull.  The Empirical Distribution Function with Arbitrarily Grouped, Censored, and Truncated Data , 1976 .

[3]  Yi-Ching Hsieh,et al.  Clonal evolution and lack of cytogenetic response are adverse prognostic factors for hematologic relapse of chronic phase CML patients treated with imatinib mesylate. , 2004, Blood.

[4]  D. Finkelstein,et al.  A proportional hazards model for interval-censored failure time data. , 1986, Biometrics.

[5]  D. Rubin.  INFERENCE AND MISSING DATA , 1975 .

[6]  David R. Cox,et al.  Regression models and life tables (with discussion , 1972 .

[7]  Jong S. Kim,et al.  Maximum likelihood estimation for the proportional hazards model with partly interval‐censored data , 2003 .

[8]  Karl E. Peace,et al.  Size and Power Assessments of Tests of Hypotheses on Survival Parameters , 1978 .

[9]  D. Cox.  Regression Models and Life-Tables , 1972 .

[10]  Susan O'Brien,et al.  Survival benefit with imatinib mesylate therapy in patients with accelerated‐phase chronic myelogenous leukemia—Comparison with historic experience , 2005, Cancer.

[11]  T. Cai,et al.  Hazard Regression for Interval‐Censored Data with Penalized Spline , 2003, Biometrics.

[12]  Wei Pan.  Extending the Iterative Convex Minorant Algorithm to the Cox Model for Interval-Censored Data , 1999 .

[13]  James M. Robins,et al.  Asymptotic Distribution of P Values in Composite Null Models , 2000 .

[14]  Tommi Härkänen.  BITE: A Bayesian Intensity Estimator , 2003, Comput. Stat..

[15]  Jong S. Kim,et al.  EFFICIENT ESTIMATION FOR THE PROPORTIONAL HAZARDS MODEL WITH LEFT-TRUNCATED AND "CASE 1" INTERVAL-CENSORED DATA , 2003 .

[16]  Wei Pan,et al.  Computation of the NPMLE of distribution functions for interval censored and truncated data with applications to the Cox model , 1998 .

[17]  Jian Huang,et al.  Efficient estimation for the proportional hazards model with interval censoring , 1996 .

[18]  D. Rubin.  Multiple Imputation After 18+ Years , 1996 .

[19]  J. Wellner,et al.  Information Bounds and Nonparametric Maximum Likelihood Estimation , 1992 .