ACAM-J-MARS: Multi-Atlas Riemannian Spectral Representations for Decoding Advanced Meditation States from fMRI
Organizations: Meditation Research Program, Psychiatry, MGH and Harvard Medical School, Boston, MA, USA · Martinos Center for Biomedical Imaging, Radiology, MGH and Harvard Medical School, Boston, MA, USA · College of Engineering, Northeastern University, Boston, MA, USA
Abstract
Inferring jhana-type advanced concentration absorption meditation (ACAM-J) states from fMRI data is challenging due to the spatially distributed nature of meditation-related effects, substantial inter-individual variability that can exceed within-subject state differences, and limited sample sizes. Conventional machine-learning approaches often formulate this problem as unimodal, Euclidean, atlas-based classification, thereby neglecting the non-Euclidean geometry of functional connectivity matrices and underrepresenting temporal information in functional magnetic resonance imaging (fMRI); these limitations can lead to suboptimal classification performance. We propose ACAM-J-MARS (Multi-Atlas Riemannian Spectral representation), which computes functional connectivity using multiple atlases, represents symmetric positive definite (SPD) connectivity matrices on a Riemannian manifold, projects them into a log-Euclidean tangent space, and augments connectivity-derived features with low-frequency BOLD spectral energy to better resolve subtle state differences under heterogeneous variability. Evaluated on 204 segments from 22 practitioners, ACAM-J-MARS improves performance across multiple backbone models, achieving a balanced accuracy of 0.757 for early versus late ACAM-J classification and 0.545 for three-class classification, corresponding to a 10.8-point absolute improvement over the strongest baseline. These findings suggest that the proposed representation improves the decoding of advanced meditative brain states and may enable downstream applications such as closed-loop neurofeedback and real-time psychophysiological brain-computer interfaces.
Figures & tables
| Representation | Binary (early vs. late) | Ternary (3-level) | ||||||||||
| Method | Atlas | Geometry | Spectral | Classifier | Acc. | BA | F1 | Acc. | BA | F1 | ||
| Unimodal Euclidean connectivity | ||||||||||||
| S200 FC | S200 | Euc. | ✗ | NCM | 0.576 | 0.579 | 0.574 | – | – | – | ||
| S200 FC | S200 | Euc. | ✗ | Ridge | 0.668 | 0.649 | 0.645 | – | – | – | ||
| All-atlas FC | All | Euc. | ✗ | Ridge | 0.659 | 0.641 | 0.637 | – | – | – | ||
| Partial geometric representations | ||||||||||||
Neural Decoding
Neural decoding maps brain activity or error-correction syndromes to stimuli, behavior, language, or corrective actions.
See all 187 papers