|
10 | 10 | colors = lines(3); |
11 | 11 |
|
12 | 12 | for im = 1:length(modelNrs) |
13 | | - |
| 13 | + figure('PaperUnits','centimeters','PaperSize',[21 29.7]) |
| 14 | + subplot(2,1,1) |
14 | 15 | modelNr = modelNrs{im}; |
15 | 16 | model_name = strsplit(modelNrs{im},'_'); |
16 | | - model_name = model_name{2}; |
17 | | - if strcmp(model_name,'Reelin') |
18 | | - model_name = 'Hass'; |
| 17 | + |
| 18 | + if strcmp(model_name{2},'Reelin') |
| 19 | + model_name{2} = 'Hass'; |
19 | 20 | end |
20 | | - if strcmp(model_name,'TGFb') |
21 | | - model_name = 'Lucarelli'; |
| 21 | + if strcmp(model_name{2},'TGFb') |
| 22 | + model_name{2} = 'Lucarelli'; |
22 | 23 | end |
23 | 24 |
|
24 | | - figure('Name',['Waterfall ' model_name],... |
25 | | - 'PaperUnits','centimeters','PaperSize',[21 29.7]) |
26 | | - subplot(2,1,1) |
27 | | - |
28 | 25 | min_model = min([ar_tmp.(modelNr).fmin_ip_log.chi2s,... |
29 | 26 | ar_tmp.(modelNr).fmin_trust_FindInputs_log.chi2s,ar_tmp.(modelNr).FindInputs_logFit.chi2s]); |
30 | 27 | s1=scatter(1:1000,sort(ar_tmp.(modelNr).FindInputs_logFit.chi2s)-min_model+1,... |
|
49 | 46 |
|
50 | 47 | set(gcf,'Color','w') |
51 | 48 | xlabel('fit number') |
52 | | - ylabel('final likelihood value') |
53 | | - title([model_name]) |
| 49 | + ylabel('final likelihood value + const.') |
| 50 | + title([model_name{2}]) |
54 | 51 | box off |
55 | 52 |
|
56 | 53 | set(gca,'FontSize',14) |
57 | 54 | axis manual |
58 | | - if ismember(model_name,{'Beer','Crauste'}) |
| 55 | + if ismember(model_name{2},{'Beer','Crauste'}) |
59 | 56 | legend(legend_string,'Location','southeast') |
60 | 57 | else |
61 | 58 | legend(legend_string,'Location','northwest') |
|
84 | 81 | end |
85 | 82 | set(gcf,'Color','w') |
86 | 83 | xlabel('optimization time [s]') |
87 | | - ylabel('final likelihood value') |
| 84 | + ylabel('final likelihood value + const.') |
88 | 85 | axis manual |
89 | 86 | x_limits = get(gca,'XLim'); |
90 | 87 | y_limits = get(gca,'YLim'); |
|
93 | 90 | text(10.^(log10(x_limits(1)) - (log10(x_limits(2)) - log10(x_limits(1)))*0.15),... |
94 | 91 | y_limits(2)+(y_limits(2)-y_limits(1))*0.2,'B','FontSize',16) |
95 | 92 |
|
96 | | -% if ~(exist('Supplement_figures')) |
97 | | -% mkdir('Supplement_figures') |
98 | | -% end |
99 | | -% print(gcf,['Supplement_figures/' model_name],'-dpdf','-fillpage') |
| 93 | + print(gcf,['Supplement_figures/' model_name{2}],'-dpdf','-fillpage') |
100 | 94 | end |
101 | | -%close all |
| 95 | +close all |
0 commit comments