Compares the measured 40Ar/36Ar ratio of an air shot on a given
detector with the atmospheric ratio.
Usage
fractionation(fname, detector, MS = "ARGUS-VI", PH = FALSE)
Arguments
fname
a .csv file with the air shot data
detector
the name of the ion detector
MS
the type of mass spectrometer
PH
TRUE if the data were recorded in 'peak hopping' mode,
FALSE if recorded in multicollector mode.
Value
an object of class logratios
Examples
data(Melbourne)
fd37file <- system.file("AirL2.csv",package="ArArRedux")
fd40file <- system.file("AirH1.csv",package="ArArRedux")
fract <- list(fractionation(fd37file,"L2",PH=TRUE),
fractionation(fd40file,"H1",PH=FALSE))
if (isTRUE(all.equal(Melbourne$fract,fract))){
print("We just re-created the fractionation correction for the Melbourne dataset")
}
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(ArArRedux)
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/ArArRedux/fractionation.Rd_%03d_medium.png", width=480, height=480)
> ### Name: fractionation
> ### Title: Compute the mass fractionation correction
> ### Aliases: fractionation
>
> ### ** Examples
>
> data(Melbourne)
> fd37file <- system.file("AirL2.csv",package="ArArRedux")
> fd40file <- system.file("AirH1.csv",package="ArArRedux")
> fract <- list(fractionation(fd37file,"L2",PH=TRUE),
+ fractionation(fd40file,"H1",PH=FALSE))
> if (isTRUE(all.equal(Melbourne$fract,fract))){
+ print("We just re-created the fractionation correction for the Melbourne dataset")
+ }
>
>
>
>
>
> dev.off()
null device
1
>