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.

Bayesian Inference for Generalized Linear Mixed Models with Predictors Subject to Detection Limits: an Approach That Leverages Information from Auxiliary Variables.

Bayesian variable selection in joint modeling of longitudinal data and interval-censored failure time data

Causal Inference using Multivariate Generalized Linear Mixed-Effects Models with Longitudinal Data

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

A Multi-parameter regression model for interval censored survival data

A Bayesian joint longitudinal-survival model with a latent stochastic process for intensive longitudinal data

A generalization of moderated statistics to data adaptive semiparametric estimation in high-dimensional biology

Detecting Disease-Associated Genomic Outcomes Using Constrained Mixture of Bayesian Hierarchical Models for Paired Data.

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

A Bayesian joint model of multiple longitudinal and categorical outcomes with application to multiple myeloma using permutation-based variable importance

Penalized Bayesian forward continuation ratio model with application to high-dimensional data with discrete survival outcomes

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

Joint model for interval-censored semi-competing events and longitudinal data with subject-specific within and between visits variabilities

Bayesian Estimation of Two-Part Joint Models for a Longitudinal Semicontinuous Biomarker and a Terminal Event with R-INLA: Interests for Cancer Clinical Trial Evaluation

A Bayesian hierarchical model for prediction of latent health states from multiple data sources with application to active surveillance of prostate cancer

Addressing Multiple Detection Limits with Semiparametric Cumulative Probability Models

Bayesian monotonic errors-in-variables models with applications to pathogen susceptibility testing

Bayesian analysis of longitudinal studies with treatment by indication

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

Semiparametric estimation of the nonmixture cure model with auxiliary survival information

Multiple-model-based robust estimation of causal treatment effect on a binary outcome with integrated information from secondary outcomes