Dear all, 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.U.name <http://dcm.u.name/> = {'null'}; 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'); Any help will be highly appreciated best Rabia