Hierarchical Models for Tumor Xenograft Experiments in Drug Development

Abstract In cancer drug development, demonstrated anticancer activity in animal models is an important step to bring a promising compound to clinic. Proper design and analysis of experiments using laboratory animals have received increasing attention recently. These experiments involve informatively censored longitudinal data with small samples. The problem is further complicated because of order constraints due to the intrinsic growth of control tumors without treatment. This article proposes a Bayesian hierarchical model to analyze informatively censored longitudinal data while accounting for the parameter constraints and providing valid small sample inference. We adopt a noniterative sampling approach, the inverse Bayes formulae (IBF) sampler, to generate independent posterior samples, which avoids convergence problems associated with Markov chain Monte-Carlo methods. To effectively deal with the restricted parameter problem, we use a linear transformation to simplify the constraints and exploit the IBF method to generate random samples from truncated multivariate normal distributions. Because diffuse priors are used, the posterior modes approximate the maximum likelihood estimates well, and the hierarchical model can be considered as an extended mixed-effects model. A real xenograft experiment on a new treatment is analyzed by using the proposed method.

[1]  Xiao-Li Meng,et al.  Maximum likelihood estimation via the ECM algorithm: A general framework , 1993 .

[2]  D A Follmann,et al.  Use of Summary Measures to Adjust for Informative Missingness in Repeated Measures Data with Random Effects , 1999, Biometrics.

[3]  Alan E. Gelfand,et al.  Bayesian statistics without tears: A sampling-resampling perspective , 1992 .

[4]  J. Gaddum Probit Analysis , 1948, Nature.

[5]  Xiao-Li Meng,et al.  Fast EM‐type implementations for mixed effects models , 1998 .

[6]  M D Schluchter,et al.  Methods for the analysis of informatively censored longitudinal data. , 1992, Statistics in medicine.

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

[8]  C. Robert Simulation of truncated normal variables , 2009, 0907.4010.

[9]  Roderick J. A. Little,et al.  Modeling the Drop-Out Mechanism in Repeated-Measures Studies , 1995 .

[10]  M. Kenward,et al.  Informative dropout in longitudinal data analysis (with discussion) , 1994 .

[11]  Hong-Bin Fang,et al.  Small‐Sample Inference for Incomplete Longitudinal Data with Truncation and Censoring in Tumor Xenograft Models , 2002, Biometrics.

[12]  D. Rubin Using the SIR algorithm to simulate posterior distributions , 1988 .

[13]  M. Kenward,et al.  Informative Drop‐Out in Longitudinal Data Analysis , 1994 .

[14]  J. Hughes,et al.  Mixed Effects Models with Censored Data with Application to HIV RNA Levels , 1999, Biometrics.

[15]  M. Tan,et al.  Antitumor activity of temozolomide combined with irinotecan is partly independent of O6-methylguanine-DNA methyltransferase and mismatch repair phenotypes in xenograft models. , 2000, Clinical cancer research : an official journal of the American Association for Cancer Research.

[16]  Hong-Bin Fang,et al.  Experimental design and sample size determination for testing synergism in drug combination studies based on uniform measures , 2003, Statistics in medicine.

[17]  N M Laird,et al.  Mixture models for the joint distribution of repeated measures and event times. , 1997, Statistics in medicine.

[18]  G. Tian,et al.  A NONITERATIVE SAMPLING METHOD FOR COMPUTING POSTERIORS IN THE STRUCTURE OF EM-TYPE ALGORITHMS , 2003 .

[19]  J. A. Breslaw Random Sampling from a Truncated Multivariate Normal Distribution , 1994 .

[20]  J. Kalbfleisch Statistical Inference Under Order Restrictions , 1975 .

[21]  Ming Tan,et al.  BAYESIAN HIERARCHICAL MODELS FOR MULTI-LEVEL REPEATED ORDINAL DATA USING WinBUGS , 2002, Journal of biopharmaceutical statistics.

[22]  Douglas G Altman,et al.  Guidelines for the design and statistical analysis of experiments using laboratory animals. , 2002, ILAR journal.

[23]  D. Oakes Direct calculation of the information matrix via the EM , 1999 .