Learn R Programming

IsoplotR (version 0.5)

age: Calculate isotopic ages

Description

Calculates U-Pb ages and propagates their analytical uncertainties. Evaluates the equivalence of multiple ($^{206}$Pb/$^{238}$U-$^{207}$Pb/$^{235}$U or $^{207}$Pb/$^{206}$Pb-$^{206}$Pb/$^{238}$U) compositions, computes the weighted mean isotopic composition and the corresponding concordia age using the method of maximum likelihood, computes the mswd of equivalence and concordance and their respective Chi-squared p-values. Performs linear regression of U-Pb data on Wetherill and Tera-Wasserburg concordia diagrams. Computes the upper and lower intercept ages (for Wetherill) or the lower intercept age and the $^{207}$Pb/$^{206}$Pb intercept (for Tera-Wasserburg), taking into account error correlations and decay constant uncertainties.

Usage

age(x, ...)
"age"(x, method = "Pb206U238", dcu = TRUE, ...)
"age"(x, concordia = 1, wetherill = TRUE, dcu = TRUE, i = NA, ...)
"age"(x, ...)
"age"(x, isochron = FALSE, dcu = TRUE, i = NA, ...)

Arguments

x
a scalar containing an isotopic ratio, a two element vector containing an isotopic ratio and its standard error, or an object of class UPb or detritals.
...
optional arguments
method
one of either 'Pb206U238', 'Pb207U235', or 'Pb207Pb206'
dcu
propagate the decay constant uncertainties?
concordia
scalar flag indicating whether each U-Pb analysis should be considered separately (concordia=1), a concordia age should be calculated from all U-Pb analyses together (concordia=2), or a discordia line should be fit through all the U-Pb analyses (concordia=2).
wetherill
boolean flag to indicate whether the data should be evaluated in Wetherill (TRUE) or Tera-Wasserburg (FALSE) space. This option is only used when concordia=2
i
(optional) index of a particular aliquot
isochron
boolean flag indicating whether each Ar-Ar analysis should be considered separately (isochron=FALSE) or an isochron age should be calculated from all Ar-Ar analyses together (isochron=TRUE).

Value

if x is a scalar or a vector, returns the age using the geochronometer given by method and its standard error.if x has class UPb and concordia=1, returns a table with the following columns: `t.75', `err[t.75]', `t.68', `err[t.68]', `t.76',`err[t.76]', `t.conc', `err[t.conc]', containing the 207Pb/235U-age and standard error, the $^{206}$Pb/$^{238}$U-age and standard error, the $^{207}$Pb/$^{206}$ Pb-age and standard error, and the concordia age and standard error, respectively.if x has class UPb and concordia=2, returns a list with the following items:
x
a named vector with the (weighted mean) U-Pb composition
cov
the covariance matrix of the (mean) U-Pb composition
age
the concordia age (in Ma)
age.err
the standard error of the concordia age
mswd
a list with two items (equivalence and concordance) containing the MSWD (Mean of the Squared Weighted Deviates, a.k.a the reduced Chi-squared statistic outside of geochronology) of isotopic equivalence and age concordance, respectively.
p.value
a list with two items (equivalence and concordance) containing the p-value of the Chi-square test for isotopic equivalence and age concordance, respectively.
if x has class UPb and concordia=3, returns a list with the following items:
x
a two element vector with the upper and lower intercept ages (if wetherill==TRUE) or the lower intercept age and $^{207}$Pb/$^{206}$Pb intercept (for Tera-Wasserburg)
cov
the covariance matrix of the elements in x

Examples

Run this code
data(examples)
print(age(examples$UPb))
print(age(examples$UPb,concordia=1))
print(age(examples$UPb,concordia=2))

Run the code above in your browser using DataLab