tadd more parameters to output table - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
(HTM) git clone git://src.adamsgaard.dk/cosmo
(DIR) Log
(DIR) Files
(DIR) Refs
(DIR) README
(DIR) LICENSE
---
(DIR) commit f7d622afe2edbcfc984dc7257f309df1521c4f56
(DIR) parent f12d359f541763b09dafe405640865f8795edb02
(HTM) Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Thu, 19 Nov 2015 10:11:50 +0100
add more parameters to output table
Diffstat:
M matlab/generate_plots.m | 154 +++++++++++++++++++++++++++++--
1 file changed, 147 insertions(+), 7 deletions(-)
---
(DIR) diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -729,33 +729,173 @@ html = [html, ...
' <td> </td>\n'...
' <td align="center">25%%</td>\n'];
for i=1:Nwalkers
- html = [html, ' <td>',num2str(epsilon_int_25(i),3),'</td>\n'];
+ html = [html, ' <td>',...
+ num2str(epsilon_int_25(i),3),'</td>\n'];
end
-html = [html, ' <td>', num2str(sum(epsilon_int_25)/Nwalkers,3),'</td>\n'...
+html = [html, ' <td>',...
+ num2str(sum(epsilon_int_25)/Nwalkers,3),'</td>\n'...
' </tr>\n'...
' <tr>\n'...
' <td align="center">ε<sub>int</sub> [m/Myr]</td>\n'...
' <td align="center">50%%</td>\n'];
for i=1:Nwalkers
- html = [html, ' <td>',num2str(epsilon_int_50(i),3),'</td>\n'];
+ html = [html, ' <td>',...
+ num2str(epsilon_int_50(i),3),'</td>\n'];
end
-html = [html, ' <td>', num2str(sum(epsilon_int_50)/Nwalkers,3),'</td>\n'...
+html = [html, ' <td>',...
+ num2str(sum(epsilon_int_50)/Nwalkers,3),'</td>\n'...
' </tr>\n'...
' <tr style="border-bottom:1px solid #D0D0D0">\n'...
' <td> </td>\n'...
' <td align="center">75%%</td>\n'];
for i=1:Nwalkers
- html = [html, ' <td>',num2str(epsilon_int_75(i),3),'</td>\n'];
+ html = [html, ' <td>',...
+ num2str(epsilon_int_75(i),3),'</td>\n'];
end
-html = [html, ' <td>', num2str(sum(epsilon_int_75)/Nwalkers,3),'</td>\n'...
+html = [html, ' <td>',...
+ num2str(sum(epsilon_int_75)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'];
+
+
+% epsilon_gla
+html = [html, ...
+ ' <th data-field="avg">Average</th>\n'...
+ ' </tr>\n'...
+ ' <thead>\n'...
+ ' <tbody>\n'...
+ ' <tr>\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">25%%</td>\n'];
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(epsilon_gla_25(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(epsilon_gla_25)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td align="center">ε<sub>gla</sub> [m/Myr]</td>\n'...
+ ' <td align="center">50%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(epsilon_gla_50(i),3),'</td>\n'];
+end
+
+
+html = [html, ' <td>',...
+ num2str(sum(epsilon_gla_50)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">75%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(epsilon_gla_75(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(epsilon_gla_75)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'];
+
+
+% record_threshold
+html = [html, ...
+ ' <th data-field="avg">Average</th>\n'...
' </tr>\n'...
- ' </tbody>\n'...
+ ' <thead>\n'...
+ ' <tbody>\n'...
+ ' <tr>\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">25%%</td>\n'];
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(record_threshold_25(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(record_threshold_25)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td align="center">δ<sup>18</sup>O<sub>threshold</sub> [‰]</td>\n'...
+ ' <td align="center">50%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(record_threshold_50(i),3),'</td>\n'];
+end
+
+
+html = [html, ' <td>',...
+ num2str(sum(record_threshold_50)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">75%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(record_threshold_75(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(record_threshold_75)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'];
+
+% E
+html = [html, ...
+ ' <th data-field="avg">Average</th>\n'...
+ ' </tr>\n'...
+ ' <thead>\n'...
+ ' <tbody>\n'...
+ ' <tr>\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">25%%</td>\n'];
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(E_25(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(E_25)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr>\n'...
+ ' <td align="center">E [m]</td>\n'...
+ ' <td align="center">50%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(E_50(i),3),'</td>\n'];
+end
+
+
+html = [html, ' <td>',...
+ num2str(sum(E_50)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'...
+ ' <tr style="border-bottom:1px solid #D0D0D0">\n'...
+ ' <td> </td>\n'...
+ ' <td align="center">75%%</td>\n'];
+
+for i=1:Nwalkers
+ html = [html, ' <td>',...
+ num2str(E_75(i),3),'</td>\n'];
+end
+
+html = [html, ' <td>',...
+ num2str(sum(E_75)/Nwalkers,3),'</td>\n'...
+ ' </tr>\n'];
+
+
+% footer
+html = [html, ' </tbody>\n'...
'</table>\n'...
];
fileID = fopen(filename,'w');