Dear all,
I am trying to estimate smoothness for each subject with the following steps (https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind1505&L=spm&D=0&P=558407), however, output values are too high.
1. calculate the square root of the ResMs.img:
spm_imcalc_ui([pwd,'\ResMS.img'],[pwd,'\sqrt_ResMS.img'],'sqrt(i1)');
2. estimate smoothness:
FWHM=spm_est_smoothness([pwd,'\sqrt_ResMS.img'],[pwd,'\mask.img'],[624,1]); % two runs, 312 scans in each, this makes n=624
%mask.img was generated at the first-level analysis for a given subject
The output of these steps are pretty high (as compared to the values in SPM.xVol.FWHM = 4.5, 4.5, 4.4):
FWHM =
58.0763 67.4205 60.5713.
In addition, according to a previous thread (https://www.jiscmail.ac.uk/cgi-bin/webadmin?A2=ind1506&L=spm&P=R1499&1=spm&9=A&J=on&d=No+Match%3BMatch%3BMatches&z=4), using ResMS.img will result in larger FWHM estimates, however, this seems not the case here:
FWHM =
21.2398 24.9775 22.9709
using: FWHM=spm_est_smoothness([pwd,'\ResMS.img'],[pwd,'\mask.img'],[624,1]).
Did I use spm_imcalc_ui or spm_est_smoothness in a wrong way? Any comments will be appreciated!
Thanks,
Chunliang