https://github.com/Data2Dynamics/d2d
Tip revision: d72d4c223fb8adc0e3e44a1d3d1eafce18b995ae authored by Andreas Raue on 06 March 2015, 00:17:07 UTC
Increase C version code, after Joep's submits today to prevent compatibility problems
Increase C version code, after Joep's submits today to prevent compatibility problems
Tip revision: d72d4c2
arPlotScan.m
% plot scan of likelihood
function arPlotScan(jks, savetofile)
global ar
if(~exist('jks','var') || isempty(jks))
jks = find(~cellfun(@isempty, ar.scan.ps));
end
if(~exist('savetofile','var'))
savetofile = false;
end
sumples = 0;
for j=jks
if(~isempty(ar.scan.ps{j}))
sumples = sumples + 1;
end
end
[nrows, ncols] = arNtoColsAndRows(sumples);
h = myRaiseFigure('likelihood scan');
clf
set(h, 'Color', [1 1 1]);
count = 1;
minchi2 = Inf;
for jk=jks
if(~isempty(ar.scan.ps{jk}))
minchi2 = min([minchi2 min(ar.scan.chi2s{jk})]);
end
end
chi2curr = ar.chi2fit;
if(ar.config.fiterrors == 1)
chi2curr = 2*ar.ndata*log(sqrt(2*pi)) + chi2curr;
end
for jk=jks
if(~isempty(ar.scan.ps{jk}))
g = subplot(nrows,ncols,count);
ps = ar.scan.ps{jk};
chi2s = ar.scan.chi2s{jk};
constrs = ar.scan.constrs{jk};
if(ar.config.fiterrors == 1)
chi2s = 2*ar.ndata*log(sqrt(2*pi)) + chi2s;
end
if(ar.ndata>0)
plot(ps, chi2s, 'k-', 'LineWidth', 1)
hold on
end
if(ar.nconstr>0)
plot(ps, constrs, 'r--', 'LineWidth', 1)
end
hold off
ax1 = g;
% optimum
line(ar.p(jk), chi2curr, 'Marker', '*', 'Color', [.5 .5 .5], 'LineWidth', 1, 'MarkerSize', 8, ...
'Parent', ax1)
xlabel(ax1, ['log_{10}(' arNameTrafo(ar.pLabel{jk}) ')'])
arSpacedAxisLimits(g)
if(mod(count-1,ncols)==0)
if(ar.config.fiterrors == 1)
ylabel(ax1, '-2*log(L)');
else
ylabel(ax1, '\chi^2');
end
end
count = count + 1;
end
end
% save
if(savetofile && exist(pleGlobals.savePath, 'dir'))
pleGlobals.figPathMulti{jk} = [pleGlobals.savePath '/multi_plot'];
save([pleGlobals.savePath '/results.mat'], 'pleGlobals');
saveas(gcf, [pleGlobals.savePath '/multi_plot'], 'fig')
print('-depsc2', [pleGlobals.savePath '/multi_plot']);
end
function h = myRaiseFigure(figname)
global pleGlobals
openfigs = get(0,'Children');
figcolor = [1 1 1];
if(isfield(pleGlobals, 'fighandel_multi') && ~isempty(pleGlobals.fighandel_multi) && ...
pleGlobals.fighandel_multi ~= 0 && ...
sum(pleGlobals.fighandel_multi==openfigs)>0 && ...
strcmp(get(pleGlobals.fighandel_multi, 'Name'), figname))
h = pleGlobals.fighandel_multi;
figure(h);
else
h = figure('Name', figname, 'NumberTitle','off', ...
'Units', 'normalized', 'Position', ...
[0.1 0.1 0.6 0.8]);
set(h,'Color', figcolor);
pleGlobals.fighandel_multi = h;
end