diff options
Diffstat (limited to 'R/D24_2014.R')
| -rw-r--r-- | R/D24_2014.R | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/R/D24_2014.R b/R/D24_2014.R index aaed029c..8a077ff1 100644 --- a/R/D24_2014.R +++ b/R/D24_2014.R @@ -7,6 +7,9 @@ #' context of pesticide registrations, as the use of the data may be #' constrained by data protection regulations. #' +#' The acronyms used for the 2,4-D transformation products are DCP for +#' 2,4-dichlorophenol and DCA for 2,4-dichloroanisole. +#' #' Data for the first dataset are from p. 685. Data for the other four #' datasets were used in the preprocessed versions given in the kinetics #' section (p. 761ff.), with the exception of residues smaller than 1 for DCP @@ -33,5 +36,25 @@ #' DCP = mkinsub("SFO", to = "DCA"), #' DCA = mkinsub("SFO")) #' print(m_D24_2) +#' D24_2014_data <- lapply(D24_2014$ds, function(x) x$data) +#' names(D24_2014_data) <- sapply(D24_2014$ds, function(x) x$title) +#' f_D24_2014_const <- mmkin( +#' models = list( +#' "SFO-SFO-SFO" = m_D24, +#' "DFOP-SFO-SFO" = m_D24_2), +#' data = D24_2014_data, +#' quiet = TRUE) +#' print(f_D24_2014_const) +#' plot(f_D24_2014_const[, 3]) +#' f_D24_2014_tc <- update(f_D24_2014_const, error_model = "tc") +#' print(f_D24_2014_tc) +#' plot(f_D24_2014_tc[, 3]) +#' # For dataset 3, the best fit is obtained using constant +#' # variance and the DFOP-SFO-SFO model +#' AIC( +#' f_D24_2014_const[["SFO-SFO-SFO", 3]], +#' f_D24_2014_const[["DFOP-SFO-SFO", 3]], +#' f_D24_2014_tc[["SFO-SFO-SFO", 3]], +#' f_D24_2014_tc[["DFOP-SFO-SFO", 3]]) #' } "D24_2014" |