PennLINC / xcpEngine

Official public repository for the XCP Engine. This tool is deprecated in favor of XCP-D and ASLPrep.
MIT License
66 stars 42 forks source link

Do I need multiple dsn files for different TRs ? #429

Open Figsandcitrus opened 3 years ago

Figsandcitrus commented 3 years ago

Describe the bug We use a modified fc-36p_spkreg.dsn with fmriprep20.2 results as input. The only varation from original dsn is that I set confound2_censor[2]=0 to avoid scrubbing (see pasted portion below). My question is: Do I have to change the "confound2_framewise [2]" value based on the TR? I am not sure if I can tell if the spike regressor considers the TR for detecting spikes, otherwise the spike regression criteria would be more stringent for TR=2s (.082=.16mm) than TR=3s (.083=.24).

I read #214 but I am not sure if it is suggesting to change fds, dv, or both for TR other than 3. and if fc-36p_spkreg.dsn scheme is much affected by this compared to despiking.

Design File

...
###################################################################
# ANALYSIS VARIABLES tg: changed analysis name and design
###################################################################

analysis=fc_tgholipour
design=/mnt/xcpEngine/xcp_design_files/fc-36p_spkreg_tg.dsn
sequence=anatomical
standard=MNI%2x2x2

###################################################################
# PIPELINE
###################################################################

pipeline=prestats,confound2,regress,fcon,reho,alff,roiquant,norm,qcfc

###################################################################
# 1 PRESTATS
###################################################################

prestats_rerun[1]=1
prestats_cleanup[1]=1
prestats_process[1]=FMP

###################################################################
# 2 CONFOUND2 tg:change confound2_censor[2]= from 1 to 0
###################################################################

confound2_rps[2]=1
confound2_rms[2]=0
confound2_wm[2]=1
confound2_csf[2]=1
confound2_gsr[2]=1
confound2_acompcor[2]=0
confound2_tcompcor[2]=0
confound2_aroma[2]=0
confound2_past[2]=0
confound2_dx[2]=1
confound2_sq[2]=2
confound2_custom[2]=
confound2_censor[2]=0
confound2_censor_contig[2]=0
confound2_framewise[2]=fds:0.08,dv:2
confound2_rerun[2]=1
confound2_cleanup[2]=1
...

Error message None

Runtime Information Singluarity

Additional context Is qcfc utility supposed to help with these problems?

Figsandcitrus commented 3 years ago

@a3sha2 could you please comment in this?

a3sha2 commented 3 years ago

I am sorry @Figsandcitrus, for late response. Yes you are right, you have to change fds, fds:0.083 is based on 3 seconds TR to obtained fd threshold of 0.25 mm. With TR of 2s, it will be fds:0.125 (that is 2*.125= 0.25mm).

let me know it is clear or not