Hi List!
I have generated some data for use in a functional connectivity analysis, by temporally filtering a dataset in matlab, then writing analyse files using spm_write_vol/plane.
When I run an analysis on this data, the ist pass goes okay, then the REML estimation. Then the analysis crashes with the following error:
----------------------------
??? Index exceeds matrix dimensions.
Error in ==> spdiags at 114
a((len(k)+1):len(k+1),:) = [i i+d(k) B(i+(m>=n)*d(k),k)];
Error in ==> spm_spm at 420
s = spdiags(1./sqrt(diag(s)),0,nScan,nScan);
Error in ==> spm_spm at 827
SPM = spm_spm(SPM);
--------------------------------
I reran the same analysis on the unfiltered data and it ran okay.
The only difference I could detect was that when spm_svd runs on xVi.V on line 419 in spm_spm.m, the resulting "s" matrix is not nScans x nScans in size, leading to the crash.
I don't see an obvious way to check xVi.V in each case, or have i corrupted my data in the filtering?
any help much appreciated!
Tony
***************************************
Tony Waites, PhD
MR physicist
MR section
University Hospital in Lund
22185 Lund
SWEDEN
tel work: +46 46 17 7027
fax : +46 46 12 1837
email : [log in to unmask]
|