IsoplotR (version 0.17)

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, ...)
# S3 method for default
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  colmap = "heat.colors", col = "black", ...)
# S3 method for detritals
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  colmap = "heat.colors", ...)
# S3 method for UPb
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", type = 4, cutoff.76 = 1100, cutoff.disc = c(-15, 5),
  common.Pb = 0, ...)
# S3 method for PbPb
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = FALSE, ...)
# S3 method for ArAr
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = FALSE, ...)
# S3 method for ThU
cad(x, pch = NA, verticals = TRUE, xlab = "age [ka]",
  col = "black", i2i = FALSE, ...)
# S3 method for ReOs
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = TRUE, ...)
# S3 method for SmNd
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = TRUE, ...)
# S3 method for RbSr
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = TRUE, ...)
# S3 method for LuHf
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", i2i = TRUE, ...)
# S3 method for UThHe
cad(x, pch = NA, verticals = TRUE, xlab = "age [Ma]",
  col = "black", ...)
# S3 method for fissiontracks
cad(x, pch = NA, verticals = TRUE,
  xlab = "age [Ma]", col = "black", ...)

Arguments

x

a numerical vector OR an object of class UPb, PbPb, ArAr, UThHe, fissiontracks, ReOs, RbSr, SmNd, LuHf, ThU 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.

common.Pb

apply a common lead correction using one of three methods:

1: use the isochron intercept as the initial Pb-composition

2: use the Stacey-Kramer two-stage model to infer the initial Pb-composition

3: use the Pb-composition stored in settings('iratio','Pb206Pb204') and settings('iratio','Pb207Pb204')

i2i

`isochron to intercept': calculates the initial (aka `inherited', `excess', or `common') \(^{40}\)Ar/\(^{36}\)Ar, \(^{207}\)Pb/\(^{204}\)Pb, \(^{87}\)Sr/\(^{86}\)Sr, \(^{143}\)Nd/\(^{144}\)Nd, \(^{187}\)Os/\(^{188}\)Os or \(^{176}\)Hf/\(^{177}\)Hf ratio from an isochron fit. Setting i2i to FALSE uses the default values stored in settings('iratio',...) or zero (for the Pb-Pb method). When applied to data of class ThU, setting i2i to TRUE applies a detrital Th-correction.

References

Vermeesch, P., 2007. Quantitative geomorphology of the White Mountains (California) using detrital apatite fission track thermochronology. Journal of Geophysical Research: Earth Surface, 112(F3).

Examples

# NOT RUN {
data(examples)
cad(examples$DZ,verticals=FALSE,pch=20)
# }