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 = "purple", 
     psize = 1.5, outfile = NULL)

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

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: numeric(optional): desired lower limit on x-axis
to: numeric(optional): desired upper limit on x-axis
step: numeric(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: numeric(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

References

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.

Examples

Run this code

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