HI - this is an implementation of the variance normalisation used by MELODIC -
see Beckmann TMI 2004.
The idea is to remove the 30 strongest PCA components (temporarily) and use the
remaining weakest components (ie residual "noise") to compute the voxelwise
scalings for the variance
Hi HCP!
As I was going through the scripts for generating HCP 820-subject PTN
release, there is a step in script "subproc_DO_2_groupPCA.m" that seems
quite confusing to me. The part of code is pasted as follows:
*BO=ciftiopen(f{r(1)},WBC); grot=demean(double(BO.cdata)'); clear BO.cdata;*