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')