dbED: De distribution summarization

In numOSL: Numeric Routines for Optically Stimulated Luminescence Dating

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)

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

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

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; EDdata

Examples

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

numOSL documentation built on Sept. 18, 2023, 9:07 a.m.