Thanks for creating this toolbox. I was thrilled to find it. 

I am not very good with Matlab and am new to this form of analysis so apologies if this is a stupid question.

I am running the first level analysis on a bunch of participants and have just adapted your anadef script to do this. However, for many of the subjects, I am finding that I get the following error:

Retrieving data from ROI 17 using summary function mean ...
Summarizing data with summary function: mean
Initialising parameters : ...done
Covariance estimate : ...voxelwise
Model : estimation
Temporal non-sphericity (over voxels) : ...estimation
ReML Block - 1
ReML Iteration : 1 ...1.134073e+01
ReML Iteration : 2 ...2.286396e+04
ReML Iteration : 3 ...5.594243e+00
ReML Iteration : 4 ...6.116961e+03
ReML Iteration : 5 ...3.742198e+00
ReML Iteration : 6 ...8.546006e+00
ReML Iteration : 7 ...6.250731e-01
ReML Iteration : 8 ...4.084047e-01
ReML Iteration : 9 ...8.420834e-02
ReML Iteration : 10 ...3.495538e-02
ReML Iteration : 11 ...9.900201e-03
Initialising parameters : ...computingIndex exceeds matrix dimensions.

Error in spdiags (line 102)
a((len(k)+1):len(k+1),:) = [i i+d(k) B(i+(m>=n)*d(k),k)];
Error in pr_estimate (line 128)
s = spdiags(1./sqrt(diag(s)),0,nScan,nScan);
Error in pr_estimate (line 345)
SPM = pr_estimate(SPM,marsY);
Error in mardo_5/estimate (line 72)
SPM = pr_estimate(SPM, marsY);
Error in BASCO>pushbuttonmodelspecest_Callback (line 424)
E = estimate(D,Y); % estimate model based on ROI summary
Error in gui_mainfcn (line 95)
feval(varargin{:});
Error in BASCO (line 20)
gui_mainfcn(gui_State, varargin{:});
Error in
matlab.graphics.internal.figfile.FigFile/read>@(hObject,eventdata)BASCO('pushbuttonmodelspecest_Callback',hObject,eventdata,guidata(hObject))
Error while evaluating UIControl Callback

It is difficult to trouble shoot without being able to see what SPM is doing and I wondered if you have any ideas? I have different trial numbers for each condition in each case as it is a correct vs. incorrect trial analysis. It also seems to be related to different ROIs in each case .

Thanks in advance

I got the same error! Any help would be appreciated!

Hi,

this error was reported several times in the marsbar forum. Matthew Brett suggested to use K. Worsleys fmristat instead of the autocorrelation implementation in SPM. For that you must insert the following code in the file basco_checkestimation on line 16.

D = autocorr(D, 'fmristat', 2);

I have also noticed that sometimes I get the error "Check number of ROIs in txt-file"
I changed line 225 in BASCO.m from
if length(thenames)~=ROINum
to
if size(thenames,1)~=ROINum

and it worked.

Best
Ralf

Thanks a lot for the code snippet, Ralf. You are awesome!


also the initial source of info
http://marsbar.sourceforge.net/faq.html#...


There is also a function called autocorr from econ toolbox. Type "which autocorr" to see which autocorr is working.
if necessary, one can copy @mardo_2/autocorr.m 

Originally posted by Ralf Veit:

Hi,
I made the changes to the file basco_checkestimation on line 16 and to line 225 in BASCO.m but I still receive the following error message for some subjects:

Retrieving data from ROI 57 using summary function mean ...
Summarizing data with summary function: mean
Re-estimating covariance
Re-estimating whitening filter
Initialising parameters : ...done
Covariance estimate : ...voxelwise
Model : estimation
Temporal non-sphericity (over voxels) : ...estimation
ReML Block - 1
ReML Iteration : 1 ...3.410752e+25
ReML Iteration : 2 ...2.828835e+03
ReML Iteration : 3 ...4.421080e+01
ReML Iteration : 4 ...1.613961e+06
ReML Iteration : 5 ...4.180238e+01
ReML Iteration : 6 ...4.352357e+06
ReML Iteration : 7 ...4.636968e+01
ReML Iteration : 8 ...4.689308e+01
ReML Iteration : 9 ...1.634696e+04
ReML Iteration : 10 ...5.011224e+01
ReML Iteration : 11 ...3.862712e+03
ReML Iteration : 12 ...8.503537e+01
ReML Iteration : 13 ...5.833874e+02
ReML Iteration : 14 ...8.889501e+02
ReML Iteration : 15 ...1.863922e+02
ReML Iteration : 16 ...4.274574e+02
ReML Iteration : 17 ...1.672121e+03
ReML Iteration : 18 ...5.110463e+05
ReML Iteration : 19 ...4.828098e+01
ReML Iteration : 20 ...4.704126e+01
ReML Iteration : 21 ...6.273130e+05
ReML Iteration : 22 ...4.188783e+01
ReML Iteration : 23 ...7.261460e+03
ReML Iteration : 24 ...5.100982e+01
ReML Iteration : 25 ...8.972958e+03
ReML Iteration : 26 ...1.505456e+06
ReML Iteration : 27 ...1.283097e+02
ReML Iteration : 28 ...2.008960e+03
ReML Iteration : 29 ...2.246124e+03
ReML Iteration : 30 ...3.719424e+03
ReML Iteration : 31 ...7.797202e+06
ReML Iteration : 32 ...6.124269e-03
ReML Block - 2
ReML Iteration : 1 ...6.028477e+24
ReML Iteration : 2 ...6.179355e+02
ReML Iteration : 3 ...9.867692e+02
ReML Iteration : 4 ...2.847287e+02
ReML Iteration : 5 ...3.228575e+02
ReML Iteration : 6 ...2.212133e+02
ReML Iteration : 7 ...3.051855e+04
ReML Iteration : 8 ...1.003979e+01
ReML Iteration : 9 ...1.422214e+03
ReML Iteration : 10 ...5.317306e+01
ReML Iteration : 11 ...6.749793e+00
ReML Iteration : 12 ...1.624683e+06
ReML Iteration : 13 ...1.148070e+05
ReML Iteration : 14 ...1.237820e+11
ReML Iteration : 15 ...3.307808e+10
ReML Iteration : 16 ...2.169239e+01
ReML Iteration : 17 ...2.275552e-01
ReML Iteration : 18 ...2.411023e-02
ReML Iteration : 19 ...8.187992e-03
ReML Block - 3
ReML Iteration : 1 ...8.628086e+24
ReML Iteration : 2 ...1.992343e-01
ReML Iteration : 3 ...1.852284e-01
ReML Iteration : 4 ...2.820160e-01
ReML Iteration : 5 ...1.189589e-01
ReML Iteration : 6 ...3.375698e-01
ReML Iteration : 7 ...9.604632e-01
ReML Iteration : 8 ...3.161417e-01
ReML Iteration : 9 ...4.307677e-01
ReML Iteration : 10 ...2.068754e-01
ReML Iteration : 11 ...3.447115e-01
ReML Iteration : 12 ...7.117137e-01
ReML Iteration : 13 ...2.499299e-01
ReML Iteration : 14 ...2.532627e-02
ReML Iteration : 15 ...1.525531e-01
ReML Iteration : 16 ...5.933029e-01
ReML Iteration : 17 ...6.832739e-01
ReML Iteration : 18 ...2.374132e+00
ReML Iteration : 19 ...3.896834e-01
ReML Iteration : 20 ...1.498764e-02
ReML Iteration : 21 ...5.230396e-01
ReML Iteration : 22 ...6.657243e-01
ReML Iteration : 23 ...1.587401e+00
ReML Iteration : 24 ...2.196471e-02
ReML Iteration : 25 ...3.874127e-01
ReML Iteration : 26 ...9.265903e-02
ReML Iteration : 27 ...6.606886e-01
ReML Iteration : 28 ...1.192551e+00
ReML Iteration : 29 ...9.113033e-01
ReML Iteration : 30 ...5.524773e-02
ReML Iteration : 31 ...7.925981e-01
ReML Iteration : 32 ...6.872270e-01
ReML Iteration : 33 ...5.076939e-01
ReML Iteration : 34 ...1.386557e-01
ReML Iteration : 35 ...5.097154e-01
ReML Iteration : 36 ...2.239037e-02
ReML Iteration : 37 ...8.979627e-02
ReML Iteration : 38 ...7.340414e-02
ReML Iteration : 39 ...1.083122e-01
ReML Iteration : 40 ...2.579397e-01
ReML Iteration : 41 ...2.568578e-01
ReML Iteration : 42 ...1.010973e+00
ReML Iteration : 43 ...3.614601e+00
ReML Iteration : 44 ...1.057727e+00
ReML Iteration : 45 ...1.504028e-01
ReML Iteration : 46 ...2.351308e-02
ReML Iteration : 47 ...1.543476e-01
ReML Iteration : 48 ...3.993500e-01
ReML Iteration : 49 ...2.469874e-01
ReML Iteration : 50 ...4.375068e-01
ReML Iteration : 51 ...2.231699e-01
ReML Iteration : 52 ...4.365207e-01
ReML Iteration : 53 ...7.082969e-01
ReML Iteration : 54 ...4.793554e-01
ReML Iteration : 55 ...9.630564e-02
ReML Iteration : 56 ...3.928655e-01
ReML Iteration : 57 ...1.110299e+00
ReML Iteration : 58 ...1.810996e+00
ReML Iteration : 59 ...1.407834e+00
ReML Iteration : 60 ...1.211868e-01
ReML Iteration : 61 ...3.618090e-01
ReML Iteration : 62 ...2.167922e-01
ReML Iteration : 63 ...1.829130e-01
ReML Iteration : 64 ...8.644111e-01
ReML Block - 4
ReML Iteration : 1 ...6.866919e+24
ReML Iteration : 2 ...5.011138e+00
ReML Iteration : 3 ...2.151095e-01
ReML Iteration : 4 ...5.406613e-02
ReML Iteration : 5 ...8.375658e-01
ReML Iteration : 6 ...1.334521e+00
ReML Iteration : 7 ...7.013734e-02
ReML Iteration : 8 ...1.285255e+00
ReML Iteration : 9 ...1.387347e+00
ReML Iteration : 10 ...1.316738e+00
ReML Iteration : 11 ...4.913568e-02
ReML Iteration : 12 ...4.446727e+00
ReML Iteration : 13 ...3.673783e-01
ReML Iteration : 14 ...5.815593e-01
ReML Iteration : 15 ...6.456046e-01
ReML Iteration : 16 ...8.945228e-01
ReML Iteration : 17 ...1.561681e+00
ReML Iteration : 18 ...4.379552e-01
ReML Iteration : 19 ...4.910517e-03
ReML Block - 5
ReML Iteration : 1 ...2.162717e+24
ReML Iteration : 2 ...8.876681e+03
ReML Iteration : 3 ...2.196735e+01
ReML Iteration : 4 ...3.817622e+01
ReML Iteration : 5 ...1.270137e+05
ReML Iteration : 6 ...4.325122e+01
ReML Iteration : 7 ...5.862207e+04
ReML Iteration : 8 ...3.096524e+02
ReML Iteration : 9 ...1.539215e+03
ReML Iteration : 10 ...7.389998e+05
ReML Iteration : 11 ...5.132941e+02
ReML Iteration : 12 ...2.444608e+02
ReML Iteration : 13 ...1.067493e+04
ReML Iteration : 14 ...2.140501e+06
ReML Iteration : 15 ...4.237760e-03
Initialising parameters : ...computingIndex exceeds matrix dimensions.

Error in spdiags (line 114)
a((len(k)+1):len(k+1),:) = [i i+d(k) B(i+(m>=n)*d(k),k)];
Error in pr_estimate (line 128)
s = spdiags(1./sqrt(diag(s)),0,nScan,nScan);
Error in pr_estimate (line 345)
SPM = pr_estimate(SPM,marsY);
Error in mardo_5/estimate (line 72)
SPM = pr_estimate(SPM, marsY);
Error in BASCO>pushbuttonmodelspecest_Callback (line 423)
E = estimate(D,Y); % estimate model based on ROI
summary
Error in gui_mainfcn (line 96)
feval(varargin{:});
Error in BASCO (line 20)
gui_mainfcn(gui_State, varargin{:});
Error in
@(hObject,eventdata)BASCO('pushbuttonmodelspecest_Callback',hObject,eventdata,guidata(hObject))
Error while evaluating uicontrol Callback


Thank you in advance for your help!