IsoplotR (version 0.8)

cad: Plot continuous data as cumulative age distributions

Description

Plot a dataset as a Cumulative Age Distribution (CAD), also known as a `empirical cumulative distribution function'.

Usage

cad(x, ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", colmap = "heat.colors", col = "black", ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", colmap = "heat.colors", ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", col = "black", type = 4, cutoff.76 = 1100, cutoff.disc = c(-15, 5), ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", col = "black", ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", col = "black", ...)
"cad"(x, pch = NA, verticals = TRUE, xlab = "age [Ma]", col = "black", ...)

Arguments

x

a numerical vector OR an object of class UPb, ArAr, UThHe, fissiontracks or detritals

...

optional arguments to the generic plot function

pch

plot character to mark the beginning of each CAD step

verticals

logical flag indicating if the horizontal lines of the CAD should be connected by vertical lines

xlab

x-axis label

colmap

an optional string with the name of one of R's built-in colour palettes (e.g., heat.colors, terrain.colors, topo.colors, cm.colors), which are to be used for plotting data of class detritals.

col

colour to give to single sample datasets (not applicable if x has class detritals)

type

scalar indicating whether to plot the $^{207}$Pb/$^{235}$U age (type=1), the $^{206}$Pb/$^{238}$U age (type=2), the $^{207}$Pb/$^{206}$Pb age (type=3), the $^{207}$Pb/$^{206}$Pb-$^{206}$Pb/$^{238}$U age (type=4), or the (Wetherill) concordia age (type=5)

cutoff.76

the age (in Ma) below which the $^{206}$Pb/$^{238}$U-age and above which the $^{207}$Pb/$^{206}$Pb-age is used. This parameter is only used if type=4.

cutoff.disc

two element vector with the maximum and minimum percentage discordance allowed between the $^{207}$Pb/$^{235}$U and $^{206}$Pb/$^{238}$U age (if $^{206}$Pb/$^{238}$U < cutoff.76) or between the $^{206}$Pb/$^{238}$U and $^{207}$Pb/$^{206}$Pb age (if $^{206}$Pb/$^{238}$U > cutoff.76). Set cutoff.disc=NA if you do not want to use this filter.

Examples

data(examples)
cad(examples$DZ,verticals=FALSE,pch=20)