Hi,
>nope, they should do the same thing. just to check: have a look at
>report.log in the ica folder and at report.com in the feat folder - the
>mcflirt command line call should be the same.
They're almost the same:
FEAT:
mcflirt -in func_BET_del4 -out prefiltered_func_data_mcf -plots -refvol 61 -rmsrel -
rmsabs
Melodic:
mcflirt -in func_BET_del4 -out prefiltered_func_data_mcf -plots -refvol 61
What are the -rmsrel and -rmsabs options? I couldn't find any reference to them in
the manual.
>You can also have a look at the ica 00index.html file: how close are
>the estimated evidence for dim23 and dim 24?
The "eigenspectrum analysis" graphs look nearly identical.
I did the following experiment on 6 different data sets:
1. take func_data and run bet
2. delete first 4 volumes using avwsplit and avwmerge
then one of three operations:
3a. command line mcflirt (default options) and Melodic (MCFLIRT, BET & STC = "off",
delete 0 vols).
3b. Melodic (MCFLIRT="on", BET & STC = "off", delete 0 vols).
3c. FEAT (prestats only, MCFLIRT="on", BET & STC = "off", FWHM=0, highpass="off",
delete 0 vols) and Melodic (MCFLIRT, BET, & STC = "off", delete 0 vols).
3a and 3b gave me the same answer every time. 3c gave approximately the same
number of ICs, sometimes 1 more or 1 less or equal. However the components
themselves were different. For example, in one data set, a component which I
believe was visual cortex activation was component 4 in 3a/3b but component 18 in
3c. The timecourses of all the components were different (although the amount of
difference varied by component). Also, some components were present in one
analysis but not in the other.
The only other thing I can think of is that the "High pass filter cutoff" = 100s in my
FEAT analysis. But it shouldn't be applied since I'm only doing prestats, and in any
case, I turned temporal filtering off. What am I doing wrong?
thanks,
jack
|