A Two-Part Mixture Model for Longitudinal Adverse Event Severity Data

We fit a mixed effects logistic regression model to longitudinal adverse event (AE) severity data (four-point ordered categorical response) to describe the dose-AE severity response for an investigational drug. The distribution of the predicted interindividual random effects (Bayes predictions) was extremely bimodal. This extreme bimodality indicated that biased parameter estimates and poor predictive performance were likely. The distribution's primary mode was composed of patients that did not experience an AE. Moreover, the Bayes predictions of these non-AE patients were nearly degenerative; i.e., the predictions were nearly identical. To resolve this extreme bimodality we propose using a two-part mixture modeling approach. The first part models the incidence of AE's, and the second part models the severity grade given the patient had an AE. Unconditional probability predictions are calculated by mixing the incidence and severity model probability predictions. We also report results of simulation studies, which assess the predictive and statistical (bias and precision) performance of our approach.

[1]  J. Hinde,et al.  A Score Test for Testing a Zero‐Inflated Poisson Regression Model Against Zero‐Inflated Negative Binomial Alternatives , 2001, Biometrics.

[2]  V. Farewell,et al.  The use of mixture models for the analysis of survival data with long-term survivors. , 1982, Biometrics.

[3]  L B Sheiner,et al.  Quantitative characterization of therapeutic index: Application of mixed‐effects modeling to evaluate oxybutynin dose–efficacy and dose–side effect relationships , 1999, Clinical pharmacology and therapeutics.

[4]  Diane Lambert,et al.  Zero-inflacted Poisson regression, with an application to defects in manufacturing , 1992 .

[5]  D R Stanski,et al.  Population pharmacodynamic model for ketorolac analgesia , 1996, Clinical pharmacology and therapeutics.

[6]  M. Hutmacher,et al.  Efficient Screening of Covariates in Population Models Using Wald's Approximation to the Likelihood Ratio Test , 2001, Journal of Pharmacokinetics and Pharmacodynamics.

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

[8]  J G Ibrahim,et al.  Maximum Likelihood Methods for Cure Rate Models with Missing Covariates , 2001, Biometrics.

[9]  G. Schwarz Estimating the Dimension of a Model , 1978 .

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

[11]  Robert V. Hogg,et al.  Introduction to Mathematical Statistics. , 1966 .

[12]  Adrian Dunne,et al.  Analysis of Nonrandomly Censored Ordered Categorical Longitudinal Data from Analgesic Trials , 1997 .

[13]  L B Sheiner,et al.  A new approach to the analysis of analgesic drug trials, illustrated with bromfenac data , 1994, Clinical pharmacology and therapeutics.

[14]  David C. Heilbron,et al.  Zero-Altered and other Regression Models for Count Data with Added Zeros , 1994 .

[15]  Gordon Graham,et al.  Determination of an Optimal Dosage Regimen Using a Bayesian Decision Analysis of Efficacy and Adverse Effect Data , 2002, Journal of Pharmacokinetics and Pharmacodynamics.