Go Multivariate: Recommendations on Bayesian Multilevel Hidden Markov Models with Categorical Data.

The multilevel hidden Markov model (MHMM) is a promising method to investigate intense longitudinal data obtained within the social and behavioral sciences. The MHMM quantifies information on the latent dynamics of behavior over time. In addition, heterogeneity between individuals is accommodated with the inclusion of individual-specific random effects, facilitating the study of individual differences in dynamics. However, the performance of the MHMM has not been sufficiently explored. We performed an extensive simulation to assess the effect of the number of dependent variables (1-8), number of individuals (5-90), and number of observations per individual (100-1600) on the estimation performance of a Bayesian MHMM with categorical data including various levels of state distinctiveness and separation. We found that using multivariate data generally alleviates the sample size needed and improves the stability of the results. Moreover, including variables only consisting of random noise was generally not detrimental to model performance. Regarding the estimation of group-level parameters, the number of individuals and observations largely compensate for each other. However, only the former drives the estimation of between-individual variability. We conclude with guidelines on the sample size necessary based on the level of state distinctiveness and separation and study objectives of the researcher.

Bayesian multilevel hidden Markov models identify stable state dynamics in longitudinal recordings from macaque primary motor cortex.

Neural populations, rather than single neurons, may be the fundamental unit of cortical computation. Analysing chronically recorded neural population activity is challenging not only because of the high dimensionality of activity but also because of changes in the signal that may or may not be due to neural plasticity. Hidden Markov models (HMMs) are a promising technique for analysing such data in terms of discrete latent states, but previous approaches have not considered the statistical properties of neural spiking data, have not been adaptable to longitudinal data, or have not modelled condition-specific differences. We present a multilevel Bayesian HMM addresses these shortcomings by incorporating multivariate Poisson log-normal emission probability distributions, multilevel parameter estimation and trial-specific condition covariates. We applied this framework to multi-unit neural spiking data recorded using chronically implanted multi-electrode arrays from macaque primary motor cortex during a cued reaching, grasping and placing task. We show that, in line with previous work, the model identifies latent neural population states which are tightly linked to behavioural events, despite the model being trained without any information about event timing. The association between these states and corresponding behaviour is consistent across multiple days of recording. Notably, this consistency is not observed in the case of a single-level HMM, which fails to generalise across distinct recording sessions. The utility and stability of this approach is demonstrated using a previously learned task, but this multilevel Bayesian HMM framework would be especially suited for future studies of long-term plasticity in neural populations.

Evidence for mood instability in patients with bipolar disorder: Applying multilevel hidden Markov modeling to intensive longitudinal ecological momentary assessment data.

Bipolar disorder (BD) is a chronic psychiatric condition characterized by large episodic changes in mood and energy. Recently, BD has been proposed to be conceptualized as chronic cyclical mood instability, as opposed to the traditional view of alternating discrete episodes with stable periods in-between. Recognizing this mood instability may improve care and call for high-frequency measures coupled with advanced statistical models. To uncover empirically derived mood states, a multilevel hidden Markov model (HMM) was applied to 4-month ecological momentary assessment data in 20 patients with BD, yielding ∼9,820 assessments in total. Ecological momentary assessment data comprised self-report questionnaires (5 × daily) measuring manic and depressive constructs. Manic and depressive symptoms were also assessed weekly using the Altman Self-Rating Mania Scale and the Quick Inventory for Depressive Symptomatology Self-Report. Alignment between HMM-uncovered momentary mood states and weekly questionnaires was assessed with a multilevel linear model. HMM uncovered four mood states: neutral, elevated, mixed, and lowered, which aligned with weekly symptom scores. On average, patients remained < 25 hr in one state. In almost half of the patients, mood instability was observed. Switching between mood states, three patterns were identified: patients switching predominantly between (a) neutral and lowered states, (b) neutral and elevated states, and (c) mixed, elevated, and lowered states. In all, elevated and lowered mood states were interspersed by mixed states. The results indicate that chronic mood instability is a key feature of BD, even in "relatively" euthymic periods. This should be considered in theoretical and clinical conceptualizations of the disorder. (PsycInfo Database Record (c) 2024 APA, all rights reserved).