Thank you very much for your answer!
According to your idea, I tried to change DCM.b and DCM.c in the batch script so that only the first (direct input) and second (modulatory input) condition and the others are set to zero (the model contains 4 regions):
DCM.b = zeros(4,4,6); DCM.b(3,1,2) = 1;
DCM.c = [1 0 0 0 0 0 ; 0 0 0 0 0 0; 0 0 0 0 0 0; 0 0 0 0 0 0];
However, DCM.U still contains all the conditions. Moreover, the DCM.mat cannot be estimated anymore:
??? Error using ==> inlineeval at 15
Error in inline expression ==> spm_int(P,M,U)
??? Index exceeds matrix dimensions.
Error in ==> inline.feval at 36
INLINE_OUT_ = inlineeval(INLINE_INPUTS_, INLINE_OBJ_.inputExpr, INLINE_OBJ_.expr);
Error in ==> spm_diff at 80
f0 = feval(f,x{:});
Error in ==> spm_nlsi_GN at 283
[dfdp f] = spm_diff(IS,Ep,M,U,1,{V});
Error in ==> spm_dcm_estimate at 184
[Ep,Cp,Eh,F] = spm_nlsi_GN(M,U,Y);
Error in ==> spm_dcm_ui at 135
spm_dcm_estimate(P{i});
??? Error while evaluating uicontrol Callback
Sorry for all these basic questions, unfortunately I am not very familiar with matlab….
|