Extrapolation to 'time zero'

Description

This function extrapolates time resolved mass spectrometer data to t=0. When fed with multicollector data, it forms the ratios of the raw signals, forms their logs and performs linear regression to t=0 When fed with single collector data, the function first takes their logs and extrapolates them to t=0 before taking ratios, unless denmass=NULL, in which case the logs of the raw signals are extrapolated.

Usage

fitlogratios(x, ...)

## Default S3 method:
fitlogratios(x, ...)

## S3 method for class 'timeresolved'
fitlogratios(x, denmass, ...)

## S3 method for class 'PHdata'
fitlogratios(x, denmass = NULL, ...)

Arguments

Value

an object of class logratios

Examples

samplefile <- system.file("Samples.csv",package="ArArRedux")
masses <- c("Ar37","Ar38","Ar39","Ar40","Ar36")
m <- loaddata(samplefile,masses) # samples and J-standards
blanklabel <- "EXB#"
l <- fitlogratios(blankcorr(m,blanklabel),"Ar40")
plotcorr(l)

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/fitlogratios.Rd_%03d_medium.png", width=480, height=480)
> ### Name: fitlogratios
> ### Title: Extrapolation to 'time zero'
> ### Aliases: fitlogratios fitlogratios.PHdata fitlogratios.default
> ###   fitlogratios.timeresolved
> 
> ### ** Examples
> 
> samplefile <- system.file("Samples.csv",package="ArArRedux")
> masses <- c("Ar37","Ar38","Ar39","Ar40","Ar36")
> m <- loaddata(samplefile,masses) # samples and J-standards
> blanklabel <- "EXB#"
> l <- fitlogratios(blankcorr(m,blanklabel),"Ar40")
> plotcorr(l)
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>