Bayesian Semi- and Non-parametric Models for Longitudinal Data with Multiple Membership Effects in R.

We introduce growcurves for R that performs analysis of repeated measures multiple membership (MM) data. This data structure arises in studies under which an intervention is delivered to each subject through the subject's participation in a set of multiple elements that characterize the intervention. In our motivating study design under which subjects receive a group cognitive behavioral therapy (CBT) treatment, an element is a group CBT session and each subject attends multiple sessions that, together, comprise the treatment. The sets of elements, or group CBT sessions, attended by subjects will partly overlap with some of those from other subjects to induce a dependence in their responses. The growcurves package offers two alternative sets of hierarchical models: 1. Separate terms are specified for multivariate subject and MM element random effects, where the subject effects are modeled under a Dirichlet process prior to produce a semi-parametric construction; 2. A single term is employed to model joint subject-by-MM effects. A fully non-parametric dependent Dirichlet process formulation allows exploration of differences in subject responses across different MM elements. This model allows for borrowing information among subjects who express similar longitudinal trajectories for flexible estimation. growcurves deploys "estimation" functions to perform posterior sampling under a suite of prior options. An accompanying set of "plot" functions allow the user to readily extract by-subject growth curves. The design approach intends to anticipate inferential goals with tools that fully extract information from repeated measures data. Computational efficiency is achieved by performing the sampling for estimation functions using compiled C++.

[1]  Susan M Paddock,et al.  Bayesian Non-Parametric Hierarchical Modeling for Multiple Membership Data in Grouped Attendance Interventions. , 2013, The annals of applied statistics.

[2]  Peter Congdon,et al.  Bayesian Models for Categorical Data , 2005 .

[3]  S. Paddock,et al.  Bayesian hierarchical semiparametric modelling of longitudinal post‐treatment outcomes from open enrolment therapy groups , 2013, Journal of the Royal Statistical Society. Series A,.

[4]  J. R. Lockwood,et al.  A Model for Teacher Effects From Longitudinal Data Without Assuming Vertical Scaling , 2010 .

[5]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[6]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[7]  J. Ibrahim,et al.  A Bayesian semiparametric joint hierarchical model for longitudinal and survival data. , 2003, Biometrics.

[8]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[9]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[10]  Conrad Sanderson,et al.  RcppArmadillo: Accelerating R with high-performance C++ linear algebra , 2014, Comput. Stat. Data Anal..

[11]  Harvey Goldstein,et al.  Multiple membership multiple classification (MMMC) models , 2001 .

[12]  P. W. Lane,et al.  Analysis of covariance and standardization as instances of prediction. , 1982, Biometrics.

[13]  B. Carlin,et al.  Frailty modeling for spatially correlated survival data, with application to infant mortality in Minnesota. , 2003, Biostatistics.

[14]  B. Carlin,et al.  On the Precision of the Conditionally Autoregressive Prior in Spatial Models , 2003, Biometrics.

[15]  Hadley Wickham,et al.  Reshaping Data with the reshape Package , 2007 .

[16]  Hadley Wickham,et al.  A Layered Grammar of Graphics , 2010 .

[17]  A. Zeileis,et al.  Extended Model Formulas in R : Multiple Parts and Multiple Responses , 2010 .

[18]  T. Ferguson A Bayesian Analysis of Some Nonparametric Problems , 1973 .

[19]  M. Escobar,et al.  Bayesian Density Estimation and Inference Using Mixtures , 1995 .

[20]  W. Fals-Stewart,et al.  Analytic complexities associated with group therapy in substance abuse treatment research: problems, recommendations, and future directions. , 2006, Experimental and clinical psychopharmacology.

[21]  K. Watkins,et al.  An effectiveness trial of group cognitive behavioral therapy for patients with persistent depressive symptoms in substance abuse treatment. , 2011, Archives of general psychiatry.

[22]  D. Blackwell,et al.  Ferguson Distributions Via Polya Urn Schemes , 1973 .

[23]  D. McCaffrey,et al.  ANALYSIS OF ROLLING GROUP THERAPY DATA USING CONDITIONALLY AUTOREGRESSIVE PRIORS. , 2011, The annals of applied statistics.

[24]  Dirk Eddelbuettel,et al.  Rcpp: Seamless R and C++ Integration , 2011 .

[25]  Harvey Goldstein,et al.  Multilevel Modeling of Educational Data With Cross-Classification and Missing Identification for Units , 1998 .

[26]  Peter Müller,et al.  DPpackage: Bayesian Semi- and Nonparametric Modeling in R , 2011 .

[27]  P. Müller,et al.  Semiparametric Bayesian Inference for Multilevel Repeated Measurement Data , 2007, Biometrics.