Luminescence (version 0.1.8)

plot_DeDistribution: Plot De distribution with a kernel density estimate (KDE)

Description

Plot De values in ascending order and combine the plot with a kernel density function.

Usage

plot_DeDistribution(values,main=expression(paste(D[e], " Distribution")),
                    distribution.parameter="",mtext="", bw="nrd0",
                    xlab=expression(paste(D[e], " [Gy]")),
                    cex.global=1)

Arguments

values

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

main

character (with default): title of the plot

distribution.parameter

character (optional): plot additional distribution parameter, i.e. mean, median or sd. Example: distribution.parameter=c("median","mean","sd"). Mean and median are potted as a line, the sd is shown as

mtext

character (optional): further sample information (works as mtext).

bw

character (with default): bin-width, choose a numeric value for manual setting.

xlab

character (with default): label for x-axis unit. Possible values are "Gy", "s", "ka" or a user defined value.

cex.global

global scaling factor.

Value

• Returns a plot.

Details

For details on the calculation of the bin-width (parameter bw) see density.

References

Berger, G.W., 2010. An alternate form of probability-distribution plot for De values. Ancient TL 28, pp. 11-21.

Berger, G.W., 2011. Response to Galbraith. Ancient TL 29, pp. 48-50.

Galbraith, R.F., 2011. Some comments arising from Berger (2010). Ancient TL 29, pp. 41-47.

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

density, plot

Examples

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