A Hierarchical Bayesian Mixture Model Approach for Analysis of Resting-State Functional Brain Connectivity: An Alternative to Thresholding

This article proposes a Bayesian hierarchical mixture model to analyze functional brain connectivity where mixture components represent “positively connected” and “non-connected” brain regions. Such an approach provides a data-informed separation of reliable and spurious connections in contrast to arbitrary thresholding of a connectivity matrix. The hierarchical structure of the model allows simultaneous inferences for the entire population as well as for each individual subject. A new connectivity measure, the posterior probability of a given pair of brain regions of a specific subject to be connected given the observed correlation of regions' activity, can be computed from the model fit. The posterior probability reflects the connectivity of a pair of regions relative to the overall connectivity pattern of an individual, which is overlooked in traditional correlation analyses. This article demonstrates that using the posterior probability might diminish the effect of spurious connections on inferences, which is present when a correlation is used as a connectivity measure. In addition, simulation analyses reveal that the sparsification of the connectivity matrix using the posterior probabilities might outperform the absolute thresholding based on correlations. Therefore, we suggest that posterior probability might be a beneficial measure of connectivity compared with the correlation. The applicability of the introduced method is exemplified by a study of functional resting-state brain connectivity in older adults.