numOSL (version 1.8)

dbED: De distribution summarization

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): a two-column matrix, equivalent dose values and associated errors are stored in the first and second column, respectively

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

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

to

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

step

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

nbin

integer(with default): desired number of intervals of a histogram

pcolor

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

psize

double(with default): size of a data point

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")