A 3‐level Bayesian mixed effects location scale model with an application to ecological momentary assessment data

Abstract

Ecological momentary assessment studies usually produce intensively measured longitudinal data with large numbers of observations per unit, and research interest is often centered around understanding the changes in variation of people's thoughts, emotions and behaviors. Hedeker et al developed a 2-level mixed effects location scale model that allows observed covariates as well as unobserved variables to influence both the mean and the within-subjects variance, for a 2-level data structure where observations are nested within subjects. In some ecological momentary assessment studies, subjects are measured at multiple waves, and within each wave, subjects are measured over time. Li and Hedeker extended the original 2-level model to a 3-level data structure where observations are nested within days and days are then nested within subjects, by including a random location and scale intercept at the intermediate wave level. However, the 3-level random intercept model assumes constant response change rate for both the mean and variance. To account for changes in variance across waves, as well as clustering attributable to waves, we propose a more comprehensive location scale model that allows subject heterogeneity at baseline as well as across different waves, for a 3-level data structure where observations are nested within waves and waves are then further nested within subjects. The model parameters are estimated using Markov chain Monte Carlo methods. We provide details on the Bayesian estimation approach and demonstrate how the Stan statistical software can be used to sample from the desired distributions and achieve consistent estimates. The proposed model is validated via a series of simulation studies. Data from an adolescent smoking study are analyzed to demonstrate this approach. The analyses clearly favor the proposed model and show significant subject heterogeneity at baseline as well as change over time, for both mood mean and variance. The proposed 3-level location scale model can be widely applied to areas of research where the interest lies in the consistency in addition to the mean level of the responses.

From Nuisance to Novel Research Questions: Using Multilevel Models to Predict Heterogeneous Variances

“…MELS models are a relatively recent development in the long history of different approaches to modeling heterogeneous variances (e.g., Aitkin, 1987;Bryk & Raudenbush, 1988;Culpepper, 2010;Goldstein, 2011;Harvey, 1976;Leonard, 1975;Lindley, 1971;Pinheiro & Bates, 2000;Raudenbush, 1988). Researchers have discussed the development of these models in statistical journals and texts (e.g., Goldstein, 2011;Hedeker, Mermelstein, & Demirtas, 2008;Lin, Mermelstein, & Hedeker, 2018a, 2018bWalters, Hoffman, & Templin, 2018) and applied these models (most often) in fields where the collection of intensive longitudinal data is more common (e.g., medicine; Pugach, Hedeker, Richmond, Sokolovsky, & Mermelstein, 2014). For example, Watts, Walters, Hoffman, and Templin (2016) examined whether time-invariant (Level 2 predictors, e.g., gender, age, Alzheimer's disease status) as well as time-varying predictors (Level 1 predictors, e.g., day monitor worn) were associated with individual differences in mean level (location side) as well as intraindividual variability (scale side) of physical activity.…”

Current State and Future Directions of Technology-Based Ecological Momentary Assessment and Intervention for Major Depressive Disorder: A Systematic Review

“…Because the data are collected during everyday life and in naturalistic environments, it becomes hard or even impossible to have complete control over the setting, and, therefore, it is not possible to rule out the role of confounding variables. Nevertheless, due to the implementation of novel statistical procedures, a balance between research necessities and clinical utility could be achieved [90]. If this were the case in the near future, EMAs and EMIs would undoubtedly transform the field of mental health, greatly contributing to the bridging of science and practice [91,92].…”

Unifying individual differences in personality, predictability and plasticity: A practical guide

“…In the social and medical sciences, DHGLMs are also known as location-scale regression models (with 'location' indicating the mean, and 'scale' indicating the variance; e.g. Lin et al, 2018;Rast et al, 2012).…”