Statistical flattening of MEG beamformer images

We propose a method of correction for multiple comparisons in MEG beamformer based Statistical Parametric Maps (SPMs). We introduce a modification to the minimum-variance beamformer, in which beamformer weights and SPMs of source-power change are computed in distinct steps. This approach allows the calculation of image smoothness based on the computed weights alone. In the first instance we estimate image smoothness by looking at local spatial correlations in residual images generated using random data; we then go on to show how the smoothness of the SPM can be obtained analytically by measuring the correlations between the adjacent weight vectors. In simulations we show that the smoothness of the SPM is highly inhomogeneous and depends on the source strength. We show that, for the minimum variance beamformer, knowledge of image smoothness is sufficient to allow for correction of the multiple comparison problem. Per-voxel threshold estimates, based on the voxels extent (or cluster size) in flattened space, provide accurate corrected false positive error rates for these highly inhomogeneously smooth images.