Modelling longitudinal semicontinuous emesis volume data with serial correlation in an acupuncture clinical trial

In longitudinal studies, we are often interested in modelling repeated assessments of volume over time. Our motivating example is an acupuncture clinical trial in which we compare the effects of active acupuncture, sham acupuncture and standard medical care on chemotherapy-induced nausea in patients being treated for advanced stage breast cancer. An important end point for this study was the daily measurement of the volume of emesis over a 14-day follow-up period. The repeated volume data contained many 0s, had apparent serial correlation and had missing observations, making analysis challenging. The paper proposes a two-part latent process model for analysing the emesis volume data which addresses these challenges. We propose a Monte Carlo EM algorithm for parameter estimation and we use this methodology to show the beneficial effects of acupuncture on reducing the volume of emesis in women being treated for breast cancer with chemotherapy. Through simulations, we demonstrate the importance of correctly modelling the serial correlation for making conditional inference. Further, we show that the correct model for the correlation structure is less important for making correct inference on marginal means. Copyright 2005 Royal Statistical Society.

[1]  Gary K Grunwald,et al.  Analysis of repeated measures data with clumping at zero , 2002, Statistical methods in medical research.

[2]  Peter A Lachenbruch,et al.  Analysis of data with excess zeros , 2002, Statistical methods in medical research.

[3]  Joseph L Schafer,et al.  A Two-Part Random-Effects Model for Semicontinuous Longitudinal Data , 2001 .

[4]  Hulin Wu,et al.  Design of viral dynamic studies for efficiently assessing potency of anti-HIV therapies in AIDS Clinical Trials , 2002 .

[5]  J Rochon,et al.  Analyzing bivariate repeated measures for discrete and continuous outcome variables. , 1996, Biometrics.

[6]  P. Shekelle,et al.  Electroacupuncture for control of myeloablative chemotherapy-induced emesis: A randomized controlled trial. , 2000, JAMA.

[7]  M. Tan,et al.  Robustness of the Latent Variable Model for Correlated Binary Data , 1999, Biometrics.

[8]  C. Morris,et al.  A Comparison of Alternative Models for the Demand for Medical Care , 1983 .

[9]  W. Tu,et al.  Comparison of Several Independent Population Means When Their Samples Contain Log‐Normal and Possibly Zero Observations , 1999, Biometrics.

[10]  J K Lindsey,et al.  On the appropriateness of marginal models for repeated measurements in clinical trials. , 1998, Statistics in medicine.

[11]  S. Zeger,et al.  Longitudinal data analysis using generalized linear models , 1986 .

[12]  W. Tu,et al.  A Wald test comparing medical costs based on log-normal distributions with zero valued costs. , 1999, Statistics in medicine.

[13]  P. Heagerty,et al.  Misspecified maximum likelihood estimates and generalised linear mixed models , 2001 .

[14]  Raymond J. Carroll,et al.  Estimation and comparison of changes in the presence of informative right censoring by modeling the censoring process , 1988 .

[15]  C. McCulloch.  Maximum Likelihood Algorithms for Generalized Linear Mixed Models , 1997 .

[16]  J. Booth,et al.  Maximizing generalized linear mixed model likelihoods with an automated Monte Carlo EM algorithm , 1999 .

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

[18]  G. C. Wei,et al.  A Monte Carlo Implementation of the EM Algorithm and the Poor Man's Data Augmentation Algorithms , 1990 .

[19]  C. Geyer,et al.  Constrained Monte Carlo Maximum Likelihood for Dependent Data , 1992 .

[20]  Paul S Albert,et al.  A Latent Autoregressive Model for Longitudinal Binary Data Subject to Informative Missingness , 2002, Biometrics.