Estimating time-varying effects for overdispersed recurrent events data with treatment switching.

In the analysis of multivariate event times, frailty models assuming time-independent regression coefficients are often considered, mainly due to their mathematical convenience. In practice, regression coefficients are often time dependent and the temporal effects are of clinical interest. Motivated by a phase III clinical trial in multiple sclerosis, we develop a semiparametric frailty modelling approach to estimate time-varying effects for overdispersed recurrent events data with treatment switching. The proposed model incorporates the treatment switching time in the time-varying coefficients. Theoretical properties of the proposed model are established and an efficient expectation-maximization algorithm is derived to obtain the maximum likelihood estimates. Simulation studies evaluate the numerical performance of the proposed model under various temporal treatment effect curves. The ideas in this paper can also be used for time-varying coefficient frailty models without treatment switching as well as for alternative models when the proportional hazard assumption is violated. A multiple sclerosis dataset is analysed to illustrate our methodology.

[1]  Carl de Boor,et al.  A Practical Guide to Splines , 1978, Applied Mathematical Sciences.

[2]  L. Schumaker.  Spline Functions: Basic Theory , 1981 .

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

[4]  D. Clayton,et al.  Multivariate generalizations of the proportional hazards model , 1985 .

[5]  L. J. Wei,et al.  Regression analysis of multivariate incomplete failure time data by modeling marginal distributions , 1989 .

[6]  A. Karr,et al.  Nonparametric Survival Analysis with Time-Dependent Covariate Effects: A Penalized Partial Likelihood Approach , 1990 .

[7]  J. Robins,et al.  Correcting for non-compliance in randomized trials using rank preserving structural failure time models , 1991 .

[8]  Robert Gray,et al.  Flexible Methods for Analyzing Survival Data Using Splines, with Applications to Breast Cancer Prognosis , 1992 .

[9]  Richard D. Gill,et al.  A counting process approach to maximum likelihood estimation in frailty models , 1992 .

[10]  J P Klein,et al.  Semiparametric estimation of random effects using the Cox model based on the EM algorithm. , 1992, Biometrics.

[11]  Susan A. Murphy,et al.  Testing for a Time Dependent Coefficient in Cox's Regression Model , 1993 .

[12]  Jianwen Cai,et al.  Some Graphical Displays and Marginal Regression Analyses for Recurrent Failure Times and Time Dependent Covariates , 1993 .

[13]  Susan A. Murphy,et al.  Consistency in a Proportional Hazards Model Incorporating a Random Effect , 1994 .

[14]  Susan A. Murphy,et al.  Asymptotic Theory for the Frailty Model , 1995 .

[15]  Jianwen Cai,et al.  Estimating equations for hazard ratio parameters based on correlated failure time data , 1995 .

[16]  T. Louis,et al.  Survival analysis using a scale change random effects model , 1995 .

[17]  R. Prentice,et al.  Regression Estimation Using Multivariate Failure Time Data and a Common Baseline Hazard Function Model , 1997, Lifetime data analysis.

[18]  Erik T. Parner,et al.  Asymptotic theory for the correlated gamma-frailty model , 1998 .

[19]  D. Oakes,et al.  Frailty Models and Rank Tests , 1998, Lifetime data analysis.

[20]  Danyu Lin,et al.  Marginal Regression Models for Multivariate Failure Time Data , 1998 .

[21]  Zhiliang Ying,et al.  Semiparametric regression for the mean and rate functions of recurrent events , 2000 .

[22]  Jianqing Fan,et al.  Efficient Estimation and Inferences for Varying-Coefficient Models , 2000 .

[23]  Ib M. Skovgaard,et al.  Efficient Estimation of Fixed and Time‐varying Covariate Effects in Multiplicative Intensity Models , 2002 .

[24]  J. Whitehead,et al.  Estimating a treatment effect in survival studies in which patients switch treatment , 2002, Statistics in medicine.

[25]  Jianqing Fan,et al.  Variable Selection for Cox's proportional Hazards Model and Frailty Model , 2002 .

[26]  Zongwu Cai,et al.  Local Linear Estimation for Time-Dependent Coefficients in Cox's Regression Models , 2003 .

[27]  B. Svennerholm,et al.  Acyclovir treatment of relapsing-remitting multiple sclerosis , 1996, Journal of Neurology.

[28]  D. Gilden.  Infectious causes of multiple sclerosis , 2005, The Lancet Neurology.

[29]  L. J. Wei,et al.  On the Cox Model With Time-Varying Regression Coefficients , 2005 .

[30]  S. Chow,et al.  Statistical inference for cancer trials with treatment switching , 2005, Statistics in medicine.

[31]  D. Zeng,et al.  Semiparametric Transformation Models With Random Effects for Recurrent Events , 2007 .

[32]  Jianwen Cai,et al.  Hazard models with varying coefficients for multivariate failure time data , 2007, 0708.0519.

[33]  Jianwen Cai,et al.  Partially linear hazard regression with varying coefficients for multivariate survival data , 2008 .

[34]  L. Meyerson,et al.  Statistical methods for the analysis of relapse data in MS clinical trials , 2009, Journal of the Neurological Sciences.

[35]  Norbert Benda,et al.  The impact of dropouts on the analysis of dose‐finding studies with recurrent event data , 2010, Statistics in medicine.

[36]  Heinz Schmidli,et al.  Blinded sample size reestimation with count data: Methods and applications in multiple sclerosis , 2010, Statistics in medicine.

[37]  Ludwig Kappos,et al.  A placebo-controlled trial of oral fingolimod in relapsing multiple sclerosis. , 2010, The New England journal of medicine.

[38]  T Friede,et al.  Blinded Sample Size Reestimation with Negative Binomial Counts in Superiority and Non-inferiority Trials , 2010, Methods of Information in Medicine.

[39]  Torben Martinussen,et al.  The Aalen additive gamma frailty hazards model , 2011 .

[40]  Jianwen Cai,et al.  Additive Mixed Effect Model for Clustered Failure Time Data , 2011, Biometrics.

[41]  Joseph G Ibrahim,et al.  Estimating treatment effects with treatment switching via semicompeting risks models: an application to a colorectal cancer study. , 2012, Biometrika.

[42]  Richard Nicholas,et al.  Time-patterns of annualized relapse rates in randomized placebo-controlled clinical trials in relapsing multiple sclerosis: A systematic review and meta-analysis , 2012, Multiple sclerosis.

[43]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .