dbED: De distribution summarization

Description Usage Arguments Value References See Also Examples

View source: R/dbED.R

Description

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

Usage

1
2
3
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

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

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

numOSL documentation built on May 1, 2019, 8:45 p.m.