How to caculate how many volumns were scrubbed for each subject using DPARSFA?

Submitted by Yangmei Luo on Tue, 07/07/2015 - 01:49
Hi DPARSFers,
 
I used the "scrubbing"  module in DPARSF  to clean my data. Is there any method to caculate how many volumns were scrubbed for each participant? The DPARSF seems not to report this data. 
 
Thank you in advance for your kind help. Any help is appreciated. 
 
Yangmei Luo

Forums:

YAN Chao-Gan

Tue, 07/07/2015 - 01:58

Hi Yangmei,

You can find "Number of FD_Power>0.2" in HeadMotion.tsv under RealignParameter.

If you use other criterion to scrub the time series, you can use the file FD_Power_SubXXX.txt or FD_Jenkinson_SubXXX.txt under RealignParameter/SubXXX to count the number of time points scrubbed.

Best,

Chao-Gan

Yangmei Luo

Tue, 07/07/2015 - 02:14

Hi Chao-Gan,

Thank you for your quick response. I used the criterion " FD_Jenkinson > 0.5". So,  I need to caculate the number for each subject one by one, although it is a little bit time-comsuming. 

I am greatly appreciated your kind help. 

Best, 

Yangmei

YAN Chao-Gan

Tue, 07/07/2015 - 02:24

Hi Yangmei,

Make a stop point on the beginning (e.g., Line 40) of DPARSFA_run.m

They run DPARSFA, and run the following scripts when matlab stop at debugging:

for i=1:AutoDataProcessParameter.SubjectNum

    cd([AutoDataProcessParameter.DataProcessDir,filesep,'RealignParameter',filesep,AutoDataProcessParameter.SubjectID{i}]);
    FD_Jenkinson=load(['FD_Jenkinson_',AutoDataProcessParameter.SubjectID{i},'.txt']);
    NumberFD_Jenkinson_05(i,1) = length(find(FD_Jenkinson>0.5));
end
 
Best,
Chao-Gan

Yangmei Luo

Tue, 07/07/2015 - 07:24

Hi Chao-Gan, 

Thanks a lot! I got it!