Learn R Programming

IsoplotR (version 0.11)

age: Calculate isotopic ages

Description

Calculates ages and propagates their analytical uncertainties.

Usage

age(x, ...)

# S3 method for default age(x, method = "U238-Pb206", exterr = TRUE, J = c(NA, NA), zeta = c(NA, NA), rhoD = c(NA, NA), ...)

# S3 method for UPb age(x, type = 1, wetherill = TRUE, exterr = TRUE, i = NA, sigdig = NA, ...)

# S3 method for ArAr age(x, isochron = FALSE, i2i = TRUE, exterr = TRUE, i = NA, sigdig = NA, ...)

# S3 method for UThHe age(x, central = FALSE, i = NA, sigdig = NA, ...)

# S3 method for fissiontracks age(x, central = FALSE, i = NA, sigdig = NA, exterr = TRUE, ...)

# S3 method for ReOs age(x, isochron = TRUE, i2i = TRUE, exterr = TRUE, i = NA, sigdig = NA, ...)

# S3 method for SmNd age(x, isochron = TRUE, i2i = TRUE, exterr = TRUE, i = NA, sigdig = NA, ...)

# S3 method for RbSr age(x, isochron = TRUE, i2i = TRUE, exterr = TRUE, i = NA, sigdig = NA, ...)

Arguments

x
can be:
  • a scalar containing an isotopic ratio,

  • a two element vector containing an isotopic ratio and its standard error, or the spontaneous and induced track densities Ns and Ni (if method='fissiontracks'),

  • a four element vector containing Ar40Ar39, s[Ar40Ar39], J, s[J],

  • a six element vector containing U, s[U], Th, s[Th], He and s[He],

  • an eight element vector containing U, s[U], Th, s[Th], He, s[He], Sm and s[Sm]

  • a six element vector containing Rb, s[Rb], Sr, s[Sr], Sr87Sr86, and s[Sr87Sr86]

  • a six element vector containing Re, s[Re], Os, s[Os], Os187Os188, and s[Os187Os188]

  • a six element vector containing Sm, s[Sm], Nd, s[Nd], Nd143Nd144, and s[Nd144Nd143]

OR

  • an object of class UPb, ArAr, RbSr, SmNd, ReOs UThHe or fissiontracks.

...
additional arguments
method
one of either 'U238-Pb206', 'U235-Pb207', 'Pb207-Pb206', 'Ar-Ar', 'Re-Os', 'Sm-Nd', 'Rb-Sr', 'U-Th-He' or 'fissiontracks'
exterr
propagate the external (decay constant and calibration factor) uncertainties?
J
two-element vector with the J-factor and its standard error.
zeta
two-element vector with the zeta-factor and its standard error.
rhoD
two-element vector with the track density of the dosimeter glass and its standard error.
type
scalar flag indicating whether each U-Pb analysis should be considered separately (type=1), a concordia age should be calculated from all U-Pb analyses together (type=2), or a discordia line should be fit through all the U-Pb analyses (type=3).
wetherill
logical flag to indicate whether the data should be evaluated in Wetherill (TRUE) or Tera-Wasserburg (FALSE) space. This option is only used when type=2
i
(optional) index of a particular aliquot
sigdig
number of significant digits for the uncertainty estimate (only used if type=1, isochron=FALSE or central=FALSE).
isochron
logical 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).
i2i
`isochron to intercept': calculates the initial (aka `inherited', `excess', or `common') \(^{40}\)Ar/\(^{36}\)Ar or \(^{187}\)Os/\(^{188}\)Os ratio from an isochron fit. Setting i2i to FALSE uses the default values stored in settings('iratio',...)
central
logical flag indicating whether each U-Th-He analysis should be considered separately (central=FALSE) or a central age should be calculated from all U-Th-He analyses together (central=TRUE).

Value

  1. if x is a scalar or a vector, returns the age using the geochronometer given by method and its standard error.
  2. if x has class UPb and type=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 \(^{207}\)Pb/\(^{235}\)U-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.
  3. if x has class UPb and type=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.
  4. if x has class UPb and type=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
  5. if x has class ArAr, RbSr, SmNd or ReOs and isochron=FALSE, returns a table of Ar-Ar, Rb-Sr, Sm-Nd or Re-Os ages and standard errors.
  6. if x has class ArAr, RbSr, SmNd or ReOs and isochron=TRUE, returns a list with the following items:
    a
    the intercept of the straight line fit and its standard error.
    b
    the slope of the fit and its standard error.
    y0
    the atmospheric \(^{40}\)Ar/\(^{36}\)Ar or initial \(^{187}\)Os/\(^{188}\)Os, \(^{87}\)Sr/\(^{86}\)Sr, or \(^{143}\)Nd/\(^{144}\)Nd ratio and its standard error.
    age
    the \(^{40}\)Ar/\(^{39}\)Ar, \(^{187}\)Os/\(^{187}\)Re, \(^{87}\)Sr/\(^{86}\)Sr, or \(^{143}\)Nd/\(^{144}\)Nd age and its standard error.
  7. if x has class UThHe and central=FALSE, returns a table of U-Th-He ages and standard errors.
  8. if x has class UThHe and central=TRUE, returns a list with the following items:
    uvw
    a three-element list with the weighted mean log[U/He], log[Th/He] and log[Sm/He] compositions.
    covmat
    a 3x3 covariance matrix for uvw
    mswd
    the reduced Chi-square value for the log[U/He]-log[Th/He] compositions.
    p.value
    the p-value of concordance between the log[U/He]-log[Th/He] compositions.
    age
    two-element vector with the central age and its standard error.
  9. if x has class fissiontracks and central=FALSE, returns a table of fission track ages and standard errors.
  10. if x has class fissiontracks and central=TRUE, returns a list with the following items:
    mswd
    the reduced Chi-square value for the fission track ages.
    p.value
    the p-value of concordance between the fission track ages.
    age
    a two-element vector with the central age and its standard error.
    disp
    the (over)dispersion of the single grain ages beyond the formal analytical uncertainties.

Examples

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

Run the code above in your browser using DataLab