calc_ratios: Removing outliers and summarising Pb isotope ratio
In andreabz/pbratios: Processing Pb isotope ratio data

Description Usage Arguments Value See Also Examples

The function removes replicates outside the median ± 2.5 mad range from raw Pb isotope ratios and calculates the average and standard error for each sample.

1	calc_ratios(data)

data

A data.table produced by the extract.data function, containing raw Pb isotope data. The dataset must have a column "sample" containing the sample name and the columns Pb208207, Pb206207, Pb208206, Pb207206 containing the measured Pb isotope ratio values for 208Pb/207Pb, 206Pb/207Pb, 208Pb/206Pb and 207Pb/206Pb, respectively.

A data.table with average Pb isotope ratio values for each sample. The number of replicates outside the median ± 2.5 mad range for each sample are printed on screen. Additionally, the standard errors for the measured Pb isotope ratios are stored in Pb20x20y.se columns.

extract_data

1
2
3

file.long <- system.file("extdata", "spmnya_2012.csv", package = "pbratios")
spmnya.2012 <- extract_data(file.long, report = "long")
spmnya.2012avg <- calc_ratios(spmnya.2012)