Commit e464d2e4 authored by Marco Ratto's avatar Marco Ratto
Browse files

make computation of log(det(sigma)) robust numerically

parent 6d0c4133
......@@ -69,7 +69,7 @@ parameter_names =;
save([M_.dname filesep 'Output' filesep M_.fname '_mean.mat'],'xparam1','hh','parameter_names','SIGMA');
fprintf('Estimation::marginal density: I''m computing the posterior log marginal density (modified harmonic mean)... ');
logdetSIGMA = log(det(SIGMA));
logdetSIGMA = sum(log(eig(SIGMA)));
invSIGMA = hh;
marginal = zeros(9,2);
linee = 0;
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment