Efficient computation of subset selection probabilities with application to Cox regression

An efficient method is presented to compute the probabilityof selection of a specified subset from the set of all subsetsof a fixed size where the subsets are taken from a populationwhose units have varying individual probabilities of selection.The problem is motivated by the computation of the exact marginallikelihood for the Cox proportional hazards model.     

Polynomial spline estimation and inference of proportional hazards regression models with flexible relative risk form

The Cox proportional hazards model usually assumes an exponential form for the dependence of the hazard function on covariate variables. However, in practice this assumption may be violated and other relative risk forms may be more appropriate. In this article, we consider the proportional hazards model with an unknown relative risk form. Issues in model interpretation are addressed. We propose a method to estimate the relative risk form and the regression parameters simultaneously by first approximating the logarithm of the relative risk form by a spline, and then employing the maximum partial likelihood estimation. An iterative alternating optimization procedure is developed for efficient implementation. Statistical inference of the regression coefficients and of the relative risk form based on parametric asymptotic theory is discussed. The proposed methods are illustrated using simulation and an application to the Veteran's Administration lung cancer data.     

Hazard regression for interval-censored data with penalized spline

This article introduces a new approach for estimating the hazard function for possibly interval- and right-censored survival data. We weakly parameterize the log-hazard function with a piecewise-linear spline and provide a smoothed estimate of the hazard function by maximizing the penalized likelihood through a mixed model-based approach. We also provide a method to estimate the amount of smoothing from the data. We illustrate our approach with two well-known interval-censored data sets. Extensive numerical studies are conducted to evaluate the efficacy of the new procedure.     

Additive hazard regression with auxiliary covariates

We consider the additive hazard model when some of the truecovariates are measured only on a randomly selected validationset whereas auxiliary covariates are observed for all studysubjects. An updated pseudoscore estimation approach is proposedfor the parameters of the additive hazard model. It allows oneto fit the model with auxiliary covariates, while leaving thebaseline hazard unspecified. Asymptotic properties of the proposedestimators are established, and consistent standard error estimatorsare developed. Simulations demonstrate that the asymptotic approximationsof the proposed estimates are adequate for practical use. Areal example is used to illustrate the performance of the proposedmethod.     

Estimating regression parameters and degree of dependence for multivariate failure time data

Multivariate failure time data are frequently encountered in longitudinal studies when subjects may experience several events or when there is a grouping of individuals into a cluster. To take into account the dependence of the failure times within the unit (the individual or the cluster) as well as censoring, two multivariate generalizations of the Cox proportional hazards model are commonly used. The marginal hazard model is used when the purpose is to estimate mean regression parameters, while the frailty model is retained when the purpose is to assess the degree of dependence within the unit. We propose a new approach based on the combination of the two aforementioned models to estimate both these quantities. This two-step estimation procedure is quicker and more simple to implement than the EM algorithm used in frailty models estimation. Simulation results are provided to illustrate robustness, consistency, and large-sample properties of estimators. Finally, this method is exemplified on a diabetic retinopathy study in order to assess the effect of photocoagulation in delaying the onset of blindness as well as the dependence between the two eyes blindness times of a patient.     

"Smooth" semiparametric regression analysis for arbitrarily censored time-to-event data

Summary .   A general framework for regression analysis of time-to-event data subject to arbitrary patterns of censoring is proposed. The approach is relevant when the analyst is willing to assume that distributions governing model components that are ordinarily left unspecified in popular semiparametric regression models, such as the baseline hazard function in the proportional hazards model, have densities satisfying mild "smoothness" conditions. Densities are approximated by a truncated series expansion that, for fixed degree of truncation, results in a "parametric" representation, which makes likelihood-based inference coupled with adaptive choice of the degree of truncation, and hence flexibility of the model, computationally and conceptually straightforward with data subject to any pattern of censoring. The formulation allows popular models, such as the proportional hazards, proportional odds, and accelerated failure time models, to be placed in a common framework; provides a principled basis for choosing among them; and renders useful extensions of the models straightforward. The utility and performance of the methods are demonstrated via simulations and by application to data from time-to-event studies.     

Competing risks regression for stratified data

For competing risks data, the Fine-Gray proportional hazards model for subdistribution has gained popularity for its convenience in directly assessing the effect of covariates on the cumulative incidence function. However, in many important applications, proportional hazards may not be satisfied, including multicenter clinical trials, where the baseline subdistribution hazards may not be common due to varying patient populations. In this article, we consider a stratified competing risks regression, to allow the baseline hazard to vary across levels of the stratification covariate. According to the relative size of the number of strata and strata sizes, two stratification regimes are considered. Using partial likelihood and weighting techniques, we obtain consistent estimators of regression parameters. The corresponding asymptotic properties and resulting inferences are provided for the two regimes separately. Data from a breast cancer clinical trial and from a bone marrow transplantation registry illustrate the potential utility of the stratified Fine-Gray model.     

Shared frailty models for recurrent events and a terminal event

There has been an increasing interest in the analysis of recurrent event data (Cook and Lawless, 2002, Statistical Methods in Medical Research 11, 141-166). In many situations, a terminating event such as death can happen during the follow-up period to preclude further occurrence of the recurrent events. Furthermore, the death time may be dependent on the recurrent event history. In this article we consider frailty proportional hazards models for the recurrent and terminal event processes. The dependence is modeled by conditioning on a shared frailty that is included in both hazard functions. Covariate effects can be taken into account in the model as well. Maximum likelihood estimation and inference are carried out through a Monte Carlo EM algorithm with Metropolis-Hastings sampler in the E-step. An analysis of hospitalization and death data for waitlisted dialysis patients is presented to illustrate the proposed methods. Methods to check the validity of the proposed model are also demonstrated. This model avoids the difficulties encountered in alternative approaches which attempt to specify a dependent joint distribution with marginal proportional hazards and yields an estimate of the degree of dependence.     

Regression analysis of multivariate interval-censored failure time data with application to tumorigenicity experiments

This paper discusses multivariate interval-censored failure time data that occur when there exist several correlated survival times of interest and only interval-censored data are available for each survival time. Such data occur in many fields. One is tumorigenicity experiments, which usually concern different types of tumors, tumors occurring in different locations of animals, or together. For regression analysis of such data, we develop a marginal inference approach using the additive hazards model and apply it to a set of bivariate interval-censored data arising from a tumorigenicity experiment. Simulation studies are conducted for the evaluation of the presented approach and suggest that the approach performs well for practical situations.