Maybe the following scripts will help?
load('SPM.mat');
% contrasts
%---------------------------------------------------------------------------
nContrast = length(SPM.xCon);
if nContrast == 0
SPM.xCon = spm_FcUtil('Set',cname,FT,'c',c(:),SPM.xX.xKXs);
else
SPM.xCon(nContrast + 1) = spm_FcUtil('Set',cname,FT,'c',c(:),SPM.xX.xKXs);
end
index = nContrast + 1;
% and evaluate
%---------------------------------------------------------------------------
spm_contrasts(SPM, nContrast + 1);
--
----------------------------------------------------------------
Xu Cui
Department of Psychiatry and Behavioral Sciences
Stanford University
Blog:
http://www.alivelearn.net
----------------------------------------------------------------