I am trying to use stochastic DCM (I only want to use stochastic, and not spectral DCM at this stage). However I m uncertain about certain issues. Are the options I am writing below are OK ? I am not sure, what should I put in options.analysis.
DCM.Y.Q = spm_Ce(ones(1,n)*v);
DCM.U.u = zeros(v,1);
DCM.a = ones(n,n);
DCM.b = zeros(n,n,0);
DCM.c = zeros(n,0);
DCM.d = zeros(n,n,0);
DCM.TE = 0.04;
DCM.delays = repmat(RT,DCM.n,1);
DCM.options.nonlinear = 0;
DCM.options.two_state = 0;
DCM.options.stochastic = 1;
options.analysis = 'CSD';
% DCM.options.nograph = 1;
str = sprintf('DCM_DMN');
DCM.name = str;
save(fullfile(glmdir,str),'DCM');