Dear Karen,
when you are looking at results, you can type the following at the
MATLAB prompt:
[N,Z,M,A,XYZ] = spm_max(xSPM.Z,xSPM.XYZ);
res = zeros(SPM.nscan,numel(Z));
for i=1:numel(Z)
Y = spm_get_data(SPM.xY.VY,M(:,i));
KWY = spm_filter(SPM.xX.K,SPM.xX.W*Y);
res(:,i) = spm_sp('r',SPM.xX.xKXs,KWY);
end
it will be looping over all local maxima of the excursion set and
compute the residuals time series in variable 'res'.
Best regards,
Guillaume.
On 27/08/10 14:30, Karen wrote:
> Hi all,
>
> How can I obtain residuals for each peak voxel?
>
> Thanks,
> -k
--
Guillaume Flandin, PhD
Wellcome Trust Centre for Neuroimaging
University College London
12 Queen Square
London WC1N 3BG
|