From c58ccd73951b2000a7a254fb36bbd9f0733db6cd Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Mon, 12 May 2025 22:16:10 +0200 Subject: Check and test locally --- docs/articles/prebuilt/2022_dmta_parent.html | 1039 +++++++++----------------- 1 file changed, 357 insertions(+), 682 deletions(-) (limited to 'docs/articles/prebuilt/2022_dmta_parent.html') diff --git a/docs/articles/prebuilt/2022_dmta_parent.html b/docs/articles/prebuilt/2022_dmta_parent.html index 14f12c0e..10035c99 100644 --- a/docs/articles/prebuilt/2022_dmta_parent.html +++ b/docs/articles/prebuilt/2022_dmta_parent.html @@ -1,396 +1,100 @@ - - - + - - - - - - - - - -Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - + + + + +Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P • mkin + + + + + - + + + +

Testing hierarchical parent degradation -kinetics with residue data on dimethenamid and dimethenamid-P

Johannes Ranke

Last change on 5 January 2023, last compiled on 13 -Februar 2025

Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P

Johannes +Ranke

+ +

Last change on 5 January +2023, last compiled on 12 Mai 2025

+ + Source: vignettes/prebuilt/2022_dmta_parent.rmd +

2022_dmta_parent.rmd

- -

-true -

- -

Introduction

+ + +

Introduction +

The purpose of this document is to demonstrate how nonlinear hierarchical models (NLHM) based on the parent degradation models SFO, FOMC, DFOP and HS can be fitted with the mkin package.

@@ -398,7 +102,7 @@ FOMC, DFOP and HS can be fitted with the mkin package.

173340 (Application of nonlinear hierarchical models to the kinetic evaluation of chemical degradation data) of the German Environment Agency carried out in 2022 and 2023.

The mkin package is used in version 1.2.9. It contains the test data +

The mkin package is used in version 1.2.10. It contains the test data and the functions used in the evaluations. The saemix package is used as a backend for fitting the NLHM, but is also loaded to make the convergence plot function available.

@@ -406,19 +110,21 @@ make the convergence plot function available.

also provides the kable function that is used to improve the display of tabular data in R markdown documents. For parallel processing, the parallel package is used.

library(mkin)
-library(knitr)
-library(saemix)
-library(parallel)
-n_cores <- detectCores()
-if (Sys.info()["sysname"] == "Windows") {
-  cl <- makePSOCKcluster(n_cores)
-} else {
-  cl <- makeForkCluster(n_cores)
-}

+library(mkin)
+library(knitr)
+library(saemix)
+library(parallel)
+n_cores <- detectCores()
+if (Sys.info()["sysname"] == "Windows") {
+  cl <- makePSOCKcluster(n_cores)
+} else {
+  cl <- makeForkCluster(n_cores)
+}

Data

Data +

The test data are available in the mkin package as an object of class mkindsg (mkin dataset group) under the identifier dimethenamid_2018. The following preprocessing steps are @@ -440,37 +146,37 @@ and Elliot 2) are combined, resulting in dimethenamid (DMTA) data from six soils.

The following commented R code performs this preprocessing.

# Apply a function to each of the seven datasets in the mkindsg object to create a list
-dmta_ds <- lapply(1:7, function(i) {
-  ds_i <- dimethenamid_2018$ds[[i]]$data                     # Get a dataset
-  ds_i[ds_i$name == "DMTAP", "name"] <-  "DMTA"              # Rename DMTAP to DMTA
-  ds_i <- subset(ds_i, name == "DMTA", c("name", "time", "value")) # Select data
-  ds_i$time <- ds_i$time * dimethenamid_2018$f_time_norm[i]  # Normalise time
-  ds_i                                                       # Return the dataset
-})
-
-# Use dataset titles as names for the list elements
-names(dmta_ds) <- sapply(dimethenamid_2018$ds, function(ds) ds$title)
-
-# Combine data for Elliot soil to obtain a named list with six elements
-dmta_ds[["Elliot"]] <- rbind(dmta_ds[["Elliot 1"]], dmta_ds[["Elliot 2"]]) #
-dmta_ds[["Elliot 1"]] <- NULL
-dmta_ds[["Elliot 2"]] <- NULL

+# Apply a function to each of the seven datasets in the mkindsg object to create a list
+dmta_ds <- lapply(1:7, function(i) {
+  ds_i <- dimethenamid_2018$ds[[i]]$data                     # Get a dataset
+  ds_i[ds_i$name == "DMTAP", "name"] <-  "DMTA"              # Rename DMTAP to DMTA
+  ds_i <- subset(ds_i, name == "DMTA", c("name", "time", "value")) # Select data
+  ds_i$time <- ds_i$time * dimethenamid_2018$f_time_norm[i]  # Normalise time
+  ds_i                                                       # Return the dataset
+})
+
+# Use dataset titles as names for the list elements
+names(dmta_ds) <- sapply(dimethenamid_2018$ds, function(ds) ds$title)
+
+# Combine data for Elliot soil to obtain a named list with six elements
+dmta_ds[["Elliot"]] <- rbind(dmta_ds[["Elliot 1"]], dmta_ds[["Elliot 2"]]) #
+dmta_ds[["Elliot 1"]] <- NULL
+dmta_ds[["Elliot 2"]] <- NULL

The following tables show the 6 datasets.

for (ds_name in names(dmta_ds)) {
-    print(kable(mkin_long_to_wide(dmta_ds[[ds_name]]),
-      caption = paste("Dataset", ds_name),
-      label = paste0("tab:", ds_name), booktabs = TRUE))
-    cat("\n\\clearpage\n")
-}

- +

+for (ds_name in names(dmta_ds)) {
+    print(kable(mkin_long_to_wide(dmta_ds[[ds_name]]),
+      caption = paste("Dataset", ds_name),
+      label = paste0("tab:", ds_name), booktabs = TRUE))
+    cat("\n\\clearpage\n")
+}

- - + - - + @@ -502,14 +208,12 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset Calke
time	DMTA
0

- +

- - + - - + @@ -577,14 +281,12 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset Borstel
time	DMTA
0.000000

- +

- - + - - + @@ -732,14 +434,12 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset Flaach
time	DMTA
0.0000000

- +

- - + - - + @@ -839,14 +539,12 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset BBA 2.2
time	DMTA
0.0000000

- +

- - + - - + @@ -946,14 +644,12 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset BBA 2.3
time	DMTA
0.0000000

- +

- - + - - + @@ -1150,24 +846,25 @@ dmta_ds[["Elliot 2"]] <- NULL

Dataset Elliot
time	DMTA
0.000000

Separate evaluations

Separate evaluations +

In order to obtain suitable starting parameters for the NLHM fits, separate fits of the four models to the data for each soil are generated using the mmkin function from the mkin package. In a first step, constant variance is assumed. Convergence is checked with the status function.

deg_mods <- c("SFO", "FOMC", "DFOP", "HS")
-f_sep_const <- mmkin(
-  deg_mods,
-  dmta_ds,
-  error_model = "const",
-  quiet = TRUE)
-
-status(f_sep_const) |> kable()

- - - +

+deg_mods <- c("SFO", "FOMC", "DFOP", "HS")
+f_sep_const <- mmkin(
+  deg_mods,
+  dmta_ds,
+  error_model = "const",
+  quiet = TRUE)
+
+status(f_sep_const) |> kable()

+ @@ -1175,8 +872,7 @@ status(f_sep_const) |> kable() - - + @@ -1221,11 +917,11 @@ to converge. All separate fits with constant variance converged, with the sole exception of the HS fit to the BBA 2.2 data. To prepare for fitting NLHM using the two-component error model, the separate fits are updated assuming two-component error.

f_sep_tc <- update(f_sep_const, error_model = "tc")
-status(f_sep_tc) |> kable()

	Calke	Borstel	BBA 2.2	BBA 2.3	Elliot
SFO

- - +

+f_sep_tc <- update(f_sep_const, error_model = "tc")
+status(f_sep_tc) |> kable()

+ @@ -1233,8 +929,7 @@ status(f_sep_tc) |> kable() - - + @@ -1278,8 +973,9 @@ status(f_sep_tc) |> kable() converge with constant variance did converge, but other non-SFO fits failed to converge.

Hierarchichal model fits

Hierarchichal model fits +

The following code fits eight versions of hierarchical models to the data, using SFO, FOMC, DFOP and HS for the parent compound, and using either constant variance or two-component error for the error model. The @@ -1291,18 +987,18 @@ all eight versions in parallel (given a sufficient number of computing cores being available) to save execution time.

Convergence plots and summaries for these fits are shown in the appendix.

f_saem <- mhmkin(list(f_sep_const, f_sep_tc), transformations = "saemix")

+f_saem <- mhmkin(list(f_sep_const, f_sep_tc), transformations = "saemix")

The output of the status function shows that all fits terminated successfully.

status(f_saem) |> kable()

	Calke	Borstel	BBA 2.2	BBA 2.3	Elliot
SFO

- - +

+status(f_saem) |> kable()

+ - - + @@ -1328,17 +1024,16 @@ terminated successfully.

	const	tc
SFO

The AIC and BIC values show that the biphasic models DFOP and HS give the best fits.

anova(f_saem) |> kable(digits = 1)

- - - +

+anova(f_saem) |> kable(digits = 1)

+ - - + @@ -1404,22 +1099,22 @@ it shows the lowest AIC and BIC values in the first set of fits when combined with the two-component error model. Therefore, the DFOP model was selected for further refinements of the fits with the aim to make the model fully identifiable.

Parameter identifiability based on the Fisher Information -Matrix

Parameter identifiability based on the Fisher Information +Matrix +

Using the illparms function, ill-defined statistical model parameters such as standard deviations of the degradation parameters in the population and error model parameters can be found.

illparms(f_saem) |> kable()

	npar	AIC	BIC	Lik
SFO const

- - +

+illparms(f_saem) |> kable()

+ - - + @@ -1450,32 +1145,34 @@ be obtained for this standard deviation when using different starting values.

The thus identified overparameterisation is addressed by removing the random effect for k2 from the parameter model.

f_saem_dfop_tc_no_ranef_k2 <- update(f_saem[["DFOP", "tc"]],
-  no_random_effect = "k2")

+f_saem_dfop_tc_no_ranef_k2 <- update(f_saem[["DFOP", "tc"]],
+  no_random_effect = "k2")

For the resulting fit, it is checked whether there are still ill-defined parameters,

illparms(f_saem_dfop_tc_no_ranef_k2)

+illparms(f_saem_dfop_tc_no_ranef_k2)

which is not the case. Below, the refined model is compared with the previous best model. The model without random effect for k2 is a reduced version of the previous model. Therefore, the models are nested and can be compared using the likelihood ratio test. This is achieved with the argument test = TRUE to the anova function.

anova(f_saem[["DFOP", "tc"]], f_saem_dfop_tc_no_ranef_k2, test = TRUE) |>
-  kable(format.args = list(digits = 4))

	const	tc
SFO

+anova(f_saem[["DFOP", "tc"]], f_saem_dfop_tc_no_ranef_k2, test = TRUE) |>
+  kable(format.args = list(digits = 4))

--------++++++++ - - + @@ -1484,8 +1181,7 @@ achieved with the argument test = TRUE to the - - + @@ -1513,35 +1209,35 @@ achieved with the argument test = TRUE to the random effect for k2. The p value of the likelihood ratio test is much greater than 0.05, indicating that the model with the higher likelihood (here the model with random effects for all -degradation parameters f_saem[["DFOP", "tc"]]) does not fit +degradation parameters f_saem[["DFOP", "tc"]]) does not fit significantly better than the model with the lower likelihood (the reduced model f_saem_dfop_tc_no_ranef_k2).

Therefore, AIC, BIC and likelihood ratio test suggest the use of the reduced model.

The convergence of the fit is checked visually.

+ Convergence plot for the NLHM DFOP fit with two-component error and without a random effect on 'k2'

Convergence plot for the NLHM DFOP fit with two-component error and without a random effect on ‘k2’

All parameters appear to have converged to a satisfactory degree. The final fit is plotted using the plot method from the mkin package.

plot(f_saem_dfop_tc_no_ranef_k2)

+plot(f_saem_dfop_tc_no_ranef_k2)

+ Plot of the final NLHM DFOP fit

Plot of the final NLHM DFOP fit

Finally, a summary report of the fit is produced.

summary(f_saem_dfop_tc_no_ranef_k2)

+summary(f_saem_dfop_tc_no_ranef_k2)

saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:33 2025 
-Date of summary: Thu Feb 13 16:33:34 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:46 2025 
+Date of summary: Mon May 12 21:06:46 2025 
 
 Equations:
 d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -1553,7 +1249,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 3.778 s
+Fitted in 3.953 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Two-component variance function 
@@ -1615,8 +1311,9 @@ Estimated disappearance times:
       DT50  DT90 DT50back DT50_k1 DT50_k2
 DMTA 11.45 41.32    12.44   10.82   81.85

Alternative check of parameter identifiability

Alternative check of parameter identifiability +

The parameter check used in the illparms function is based on a quadratic approximation of the likelihood surface near its optimum, which is calculated using the Fisher Information Matrix (FIM). @@ -1626,12 +1323,13 @@ approach has recently been implemented in mkin.

of the saem algorithm with different parameter combinations, sampled from the range of the parameters obtained for the individual datasets fitted separately using nonlinear regression.

f_saem_dfop_tc_multi <- multistart(f_saem[["DFOP", "tc"]], n = 50, cores = 15)

par(mar = c(6.1, 4.1, 2.1, 2.1))
-parplot(f_saem_dfop_tc_multi, lpos = "bottomright", ylim = c(0.3, 10), las = 2)

+f_saem_dfop_tc_multi <- multistart(f_saem[["DFOP", "tc"]], n = 50, cores = 15)

+par(mar = c(6.1, 4.1, 2.1, 2.1))
+parplot(f_saem_dfop_tc_multi, lpos = "bottomright", ylim = c(0.3, 10), las = 2)

+ Scaled parameters from the multistart runs, full model

Scaled parameters from the multistart runs, full model

@@ -1642,34 +1340,36 @@ parameter g.

The parameter boxplots of the multistart runs with the reduced model shown below indicate that all runs give similar results, regardless of the starting parameters.

f_saem_dfop_tc_no_ranef_k2_multi <- multistart(f_saem_dfop_tc_no_ranef_k2,
-  n = 50, cores = 15)

par(mar = c(6.1, 4.1, 2.1, 2.1))
-parplot(f_saem_dfop_tc_no_ranef_k2_multi, ylim = c(0.5, 2), las = 2,
-  lpos = "bottomright")

+f_saem_dfop_tc_no_ranef_k2_multi <- multistart(f_saem_dfop_tc_no_ranef_k2,
+  n = 50, cores = 15)

+par(mar = c(6.1, 4.1, 2.1, 2.1))
+parplot(f_saem_dfop_tc_no_ranef_k2_multi, ylim = c(0.5, 2), las = 2,
+  lpos = "bottomright")

+ Scaled parameters from the multistart runs, reduced model

Scaled parameters from the multistart runs, reduced model

When only the parameters of the top 25% of the fits are shown (based on a feature introduced in mkin 1.2.2 currently under development), the scatter is even less as shown below.

par(mar = c(6.1, 4.1, 2.1, 2.1))
-parplot(f_saem_dfop_tc_no_ranef_k2_multi, ylim = c(0.5, 2), las = 2, llquant = 0.25,
-  lpos = "bottomright")

+par(mar = c(6.1, 4.1, 2.1, 2.1))
+parplot(f_saem_dfop_tc_no_ranef_k2_multi, ylim = c(0.5, 2), las = 2, llquant = 0.25,
+  lpos = "bottomright")

- $Scaled parameters from the multistart runs, reduced model, fits with the top 25\% likelihood values$ -

+ $Scaled parameters from the multistart runs, reduced model, fits with the top 25\% likelihood values$

Scaled parameters from the multistart runs, reduced model, fits with the top 25% likelihood values

Conclusions

Conclusions +

Fitting the four parent degradation models SFO, FOMC, DFOP and HS as part of hierarchical model fits with two different error models and normal distributions of the transformed degradation parameters works @@ -1681,35 +1381,39 @@ was fully identifiable and showed the lowest values for the model selection criteria AIC and BIC. The reliability of the identification of all model parameters was confirmed using multiple starting values.

Acknowledgements

Acknowledgements +

The helpful comments by Janina Wöltjen of the German Environment Agency are gratefully acknowledged.

References

References +

Duchesne, Ronan, Anissa Guillemin, Olivier Gandrillon, and Fabien Crauste. 2021. “Practical Identifiability in the Frame of Nonlinear Mixed Effects Models: The Example of the in Vitro -Erythropoiesis.” BMC Bioinformatics 22 (478). https://doi.org/10.1186/s12859-021-04373-4. +Erythropoiesis.” BMC Bioinformatics 22 (478). https://doi.org/10.1186/s12859-021-04373-4.

Appendix

Hierarchical model fit listings

Appendix +

Hierarchical model fit listings +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:26 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:38 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - k_DMTA * DMTA
@@ -1719,7 +1423,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 0.792 s
+Fitted in 0.793 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Constant variance 
@@ -1771,17 +1475,17 @@ Estimated disappearance times:
       DT50  DT90
 DMTA 12.24 40.65
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:28 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:40 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - k_DMTA * DMTA
@@ -1791,7 +1495,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 2.245 s
+Fitted in 2.462 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Two-component variance function 
@@ -1845,17 +1549,17 @@ Estimated disappearance times:
       DT50  DT90
 DMTA 12.24 40.65
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:27 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:39 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - (alpha/beta) * 1/((time/beta) + 1) * DMTA
@@ -1865,7 +1569,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 1.409 s
+Fitted in 1.408 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Constant variance 
@@ -1922,17 +1626,17 @@ Estimated disappearance times:
       DT50  DT90 DT50back
 DMTA 11.41 42.53     12.8
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:28 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:40 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - (alpha/beta) * 1/((time/beta) + 1) * DMTA
@@ -1942,7 +1646,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 2.811 s
+Fitted in 2.735 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Two-component variance function 
@@ -2001,17 +1705,17 @@ Estimated disappearance times:
       DT50  DT90 DT50back
 DMTA 11.05 42.81    12.89
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:27 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:39 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -2023,7 +1727,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 1.638 s
+Fitted in 1.769 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Constant variance 
@@ -2085,17 +1789,17 @@ Estimated disappearance times:
       DT50 DT90 DT50back DT50_k1 DT50_k2
 DMTA 11.79 42.8    12.88   11.09   76.46
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:28 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:41 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -2107,7 +1811,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 3.024 s
+Fitted in 3.069 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Two-component variance function 
@@ -2171,17 +1875,17 @@ Estimated disappearance times:
       DT50  DT90 DT50back DT50_k1 DT50_k2
 DMTA 11.39 41.36    12.45   10.74   83.48
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:28 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:39 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - ifelse(time <= tb, k1, k2) * DMTA
@@ -2191,7 +1895,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 2.301 s
+Fitted in 1.996 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Constant variance 
@@ -2253,17 +1957,17 @@ Estimated disappearance times:
       DT50  DT90 DT50back DT50_k1 DT50_k2
 DMTA 8.801 67.91    20.44   8.801   31.13
 
-

+ +


 saemix version used for fitting:      3.3 
-mkin version used for pre-fitting:  1.2.9 
-R version used for fitting:         4.4.2 
-Date of fit:     Thu Feb 13 16:33:29 2025 
-Date of summary: Thu Feb 13 16:34:39 2025 
+mkin version used for pre-fitting:  1.2.10 
+R version used for fitting:         4.5.0 
+Date of fit:     Mon May 12 21:06:41 2025 
+Date of summary: Mon May 12 21:07:51 2025 
 
 Equations:
 d_DMTA/dt = - ifelse(time <= tb, k1, k2) * DMTA
@@ -2273,7 +1977,7 @@ Data:
 
 Model predictions using solution type analytical 
 
-Fitted in 3.264 s
+Fitted in 3.389 s
 Using 300, 100 iterations and 9 chains
 
 Variance model: Two-component variance function 
@@ -2337,69 +2041,63 @@ Estimated disappearance times:
       DT50  DT90 DT50back DT50_k1 DT50_k2
 DMTA 11.84 51.71    15.57   11.84   29.66
 
-

+ +

Hierarchical model convergence plots

Hierarchical model convergence plots +

+ Convergence plot for the NLHM SFO fit with constant variance

Convergence plot for the NLHM SFO fit with constant variance

+ Convergence plot for the NLHM SFO fit with two-component error

Convergence plot for the NLHM SFO fit with two-component error

+ Convergence plot for the NLHM FOMC fit with constant variance

Convergence plot for the NLHM FOMC fit with constant variance

+ Convergence plot for the NLHM FOMC fit with two-component error

Convergence plot for the NLHM FOMC fit with two-component error

+ Convergence plot for the NLHM DFOP fit with constant variance

Convergence plot for the NLHM DFOP fit with constant variance

+ Convergence plot for the NLHM DFOP fit with two-component error

Convergence plot for the NLHM DFOP fit with two-component error

+ Convergence plot for the NLHM HS fit with constant variance

Convergence plot for the NLHM HS fit with constant variance

+ Convergence plot for the NLHM HS fit with two-component error

Convergence plot for the NLHM HS fit with two-component error

Session info

R version 4.4.2 (2024-10-31)
+
+Session info
+
+R version 4.5.0 (2025-04-11)
 Platform: x86_64-pc-linux-gnu
 Running under: Debian GNU/Linux 12 (bookworm)
 
 Matrix products: default
 BLAS:   /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.11.0 
-LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.11.0
+LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.11.0  LAPACK version 3.11.0
 
 locale:
  [1] LC_CTYPE=de_DE.UTF-8       LC_NUMERIC=C              
@@ -2417,73 +2115,50 @@ attached base packages:
 [8] base     
 
 other attached packages:
-[1] saemix_3.3      npde_3.5        knitr_1.49      mkin_1.2.9     
-[5] rmarkdown_2.29  nvimcom_0.9-167
+[1] rmarkdown_2.29  nvimcom_0.9-167 saemix_3.3      npde_3.5       
+[5] knitr_1.49      mkin_1.2.10    
 
 loaded via a namespace (and not attached):
- [1] jsonlite_1.8.9    gtable_0.3.6      dplyr_1.1.4       compiler_4.4.2   
- [5] tidyselect_1.2.1  colorout_1.3-2    tinytex_0.54      gridExtra_2.3    
- [9] jquerylib_0.1.4   scales_1.3.0      yaml_2.3.10       fastmap_1.2.0    
-[13] lattice_0.22-6    ggplot2_3.5.1     R6_2.5.1          generics_0.1.3   
-[17] lmtest_0.9-40     MASS_7.3-61       tibble_3.2.1      munsell_0.5.1    
-[21] bslib_0.8.0       pillar_1.9.0      rlang_1.1.4       utf8_1.2.4       
-[25] cachem_1.1.0      xfun_0.49         sass_0.4.9        cli_3.6.3        
-[29] magrittr_2.0.3    digest_0.6.37     grid_4.4.2        mclust_6.1.1     
-[33] lifecycle_1.0.4   nlme_3.1-166      vctrs_0.6.5       evaluate_1.0.1   
-[37] glue_1.8.0        codetools_0.2-20  zoo_1.8-12        fansi_1.0.6      
-[41] colorspace_2.1-1  tools_4.4.2       pkgconfig_2.0.3   htmltools_0.5.8.1
+ [1] gtable_0.3.6      jsonlite_1.9.0    dplyr_1.1.4       compiler_4.5.0   
+ [5] tidyselect_1.2.1  gridExtra_2.3     jquerylib_0.1.4   systemfonts_1.2.1
+ [9] scales_1.3.0      textshaping_1.0.0 yaml_2.3.10       fastmap_1.2.0    
+[13] lattice_0.22-6    ggplot2_3.5.1     R6_2.6.1          generics_0.1.3   
+[17] lmtest_0.9-40     MASS_7.3-65       htmlwidgets_1.6.4 tibble_3.2.1     
+[21] desc_1.4.3        munsell_0.5.1     bslib_0.9.0       pillar_1.10.1    
+[25] rlang_1.1.5       cachem_1.1.0      xfun_0.51         fs_1.6.5         
+[29] sass_0.4.9        cli_3.6.4         pkgdown_2.1.1     magrittr_2.0.3   
+[33] digest_0.6.37     grid_4.5.0        mclust_6.1.1      lifecycle_1.0.4  
+[37] nlme_3.1-168      vctrs_0.6.5       evaluate_1.0.3    glue_1.8.0       
+[41] codetools_0.2-20  ragg_1.3.3        zoo_1.8-13        colorspace_2.1-1 
+[45] tools_4.5.0       pkgconfig_2.0.3   htmltools_0.5.8.1

Hardware info

Hardware info +

CPU model: AMD Ryzen 9 7950X 16-Core Processor

MemTotal:       64927788 kB

MemTotal:       64927780 kB

- - - - +

- - - + - - - - - + -- cgit v1.2.1

Hierarchical mkin fit of the SFO model with error model const
	npar	AIC	Chisq	Df	Pr(>Chisq)
f_saem_dfop_tc_no_ranef_k2

Testing hierarchical parent degradation -kinetics with residue data on dimethenamid and dimethenamid-P

Johannes Ranke

Last change on 5 January 2023, last compiled on 13 -Februar 2025

Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P

Johannes +Ranke

Last change on 5 January +2023, last compiled on 12 Mai 2025

Introduction

Introduction +

Data

Data +

Separate evaluations

Separate evaluations +

Hierarchichal model fits

Hierarchichal model fits +

Parameter identifiability based on the Fisher Information -Matrix

Parameter identifiability based on the Fisher Information +Matrix +

Alternative check of parameter identifiability

Alternative check of parameter identifiability +

Conclusions

Conclusions +

Acknowledgements

Acknowledgements +

References

References +

Appendix

Hierarchical model fit listings

Appendix +

Hierarchical model fit listings +

Hierarchical model convergence plots

Hierarchical model convergence plots +

Session info

Session info +

Hardware info

Hardware info +