Dear Users,
I’m trying to create a script to do Multiple ROIs on Multiple
Subjects and create stats table for each ROI/ Subject. I’m using SPM5. I
keep getting the following error and I’m not sure if I’m somehow over
writing my mars estimated design, I hope someone can point me in the right direction.
My script is also posted below, in case I’m doing something blindingly
obvious wrong…
Thank you in advance,
Jac
ERROR
??? Reference to
non-existent field 'betas'.
Error in ==> pr_stat_compute at 39
betas = SPM.betas;
Error in ==> mardo_2.compute_contrasts at 43
[marsS.con marsS.stat, marsS.P, marsS.Pc] = ...
Error in ==> mardo.stat_table at 29
[marsS] = compute_contrasts(marsD, Ic);
Error in ==> v_batch_mars at 116
[rep_strs, marsS, marsD, changef] = stat_table(marsD, 1:length(xCon))
SCRIPT
clear all
% Start marsbar to make sure spm_get works
marsbar('on')
SubsSPMs={'C:\xxx\results\SPM.mat';...
'C:\yyy\results\SPM.mat’;...
};
r1={'C:\xxx\MT+_roi.mat'
'C:\xxx\MT_roi.mat';...
};
r2={'C:\xxx\MT+_roi.mat'};
'C:\xxx\MT_roi.mat';...
};
for
z=1:length(SubsSPMs)
spm_name =SubsSPMs{z};
disp (['Subject ' ,num2str(z), ' ...loading SPM'])
% Make marsbar design object
marsD = mardo(spm_name);
rnam=strcat('r',num2str(z));
eval ([ 'rr = ' rnam ';'])
for i=1:length(rr)
roi_file = rr{i}
% Make marsbar ROI object
R = maroi(roi_file);
% Fetch data into marsbar data object
Y = get_marsy(R, marsD, 'mean');
% Get contrasts from original design
xCon = get_contrasts(marsD);
% Estimate design on ROI data
E = estimate(marsD, Y);
% Put contrasts from original design back into design object
E = set_contrasts(E, xCon);
% get design betas
b = betas(E);
% get stats and stuff for all contrasts into statistics
structure
marsS = compute_contrasts(E, 1:length(xCon));
% create stats table. IC - indices for constrast to be
displayed MarsD -
% Marsbar design structure
%[rep_strs, marsS, marsD, changef] = stat_table(marsD, Ic)
[rep_strs, marsS, marsD, changef] = stat_table(marsD,
1:length(xCon))
end
end
RSS Feed