Interval estimation of random effects in proportional hazards models with frailties

Semi-parametric frailty models are widely used to analyze clustered survival data. In this article, we propose the use of the hierarchical likelihood interval for individual frailties. We study the relationship between hierarchical likelihood, empirical Bayesian, and fully Bayesian intervals for frailties. We show that our proposed interval can be interpreted as a frequentist confidence interval and Bayesian credible interval under a uniform prior. We also propose an adjustment of the proposed interval to avoid null intervals. Simulation studies show that the proposed interval preserves the nominal confidence level. The procedure is illustrated using data from a multicenter lung cancer clinical trial.

[1]  Robert Tibshirani,et al.  An Introduction to the Bootstrap , 1994 .

[2]  J. Booth,et al.  Standard Errors of Prediction in Generalized Linear Mixed Models , 1998 .

[3]  D. Ettinger,et al.  A randomized comparison of standard chemotherapy versus alternating chemotherapy and maintenance versus no maintenance therapy for extensive-stage small-cell lung cancer: a phase III study of the Eastern Cooperative Oncology Group. , 1990, Journal of clinical oncology : official journal of the American Society of Clinical Oncology.

[4]  Huilin Li,et al.  An adjusted maximum likelihood method for solving small area estimation problems , 2010, J. Multivar. Anal..

[5]  Youngjo Lee,et al.  Model selection for multi‐component frailty models , 2007, Statistics in medicine.

[6]  F. Vaida,et al.  Proportional hazards model with random effects. , 2000, Statistics in medicine.

[7]  R. Gray.  A Bayesian analysis of institutional effects in a multicenter cancer clinical trial. , 1994, Biometrics.

[8]  Bias Reduction of Likelihood Estimators in Semiparametric Frailty Models , 2010 .

[9]  M. Plummer.  Penalized loss functions for Bayesian model comparison. , 2008, Biostatistics.

[10]  C. R. Henderson,et al.  Best linear unbiased estimation and prediction under a selection model. , 1975, Biometrics.

[11]  Philip Hougaard,et al.  Analysis of Multivariate Survival Data , 2001 .

[12]  J. Nelder,et al.  Double hierarchical generalized linear models (with discussion) , 2006 .

[13]  T. Louis,et al.  Empirical Bayes Confidence Intervals Based on Bootstrap Samples , 1987 .

[14]  J. Pignon,et al.  Investigating trial and treatment heterogeneity in an individual patient data meta‐analysis of survival data by means of the penalized maximum likelihood approach , 2008, Statistics in medicine.

[15]  C. Morris.  Parametric Empirical Bayes Inference: Theory and Applications , 1983 .

[16]  R. Kass,et al.  Approximate Bayesian Inference in Conditionally Independent Hierarchical Models (Parametric Empirical Bayes Models) , 1989 .

[17]  Bradley P. Carlin,et al.  Approaches for Empirical Bayes Confidence Intervals , 1990 .

[18]  John A. Nelder,et al.  Likelihood Inference for Models with Unobservables: Another View , 2009 .

[19]  Youngjo Lee,et al.  Prediction interval for disease mapping using hierarchical likelihood , 2011, Comput. Stat..

[20]  Youngjo Lee,et al.  Orthodox BLUP versus h-likelihood methods for inferences about  random effects in Tweedie mixed models , 2010, Stat. Comput..

[21]  J. Nelder,et al.  Generalized Linear Models with Random Effects: Unified Analysis via H-likelihood , 2006 .

[22]  Jiming Jiang,et al.  Mixed model prediction and small area estimation , 2006 .

[23]  J. Palmgren,et al.  Estimation of Multivariate Frailty Models Using Penalized Partial Likelihood , 2000, Biometrics.

[24]  A. Rukhin.  Bayes and Empirical Bayes Methods for Data Analysis , 1997 .

[25]  P. Janssen,et al.  A Bayesian approach to jointly estimate centre and treatment by centre heterogeneity in a proportional hazards model , 2005, Statistics in medicine.

[26]  C Morris,et al.  Fitting Weibull duration models with random effects , 1995, Lifetime data analysis.

[27]  Catherine Legrand,et al.  Baseline and treatment effect heterogeneity for survival times between centers using a random effects accelerated failure time model with flexible error distribution , 2007, Statistics in medicine.

[28]  Y. Benjamini,et al.  False Discovery Rate–Adjusted Multiple Confidence Intervals for Selected Parameters , 2005 .

[29]  Youngjo Lee,et al.  Hierarchical likelihood approach for frailty models , 2001 .

[30]  R. Sylvester,et al.  Frailty modelling for survival data from multi‐centre clinical trials , 2011, Statistics in medicine.

[31]  Matthew Knuiman,et al.  A hybrid ML-EM algorithm for calculation of maximum likelihood estimates in semiparametric shared frailty models , 2002 .

[32]  J. F. Bjørnstad.  Predictive Likelihood: A Review , 1990 .

[33]  Youngjo Lee,et al.  Estimating Frailty Models via Poisson Hierarchical Generalized Linear Models , 2003 .

[34]  R. Fay,et al.  Estimates of Income for Small Places: An Application of James-Stein Procedures to Census Data , 1979 .

[35]  Youngjo Lee,et al.  Comparison of hierarchical likelihood versus orthodox best linear unbiased predictor approaches for frailty models , 2005 .