Submitted by zqbaxfgg on Mon, 04/20/2009 - 15:59
我想请问一下rest中fc.m关于去除头动干扰的源码程序
if ACovariablesDef.polort>=0,
thePolOrt =(1:sampleLength)';
thePolOrt =repmat(thePolOrt, [1, (1+ACovariablesDef.polort)]);
for x=1:(ACovariablesDef.polort+1),
thePolOrt(:, x) =thePolOrt(:, x).^(x-1) ;
end
end
% switch ACovariablesDef.polort,
% case 0,
% thePolOrt =ones(sampleLength, 1);
% case 1,
% thePolOrt =[ones(sampleLength, 1), (1:sampleLength)'];
% case 2,
% thePolOrt =[ones(sampleLength, 1), (1:sampleLength)', (1:sampleLength)'.^2];
% otherwise
% end
%Regress out the covariables
theCovariables =[thePolOrt, theCovariables];
举出的情况说明ACovariablesDef.polort可能为0,1,2,....
而在fc.gui.m中有
handles.Covariables.polort =1;
我想知道ACovariablesDef.polort是什么
如果自己编程序的话,
theCovariables具体是什么
Submitted by admin on Mon, 04/20/2009 - 18:26 Permalink
Reference from AFNI's 3dfim+ command
"fc.m" came from AFNI's 3dfim+ command, and "fc.gui.m" only considers the simplest configuration (a linear trend covariable including constant) for most common users. They should generate the same result as AFNI's with "3dfim+ -polort=1".
As for you, you could read the codes and I think you could also re-script it to meet your situation, even re-write the codes.
And if you rewrite some codes for a specific situation, sharing your codes by this site would be appreciated.
Add Matlab codes with some colors is simple, just change the defalt editor of this site to "Switch to plain text editor", then surround your codes with <blockcode language='Matlab'>
...
</blockcode>
Please refer "http://restfmri.net/forum/history#comment-13" for more.
Submitted by zqbaxfgg on Tue, 04/21/2009 - 12:57 Permalink
呵呵,我看源代码只是想找出回归怎么做的,基本也看不太懂,只
呵呵,我看源代码只是想找出回归怎么做的,基本也看不太懂,只是想引用一些稍作修改,用以现在自己想要的。