Scrubbing function

Submitted by YukiSakai on

Dear DPARSFA experts,

 Thank you very much for the powerful tool. I have several questions.

 1. I cannot find the function to calcurate the FD of Power's method in DPARSFA folder. Where is it and what is the name of function?

 2. I can easily recognize the propotion of excluded frame by HeadMotion.tsv when the threshold is 0.5 or 0.2. In the case of other settings of threshold, how can I get that information?

Best

 Yuki

YAN Chao-Gan

Wed, 11/25/2015 - 07:51

Permalink

Hi Yuki,

Please see Line 827 in DPARSFA_run for calculating FD Power.

rpname=dir([FunSessionPrefixSet{iFunSession},'rp*']);

RP=load(rpname.name);

%Calculate FD Power (Power, J.D., Barnes, K.A., Snyder, A.Z., Schlaggar, B.L., Petersen, S.E., 2012. Spurious but systematic correlations in functional connectivity MRI networks arise from subject motion. Neuroimage 59, 2142-2154.) 
            RPDiff=diff(RP);
            RPDiff=[zeros(1,6);RPDiff];
            RPDiffSphere=RPDiff;
            RPDiffSphere(:,4:6)=RPDiffSphere(:,4:6)*50;
            FD_Power=sum(abs(RPDiffSphere),2);
            save([FunSessionPrefixSet{iFunSession},'FD_Power_',AutoDataProcessParameter.SubjectID{i},'.txt'], 'FD_Power', '-ASCII', '-DOUBLE','-TABS');
            MeanFD_Power = mean(FD_Power);
            
            NumberFD_Power_05 = length(find(FD_Power>0.5));
            PercentFD_Power_05 = length(find(FD_Power>0.5)) / length(FD_Power);
            NumberFD_Power_02 = length(find(FD_Power>0.2));
            PercentFD_Power_02 = length(find(FD_Power>0.2)) / length(FD_Power);
 
 
If you change 0.5 or 0.2 here, you can get the related information.
 
Best,
 
Chao-Gan
Forums
DPABI/DPABISurf/DPARSF