Hello,

Previously, fMRIPrep computed anatomical CompCor using a single mask generated as the union of eroded WM and CSF masks following the original Behzadi implementation. Starting in v1.4.0, fMRIPrep computes this decomposition, but also separate decompositions within the CSF and WM masks. It returns a sufficient number of components such that 50 percent of the variance in each compartment (CSF, WM, and combined) can be explained (and the newer Muschelli CompCor implementations are also possible). For larger functional datasets (e.g., multiband), it’s possible that CompCor will return hundreds of components.

Currently, you can determine the mask in which each component is computed by examining or querying the associated metadata file, `desc-confounds_regressors.json`

.

- The
`Mask`

field denotes whether the component was generated from the CSF, WM, or combined mask.
- The components will be ordered according to variance explained such that the first 5 components with
`"Mask": "CSF"`

will be the first 5 CSF confounds. If you want to make sure, you can always check the `SingularValue`

or `VarianceExplained`

fields: the first 5 components should have the largest singular values and explain the greatest fraction of variance.

Additional documentation is available here. Please let me know if you have any other questions!