Dear Guillaume
Thank you very much for your suggestion. I enabled the "write residuals"
matlabbatch{2}.spm.stats.fmri_est.write_residuals = 1;
I just want to make sure that I'm on the right path. To extract residuals, I've to go to the voxels I'm interested in in each Res_xxxx.nii (xxxx=volume) and extract Y. so, here is the script:
V = spm_vol('Res_0600.nii');
[Y,XYZ] = spm_read_vols(V);
[Y] = spm_get_data(V,XYZ,true); % Y = residuals for all voxels in volume 600, right?
is this correct???
Thanks for your help
Maaly
|