MBB-team / VBA-toolbox

The VBA toolbox
GNU General Public License v3.0
129 stars 67 forks source link

VBA_bDCM_lesion - Does not support additional Gaussian Sources #21

Closed SRSteinkamp closed 7 years ago

SRSteinkamp commented 7 years ago

Dear all,

I found a small issue in VBA_bDCM_lesion in the two simulateNLSS definitions (line 18 and 37, or more specifically 25 and 44):

[yp,~,~,~,er] = simulateNLSS(...
        out.dim.n_t,...
        out.options.f_fname,...
        out.options.g_fname,...
        posterior.muTheta,...
        posterior.muPhi,...
        out.u,...
        Inf, Inf,out.options,...
        posterior.muX0);

        rnormal.y=yp-er;

The second Inf causes an error as simulateNLSS expects a sigma for each Gaussian source. I am really unsure about the parameters that should be fed into the function, but something like posterior.a_sigma * Inf could be an easy fix.

Best, Simon