From 1519d761ebb8c10350402b294907102d960b15a2 Mon Sep 17 00:00:00 2001
From: Johannes Ranke <johannes.ranke@jrwb.de>
Date: Fri, 28 Nov 2025 09:47:18 +0100
Subject: Add two-component error model to 2,4-D example

Also, explain the names of the metabolites.
---
 man/D24_2014.Rd | 19 ++++++++++++++++---
 1 file changed, 16 insertions(+), 3 deletions(-)

(limited to 'man')

diff --git a/man/D24_2014.Rd b/man/D24_2014.Rd
index 59ea5a9e..8a97dc47 100644
--- a/man/D24_2014.Rd
+++ b/man/D24_2014.Rd
@@ -25,6 +25,9 @@ context of pesticide  registrations, as the use of the data may be
 constrained by data protection regulations.
 }
 \details{
+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
@@ -48,14 +51,24 @@ m_D24_2 = mkinmod(D24 = mkinsub("DFOP", to = "DCP"),
 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 <- mmkin(
+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)
-plot(f_D24_2014[, 3])
+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]])
 }
 }
 \keyword{datasets}
-- 
cgit v1.2.3