dbED: Summarizing the distribution of equivalent dose values

Description

Calculating statistical parameters (skewness, kurtosis, quantiles) for a number of equivalent dose values.

Usage

dbED(EDdata, plot = TRUE, typ = c("pdf","hist","scatter"), 
     from = NULL, to = NULL, step = 0.1, nbin = 15, 
     pcolor = "grey", psize = 1.5)

Arguments

EDdata

matrix(required): two columns, i.e., equivalent dose values and their standard errors

plot

logical(with default): draw a plot or not

typ

character(with default): type of plot, "pdf" means a probability density plot, "hist" means a histogram, and "scatter" means a scatterplot, default typ="pdf"

from

numeric(optional): the desired lower limit on x-axis

numeric(optional): the desired upper limit on x-axis

step

numeric(with default): a step-size used for constructing the probability density plot. Smaller step value gives smoother density curve

nbin

numeric(with default): the desired number of intervals of a histogram

pcolor

character(with default): the color of a data point, input colors() to see more available colors

psize

numeric(with default): the size of a data point

Value

Return a list that contains following elements:
weight.EDthe weigthed mean of equivalent dose values and its standard error
skewnessthe weighted skewness of equivalent dose values and its standard error
kurtosisthe kurtosis of equivalent dose values and its standard error
quantile.EDthe quantiles of ED values

References

Bailey RM, Arnold LJ, 2006. Statistical modelling of single grain quartz De distributions and an assessment of procedures for estimating burial dose. Quaternary Science Reviews, 25: 2475-2502.

Galbraith RF, 2010. On plotting OSL equivalent doses. Ancient TL, 28: 1-10.

Galbraith RF, Roberts RG, 2012. Statistical aspects of equivalent dose and error calculation and display in OSL dating: An overview and some recommendations. Quaternary Geochronology, 11: 1-27.

Examples

Run this code

data(EDdata)
 dbED(EDdata$gl11,typ="pdf")