age2ratio: Predict isotopic ratios from ages

Description Usage Arguments Value Examples

View source: R/age.R

Description

Groups a set of functions that take one (or more) ages (and their uncertainties) as input and produces the U–Pb, Th–Pb, Pb–Pb, Ar–Ar, K–Ca, Rb–Sr, Sm–Nd, Lu–Hf, Re–Os, concordia or Stacey-Kramers ratios as output.

Usage

1
age2ratio(tt, st = 0, ratio = "Pb206U238", exterr = FALSE, d = diseq(), J, sJ)

Arguments

tt

a scalar or (except when ratio = 'Wetherill', 'Tera-Wasserburg' or 'U-Th-Pb') vector of ages.

st

a scalar or (except when ratio = 'Wetherill', 'Tera-Wasserburg' or 'U-Th-Pb') vector with the standard error(s) of tt. Not used when ratio = 'Stacey-Kramers'.

ratio

one of 'Pb207U235', 'U238Pb206', 'Pb207Pb206', 'Pb208Th232', 'Wetherill', 'Tera-Wasserburg', 'U-Th-Pb', 'Ar40Ar39', 'Ca40K40', 'Hf176Lu176', 'Sr87Rb87', 'Os187Re187', 'Nd143Sm147' or 'Stacey-Kramers'.

exterr

logical. If TRUE, propagates decay constant uncertainties into st. Not used when ratio = 'Stacey-Kramers'.

d

an object of class diseq.

J

the J-factor of the Ar–Ar system (only used if ratio is 'Ar40Ar39').

sJ

the standard error of J (only used if ratio is 'Ar40Ar39').

Value

If ratio is 'Pb207U235', 'U238Pb206', 'Pb207Pb206', 'Pb208Th232', 'Ar40Ar39', 'Ca40K40', 'Hf176Lu176', 'Sr87Rb87', 'Os187Re187' or 'Nd143Sm147': either a two-element vector or a two-column matrix with the predicted isotopic ratio(s) and its/their standard error(s).

If ratio is 'Wetherill', 'Tera-Wasserburg' or 'U-Th-Pb': a two-element list containing

x: the concordia ratios

cov: the covariance matrix of the concordia ratios

If ratio is 'Stacey-Kramers': a three-column matrix with predicted ^{206}Pb/^{204}Pb, ^{207}Pb/^{204}Pb and ^{208}Pb/^{204}Pb ratios.

Examples

1
2
3
4
5
6
7
8
ratios <- c('Pb207U235','U238Pb206','Pb207Pb206','Pb208Th232',
            'Wetherill','Tera-Wasserburg','U-Th-Pb','Ar40Ar39',
            'Ca40K40','Hf176Lu176','Sr87Rb87','Os187Re187',
            'Nd143Sm147','Stacey-Kramers')
for (ratio in ratios){
     r <- age2ratio(tt=1000,st=10,ratio=ratio,J=1,sJ=0.1)
     print(r)
}

IsoplotR documentation built on July 10, 2021, 1:06 a.m.