An F-ratio-based method for estimating the number of active sources in MEG

Abstract

IntroductionMagnetoencephalography (MEG) is a powerful technique for studying the human brain function. However, accurately estimating the number of sources that contribute to the MEG recordings remains a challenging problem due to the low signal-to-noise ratio (SNR), the presence of correlated sources, inaccuracies in head modeling, and variations in individual anatomy.MethodsTo address these issues, our study introduces a robust method for accurately estimating the number of active sources in the brain based on the F-ratio statistical approach, which allows for a comparison between a full model with a higher number of sources and a reduced model with fewer sources. Using this approach, we developed a formal statistical procedure that sequentially increases the number of sources in the multiple dipole localization problem until all sources are found.ResultsOur results revealed that the selection of thresholds plays a critical role in determining the method's overall performance, and appropriate thresholds needed to be adjusted for the number of sources and SNR levels, while they remained largely invariant to different inter-source correlations, translational modeling inaccuracies, and different cortical anatomies. By identifying optimal thresholds and validating our F-ratio-based method in simulated, real phantom, and human MEG data, we demonstrated the superiority of our F-ratio-based method over existing state-of-the-art statistical approaches, such as the Akaike Information Criterion (AIC) and Minimum Description Length (MDL).DiscussionOverall, when tuned for optimal selection of thresholds, our method offers researchers a precise tool to estimate the true number of active brain sources and accurately model brain function.

Localization of Spatially Extended Brain Sources by Flexible Alternating Projection (Flex-AP)

“…To simulate the source time courses, we first generated two random sequences following a 1 /f frequency spectrum. Subsequently, we generated the time courses of these two sources by mixing the sequences in a manner that adhered to a moderate inter-source correlation with a coefficient of ρ = 0.5 (as referenced in [14]). Within each condition, we examined these original sources with no extent (discrete source), but also with an extent of a given neighborhood order (extended source).…”