Glow builds brain regions linked to variables via improved clustering

GLOW: General Linear Models Optimized with Ward's Method

Machine Learning

Summary

Finding brain areas connected to specific traits often uses simple models per tiny spots in brain images, which can miss broader patterns. The authors propose GLOW, which groups nearby brain spots by how well linear models fit, making regions that better capture effects. Tested on brain scans and fake signals, GLOW finds more complete areas but may include extra nearby parts. It runs slower but shows promise in detecting linked brain regions more reliably than some existing methods.

What this means in practice

  • For neuroimaging analysis teams: Produce more spatially consistent brain regions associated with variables of interest by clustering voxels minimizing GLM error.
  • For medical image analysts: Improve detection of disease-related brain regions by applying GLOW to imaging datasets with spatial hierarchy control.

Authors

Matt Higger

Abstract

Neuroimaging seeks brain regions whose imaging features correlate with variables of interest (e.g., diagnosis, age), usually via the General Linear Model (GLM). Region-Based Analysis dilutes or splits any effect that does not fill exactly one atlas region. Voxel-Based Analysis, with Cluster Extent Thresholding and Threshold-Free Cluster Enhancement built on it, fits a GLM per voxel and groups adjacent significant voxels. But neighboring voxels can share a statistic without sharing a GLM, which a per-voxel statistic cannot detect. GLOW (General Linear Models Optimized with Ward's method) instead builds a hierarchy of nested regions by greedily minimizing GLM error (provably a special case of Ward's clustering). Maximum-statistic (Westfall-Young) permutation controls family-wise error over all regions. Significance spreads to an effect's nested fragments and the regions containing it; a final pruning step seeks whole, pure effect extents. We plant synthetic effects of known extent and varying strength in Human Connectome Project (HCP) diffusion images and in white noise. Pruning delivers wholeness: GLOW returns a spatially consistent effect as one region on 57% of HCP trials, the best baseline on 35%; baselines recover comparable volume, but in many pure pieces. The cost is positive predictive value (PPV): GLOW is more sensitive (by up to 0.29), but its regions overrun the effect, mostly because no hierarchy region isolates it, not because of pruning. On HCP, GLOW's mean Dice gains 0.059 over the best baseline at weaker effects and gives back 0.005 at stronger ones. On white noise, GLOW trails by up to 0.61 Dice: its hierarchy holds the effect, but its permutation test does not certify it. Runtime grows roughly linearly with voxel count, about 4.5 times the voxel-wise methods' on the full HCP volume. The PPV loss may trace to the assumption of spatially independent residuals.