去除协变量

One is for voxel-wise, the other is for ROI-wise.

 0356 function Result =Brain4D_RegressOutCovariables(ABrain4D, ABasisFunc)

0357 %20070926, Regress some covariables out first
0358     %Result =( E - X(X'X)~X')Y
0359     [nDim1, nDim2, nDim3, nDim4]=size(ABrain4D);
0360     
0361     %Make sure the 1st dim of ABasisFunc is nDim4 long
0362     if size(ABasisFunc,1)~=nDim4, error('The length of Covariables don''t match with the volume.'); end
0363     
0364     % (1*sampleLength) A matrix * B matrix (sampleLength * VoxelCount)
0365     ABrain4D =reshape(ABrain4D, nDim1*nDim2*nDim3, nDim4)';
0366     Result =(eye(nDim4) - ABasisFunc * inv(ABasisFunc' * ABasisFunc)* ABasisFunc') * ABrain4D;
0367     %20071102 Bug fixed squeeze must not be excluded because nDim1 may be ONE !!!
0368     %Result =squeeze(reshape(Result', nDim1, nDim2, nDim3, nDim4));
0369     Result =reshape(Result', nDim1, nDim2, nDim3, nDim4);
0370 
0371 function Result =Brain1D_RegressOutCovariables(ABrain1D, ABasisFunc)
0372 %20070926, Regress some covariables out first
0373     %Result =( E - X(X'X)~X')Y
0374     %Make sure the input is a column vector
0375     % ABrain1D =reshape(ABrain1D, prod(size(ABrain1D)), 1);
0376     
0377     %Make sure the 1st dim of ABasisFunc is nDim4 long
0378     if size(ABasisFunc,1)~=length(ABrain1D), error('The length of Covariables don''t match with the volume.'); end
0379     
0380     % (1*sampleLength) A matrix * B matrix (sampleLength * VoxelCount)
0381     Result =(eye(size(ABrain1D, 1)) - ABasisFunc * inv(ABasisFunc' * ABasisFunc)* ABasisFunc') * ABrain1D;