Luminescence (version 0.1.7)

plot_Histogram: Plot a histogram with a separate error plot.

Description

Plot a predefined histogram with an accompanying error plot suggested by Rex Galbraith on the UK LED in Oxford 2010.

Usage

plot_Histogram(values,
              main="Histogram",
              mtext="",
              xlab=expression(paste(D[e]," [Gy]")), 
              cex.global=1,
              breaks="Sturges",normal_curve=FALSE)

Arguments

values

data.frame (required): De (values[,1]) and De error (values[,2])

main

character (with default): title of the plot

mtext

character (optional): further sample information (mtext)

xlab

character (with default): label x-axis (xlab, see also plot)

cex.global

numeric (with default): global scale factor.

breaks

(with default): set breakpoints for histogram. Works as in hist.

normal_curve

logical (with default): add a normal curve to the histogram. Mean and sd are calculated from the input data. More see details section.

Value

• Returns a plot.

Details

If the normal curve is added, the y-axis in the histogram will show the probability density. Notice: The added curve is not a kernel density function!

References

Galbraith, R. (2010). Statistics in OSL: Some Current Questions; Ask Rex. Oral presentation during the UK TL/OSL/ESR Meeting at the School of Geography and the Environment, University of Oxford,8-10 September 2010.

Galbraith, R. F., Roberts, R. G. (in press), Statistical aspects of equivalent dose and error calculation and display in OSL dating: An overview and some recommendations, Quaternary Geochronology doi:10.1016/j.quageo.2012.04.020.

See Also

hist, plot

Examples

data(ExampleData.DeValues)
plot_Histogram(ExampleData.DeValues, xlab=expression(paste(D[e],"[s]")))