SurfStat: A Matlab toolbox for the statistical analysis of univariate and multivariate surface and volumetric data using linear mixed effects models and random field theory
Introduction
SurfStat does statistical analysis of univariate and multivariate surface data using linear mixed effects models (fitted by ReML) and random field theory (RFT). SurfStat is intended for cortical thickness data on triangular meshes, but it will handle any triangulated surface data; the only requirement is that the triangulation scheme must be the same for all surfaces, i.e. the data must be registered to a common surface.
Section snippets
Methods
Inference uses RFT for T, F, Hotelling's T2 and Roy's maximum root statistics. Model formulae rather than design matrices specify the linear model. It is fast, because everything is loaded into memory, permitting a truly interactive analysis. Finally, off-the-shelf Matlab graphics are ready to publish.
It will also analyse volumetric data, e.g. VBM, DBM and PET data. To cope with increased memory demands, while maintaining fast execution, memory mapping is implemented.
Results
SurfStat will handle registered FreeSurfer data, but here we shall analyse 147 subjects from the ICBM data base. The scientific question is if cortical thickness is related to gender or age. Figure 1 shows data on a single subject.
Cortical thickness has been read into a single subjects(147) × vertices(81920) matrix Y, and gender and age have been read into subjects(147) × 1 vectors gender (values ‘male’ and ‘female’) and age (in years). We now convert these to terms that can be combined into a
Conclusions
SurfStat is a very flexible tool for the statistical analysis of surface or volumetric data via linear mixed effects models.