1. Introduction

Dual permutation (DP) is a method to compute the p-value of an inter-modal spatial Pearson correlation coefficient (PCC) of cross-subject mean maps. For example, the result of an fMRI study and a mean map of serotonin receptor density as measured by PET. Each modality is represented by a "stack" of subjects, ie, a map with multiple frames where each frame is a subject. The mean maps to be correlated are computed from these stacks by fitting a regression model in a GLM, then computing a contrast of regression coefficients to reduce the stack to a single univariate map. Dual permutation works by creating "surrogate" maps of each modality and then computing the correlation between the surrogates. The surrogate maps are constructed by first residualizing the stack using the GLM, then permuting the design matrix, then reanalyzing the resdiualized stack with the permuted design matrix, followed by applying the contrast to create a univariate surrogate. This is repeated many times to create a list of surrogate PCCs. The p-value is then computed at the fraction of times the surrogate PCCs exceeded the real PCC. DP should only be used in cases where the two modalities come from different subject populations.

DP has been rigorously verified. The surrogate PCC has been shown analytically to have the same variance of that of the real PCC under the null regardless of spatial autocorrelation function (SACF, including spatial non-staionarity of the SACF), spatial non-stationarity in noise power, and group level design. As a non-parametric method, it does not rely on assumptions about the distribution of the noise. DP has been shown to control the false positive rate (also known as the type I error rate) in pairs of real data configured to be null, including, fMRI, cortical thickness, cortical curvature, DTI fractional anisotropy, and serotonin transporter binding potential as measured by PET. DP works in cases where both modalities are null (known as the "Full Null") as well as when one of the modalities is null but the other is non-null (known as the "Partial Null").

mri_dualperm is the FreeSurfer program that implements DP. It can be run in several ways:

1.1. Overview

1.2. How to cite

1.3. Example

mri_dualperm