Hello, i was wondering if it makes sense to regress out the "noise component" time series from an FSL melodic spatial ICA decomposition (ie, the file created by melodic called "melodic_mix") using 3dTproject. this file contains a time series associated to each spatial ICA component. i would like to regress out some of the melodic ICA components.

typically, i use the fsl command (fsl_regfilt) for subtracting components from the BOLD signal, but i would also like to remove the global signal and do bandpass filtering in the same step, using 3dTproject (i learned recently that it was not mathematically correct to regress out these components sequentially). i'm aware that it's not good practice to combine different software packages in a single analysis pipeline, but fsl_regfilt doesn't allow for bandpass filtering in a single step, whereas 3dTproject does.

if anyone here has experience with fsl's melodic, do you think what i am proposing sounds correct? thank you,

Russell