tadd percentiles, add permil to d18O threshold label - cosmo - front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
HTML git clone git://src.adamsgaard.dk/cosmo
DIR Log
DIR Files
DIR Refs
DIR README
DIR LICENSE
---
DIR commit b07a1882fe622f24f138058c07a0c9910ced2b4d
DIR parent dab7ed7cd53163e444bc2f8d5310ac1b4a2b36d3
HTML Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date: Mon, 9 Nov 2015 14:39:59 +0100
add percentiles, add permil to d18O threshold label
Diffstat:
M matlab/generate_plots.m | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
---
DIR diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m
t@@ -229,7 +229,7 @@ for i1 = 1:M % for each model parameter
text(0.02,0.98,'c', 'Units', ...
'Normalized', 'VerticalAlignment', 'Top')
elseif i1 == 4
- xlabel('$\delta^{18}$O$_\mathrm{threshold}$', ...
+ xlabel(['$\delta^{18}$O$_\mathrm{threshold}$ [' char(8240) ']'],...
'Interpreter', 'LaTeX')
text(0.02,0.98,'d','Units', ...
'Normalized', 'VerticalAlignment', 'Top')
t@@ -338,13 +338,16 @@ for i1 = 1:M % for each model parameter
% 2nd quartile = median = 50th percentile
med = median(data);
plot([med, med], get(gca,'YLim'), 'm-')
+ ylims = get(gca,'YLim');
+ text(med, ylims(1) + (ylims(2) - ylims(1))*0.9, ...
+ [' \leftarrow ' num2str(med)]);
% 1st quartile = 25th percentile
prctile25 = prctile(data, 25);
plot([prctile25, prctile25], get(gca,'YLim'), 'm--')
% 3rd quartile = 75th percentile
- prctile75 = prctile(data, 75);
+ prctile75 = prctile(data, 75);
plot([prctile75, prctile75], get(gca,'YLim'), 'm--')
hold off
t@@ -362,7 +365,7 @@ for i1 = 1:M % for each model parameter
text(0.02,0.98,'c', 'Units', ...
'Normalized', 'VerticalAlignment', 'Top')
elseif i1 == 4
- xlabel('$\delta^{18}$O$_\mathrm{threshold}$', ...
+ xlabel(['$\delta^{18}$O$_\mathrm{threshold}$ [' char(8240) ']'],...
'Interpreter', 'LaTeX')
text(0.02,0.98,'d','Units', ...
'Normalized', 'VerticalAlignment', 'Top')