Dear FSL Experts,
Thank you in advance with advice on this problem. I am having a bit of a confusing problem. It is as such:
I have a 2x2 ANOVA experimental design for a rs-fMRI experiment (cases vs. controls ; pre-intervention vs. post-intervention). I ran MELODIC with d=30 components, and then ran a dual regression. When I check my melodic_IC.nii.gz (2 scans per case, 2 scans per control concatenated) file against reference networks (using fslcc, with -t=.204) provided from Yeo2011 and Smith 2009 they do not really correlate with any of the reference networks.
So I figured that I probably shouldn't collapse across cases and controls when I am identifying RSNs. So I took out the cases (and left pre-intervention and post intervention) and ran MELODIC again. This outputted melodic_IC_controls.nii.gz. When I use fslcc to compare this file with reference networks, I got quite a bit of correlation.
My main question is A) Is it ok to do this just to identify which components are which types of networks and B) If I use the original melodic_IC.nii.gz (2 scans per case, 2 scans per control concatenated) for the dual_regression, do the component numbers still correlate with the networks i've identified using melodic_IC_controls.nii.gz? That's my biggest problem.
Thank you so much for your help! I'm extremely new to fMRI so any help with this matter would be greatly appreciated.
Cheers,
Pranish Kantak
|