kinfit: Fit kinetic models to chemical degradation data

Description

A selection of kinetic models as defined in the FOCUS kinetics report from 2006 are fitted to a given dataframe of chemical degradation data.

Usage

kinfit(kindata, kinmodels = c("SFO"), parent.0.user = NA, parent.0.fixed = FALSE, start.SFO = list(parent.0 = NA, k = NA),  start.FOMC = list(parent.0 = NA, alpha = NA, beta = NA),  start.DFOP = list(parent.0 = NA, k1 = NA, k2 = NA, g = NA), start.HS = list(parent.0 = NA, k1 = NA, k2 = NA, tb = NA), algorithm = "default")

Arguments

kindata

A data frame containing a time variable named t and concentration data for the parent compound named parent.

kinmodels

An array of character strings which are names of the models to be fit. Possible names are SFO, FOMC, DFOP and HS.

parent.0.user

The user can give a starting estimate for parent.0 here, overriding other potential sources for starting values as specified below.

parent.0.fixed

If TRUE, the initial value for the observed value will be fixed to the value specified by parent.0.user.

start.SFO

A list of starting parameters for fitting the SFO model, containing parent.0 and k.

start.FOMC

A list of starting parameters for fitting the FOMC model, containing parent.0, alpha and beta.

start.DFOP

A list of starting parameters for fitting the DFOP model, containing parent.0, k1, k2 and g.

start.HS

A list of starting parameters for fitting the HS model, containing parent.0, k1, k2 and tb.

algorithm

The algorithm to use for the calls to nls.

Value

A list of models of class nls representing the models that were fitted successfully.

Details

Per default all starting parameters are NA and the function tries to find suitable starting parameters on its own.

References

FOCUS (2006) “Guidance Document on Estimating Persistence and Degradation Kinetics from Environmental Fate Studies on Pesticides in EU Registration” Report of the FOCUS Work Group on Degradation Kinetics, EC Document Reference Sanco/10058/2005 version 2.0, 434 pp, http://focus.jrc.ec.europa.eu/dk

Examples

Run this code

data(FOCUS_2006_A)
(kinfits <- kinfit(FOCUS_2006_A))
data(FOCUS_2006_B)
(kinfits <- kinfit(FOCUS_2006_B, kinmodels=c("SFO","FOMC")))

Run the code above in your browser using DataLab