Dear DCM-experts, I am experimenting with the estimator of DCM (spm_nlsi_GN.m) and I am wondering what is the best way to change the starting position of the EM while preserving the priors for the estimation. I implemented the following hack. It would be great if anybody could confirm that this is correct. In spm_nlsi_GN.m I added the following lines just in front of the main loop ... % EM %========================================================================== C.F = -Inf; % free energy v = 0; % log ascent rate dFdh = zeros(nh,1); dFdhh = zeros(nh,nh); %%%%%%%%%%%%%%%%%%%%%%%%% % Change start position load newStart # load the variable 'start' Ep = start; # changes the MAP p = [V'*(spm_vec(start) - spm_vec(M.pE)); Eu]; # adjusts the deviation from the prior expectations %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % end of Change start position for k = 1:64 % time %---------------------------------------------------------------------- Ti = clock; .... Thank you, Martin -- SPM for programmers http://spm.martinpyka.de