numOSL (version 2.0)

dbED: De distribution summarization

Description

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

Usage

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

Arguments

EDdata

matrix(required): a two-column matrix (i.e., equivalent dose values and associated standard errors)

plot

logical(with default): draw a plot or not

typ

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

from

double(optional): desired lower limit on x-axis

to

double(optional): desired upper limit on x-axis

step

double(optional): a step-size used for constructing the probability density plot (if typ="pdf"). Smaller step value gives smoother density curve

nbin

integer(with default): desired number of intervals for the histogram (if typ="hist")

pcolor

character(with default): color of data points, input colors() to see available colors

psize

double(with default): size of data points

outfile

character(optional): if specified, calculated probability densities (if typ="pdf") will be written to a CSV file named "outfile" and saved to the current work directory

Value

Return a list that contains the following elements:

weight.ED

weigthed mean of equivalent dose values and associated standard error

skewness

weighted skewness of equivalent dose values and associated standard error

kurtosis

kurtosis of equivalent dose values and associated standard error

quantile.ED

quantiles of equivalent dose 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(19): 2475-2502.

Galbraith RF, 2010. On plotting OSL equivalent doses. Ancient TL, 28(1): 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.

See Also

psRadialPlot; RadialPlotter

Examples

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