IsoplotR (version 4.0)

age2ratio: Predict isotopic ratios from ages

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

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

# NOT RUN {
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)
}
# }