The model parameters are obtained according to a maximum likelihood criterion employing the expectation maximization (EM) algorithm.
An averaged control subject image is computed to obtain a threshold-based mask that selects only the voxels inside the skull.