Two step Gaussian mixture model approach to characterize white matter disease based on distributional changes

Background Magnetic resonance imaging reveals macro- and microstructural correlates of neurodegeneration, which are often assessed using voxel-by-voxel t-tests for comparing mean image intensities measured by fractional anisotropy (FA) between cases and controls or regression analysis for associating mean intensity with putative risk factors. This analytic strategy focusing on mean intensity in individual voxels, however, fails to account for change in distribution of image intensities due to disease. New method We propose a method that aims to facilitate simple and clear characterization of underlying distribution. Our method consists of two steps: subject-level (Step 1) and group-level or a specific risk-level density function estimation across subjects (Step 2). Results The proposed method was demonstrated with a simulated data set and real FA data sets from two white matter tracts, where the proposed method successfully detected any departure of the FA distribution from the normal state by disease: p