cosmo

Front and backend for Markov-Chain Monte Carlo inversion of cosmogenic nuclide concentrations
git clone git://src.adamsgaard.dk/cosmo
Log | Files | Refs | README | LICENSE

commit 3a499ebc15c7a4397c2069603a58f567ded159bf
parent b07a1882fe622f24f138058c07a0c9910ced2b4d
Author: Anders Damsgaard <anders.damsgaard@geo.au.dk>
Date:   Mon,  9 Nov 2015 14:52:48 +0100

put median value in xlabel

Diffstat:
Mmatlab/generate_plots.m | 20++++++++++++--------
1 file changed, 12 insertions(+), 8 deletions(-)

diff --git a/matlab/generate_plots.m b/matlab/generate_plots.m @@ -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}$ [' char(8240) ']'],... + xlabel('$\delta^{18}$O$_\mathrm{threshold}$ [$^o/_{oo}$]', ... 'Interpreter', 'LaTeX') text(0.02,0.98,'d','Units', ... 'Normalized', 'VerticalAlignment', 'Top') @@ -338,9 +338,9 @@ 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)]); + %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); @@ -353,19 +353,23 @@ for i1 = 1:M % for each model parameter hold off if i1 == 1 - xlabel('Interglacial erosion rate [mm/yr]') + xlabel(['Interglacial erosion rate [mm/yr], median = ' ... + num2str(med) ' mm/yr']) text(0.02,0.98,'a', 'Units', ... 'Normalized', 'VerticalAlignment', 'Top') elseif i1 == 2 - xlabel('Glacial erosion rate [mm/yr]') + xlabel(['Glacial erosion rate [mm/yr], median = ' ... + num2str(med) ' mm/yr']) text(0.02,0.98,'b', 'Units', ... 'Normalized', 'VerticalAlignment', 'Top') elseif i1 == 3 - xlabel('Timing of last deglaciation [yr]') + xlabel(['Timing of last deglaciation [yr], median = ' ... + num2str(med) ' yr']) text(0.02,0.98,'c', 'Units', ... 'Normalized', 'VerticalAlignment', 'Top') elseif i1 == 4 - xlabel(['$\delta^{18}$O$_\mathrm{threshold}$ [' char(8240) ']'],... + xlabel(['$\delta^{18}$O$_\mathrm{threshold}$ [$^o/_{oo}$]'... + ', median = ' num2str(med) ' $^o/_{oo}$'], ... 'Interpreter', 'LaTeX') text(0.02,0.98,'d','Units', ... 'Normalized', 'VerticalAlignment', 'Top')