interference {ArArRedux}  R Documentation 
create a new object of class logratios
containing
the interferences from neutron reactions on Ca and K
interference(intercepts, covmat, num, den, irr, label)
intercepts 
a vector with logratios 
covmat 
the covariance matrix of the logratios 
num 
a vector of strings marking the numerator isotopes of

den 
a vector of strings marking the denominator isotopes of

irr 
an object of class 
label 
a string with a name which can be used to identify the interference data in subsequent calculations 
an object of class logratios
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 coirradiate Casalt # 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="Casalt") # assume log(39Ar/40Ar) = 4.637788 in coirradiate Kglass # with variance 7.9817e4 k < interference(intercepts=4.637788,covmat=7.9817e4, num="Ar39",den="Ar40",irr=X$irr[1], label="Kglass") ages < process(X,irr,ca=ca,k=k) summary(ages)