R Graphical Manual

Browse All

Last data update: 2014.03.03

R: define the interference corrections

interference

R Documentation

define the interference corrections

Description

create a new object of class logratios containing the interferences from neutron reactions on Ca and K

Usage

interference(intercepts, covmat, num, den, irr, label)

Arguments

`intercepts`	a vector with logratios
`covmat`	the covariance matrix of the logratios
`num`	a vector of strings marking the numerator isotopes of `intercepts`
`den`	a vector of strings marking the denominator isotopes of `intercepts`
`irr`	an object of class `irradiations`
`label`	a string with a name which can be used to identify the interference data in subsequent calculations

Value

an object of class logratios

Examples

samplefile <- system.file("Samples.csv",package="ArArRedux")
irrfile <- system.file("irradiations.csv",package="ArArRedux")
masses <- c("Ar37","Ar38","Ar39","Ar40","Ar36")
X <- read(samplefile,masses,blabel="EXB#",Jpos=c(3,15))
irr <- loadirradiations(irrfile)
# assume log(36Ar/37Ar) = log(39Ar/37Ar) = 1 in co-irradiate Ca-salt
# with variances of 0.0001 and zero covariances
ca <- interference(intercepts=c(1,1),
                   covmat=matrix(c(0.001,0,0,0.001),nrow=2),
                   num=c("Ar39","Ar36"),den=c("Ar37","Ar37"),
                   irr=X$irr[1],label="Ca-salt")
# assume log(39Ar/40Ar) = 4.637788 in co-irradiate K-glass
# with variance 7.9817e-4
k <- interference(intercepts=4.637788,covmat=7.9817e-4,
                  num="Ar39",den="Ar40",irr=X$irr[1],
                  label="K-glass")
ages <- process(X,irr,ca=ca,k=k)
summary(ages)

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/interference.Rd_%03d_medium.png", width=480, height=480)
> ### Name: interference
> ### Title: define the interference corrections
> ### Aliases: interference
> 
> ### ** Examples
> 
> samplefile <- system.file("Samples.csv",package="ArArRedux")
> irrfile <- system.file("irradiations.csv",package="ArArRedux")
> masses <- c("Ar37","Ar38","Ar39","Ar40","Ar36")
> X <- read(samplefile,masses,blabel="EXB#",Jpos=c(3,15))
> irr <- loadirradiations(irrfile)
> # assume log(36Ar/37Ar) = log(39Ar/37Ar) = 1 in co-irradiate Ca-salt
> # with variances of 0.0001 and zero covariances
> ca <- interference(intercepts=c(1,1),
+                    covmat=matrix(c(0.001,0,0,0.001),nrow=2),
+                    num=c("Ar39","Ar36"),den=c("Ar37","Ar37"),
+                    irr=X$irr[1],label="Ca-salt")
> # assume log(39Ar/40Ar) = 4.637788 in co-irradiate K-glass
> # with variance 7.9817e-4
> k <- interference(intercepts=4.637788,covmat=7.9817e-4,
+                   num="Ar39",den="Ar40",irr=X$irr[1],
+                   label="K-glass")
> ages <- process(X,irr,ca=ca,k=k)
> summary(ages)
      name      age[Ma]            s.e.[Ma]            
 [1,] "FC1_MD2" "28.1878378663059" "0.0243486389540383"
 [2,] "FC2-MD2" "28.2059707626935" "0.0238000488453725"
 [3,] "FC3-MD2" "28.2061617566445" "0.0238808633281276"
 [4,] "FC4-MD2" "28.2127749148053" "0.0241269932969587"
 [5,] "MD2-1a"  "101.273103914924" "0.1722365046112"   
 [6,] "MD2-1b"  "100.507680197963" "0.107838907602565" 
 [7,] "MD2-1c"  "99.8271292758592" "0.0901007895945283"
 [8,] "MD2-1d"  "99.7770601461817" "0.0889320487917834"
 [9,] "MD2-1e"  "99.803127628833"  "0.0874932779592161"
[10,] "MD2-1f"  "99.7593411680905" "0.0881604822499357"
[11,] "MD2-1g"  "99.7624512608441" "0.088402945869156" 
[12,] "MD2-1h"  "99.7657997135192" "0.0892932064066358"
[13,] "MD2-1i"  "99.7809850038155" "0.0885421365268459"
[14,] "MD2-1j"  "99.7963859579812" "0.0884435055661961"
[15,] "MD2-1k"  "99.8739009945625" "0.089214938016399" 
[16,] "MD2-1l"  "99.8376879696531" "0.0892360088523434"
[17,] "MD2-1m"  "99.8445421567958" "0.0909304504563129"
[18,] "MD2-1n"  "99.6802739993537" "0.0885902073258433"
[19,] "MD2-1o"  "99.4927099526701" "0.0878812462532874"
[20,] "MD2-1p"  "99.3816081556863" "0.088417374524879" 
[21,] "MD2-1r"  "99.6211030333348" "0.0944821681182069"
[22,] "MD2-1s"  "99.6385400875427" "0.0951221259029342"
[23,] "MD2-1t"  "99.7334321947907" "0.0910234167472191"
[24,] "MD2-2a"  "105.545383781371" "1.69685053245541"  
[25,] "MD2-2b"  "97.4931531999007" "0.266616333713798" 
[26,] "MD2-2c"  "97.6569716418552" "0.206099368489869" 
[27,] "MD2-2d"  "100.690549868299" "0.171232044344712" 
[28,] "MD2-2e"  "100.552754245508" "0.132448604867874" 
[29,] "MD2-2f"  "100.0546043186"   "0.0985526215893848"
[30,] "MD2-2g"  "100.036995868084" "0.0881731823043542"
[31,] "MD2-2h"  "99.7858398767599" "0.0883969007980255"
[32,] "MD2-2i"  "99.7818974442849" "0.0896261406771558"
[33,] "MD2-2j"  "99.7986673890951" "0.0877493336566643"
[34,] "MD2-2k"  "99.8329346997003" "0.0882969840945032"
[35,] "MD2-2l"  "99.8415070406617" "0.0874428921770652"
[36,] "MD2-2m"  "99.8331728076504" "0.0876848050174195"
[37,] "MD2-2n"  "99.8730792778322" "0.0880697014454322"
[38,] "MD2-2o"  "99.8757679521057" "0.08765956116161"  
[39,] "MD2-2p"  "99.8958516907824" "0.087483145489674" 
[40,] "MD2-2q"  "99.8980429765299" "0.0886526344390011"
[41,] "MD2-2r"  "99.8613246041806" "0.088877795981245" 
[42,] "MD2-2s"  "99.8877983286943" "0.0878257636232229"
[43,] "MD2-2t"  "99.8796672694624" "0.0880955189748692"
[44,] "MD2-2u"  "99.8061216673511" "0.0878529913923136"
[45,] "MD2-2v"  "99.7125959777839" "0.0876293051858243"
[46,] "MD2-2w"  "99.7417322780516" "0.0881991272630938"
[47,] "MD2-2x"  "99.6276526960413" "0.0876072633521482"
[48,] "MD2-2y"  "99.5352270990521" "0.0879587387153691"
[49,] "MD2-2z"  "99.5087956755033" "0.088644270405511" 
[50,] "MD2-2xz" "99.332575518587"  "0.0874330733602781"
[51,] "MD2-2yz" "99.3596055255094" "0.0895842547146868"
[52,] "MD2-1"   "101.641349150571" "0.1026265090852"   
[53,] "MD2-2"   "102.150589336283" "0.10103835424452"  
[54,] "MD2-3"   "100.802643459406" "0.112901997389179" 
[55,] "MD2-4"   "100.711816060455" "0.102831530360239" 
[56,] "MD2-5"   "98.4936827198769" "0.168934088907396" 
[57,] "MD2-6"   "100.355432895495" "0.0955440110700883"
[58,] "MD2-7"   "101.876328107968" "0.116688842247282" 
[59,] "MD2-8"   "101.690238866155" "0.0977111792987144"
[60,] "MD2-9"   "100.601277834395" "0.0961583533270614"
[61,] "MD2-10"  "100.212508831827" "0.0912318599415401"
      name      age[Ma]            s.e.[Ma]            
 [1,] "FC1_MD2" "28.1878378663059" "0.0243486389540383"
 [2,] "FC2-MD2" "28.2059707626935" "0.0238000488453725"
 [3,] "FC3-MD2" "28.2061617566445" "0.0238808633281276"
 [4,] "FC4-MD2" "28.2127749148053" "0.0241269932969587"
 [5,] "MD2-1a"  "101.273103914924" "0.1722365046112"   
 [6,] "MD2-1b"  "100.507680197963" "0.107838907602565" 
 [7,] "MD2-1c"  "99.8271292758592" "0.0901007895945283"
 [8,] "MD2-1d"  "99.7770601461817" "0.0889320487917834"
 [9,] "MD2-1e"  "99.803127628833"  "0.0874932779592161"
[10,] "MD2-1f"  "99.7593411680905" "0.0881604822499357"
[11,] "MD2-1g"  "99.7624512608441" "0.088402945869156" 
[12,] "MD2-1h"  "99.7657997135192" "0.0892932064066358"
[13,] "MD2-1i"  "99.7809850038155" "0.0885421365268459"
[14,] "MD2-1j"  "99.7963859579812" "0.0884435055661961"
[15,] "MD2-1k"  "99.8739009945625" "0.089214938016399" 
[16,] "MD2-1l"  "99.8376879696531" "0.0892360088523434"
[17,] "MD2-1m"  "99.8445421567958" "0.0909304504563129"
[18,] "MD2-1n"  "99.6802739993537" "0.0885902073258433"
[19,] "MD2-1o"  "99.4927099526701" "0.0878812462532874"
[20,] "MD2-1p"  "99.3816081556863" "0.088417374524879" 
[21,] "MD2-1r"  "99.6211030333348" "0.0944821681182069"
[22,] "MD2-1s"  "99.6385400875427" "0.0951221259029342"
[23,] "MD2-1t"  "99.7334321947907" "0.0910234167472191"
[24,] "MD2-2a"  "105.545383781371" "1.69685053245541"  
[25,] "MD2-2b"  "97.4931531999007" "0.266616333713798" 
[26,] "MD2-2c"  "97.6569716418552" "0.206099368489869" 
[27,] "MD2-2d"  "100.690549868299" "0.171232044344712" 
[28,] "MD2-2e"  "100.552754245508" "0.132448604867874" 
[29,] "MD2-2f"  "100.0546043186"   "0.0985526215893848"
[30,] "MD2-2g"  "100.036995868084" "0.0881731823043542"
[31,] "MD2-2h"  "99.7858398767599" "0.0883969007980255"
[32,] "MD2-2i"  "99.7818974442849" "0.0896261406771558"
[33,] "MD2-2j"  "99.7986673890951" "0.0877493336566643"
[34,] "MD2-2k"  "99.8329346997003" "0.0882969840945032"
[35,] "MD2-2l"  "99.8415070406617" "0.0874428921770652"
[36,] "MD2-2m"  "99.8331728076504" "0.0876848050174195"
[37,] "MD2-2n"  "99.8730792778322" "0.0880697014454322"
[38,] "MD2-2o"  "99.8757679521057" "0.08765956116161"  
[39,] "MD2-2p"  "99.8958516907824" "0.087483145489674" 
[40,] "MD2-2q"  "99.8980429765299" "0.0886526344390011"
[41,] "MD2-2r"  "99.8613246041806" "0.088877795981245" 
[42,] "MD2-2s"  "99.8877983286943" "0.0878257636232229"
[43,] "MD2-2t"  "99.8796672694624" "0.0880955189748692"
[44,] "MD2-2u"  "99.8061216673511" "0.0878529913923136"
[45,] "MD2-2v"  "99.7125959777839" "0.0876293051858243"
[46,] "MD2-2w"  "99.7417322780516" "0.0881991272630938"
[47,] "MD2-2x"  "99.6276526960413" "0.0876072633521482"
[48,] "MD2-2y"  "99.5352270990521" "0.0879587387153691"
[49,] "MD2-2z"  "99.5087956755033" "0.088644270405511" 
[50,] "MD2-2xz" "99.332575518587"  "0.0874330733602781"
[51,] "MD2-2yz" "99.3596055255094" "0.0895842547146868"
[52,] "MD2-1"   "101.641349150571" "0.1026265090852"   
[53,] "MD2-2"   "102.150589336283" "0.10103835424452"  
[54,] "MD2-3"   "100.802643459406" "0.112901997389179" 
[55,] "MD2-4"   "100.711816060455" "0.102831530360239" 
[56,] "MD2-5"   "98.4936827198769" "0.168934088907396" 
[57,] "MD2-6"   "100.355432895495" "0.0955440110700883"
[58,] "MD2-7"   "101.876328107968" "0.116688842247282" 
[59,] "MD2-8"   "101.690238866155" "0.0977111792987144"
[60,] "MD2-9"   "100.601277834395" "0.0961583533270614"
[61,] "MD2-10"  "100.212508831827" "0.0912318599415401"
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>