data2york: Prepare geochronological data for York regression

Description Usage Arguments Value See Also Examples

View source: R/york.R

Description

Takes geochronology data as input and produces a five-column table with the variables, their uncertainties and error correlations as output. These can subsequently be used for York regression.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
data2york(x, ...)

## Default S3 method:
data2york(x, format = 1, ...)

## S3 method for class 'UPb'
data2york(x, option = 1, tt = 0, ...)

## S3 method for class 'ArAr'
data2york(x, inverse = TRUE, ...)

## S3 method for class 'KCa'
data2york(x, inverse = FALSE, ...)

## S3 method for class 'PbPb'
data2york(x, inverse = TRUE, ...)

## S3 method for class 'PD'
data2york(x, exterr = FALSE, inverse = FALSE, ...)

## S3 method for class 'UThHe'
data2york(x, ...)

## S3 method for class 'ThU'
data2york(x, type = 2, generic = TRUE, ...)

Arguments

x

a five or six column matrix OR an object of class UPb, PbPb, ArAr, ThU, UThHe, or PD (which includes objects of class RbSr, SmNd, LuHf and ReOs), generated by the read.data(...) function

...

optional arguments

format

one of

1 or 2: X, s[X], Y, s[Y], rho; where rho is the error correlation between X and Y; or

3: X/Z, s[X/Z], Y/Z, s[Y/Z], X/Y, s[X/Y]; for which the error correlations are automatically computed from the redundancy of the three ratios.

option

one of

1: Wetherill concordia ratios: X=07/35, sX=s[07/35], Y=06/38, sY=s[06/38], rXY.

2: Tera-Wasserburg ratios: X=38/06, sX=s[38/06], Y=07/06, sY=s[07/06], rho=rXY.

3: X=04/06, sX=s[04/06], Y=38/06, sY=s[38/06], rho=rXY (only valid if format=4,5 or 6).

4: X=04/07, sX=s[04/07], Y=05/07, sY=s[35/07], rho=rXY (only valid if format=4,5 or 6).

5: U-Th-Pb concordia ratios: X=08/32, sX=s[08/32], Y=06/38, sY=s[06/38], rho=rXY (only valid if format=7 or 8).

6: X=38/06, sX=s[38/06], Y=08/06, sY=s[08/06], rho=rXY (only valid if format=7 or 8).

7: X=38/06, sX=s[38/06], Y=08/06, sY=s[08/06], rho=rXY (only valid if format=7 or 8).

tt

the age of the sample. This is only used if x has class 7 or 8, in order to calculate the inherited {}^{208}Pb/{}^{232}Th ratio.

inverse

toggles between normal and inverse isochron ratios. data2york returns five columns X, s[X], Y, s[Y] and r[X,Y].

If inverse=FALSE, then X = {}^{206}Pb/{}^{204}Pb and Y = {}^{207}Pb/{}^{204}Pb (if x has class PbPb), or X = {}^{39}Ar/{}^{36}Ar and Y = {}^{40}Ar/{}^{36}Ar (if x has class ArAr), or X = {}^{40}K/{}^{44}Ca and Y = {}^{40}Ca/{}^{44}Ca (if x has class KCa), or X = {}^{87}Rb/{}^{86}Sr and Y = {}^{87}Sr/{}^{86}Sr (if x has class RbSr), or X = {}^{147}Sm/{}^{144}Nd and Y = {}^{143}Nd/{}^{144}Nd (if x has class SmNd), or X = {}^{187}Re/{}^{188}Os and Y = {}^{187}Os/{}^{188}Os (if x has class ReOs), or X = {}^{176}Lu/{}^{177}Hf and Y = {}^{176}Hf/{}^{177}Hf (if x has class LuHf).

If inverse=TRUE, then X = {}^{204}Pb/{}^{206}Pb and Y = {}^{207}Pb/{}^{206}Pb (if x has class PbPb), or X = {}^{39}Ar/{}^{40}Ar and Y = {}^{36}Ar/{}^{40}Ar (if x has class ArAr), or X = {}^{40}K/{}^{40}Ca and Y = {}^{44}Ca/{}^{40}Ca (if x has class KCa), or X = {}^{87}Rb/{}^{87}Sr and Y = {}^{86}Sr/{}^{87}Sr (if x has class RbSr), or X = {}^{147}Sm/{}^{143}Nd and Y = {}^{144}Nd/{}^{143}Nd (if x has class SmNd), or X = {}^{187}Re/{}^{187}Os and Y = {}^{188}Os/{}^{187}Os (if x has class ReOs), or X = {}^{176}Lu/{}^{176}Hf and Y = {}^{177}Hf/{}^{176}Hf (if x has class LuHf).

exterr

If TRUE, propagates the external uncertainties (e.g. decay constants) into the output errors.

type

Return ‘Rosholt’ or ‘Osmond’ ratios?

Rosholt (type=1) returns X=8/2, sX=s[8/2], Y=0/2, sY=s[0/2], rXY.

Osmond (type=2) returns X=2/8, sX=s[2/8], Y=0/8, sY=s[0/8], rXY.

generic

If TRUE, uses the following column headers: X, sX, Y, sY, rXY.

If FALSE and type=1, uses U238Th232, errU238Th232, Th230Th232, errTh230Th232, rho

or if FALSE and type=2, uses Th232U238, errTh232U238, Th230U238, errTh230U238, rho.

Value

a five-column table that can be used as input for york-regression.

See Also

york

Examples

1
2
3
4
f <- system.file("RbSr1.csv",package="IsoplotR")
dat <- read.csv(f)
yorkdat <- data2york(dat)
fit <- york(yorkdat)

pvermees/IsoplotR documentation built on Jan. 8, 2020, 7:03 p.m.