Luminescence (version 0.2.1)

plot_Histogram: Plot a histogram with a separate error plot

Description

Plot a predefined histogram with an accompanying error plot as suggested by Rex Galbraith at 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 scaling 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.

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., 2012. Statistical aspects of equivalent dose and error calculation and display in OSL dating: An overview and some recommendations. Quaternary Geochronology, 11, pp.1-27.

See Also

hist, plot

Examples

##load data
data(ExampleData.DeValues)

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