Abstract:The Women's Health Initiative (WHI) Life and Longevity After Cancer (LILAC) study is an excellent resource for studying the quality of life following breast cancer treatment. At study entry, women were asked about new symptoms that appeared following their initial cancer treatment. In this article, we were interested in using regression modeling to estimate associations of clinical and lifestyle factors at cancer diagnosis (independent variables) with the number of new symptoms (dependent variable). Although clinical and lifestyle data were collected longitudinally, few measurements were obtained at diagnosis or at a consistent timepoint prior to diagnosis, which complicates the analysis. Furthermore, parametric count models, such as the Poisson and negative binomial, do not fit the symptom data well. Thus, motivated by the issues encountered in LILAC, we propose two Bayesian joint models for longitudinal data and a count outcome. Our two models differ according to the assumption on the outcome distribution: one uses a negative binomial (NB) distribution and the other a nonparametric rounded mixture of Gaussians (RMG). The mean of each count distribution is dependent on imputed values of continuous, binary, and ordinal variables at a time point of interest (e.g. diagnosis). To facilitate imputation, longitudinal variables are modeled jointly using a linear mixed model for a latent underlying normal random variable, and a Dirichlet process prior is assigned to the random subject-specific effects to relax distribution assumptions. In simulation studies, the RMG joint model exhibited superior power and predictive accuracy over the NB model when the data were not NB. The RMG joint model also outperformed an RMG model containing predictors imputed using the last value carried forward, which generated estimates that were biased toward the null. We used our models to examine the relationship between sleep health at diagnosis and the number of new symptoms following breast cancer treatment in LILAC.

Application of Beta-Distribution and Combined Uniform and Binomial Methods in Longitudinal Modeling of Bounded Outcome Score Data

Modeling of Correlated Cognitive Function and Functional Disability Outcomes with Bounded and Missing Data in a Longitudinal Aging Study

A robust regression model for bounded count health data

Semiparametric multivariate joint model for skewed-longitudinal and survival data: A Bayesian approach

Bayesian joint analysis of heterogeneous- and skewed-longitudinal data and a binary outcome, with application to AIDS clinical studies

A Bayesian Phase II Proof-of-concept Design for Clinical Trials with Longitudinal Endpoints.

Robust Bayesian cumulative probit linear mixed models for longitudinal ordinal data

Bayesian Nonparametric Model for Heterogeneous Treatment Effects With Zero‐Inflated Data

A threshold mixed-effects Tobit model for treatment-sensitive subgroup identification based on longitudinal measures with floor and ceiling effects and a continuous covariate

A Single Index Model for Longitudinal Outcomes to Optimize Individual Treatment Decision Rules

Use of copula to model within-study association in bivariate meta-analysis of binomial data at the aggregate level a Bayesian approach and application to surrogate endpoint evaluation

Flexible Bayesian Modeling for Longitudinal Binary and Ordinal Responses

A Bayesian beta-binomial piecewise growth mixture model for longitudinal overdispersed binomial data

Bayesian semiparametric joint modeling of a count outcome and inconveniently timed longitudinal predictors

A joint model for (un)bounded longitudinal markers, competing risks, and recurrent events using patient registry data

Random Effects Models for Analyzing Mixed Overdispersed Binomial and Normal Longitudinal Responses With Application to Kidney Function Data of Cancer Patients

A Bayesian Joint model for Longitudinal DAS28 Scores and Competing Risk Informative Drop Out in a Rheumatoid Arthritis Clinical Trial

A penalized Cox proportional hazards model with multiple time-varying exposures

A Bayesian Approach to Modeling Variance of Intensive Longitudinal Biomarker Data as a Predictor of Health Outcomes

A two-part mixed-effects model for analyzing longitudinal microbiome compositional data