Details
-
Data for the first dataset are from p. 685. Data for the other four
+
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
in the soil from Site I2, where the values given on p. 694 were used.
@@ -195,13 +197,13 @@ specific pieces of information in the comments.
#> d_DCA/dt = + f_DCP_to_DCA * k_DCP * DCP - k_DCA * DCA
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)
+
print(f_D24_2014_const)
#> <mmkin> object
#> Status of individual fits:
#>
@@ -212,8 +214,34 @@ specific pieces of information in the comments.
#>
#> Z: Observations with value of zero were removed from the data
#> OK: No warnings
-
plot(f_D24_2014[, 3])
+
plot(f_D24_2014_const[, 3])

+
f_D24_2014_tc <- update(f_D24_2014_const, error_model = "tc")
+
print(f_D24_2014_tc)
+
#> <mmkin> object
+
#> Status of individual fits:
+
#>
+
#> dataset
+
#> model Mississippi Fayette RefSol 03-G Site E1 Site I2
+
#> SFO-SFO-SFO OK Z Z Z Z
+
#> DFOP-SFO-SFO OK Z Z Z Z
+
#>
+
#> Z: Observations with value of zero were removed from the data
+
#> OK: No warnings
+
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]])
+
#> df AIC
+
#> f_D24_2014_const[["SFO-SFO-SFO", 3]] 7 220.4991
+
#> f_D24_2014_const[["DFOP-SFO-SFO", 3]] 9 179.8933
+
#> f_D24_2014_tc[["SFO-SFO-SFO", 3]] 8 200.7663
+
#> f_D24_2014_tc[["DFOP-SFO-SFO", 3]] 10 181.8933
# }
diff --git a/docs/search.json b/docs/search.json
index 79534c2e..994e0347 100644
--- a/docs/search.json
+++ b/docs/search.json
@@ -1 +1 @@
-[{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"laboratory-data-l1","dir":"Articles","previous_headings":"","what":"Laboratory Data L1","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"following code defines example dataset L1 FOCUS kinetics report, p. 284: use assumptions simple first order (SFO), case declining rate constant time (FOMC) case two different phases kinetics (DFOP). detailed discussion models, please see FOCUS kinetics report. Since mkin version 0.9-32 (July 2014), can use shorthand notation like \"SFO\" parent degradation models. following two lines fit model produce summary report model fit. covers numerical analysis given FOCUS report. plot fit obtained plot function mkinfit objects. residual plot can easily obtained comparison, FOMC model fitted well, χ2\\chi^2 error level checked. get warning default optimisation algorithm Port converge, indication model overparameterised, .e. contains many parameters ill-defined consequence. fact, due higher number parameters, lower number degrees freedom fit, χ2\\chi^2 error level actually higher FOMC model (3.6%) SFO model (3.4%). Additionally, parameters log_alpha log_beta internally fitted model excessive confidence intervals, span 25 orders magnitude (!) backtransformed scale alpha beta. Also, t-test significant difference zero indicate significant difference, p-values greater 0.1, finally, parameter correlation log_alpha log_beta 1.000, clearly indicating model overparameterised. χ2\\chi^2 error levels reported Appendix 3 Appendix 7 FOCUS kinetics report rounded integer percentages partly deviate one percentage point results calculated mkin. reason known. However, mkin gives χ2\\chi^2 error levels kinfit package calculation routines kinfit package extensively compared results obtained KinGUI software, documented kinfit package vignette. KinGUI first widely used standard package field. Also, calculation χ2\\chi^2 error levels compared KinGUII, CAKE DegKin manager project sponsored German Umweltbundesamt (Ranke 2014).","code":"library(\"mkin\", quietly = TRUE) FOCUS_2006_L1 = data.frame( t = rep(c(0, 1, 2, 3, 5, 7, 14, 21, 30), each = 2), parent = c(88.3, 91.4, 85.6, 84.5, 78.9, 77.6, 72.0, 71.9, 50.3, 59.4, 47.0, 45.1, 27.7, 27.3, 10.0, 10.4, 2.9, 4.0)) FOCUS_2006_L1_mkin <- mkin_wide_to_long(FOCUS_2006_L1) m.L1.SFO <- mkinfit(\"SFO\", FOCUS_2006_L1_mkin, quiet = TRUE) summary(m.L1.SFO) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:26 2025 ## Date of summary: Mon May 12 21:54:26 2025 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## ## Fitted using 133 model solutions performed in 0.011 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 89.85 state ## k_parent 0.10 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 89.850000 -Inf Inf ## log_k_parent -2.302585 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 93.88778 96.5589 -43.94389 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 92.470 1.28200 89.740 95.200 ## log_k_parent -2.347 0.03763 -2.428 -2.267 ## sigma 2.780 0.46330 1.792 3.767 ## ## Parameter correlation: ## parent_0 log_k_parent sigma ## parent_0 1.000e+00 6.186e-01 -1.516e-09 ## log_k_parent 6.186e-01 1.000e+00 -3.124e-09 ## sigma -1.516e-09 -3.124e-09 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 92.47000 72.13 8.824e-21 89.74000 95.2000 ## k_parent 0.09561 26.57 2.487e-14 0.08824 0.1036 ## sigma 2.78000 6.00 1.216e-05 1.79200 3.7670 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 3.424 2 7 ## parent 3.424 2 7 ## ## Estimated disappearance times: ## DT50 DT90 ## parent 7.249 24.08 ## ## Data: ## time variable observed predicted residual ## 0 parent 88.3 92.471 -4.1710 ## 0 parent 91.4 92.471 -1.0710 ## 1 parent 85.6 84.039 1.5610 ## 1 parent 84.5 84.039 0.4610 ## 2 parent 78.9 76.376 2.5241 ## 2 parent 77.6 76.376 1.2241 ## 3 parent 72.0 69.412 2.5884 ## 3 parent 71.9 69.412 2.4884 ## 5 parent 50.3 57.330 -7.0301 ## 5 parent 59.4 57.330 2.0699 ## 7 parent 47.0 47.352 -0.3515 ## 7 parent 45.1 47.352 -2.2515 ## 14 parent 27.7 24.247 3.4528 ## 14 parent 27.3 24.247 3.0528 ## 21 parent 10.0 12.416 -2.4163 ## 21 parent 10.4 12.416 -2.0163 ## 30 parent 2.9 5.251 -2.3513 ## 30 parent 4.0 5.251 -1.2513 plot(m.L1.SFO, show_errmin = TRUE, main = \"FOCUS L1 - SFO\") mkinresplot(m.L1.SFO, ylab = \"Observed\", xlab = \"Time\") m.L1.FOMC <- mkinfit(\"FOMC\", FOCUS_2006_L1_mkin, quiet=TRUE) ## Warning in mkinfit(\"FOMC\", FOCUS_2006_L1_mkin, quiet = TRUE): Optimisation did not converge: ## false convergence (8) plot(m.L1.FOMC, show_errmin = TRUE, main = \"FOCUS L1 - FOMC\") summary(m.L1.FOMC, data = FALSE) ## Warning in sqrt(diag(covar)): NaNs produced ## Warning in cov2cor(ans$covar): diag(V) had non-positive or NA entries; the ## non-finite result may be dubious ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:26 2025 ## Date of summary: Mon May 12 21:54:26 2025 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## ## Fitted using 342 model solutions performed in 0.024 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 89.85 state ## alpha 1.00 deparm ## beta 10.00 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 89.850000 -Inf Inf ## log_alpha 0.000000 -Inf Inf ## log_beta 2.302585 -Inf Inf ## ## Fixed parameter values: ## None ## ## ## Warning(s): ## Optimisation did not converge: ## false convergence (8) ## ## Results: ## ## AIC BIC logLik ## 95.88782 99.44931 -43.94391 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 92.47 1.2820 89.720 95.220 ## log_alpha 13.20 NaN NaN NaN ## log_beta 15.54 NaN NaN NaN ## sigma 2.78 0.4607 1.792 3.768 ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma ## parent_0 1.000000 NaN NaN 0.000603 ## log_alpha NaN 1 NaN NaN ## log_beta NaN NaN 1 NaN ## sigma 0.000603 NaN NaN 1.000000 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 9.247e+01 NA NA 89.720 95.220 ## alpha 5.386e+05 NA NA NA NA ## beta 5.633e+06 NA NA NA NA ## sigma 2.780e+00 NA NA 1.792 3.768 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 3.619 3 6 ## parent 3.619 3 6 ## ## Estimated disappearance times: ## DT50 DT90 DT50back ## parent 7.249 24.08 7.249"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"laboratory-data-l2","dir":"Articles","previous_headings":"","what":"Laboratory Data L2","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"following code defines example dataset L2 FOCUS kinetics report, p. 287:","code":"FOCUS_2006_L2 = data.frame( t = rep(c(0, 1, 3, 7, 14, 28), each = 2), parent = c(96.1, 91.8, 41.4, 38.7, 19.3, 22.3, 4.6, 4.6, 2.6, 1.2, 0.3, 0.6)) FOCUS_2006_L2_mkin <- mkin_wide_to_long(FOCUS_2006_L2)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"sfo-fit-for-l2","dir":"Articles","previous_headings":"Laboratory Data L2","what":"SFO fit for L2","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":", SFO model fitted result plotted. residual plot can obtained simply adding argument show_residuals plot command. χ2\\chi^2 error level 14% suggests model fit well. also obvious plots fit, included residual plot. FOCUS kinetics report, stated apparent systematic error observed residual plot measured DT90 (approximately day 5), underestimation beyond point. may add difficult judge random nature residuals just three samplings days 0, 1 3. Also, clear priori consistent underestimation approximate DT90 irrelevant. However, can rationalised fact FOCUS fate models generally implement SFO kinetics.","code":"m.L2.SFO <- mkinfit(\"SFO\", FOCUS_2006_L2_mkin, quiet=TRUE) plot(m.L2.SFO, show_residuals = TRUE, show_errmin = TRUE, main = \"FOCUS L2 - SFO\")"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"fomc-fit-for-l2","dir":"Articles","previous_headings":"Laboratory Data L2","what":"FOMC fit for L2","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"comparison, FOMC model fitted well, χ2\\chi^2 error level checked. error level χ2\\chi^2 test passes much lower case. Therefore, FOMC model provides better description data, less experimental error assumed order explain data.","code":"m.L2.FOMC <- mkinfit(\"FOMC\", FOCUS_2006_L2_mkin, quiet = TRUE) plot(m.L2.FOMC, show_residuals = TRUE, main = \"FOCUS L2 - FOMC\") summary(m.L2.FOMC, data = FALSE) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:26 2025 ## Date of summary: Mon May 12 21:54:26 2025 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## ## Fitted using 239 model solutions performed in 0.015 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 93.95 state ## alpha 1.00 deparm ## beta 10.00 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 93.950000 -Inf Inf ## log_alpha 0.000000 -Inf Inf ## log_beta 2.302585 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 61.78966 63.72928 -26.89483 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 93.7700 1.6130 90.05000 97.4900 ## log_alpha 0.3180 0.1559 -0.04149 0.6776 ## log_beta 0.2102 0.2493 -0.36460 0.7850 ## sigma 2.2760 0.4645 1.20500 3.3470 ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma ## parent_0 1.000e+00 -1.151e-01 -2.085e-01 -7.436e-09 ## log_alpha -1.151e-01 1.000e+00 9.741e-01 -1.617e-07 ## log_beta -2.085e-01 9.741e-01 1.000e+00 -1.386e-07 ## sigma -7.436e-09 -1.617e-07 -1.386e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 93.770 58.120 4.267e-12 90.0500 97.490 ## alpha 1.374 6.414 1.030e-04 0.9594 1.969 ## beta 1.234 4.012 1.942e-03 0.6945 2.192 ## sigma 2.276 4.899 5.977e-04 1.2050 3.347 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 6.205 3 3 ## parent 6.205 3 3 ## ## Estimated disappearance times: ## DT50 DT90 DT50back ## parent 0.8092 5.356 1.612"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"dfop-fit-for-l2","dir":"Articles","previous_headings":"Laboratory Data L2","what":"DFOP fit for L2","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"Fitting four parameter DFOP model reduces χ2\\chi^2 error level. , DFOP model clearly best-fit model dataset L2 based chi^2 error level criterion.","code":"m.L2.DFOP <- mkinfit(\"DFOP\", FOCUS_2006_L2_mkin, quiet = TRUE) plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, main = \"FOCUS L2 - DFOP\") summary(m.L2.DFOP, data = FALSE) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:27 2025 ## Date of summary: Mon May 12 21:54:27 2025 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * ## time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) ## * parent ## ## Model predictions using solution type analytical ## ## Fitted using 581 model solutions performed in 0.043 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 93.95 state ## k1 0.10 deparm ## k2 0.01 deparm ## g 0.50 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 93.950000 -Inf Inf ## log_k1 -2.302585 -Inf Inf ## log_k2 -4.605170 -Inf Inf ## g_qlogis 0.000000 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 52.36695 54.79148 -21.18347 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 93.950 9.998e-01 91.5900 96.3100 ## log_k1 3.113 1.849e+03 -4369.0000 4375.0000 ## log_k2 -1.088 6.285e-02 -1.2370 -0.9394 ## g_qlogis -0.399 9.946e-02 -0.6342 -0.1638 ## sigma 1.414 2.886e-01 0.7314 2.0960 ## ## Parameter correlation: ## parent_0 log_k1 log_k2 g_qlogis sigma ## parent_0 1.000e+00 6.763e-07 -8.944e-10 2.665e-01 -1.083e-09 ## log_k1 6.763e-07 1.000e+00 1.112e-04 -2.187e-04 -1.027e-05 ## log_k2 -8.944e-10 1.112e-04 1.000e+00 -7.903e-01 9.464e-09 ## g_qlogis 2.665e-01 -2.187e-04 -7.903e-01 1.000e+00 -1.532e-08 ## sigma -1.083e-09 -1.027e-05 9.464e-09 -1.532e-08 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 93.9500 9.397e+01 2.036e-12 91.5900 96.3100 ## k1 22.4900 5.533e-04 4.998e-01 0.0000 Inf ## k2 0.3369 1.591e+01 4.697e-07 0.2904 0.3909 ## g 0.4016 1.680e+01 3.238e-07 0.3466 0.4591 ## sigma 1.4140 4.899e+00 8.776e-04 0.7314 2.0960 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 2.53 4 2 ## parent 2.53 4 2 ## ## Estimated disappearance times: ## DT50 DT90 DT50back DT50_k1 DT50_k2 ## parent 0.5335 5.311 1.599 0.03083 2.058"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"laboratory-data-l3","dir":"Articles","previous_headings":"","what":"Laboratory Data L3","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"following code defines example dataset L3 FOCUS kinetics report, p. 290.","code":"FOCUS_2006_L3 = data.frame( t = c(0, 3, 7, 14, 30, 60, 91, 120), parent = c(97.8, 60, 51, 43, 35, 22, 15, 12)) FOCUS_2006_L3_mkin <- mkin_wide_to_long(FOCUS_2006_L3)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"fit-multiple-models","dir":"Articles","previous_headings":"Laboratory Data L3","what":"Fit multiple models","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"mkin version 0.9-39 (June 2015), can fit several models one datasets one call function mmkin. datasets passed list, case named list holding L3 dataset prepared . χ2\\chi^2 error level 21% well plot suggest SFO model fit well. FOMC model performs better, error level χ2\\chi^2 test passes 7%. Fitting four parameter DFOP model reduces χ2\\chi^2 error level considerably.","code":"# Only use one core here, not to offend the CRAN checks mm.L3 <- mmkin(c(\"SFO\", \"FOMC\", \"DFOP\"), cores = 1, list(\"FOCUS L3\" = FOCUS_2006_L3_mkin), quiet = TRUE) plot(mm.L3)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"accessing-mmkin-objects","dir":"Articles","previous_headings":"Laboratory Data L3","what":"Accessing mmkin objects","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"objects returned mmkin arranged like matrix, models row index datasets column index. can extract summary plot e.g. DFOP fit, using square brackets indexing result use summary plot functions working mkinfit objects. , look model plot, confidence intervals parameters correlation matrix suggest parameter estimates reliable, DFOP model can used best-fit model based χ2\\chi^2 error level criterion laboratory data L3. also example standard t-test parameter g_ilr misleading, tests significant difference zero. case, zero appears correct value parameter, confidence interval backtransformed parameter g quite narrow.","code":"summary(mm.L3[[\"DFOP\", 1]]) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:27 2025 ## Date of summary: Mon May 12 21:54:27 2025 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * ## time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) ## * parent ## ## Model predictions using solution type analytical ## ## Fitted using 376 model solutions performed in 0.024 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 97.80 state ## k1 0.10 deparm ## k2 0.01 deparm ## g 0.50 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 97.800000 -Inf Inf ## log_k1 -2.302585 -Inf Inf ## log_k2 -4.605170 -Inf Inf ## g_qlogis 0.000000 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 32.97732 33.37453 -11.48866 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 97.7500 1.01900 94.5000 101.000000 ## log_k1 -0.6612 0.10050 -0.9812 -0.341300 ## log_k2 -4.2860 0.04322 -4.4230 -4.148000 ## g_qlogis -0.1739 0.05270 -0.3416 -0.006142 ## sigma 1.0170 0.25430 0.2079 1.827000 ## ## Parameter correlation: ## parent_0 log_k1 log_k2 g_qlogis sigma ## parent_0 1.000e+00 1.732e-01 2.282e-02 4.009e-01 -9.696e-08 ## log_k1 1.732e-01 1.000e+00 4.945e-01 -5.809e-01 7.148e-07 ## log_k2 2.282e-02 4.945e-01 1.000e+00 -6.812e-01 1.022e-06 ## g_qlogis 4.009e-01 -5.809e-01 -6.812e-01 1.000e+00 -7.930e-07 ## sigma -9.696e-08 7.148e-07 1.022e-06 -7.930e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 97.75000 95.960 1.248e-06 94.50000 101.00000 ## k1 0.51620 9.947 1.081e-03 0.37490 0.71090 ## k2 0.01376 23.140 8.840e-05 0.01199 0.01579 ## g 0.45660 34.920 2.581e-05 0.41540 0.49850 ## sigma 1.01700 4.000 1.400e-02 0.20790 1.82700 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 2.225 4 4 ## parent 2.225 4 4 ## ## Estimated disappearance times: ## DT50 DT90 DT50back DT50_k1 DT50_k2 ## parent 7.464 123 37.03 1.343 50.37 ## ## Data: ## time variable observed predicted residual ## 0 parent 97.8 97.75 0.05396 ## 3 parent 60.0 60.45 -0.44933 ## 7 parent 51.0 49.44 1.56338 ## 14 parent 43.0 43.84 -0.83632 ## 30 parent 35.0 35.15 -0.14707 ## 60 parent 22.0 23.26 -1.25919 ## 91 parent 15.0 15.18 -0.18181 ## 120 parent 12.0 10.19 1.81395 plot(mm.L3[[\"DFOP\", 1]], show_errmin = TRUE)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/FOCUS_L.html","id":"laboratory-data-l4","dir":"Articles","previous_headings":"","what":"Laboratory Data L4","title":"Example evaluation of FOCUS Laboratory Data L1 to L3","text":"following code defines example dataset L4 FOCUS kinetics report, p. 293: Fits SFO FOMC models, plots summaries produced : χ2\\chi^2 error level 3.3% well plot suggest SFO model fits well. error level χ2\\chi^2 test passes slightly lower FOMC model. However, difference appears negligible.","code":"FOCUS_2006_L4 = data.frame( t = c(0, 3, 7, 14, 30, 60, 91, 120), parent = c(96.6, 96.3, 94.3, 88.8, 74.9, 59.9, 53.5, 49.0)) FOCUS_2006_L4_mkin <- mkin_wide_to_long(FOCUS_2006_L4) # Only use one core here, not to offend the CRAN checks mm.L4 <- mmkin(c(\"SFO\", \"FOMC\"), cores = 1, list(\"FOCUS L4\" = FOCUS_2006_L4_mkin), quiet = TRUE) plot(mm.L4) summary(mm.L4[[\"SFO\", 1]], data = FALSE) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:27 2025 ## Date of summary: Mon May 12 21:54:28 2025 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## ## Fitted using 142 model solutions performed in 0.01 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 96.6 state ## k_parent 0.1 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 96.600000 -Inf Inf ## log_k_parent -2.302585 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 47.12133 47.35966 -20.56067 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 96.440 1.69900 92.070 100.800 ## log_k_parent -5.030 0.07059 -5.211 -4.848 ## sigma 3.162 0.79050 1.130 5.194 ## ## Parameter correlation: ## parent_0 log_k_parent sigma ## parent_0 1.000e+00 5.938e-01 3.430e-07 ## log_k_parent 5.938e-01 1.000e+00 5.885e-07 ## sigma 3.430e-07 5.885e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 96.440000 56.77 1.604e-08 92.070000 1.008e+02 ## k_parent 0.006541 14.17 1.578e-05 0.005455 7.842e-03 ## sigma 3.162000 4.00 5.162e-03 1.130000 5.194e+00 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 3.287 2 6 ## parent 3.287 2 6 ## ## Estimated disappearance times: ## DT50 DT90 ## parent 106 352 summary(mm.L4[[\"FOMC\", 1]], data = FALSE) ## mkin version used for fitting: 1.2.10 ## R version used for fitting: 4.5.0 ## Date of fit: Mon May 12 21:54:27 2025 ## Date of summary: Mon May 12 21:54:28 2025 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## ## Fitted using 224 model solutions performed in 0.014 s ## ## Error model: Constant variance ## ## Error model algorithm: OLS ## ## Starting values for parameters to be optimised: ## value type ## parent_0 96.6 state ## alpha 1.0 deparm ## beta 10.0 deparm ## ## Starting values for the transformed parameters actually optimised: ## value lower upper ## parent_0 96.600000 -Inf Inf ## log_alpha 0.000000 -Inf Inf ## log_beta 2.302585 -Inf Inf ## ## Fixed parameter values: ## None ## ## Results: ## ## AIC BIC logLik ## 40.37255 40.69032 -16.18628 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 99.1400 1.2670 95.6300 102.7000 ## log_alpha -0.3506 0.2616 -1.0770 0.3756 ## log_beta 4.1740 0.3938 3.0810 5.2670 ## sigma 1.8300 0.4575 0.5598 3.1000 ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma ## parent_0 1.000e+00 -4.696e-01 -5.543e-01 -2.447e-07 ## log_alpha -4.696e-01 1.000e+00 9.889e-01 2.198e-08 ## log_beta -5.543e-01 9.889e-01 1.000e+00 4.923e-08 ## sigma -2.447e-07 2.198e-08 4.923e-08 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 99.1400 78.250 7.993e-08 95.6300 102.700 ## alpha 0.7042 3.823 9.365e-03 0.3407 1.456 ## beta 64.9800 2.540 3.201e-02 21.7800 193.900 ## sigma 1.8300 4.000 8.065e-03 0.5598 3.100 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df ## All data 2.029 3 5 ## parent 2.029 3 5 ## ## Estimated disappearance times: ## DT50 DT90 DT50back ## parent 108.9 1644 494.9"},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"abstract","dir":"Articles","previous_headings":"","what":"Abstract","title":"Short introduction to mkin","text":"regulatory evaluation chemical substances like plant protection products (pesticides), biocides chemicals, degradation data play important role. evaluation pesticide degradation experiments, detailed guidance developed, based nonlinear optimisation. R add-package mkin implements fitting models recommended guidance within R calculates statistical measures data series within one compartments, parent metabolites.","code":"library(\"mkin\", quietly = TRUE) # Define the kinetic model m_SFO_SFO_SFO <- mkinmod(parent = mkinsub(\"SFO\", \"M1\"), M1 = mkinsub(\"SFO\", \"M2\"), M2 = mkinsub(\"SFO\"), use_of_ff = \"max\", quiet = TRUE) # Produce model predictions using some arbitrary parameters sampling_times = c(0, 1, 3, 7, 14, 28, 60, 90, 120) d_SFO_SFO_SFO <- mkinpredict(m_SFO_SFO_SFO, c(k_parent = 0.03, f_parent_to_M1 = 0.5, k_M1 = log(2)/100, f_M1_to_M2 = 0.9, k_M2 = log(2)/50), c(parent = 100, M1 = 0, M2 = 0), sampling_times) # Generate a dataset by adding normally distributed errors with # standard deviation 3, for two replicates at each sampling time d_SFO_SFO_SFO_err <- add_err(d_SFO_SFO_SFO, reps = 2, sdfunc = function(x) 3, n = 1, seed = 123456789 ) # Fit the model to the dataset f_SFO_SFO_SFO <- mkinfit(m_SFO_SFO_SFO, d_SFO_SFO_SFO_err[[1]], quiet = TRUE) # Plot the results separately for parent and metabolites plot_sep(f_SFO_SFO_SFO, lpos = c(\"topright\", \"bottomright\", \"bottomright\"))"},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"background","dir":"Articles","previous_headings":"","what":"Background","title":"Short introduction to mkin","text":"mkin package (J. Ranke 2021) implements approach degradation kinetics recommended kinetics report provided FOrum Co-ordination pesticide fate models USe (FOCUS Work Group Degradation Kinetics 2006, 2014). covers data series describing decline one compound, data series transformation products (commonly termed metabolites) data series one compartment. possible include back reactions. Therefore, equilibrium reactions equilibrium partitioning can specified, although often leads overparameterisation model. first mkin code published 2010, commonly used tools fitting complex kinetic degradation models experimental data KinGUI (Schäfer et al. 2007), MATLAB based tool graphical user interface specifically tailored task included output proposed FOCUS Kinetics Workgroup, ModelMaker, general purpose compartment based tool providing infrastructure fitting dynamic simulation models based differential equations data. ‘mkin’ code first uploaded BerliOS development platform. taken , version control history imported R-Forge site (see e.g. initial commit 11 May 2010), code still updated. time, R package FME (Flexible Modelling Environment) (Soetaert Petzoldt 2010) already available, provided good basis developing package specifically tailored task. remaining challenge make easy possible users (including author vignette) specify system differential equations include output requested FOCUS guidance, χ2\\chi^2 error level defined guidance. Also, mkin introduced using analytical solutions parent kinetics improved optimization speed. Later, Eigenvalue based solutions introduced mkin case linear differential equations (.e. FOMC DFOP models used parent compound), greatly improving optimization speed cases. , become somehow obsolete, use compiled code described gives even faster execution times. possibility specify back-reactions biphasic model (SFORB) metabolites present mkin beginning.","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"derived-software-tools","dir":"Articles","previous_headings":"Background","what":"Derived software tools","title":"Short introduction to mkin","text":"Soon publication mkin, two derived tools published, namely KinGUII (developed Bayer Crop Science) CAKE (commissioned Tessella Syngenta), added graphical user interface (GUI), added fitting iteratively reweighted least squares (IRLS) characterisation likely parameter distributions Markov Chain Monte Carlo (MCMC) sampling. CAKE focuses smooth use experience, sacrificing flexibility model definition, originally allowing two primary metabolites parallel. current version 3.4 CAKE released May 2020 uses scheme six metabolites flexible arrangement supports biphasic modelling metabolites, support back-reactions (non-instantaneous equilibria). KinGUI offers even flexible widget specifying complex kinetic models. Back-reactions (non-instantaneous equilibria) supported early , 2014, simple first-order models specified transformation products. Starting KinGUII version 2.1, biphasic modelling metabolites also available KinGUII. graphical user interface (GUI) recently brought decent degree maturity browser based GUI named gmkin. Please see documentation page manual information. comparison scope, usability numerical results obtained tools recently published Johannes Ranke, Wöltjen, Meinecke (2018).","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"unique-features","dir":"Articles","previous_headings":"","what":"Unique features","title":"Short introduction to mkin","text":"Currently, main unique features available mkin speed increase using compiled code compiler present, parallel model fitting multicore machines using mmkin function, estimation parameter confidence intervals based transformed parameters (see ) possibility use two-component error model iteratively reweighted least squares fitting different variances variable introduced Gao et al. (2011) available mkin since version 0.9-22. release 0.9.49.5, IRLS algorithm complemented direct step-wise maximisation likelihood function, makes possible fit variance variable error model also two-component error model inspired error models developed analytical chemistry (Johannes Ranke Meinecke 2019).","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"internal-parameter-transformations","dir":"Articles","previous_headings":"","what":"Internal parameter transformations","title":"Short introduction to mkin","text":"rate constants, log transformation used, proposed Bates Watts (1988, 77, 149). Approximate intervals constructed transformed rate constants (compare Bates Watts 1988, 135), .e. logarithms. Confidence intervals rate constants obtained using appropriate backtransformation using exponential function. first version mkin allowing specifying models using formation fractions, home-made reparameterisation used order ensure sum formation fractions exceed unity. method still used current version KinGUII (v2.1 April 2014), modification allows fixing pathway sink zero. CAKE uses penalties objective function order enforce constraint. 2012, alternative reparameterisation formation fractions proposed together René Lehmann (J. Ranke Lehmann 2012), based isometric logratio transformation (ILR). aim improve validity linear approximation objective function parameter estimation procedure well subsequent calculation parameter confidence intervals. current version mkin, logit transformation used parameters bound 0 1, g parameter DFOP model.","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"confidence-intervals-based-on-transformed-parameters","dir":"Articles","previous_headings":"Internal parameter transformations","what":"Confidence intervals based on transformed parameters","title":"Short introduction to mkin","text":"first attempt providing improved parameter confidence intervals introduced mkin 2013, confidence intervals obtained FME transformed parameters simply backtransformed one one yield asymmetric confidence intervals backtransformed parameters. However, 1:1 relation rate constants model transformed parameters fitted model, parameters obtained isometric logratio transformation calculated set formation fractions quantify paths compounds formed specific parent compound, 1:1 relation exists. Therefore, parameter confidence intervals formation fractions obtained method appear valid case single transformation product, currently logit transformation used formation fraction. confidence intervals obtained backtransformation cases 1:1 relation transformed original parameter exist considered author vignette accurate obtained using re-estimation Hessian matrix backtransformation, implemented FME package.","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/mkin.html","id":"parameter-t-test-based-on-untransformed-parameters","dir":"Articles","previous_headings":"Internal parameter transformations","what":"Parameter t-test based on untransformed parameters","title":"Short introduction to mkin","text":"standard output many nonlinear regression software packages includes results test significant difference zero parameters. test also recommended check validity rate constants FOCUS guidance (FOCUS Work Group Degradation Kinetics 2014, 96ff). argued precondition test, .e. normal distribution estimator parameters, fulfilled case nonlinear regression (J. Ranke Lehmann 2015). However, test commonly used industry, consultants national authorities order decide reliability parameter estimates, based FOCUS guidance mentioned . Therefore, results one-sided t-test included summary output mkin. reasonable test significant difference transformed parameters (e.g. log(k)log(k)) zero, t-test calculated based model definition parameter transformation, .e. similar way packages apply internal parameter transformation. note included mkin output, pointing fact t-test based unjustified assumption normal distribution parameter estimators.","code":""},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"introduction","dir":"Articles > Prebuilt","previous_headings":"","what":"Introduction","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"purpose document test demonstrate nonlinear hierarchical models (NLHM) based parent degradation models SFO, FOMC, DFOP HS, serial formation two metabolites can fitted mkin package. assembled course work package 1.2 Project Number 173340 (Application nonlinear hierarchical models kinetic evaluation chemical degradation data) German Environment Agency carried 2022 2023. mkin package used version 1.2.10 currently development. newly introduced functionality used simplification excluding random effects set fits based related set fits reduced model, documentation starting parameters fit, starting parameters saem fits now listed summary. saemix package used backend fitting NLHM, also loaded make convergence plot function available. document processed knitr package, also provides kable function used improve display tabular data R markdown documents. parallel processing, parallel package used.","code":"library(mkin) library(knitr) library(saemix) library(parallel) n_cores <- detectCores() # We need to start a new cluster after defining a compiled model that is # saved as a DLL to the user directory, therefore we define a function # This is used again after defining the pathway model start_cluster <- function(n_cores) { if (Sys.info()[\"sysname\"] == \"Windows\") { ret <- makePSOCKcluster(n_cores) } else { ret <- makeForkCluster(n_cores) } return(ret) } cl <- start_cluster(n_cores)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"test-data","dir":"Articles > Prebuilt","previous_headings":"Introduction","what":"Test data","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"example data taken final addendum DAR 2014 distributed mkin package. Residue data time step normalisation factors read using function read_spreadsheet mkin package. function also performs time step normalisation. following tables show covariate data 5 datasets read spreadsheet file. Covariate data Dataset Nambsheim Dataset Tama Dataset Gross-Umstadt Dataset Sassafras Dataset Lleida","code":"data_file <- system.file( \"testdata\", \"cyantraniliprole_soil_efsa_2014.xlsx\", package = \"mkin\") cyan_ds <- read_spreadsheet(data_file, parent_only = FALSE) pH <- attr(cyan_ds, \"covariates\") kable(pH, caption = \"Covariate data\") for (ds_name in names(cyan_ds)) { print( kable(mkin_long_to_wide(cyan_ds[[ds_name]]), caption = paste(\"Dataset\", ds_name), booktabs = TRUE, row.names = FALSE)) cat(\"\\n\\\\clearpage\\n\") }"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"parent-only-evaluations","dir":"Articles > Prebuilt","previous_headings":"","what":"Parent only evaluations","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"pathway fits long run times, evaluations parent data performed first, order determine hierarchical parent degradation model random effects degradation model parameters ill-defined. fits converged successfully. almost models, random effect initial concentration parent compound ill-defined. biexponential models DFOP SFORB, random effect one additional parameter ill-defined two-component error model used. Model comparison based AIC BIC indicates two-component error model preferable parent models exception DFOP. lowest AIC BIC values obtained FOMC model, followed SFORB DFOP.","code":"cyan_sep_const <- mmkin(c(\"SFO\", \"FOMC\", \"DFOP\", \"SFORB\", \"HS\"), cyan_ds, quiet = TRUE, cores = n_cores) cyan_sep_tc <- update(cyan_sep_const, error_model = \"tc\") cyan_saem_full <- mhmkin(list(cyan_sep_const, cyan_sep_tc)) status(cyan_saem_full) |> kable() illparms(cyan_saem_full) |> kable() anova(cyan_saem_full) |> kable(digits = 1) stopCluster(cl)"},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"evaluations-with-pathway-established-previously","dir":"Articles > Prebuilt","previous_headings":"Pathway fits","what":"Evaluations with pathway established previously","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"test technical feasibility coupling relevant parent degradation models different transformation pathway models, list mkinmod models set . EU evaluation, parallel formation metabolites JCZ38 J9Z38 secondary formation metabolite JSE76 JCZ38 used. obtain suitable starting values NLHM fits, separate pathway fits performed datasets. separate fits converged successfully. biggest convergence problems seen using HS model constant variance. hierarchical pathway fits, random effects quantified corresponding parent data analyses excluded. code , output illparms function parent fits used argument no_random_effect mhmkin function. possibility introduced mkin version 1.2.2 currently development. status information individual fits shows fits completed successfully. matrix entries Fth FO indicate Fisher Information Matrix inverted fixed effects (theta) random effects (Omega), respectively. affected fits, ill-defined parameters determined using illparms function, relies Fisher Information Matrix. model comparisons suggest pathway fits using DFOP SFORB parent compound provide best fit. two parent model, successful fits shown . Plots fits parent models shown Appendix. DFOP pathway fit two-component error SFORB pathway fit two-component error closer graphical analysis Figures shows residues transformation product JCZ38 soils Tama Nambsheim observed later time points strongly systematically underestimated.","code":"if (!dir.exists(\"cyan_dlls\")) dir.create(\"cyan_dlls\") cyan_path_1 <- list( sfo_path_1 = mkinmod( cyan = mkinsub(\"SFO\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\"), quiet = TRUE, name = \"sfo_path_1\", dll_dir = \"cyan_dlls\", overwrite = TRUE), fomc_path_1 = mkinmod( cyan = mkinsub(\"FOMC\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\"), quiet = TRUE, name = \"fomc_path_1\", dll_dir = \"cyan_dlls\", overwrite = TRUE), dfop_path_1 = mkinmod( cyan = mkinsub(\"DFOP\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\"), quiet = TRUE, name = \"dfop_path_1\", dll_dir = \"cyan_dlls\", overwrite = TRUE), sforb_path_1 = mkinmod( cyan = mkinsub(\"SFORB\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\"), quiet = TRUE, name = \"sforb_path_1\", dll_dir = \"cyan_dlls\", overwrite = TRUE), hs_path_1 = mkinmod( cyan = mkinsub(\"HS\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\"), quiet = TRUE, name = \"hs_path_1\", dll_dir = \"cyan_dlls\", overwrite = TRUE) ) cl_path_1 <- start_cluster(n_cores) f_sep_1_const <- mmkin( cyan_path_1, cyan_ds, error_model = \"const\", cluster = cl_path_1, quiet = TRUE) status(f_sep_1_const) |> kable() f_sep_1_tc <- update(f_sep_1_const, error_model = \"tc\") status(f_sep_1_tc) |> kable() f_saem_1 <- mhmkin(list(f_sep_1_const, f_sep_1_tc), no_random_effect = illparms(cyan_saem_full), cluster = cl_path_1) status(f_saem_1) |> kable() illparms(f_saem_1) |> kable() anova(f_saem_1[, \"const\"]) |> kable(digits = 1) anova(f_saem_1[1:4, ]) |> kable(digits = 1) plot(f_saem_1[[\"dfop_path_1\", \"tc\"]]) plot(f_saem_1[[\"sforb_path_1\", \"tc\"]]) stopCluster(cl_path_1)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"alternative-pathway-fits","dir":"Articles > Prebuilt","previous_headings":"Pathway fits","what":"Alternative pathway fits","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"improve fit JCZ38, back-reaction JSE76 JCZ38 introduced alternative version transformation pathway, analogy back-reaction K5A78 K5A77. pairs transformation products pairs organic acid corresponding amide (Addendum 2014, p. 109). FOMC provided best fit parent, biexponential models DFOP SFORB provided best initial pathway fits, three parent models used alternative pathway fits. Using constant variance, separate fits converge exception fits Sassafras soil data. Using two-component error model, separate fits converge exception alternative pathway fit DFOP used parent Sassafras dataset. hierarchical fits alternative pathway completed successfully, exception model using FOMC parent compound constant variance error model. biphasic fits (DFOP SFORB parent compound), random effects formation fractions pathways JCZ38 JSE76, reverse pathway JSE76 JCZ38 ill-defined. variants using biexponential models DFOP SFORB parent compound two-component error model give lowest AIC BIC values plotted . Compared original pathway, AIC BIC values indicate large improvement. confirmed plots, show metabolite JCZ38 fitted much better model. FOMC pathway fit two-component error, alternative pathway DFOP pathway fit two-component error, alternative pathway SFORB pathway fit two-component error, alternative pathway","code":"cyan_path_2 <- list( fomc_path_2 = mkinmod( cyan = mkinsub(\"FOMC\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\", \"JCZ38\"), name = \"fomc_path_2\", quiet = TRUE, dll_dir = \"cyan_dlls\", overwrite = TRUE ), dfop_path_2 = mkinmod( cyan = mkinsub(\"DFOP\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\", \"JCZ38\"), name = \"dfop_path_2\", quiet = TRUE, dll_dir = \"cyan_dlls\", overwrite = TRUE ), sforb_path_2 = mkinmod( cyan = mkinsub(\"SFORB\", c(\"JCZ38\", \"J9Z38\")), JCZ38 = mkinsub(\"SFO\", \"JSE76\"), J9Z38 = mkinsub(\"SFO\"), JSE76 = mkinsub(\"SFO\", \"JCZ38\"), name = \"sforb_path_2\", quiet = TRUE, dll_dir = \"cyan_dlls\", overwrite = TRUE ) ) cl_path_2 <- start_cluster(n_cores) f_sep_2_const <- mmkin( cyan_path_2, cyan_ds, error_model = \"const\", cluster = cl_path_2, quiet = TRUE) status(f_sep_2_const) |> kable() f_sep_2_tc <- update(f_sep_2_const, error_model = \"tc\") status(f_sep_2_tc) |> kable() f_saem_2 <- mhmkin(list(f_sep_2_const, f_sep_2_tc), no_random_effect = illparms(cyan_saem_full[2:4, ]), cluster = cl_path_2) status(f_saem_2) |> kable() illparms(f_saem_2) |> kable() anova(f_saem_2[, \"tc\"]) |> kable(digits = 1) anova(f_saem_2[2:3,]) |> kable(digits = 1) plot(f_saem_2[[\"fomc_path_2\", \"tc\"]]) plot(f_saem_2[[\"dfop_path_2\", \"tc\"]]) plot(f_saem_2[[\"sforb_path_2\", \"tc\"]])"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"refinement-of-alternative-pathway-fits","dir":"Articles > Prebuilt","previous_headings":"Pathway fits","what":"Refinement of alternative pathway fits","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"ill-defined random effects identified parent fits pathway fits, excluded final evaluations . purpose, list character vectors created can indexed row column indices, contains degradation parameter names random effects excluded hierarchical fits contained f_saem_2. exception FOMC pathway fit constant variance, updated fits completed successfully. However, Fisher Information Matrix fixed effects (Fth) inverted, confidence intervals optimised parameters available. AIC BIC values best fit (DFOP pathway fit two-component error) lower previous fits alternative pathway, practical value refined evaluations limited confidence intervals obtained.","code":"no_ranef <- matrix(list(), nrow = 3, ncol = 2, dimnames = dimnames(f_saem_2)) no_ranef[[\"fomc_path_2\", \"const\"]] <- c(\"log_beta\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") no_ranef[[\"fomc_path_2\", \"tc\"]] <- c(\"cyan_0\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") no_ranef[[\"dfop_path_2\", \"const\"]] <- c(\"cyan_0\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") no_ranef[[\"dfop_path_2\", \"tc\"]] <- c(\"cyan_0\", \"log_k1\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") no_ranef[[\"sforb_path_2\", \"const\"]] <- c(\"cyan_free_0\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") no_ranef[[\"sforb_path_2\", \"tc\"]] <- c(\"cyan_free_0\", \"log_k_cyan_free_bound\", \"f_JCZ38_qlogis\", \"f_JSE76_qlogis\") clusterExport(cl_path_2, \"no_ranef\") f_saem_3 <- update(f_saem_2, no_random_effect = no_ranef, cluster = cl_path_2) status(f_saem_3) |> kable() illparms(f_saem_3) |> kable() anova(f_saem_3[, \"tc\"]) |> kable(digits = 1) anova(f_saem_3[2:3,]) |> kable(digits = 1) stopCluster(cl_path_2)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"conclusion","dir":"Articles > Prebuilt","previous_headings":"","what":"Conclusion","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"demonstrated relatively complex transformation pathway parallel formation two primary metabolites one secondary metabolite can fitted even data individual datasets quite different partly cover formation phase. run times pathway fits several hours, limiting practical feasibility iterative refinements based ill-defined parameters alternative checks parameter identifiability based multistart runs.","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"acknowledgements","dir":"Articles > Prebuilt","previous_headings":"","what":"Acknowledgements","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"helpful comments Janina Wöltjen German Environment Agency gratefully acknowledged.","code":""},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"plots-of-fits-that-were-not-refined-further","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Plots of fits that were not refined further","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"SFO pathway fit two-component error FOMC pathway fit two-component error HS pathway fit two-component error","code":"plot(f_saem_1[[\"sfo_path_1\", \"tc\"]]) plot(f_saem_1[[\"fomc_path_1\", \"tc\"]]) plot(f_saem_1[[\"sforb_path_1\", \"tc\"]])"},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"pathway-1","dir":"Articles > Prebuilt","previous_headings":"Appendix > Hierarchical fit listings","what":"Pathway 1","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"Hierarchical SFO path 1 fit constant variance Hierarchical SFO path 1 fit two-component error Hierarchical FOMC path 1 fit constant variance Hierarchical FOMC path 1 fit two-component error Hierarchical DFOP path 1 fit constant variance Hierarchical DFOP path 1 fit two-component error Hierarchical SFORB path 1 fit constant variance Hierarchical SFORB path 1 fit two-component error Hierarchical HS path 1 fit constant variance","code":"saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:35:03 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - k_cyan * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * k_cyan * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * k_cyan * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 445.08 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_cyan log_k_JCZ38 log_k_J9Z38 log_k_JSE76 95.3304 -3.8459 -3.1305 -5.0678 -5.3196 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis 0.8158 23.5335 11.8774 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_cyan log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_0 4.797 0.0000 0.000 0.000 0.0000 log_k_cyan 0.000 0.9619 0.000 0.000 0.0000 log_k_JCZ38 0.000 0.0000 2.139 0.000 0.0000 log_k_J9Z38 0.000 0.0000 0.000 1.639 0.0000 log_k_JSE76 0.000 0.0000 0.000 0.000 0.7894 f_cyan_ilr_1 0.000 0.0000 0.000 0.000 0.0000 f_cyan_ilr_2 0.000 0.0000 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.0000 0.000 0.000 0.0000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis cyan_0 0.0000 0.000 0.00 log_k_cyan 0.0000 0.000 0.00 log_k_JCZ38 0.0000 0.000 0.00 log_k_J9Z38 0.0000 0.000 0.00 log_k_JSE76 0.0000 0.000 0.00 f_cyan_ilr_1 0.7714 0.000 0.00 f_cyan_ilr_2 0.0000 9.247 0.00 f_JCZ38_qlogis 0.0000 0.000 16.61 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2693 2687 -1331 Optimised parameters: est. lower upper cyan_0 95.1279 9.354e+01 9.671e+01 log_k_cyan -3.8527 -4.367e+00 -3.338e+00 log_k_JCZ38 -3.0381 -4.187e+00 -1.889e+00 log_k_J9Z38 -5.0095 -5.623e+00 -4.396e+00 log_k_JSE76 -5.3357 -6.025e+00 -4.646e+00 f_cyan_ilr_1 0.8050 5.174e-01 1.093e+00 f_cyan_ilr_2 12.4820 -1.050e+06 1.051e+06 f_JCZ38_qlogis 1.2912 3.561e-01 2.226e+00 a.1 4.8393 NA NA SD.log_k_cyan 0.5840 NA NA SD.log_k_JCZ38 1.2740 NA NA SD.log_k_J9Z38 0.3172 NA NA SD.log_k_JSE76 0.5677 NA NA SD.f_cyan_ilr_1 0.2623 NA NA SD.f_cyan_ilr_2 1.3724 NA NA SD.f_JCZ38_qlogis 0.1464 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_cyan 0.5840 NA NA SD.log_k_JCZ38 1.2740 NA NA SD.log_k_J9Z38 0.3172 NA NA SD.log_k_JSE76 0.5677 NA NA SD.f_cyan_ilr_1 0.2623 NA NA SD.f_cyan_ilr_2 1.3724 NA NA SD.f_JCZ38_qlogis 0.1464 NA NA Variance model: est. lower upper a.1 4.839 NA NA Backtransformed parameters: est. lower upper cyan_0 95.127935 93.542456 96.713413 k_cyan 0.021221 0.012687 0.035497 k_JCZ38 0.047924 0.015189 0.151213 k_J9Z38 0.006674 0.003612 0.012332 k_JSE76 0.004817 0.002417 0.009601 f_cyan_to_JCZ38 0.757402 NA NA f_cyan_to_J9Z38 0.242597 NA NA f_JCZ38_to_JSE76 0.784347 0.588098 0.902582 Resulting formation fractions: ff cyan_JCZ38 7.574e-01 cyan_J9Z38 2.426e-01 cyan_sink 9.839e-08 JCZ38_JSE76 7.843e-01 JCZ38_sink 2.157e-01 Estimated disappearance times: DT50 DT90 cyan 32.66 108.50 JCZ38 14.46 48.05 J9Z38 103.86 345.00 JSE76 143.91 478.04 saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:34:47 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - k_cyan * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * k_cyan * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * k_cyan * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 428.46 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_cyan log_k_JCZ38 log_k_J9Z38 log_k_JSE76 96.0039 -3.8907 -3.1276 -5.0069 -4.9367 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis 0.7937 22.3422 17.8932 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_cyan log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_0 4.859 0.000 0.00 0.00 0.0000 log_k_cyan 0.000 0.962 0.00 0.00 0.0000 log_k_JCZ38 0.000 0.000 2.04 0.00 0.0000 log_k_J9Z38 0.000 0.000 0.00 1.72 0.0000 log_k_JSE76 0.000 0.000 0.00 0.00 0.9076 f_cyan_ilr_1 0.000 0.000 0.00 0.00 0.0000 f_cyan_ilr_2 0.000 0.000 0.00 0.00 0.0000 f_JCZ38_qlogis 0.000 0.000 0.00 0.00 0.0000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis cyan_0 0.0000 0.000 0.00 log_k_cyan 0.0000 0.000 0.00 log_k_JCZ38 0.0000 0.000 0.00 log_k_J9Z38 0.0000 0.000 0.00 log_k_JSE76 0.0000 0.000 0.00 f_cyan_ilr_1 0.7598 0.000 0.00 f_cyan_ilr_2 0.0000 8.939 0.00 f_JCZ38_qlogis 0.0000 0.000 14.49 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2658 2651 -1312 Optimised parameters: est. lower upper cyan_0 94.81681 NA NA log_k_cyan -3.91558 NA NA log_k_JCZ38 -3.12715 NA NA log_k_J9Z38 -5.04840 NA NA log_k_JSE76 -5.10443 NA NA f_cyan_ilr_1 0.80760 NA NA f_cyan_ilr_2 48.66960 NA NA f_JCZ38_qlogis 3.03397 NA NA a.1 3.93879 NA NA b.1 0.08057 NA NA SD.log_k_cyan 0.58921 NA NA SD.log_k_JCZ38 1.29813 NA NA SD.log_k_J9Z38 0.68372 NA NA SD.log_k_JSE76 0.35128 NA NA SD.f_cyan_ilr_1 0.38352 NA NA SD.f_cyan_ilr_2 4.98884 NA NA SD.f_JCZ38_qlogis 1.75636 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_cyan 0.5892 NA NA SD.log_k_JCZ38 1.2981 NA NA SD.log_k_J9Z38 0.6837 NA NA SD.log_k_JSE76 0.3513 NA NA SD.f_cyan_ilr_1 0.3835 NA NA SD.f_cyan_ilr_2 4.9888 NA NA SD.f_JCZ38_qlogis 1.7564 NA NA Variance model: est. lower upper a.1 3.93879 NA NA b.1 0.08057 NA NA Backtransformed parameters: est. lower upper cyan_0 94.81681 NA NA k_cyan 0.01993 NA NA k_JCZ38 0.04384 NA NA k_J9Z38 0.00642 NA NA k_JSE76 0.00607 NA NA f_cyan_to_JCZ38 0.75807 NA NA f_cyan_to_J9Z38 0.24193 NA NA f_JCZ38_to_JSE76 0.95409 NA NA Resulting formation fractions: ff cyan_JCZ38 0.75807 cyan_J9Z38 0.24193 cyan_sink 0.00000 JCZ38_JSE76 0.95409 JCZ38_sink 0.04591 Estimated disappearance times: DT50 DT90 cyan 34.78 115.54 JCZ38 15.81 52.52 J9Z38 107.97 358.68 JSE76 114.20 379.35 saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:35:21 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 462.739 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 101.2314 -3.3680 -5.1108 -5.9416 0.7144 f_cyan_ilr_2 f_JCZ38_qlogis log_alpha log_beta 7.0229 14.9234 -0.1791 2.9811 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.416 0.000 0.0 0.000 0.0000 log_k_JCZ38 0.000 2.439 0.0 0.000 0.0000 log_k_J9Z38 0.000 0.000 1.7 0.000 0.0000 log_k_JSE76 0.000 0.000 0.0 1.856 0.0000 f_cyan_ilr_1 0.000 0.000 0.0 0.000 0.7164 f_cyan_ilr_2 0.000 0.000 0.0 0.000 0.0000 f_JCZ38_qlogis 0.000 0.000 0.0 0.000 0.0000 log_alpha 0.000 0.000 0.0 0.000 0.0000 log_beta 0.000 0.000 0.0 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis log_alpha log_beta cyan_0 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 11.57 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 18.81 0.0000 0.0000 log_alpha 0.00 0.00 0.4144 0.0000 log_beta 0.00 0.00 0.0000 0.5077 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2428 2421 -1196 Optimised parameters: est. lower upper cyan_0 101.1664 98.51265 103.8202 log_k_JCZ38 -3.3883 -4.78250 -1.9941 log_k_J9Z38 -5.3087 -5.91564 -4.7017 log_k_JSE76 -6.1313 -7.30061 -4.9619 f_cyan_ilr_1 0.7456 0.43782 1.0534 f_cyan_ilr_2 0.8181 0.24956 1.3866 f_JCZ38_qlogis 2.0467 0.61165 3.4817 log_alpha -0.2391 -0.62806 0.1499 log_beta 2.8739 2.67664 3.0711 a.1 3.4160 3.17960 3.6525 SD.cyan_0 2.4355 0.40399 4.4671 SD.log_k_JCZ38 1.5654 0.57311 2.5576 SD.log_k_J9Z38 0.4645 -0.06533 0.9943 SD.log_k_JSE76 0.9841 0.10738 1.8609 SD.f_cyan_ilr_1 0.3285 0.10546 0.5515 SD.f_cyan_ilr_2 0.2276 -0.38711 0.8424 SD.f_JCZ38_qlogis 0.8340 -0.20970 1.8777 SD.log_alpha 0.4250 0.16017 0.6898 Correlation: cyan_0 l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 log_lph log_k_JCZ38 -0.0159 log_k_J9Z38 -0.0546 0.0080 log_k_JSE76 -0.0337 0.0016 0.0074 f_cyan_ilr_1 -0.0095 0.0194 -0.1573 0.0003 f_cyan_ilr_2 -0.2733 0.0799 0.3059 0.0263 0.0125 f_JCZ38_qlogis 0.0755 -0.0783 -0.0516 0.1222 -0.1155 -0.5231 log_alpha -0.0567 0.0120 0.0351 0.0189 0.0040 0.0829 -0.0502 log_beta -0.2980 0.0461 0.1382 0.0758 0.0209 0.4079 -0.2053 0.2759 Random effects: est. lower upper SD.cyan_0 2.4355 0.40399 4.4671 SD.log_k_JCZ38 1.5654 0.57311 2.5576 SD.log_k_J9Z38 0.4645 -0.06533 0.9943 SD.log_k_JSE76 0.9841 0.10738 1.8609 SD.f_cyan_ilr_1 0.3285 0.10546 0.5515 SD.f_cyan_ilr_2 0.2276 -0.38711 0.8424 SD.f_JCZ38_qlogis 0.8340 -0.20970 1.8777 SD.log_alpha 0.4250 0.16017 0.6898 Variance model: est. lower upper a.1 3.416 3.18 3.652 Backtransformed parameters: est. lower upper cyan_0 1.012e+02 9.851e+01 103.82023 k_JCZ38 3.377e-02 8.375e-03 0.13614 k_J9Z38 4.948e-03 2.697e-03 0.00908 k_JSE76 2.174e-03 6.751e-04 0.00700 f_cyan_to_JCZ38 6.389e-01 NA NA f_cyan_to_J9Z38 2.226e-01 NA NA f_JCZ38_to_JSE76 8.856e-01 6.483e-01 0.97016 alpha 7.873e-01 5.336e-01 1.16166 beta 1.771e+01 1.454e+01 21.56509 Resulting formation fractions: ff cyan_JCZ38 0.6389 cyan_J9Z38 0.2226 cyan_sink 0.1385 JCZ38_JSE76 0.8856 JCZ38_sink 0.1144 Estimated disappearance times: DT50 DT90 DT50back cyan 25.00 312.06 93.94 JCZ38 20.53 68.19 NA J9Z38 140.07 465.32 NA JSE76 318.86 1059.22 NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:35:31 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 472.067 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 101.13294 -3.32499 -5.09097 -5.93566 0.71359 f_cyan_ilr_2 f_JCZ38_qlogis log_alpha log_beta 10.30315 14.62272 -0.09633 3.10634 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.649 0.000 0.00 0.00 0.0000 log_k_JCZ38 0.000 2.319 0.00 0.00 0.0000 log_k_J9Z38 0.000 0.000 1.73 0.00 0.0000 log_k_JSE76 0.000 0.000 0.00 1.86 0.0000 f_cyan_ilr_1 0.000 0.000 0.00 0.00 0.7183 f_cyan_ilr_2 0.000 0.000 0.00 0.00 0.0000 f_JCZ38_qlogis 0.000 0.000 0.00 0.00 0.0000 log_alpha 0.000 0.000 0.00 0.00 0.0000 log_beta 0.000 0.000 0.00 0.00 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis log_alpha log_beta cyan_0 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 12.85 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 18.54 0.0000 0.0000 log_alpha 0.00 0.00 0.3142 0.0000 log_beta 0.00 0.00 0.0000 0.7333 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2424 2416 -1193 Optimised parameters: est. lower upper cyan_0 100.65667 NA NA log_k_JCZ38 -3.45782 NA NA log_k_J9Z38 -5.23476 NA NA log_k_JSE76 -5.71827 NA NA f_cyan_ilr_1 0.68389 NA NA f_cyan_ilr_2 0.61027 NA NA f_JCZ38_qlogis 116.27482 NA NA log_alpha -0.14484 NA NA log_beta 3.03220 NA NA a.1 3.11051 NA NA b.1 0.04508 NA NA SD.log_k_JCZ38 1.39961 NA NA SD.log_k_J9Z38 0.57920 NA NA SD.log_k_JSE76 0.68364 NA NA SD.f_cyan_ilr_1 0.31477 NA NA SD.f_cyan_ilr_2 0.37716 NA NA SD.f_JCZ38_qlogis 5.52695 NA NA SD.log_alpha 0.22823 NA NA SD.log_beta 0.39161 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.3996 NA NA SD.log_k_J9Z38 0.5792 NA NA SD.log_k_JSE76 0.6836 NA NA SD.f_cyan_ilr_1 0.3148 NA NA SD.f_cyan_ilr_2 0.3772 NA NA SD.f_JCZ38_qlogis 5.5270 NA NA SD.log_alpha 0.2282 NA NA SD.log_beta 0.3916 NA NA Variance model: est. lower upper a.1 3.11051 NA NA b.1 0.04508 NA NA Backtransformed parameters: est. lower upper cyan_0 1.007e+02 NA NA k_JCZ38 3.150e-02 NA NA k_J9Z38 5.328e-03 NA NA k_JSE76 3.285e-03 NA NA f_cyan_to_JCZ38 5.980e-01 NA NA f_cyan_to_J9Z38 2.273e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA alpha 8.652e-01 NA NA beta 2.074e+01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.5980 cyan_J9Z38 0.2273 cyan_sink 0.1746 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 Estimated disappearance times: DT50 DT90 DT50back cyan 25.48 276.2 83.15 JCZ38 22.01 73.1 NA J9Z38 130.09 432.2 NA JSE76 210.98 700.9 NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:36:02 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 503.671 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.0643 -3.4008 -5.0024 -5.8612 0.6855 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 g_qlogis 1.2366 13.6901 -1.8641 -4.5063 -0.6468 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 4.466 0.000 0.000 0.000 0.0000 log_k_JCZ38 0.000 2.382 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.000 1.595 0.000 0.0000 log_k_JSE76 0.000 0.000 0.000 1.245 0.0000 f_cyan_ilr_1 0.000 0.000 0.000 0.000 0.6852 f_cyan_ilr_2 0.000 0.000 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 0.0000 log_k1 0.000 0.000 0.000 0.000 0.0000 log_k2 0.000 0.000 0.000 0.000 0.0000 g_qlogis 0.000 0.000 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 g_qlogis cyan_0 0.00 0.00 0.0000 0.0000 0.000 log_k_JCZ38 0.00 0.00 0.0000 0.0000 0.000 log_k_J9Z38 0.00 0.00 0.0000 0.0000 0.000 log_k_JSE76 0.00 0.00 0.0000 0.0000 0.000 f_cyan_ilr_1 0.00 0.00 0.0000 0.0000 0.000 f_cyan_ilr_2 1.28 0.00 0.0000 0.0000 0.000 f_JCZ38_qlogis 0.00 16.08 0.0000 0.0000 0.000 log_k1 0.00 0.00 0.9866 0.0000 0.000 log_k2 0.00 0.00 0.0000 0.5953 0.000 g_qlogis 0.00 0.00 0.0000 0.0000 1.583 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2403 2395 -1182 Optimised parameters: est. lower upper cyan_0 102.5565 NA NA log_k_JCZ38 -3.4729 NA NA log_k_J9Z38 -5.1533 NA NA log_k_JSE76 -5.6669 NA NA f_cyan_ilr_1 0.6665 NA NA f_cyan_ilr_2 0.5191 NA NA f_JCZ38_qlogis 37.0113 NA NA log_k1 -1.8497 NA NA log_k2 -4.4931 NA NA g_qlogis -0.6383 NA NA a.1 3.2397 NA NA SD.log_k_JCZ38 1.4286 NA NA SD.log_k_J9Z38 0.5312 NA NA SD.log_k_JSE76 0.6627 NA NA SD.f_cyan_ilr_1 0.3013 NA NA SD.f_cyan_ilr_2 0.2980 NA NA SD.f_JCZ38_qlogis 0.1637 NA NA SD.log_k1 0.5069 NA NA SD.log_k2 0.3828 NA NA SD.g_qlogis 0.8641 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.4286 NA NA SD.log_k_J9Z38 0.5312 NA NA SD.log_k_JSE76 0.6627 NA NA SD.f_cyan_ilr_1 0.3013 NA NA SD.f_cyan_ilr_2 0.2980 NA NA SD.f_JCZ38_qlogis 0.1637 NA NA SD.log_k1 0.5069 NA NA SD.log_k2 0.3828 NA NA SD.g_qlogis 0.8641 NA NA Variance model: est. lower upper a.1 3.24 NA NA Backtransformed parameters: est. lower upper cyan_0 1.026e+02 NA NA k_JCZ38 3.103e-02 NA NA k_J9Z38 5.780e-03 NA NA k_JSE76 3.459e-03 NA NA f_cyan_to_JCZ38 5.813e-01 NA NA f_cyan_to_J9Z38 2.265e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA k1 1.573e-01 NA NA k2 1.119e-02 NA NA g 3.456e-01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.5813 cyan_J9Z38 0.2265 cyan_sink 0.1922 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 25.23 167.94 50.55 4.407 61.97 JCZ38 22.34 74.22 NA NA NA J9Z38 119.92 398.36 NA NA NA JSE76 200.41 665.76 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:38:54 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 674.81 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 101.3964 -3.3626 -4.9792 -5.8727 0.6814 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 g_qlogis 6.8713 13.6901 -1.9222 -4.5035 -0.7172 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.317 0.000 0.000 0.000 0.0000 log_k_JCZ38 0.000 2.272 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.000 1.633 0.000 0.0000 log_k_JSE76 0.000 0.000 0.000 1.271 0.0000 f_cyan_ilr_1 0.000 0.000 0.000 0.000 0.6839 f_cyan_ilr_2 0.000 0.000 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 0.0000 log_k1 0.000 0.000 0.000 0.000 0.0000 log_k2 0.000 0.000 0.000 0.000 0.0000 g_qlogis 0.000 0.000 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 g_qlogis cyan_0 0.00 0.00 0.0000 0.0000 0.000 log_k_JCZ38 0.00 0.00 0.0000 0.0000 0.000 log_k_J9Z38 0.00 0.00 0.0000 0.0000 0.000 log_k_JSE76 0.00 0.00 0.0000 0.0000 0.000 f_cyan_ilr_1 0.00 0.00 0.0000 0.0000 0.000 f_cyan_ilr_2 11.95 0.00 0.0000 0.0000 0.000 f_JCZ38_qlogis 0.00 16.08 0.0000 0.0000 0.000 log_k1 0.00 0.00 0.9496 0.0000 0.000 log_k2 0.00 0.00 0.0000 0.5846 0.000 g_qlogis 0.00 0.00 0.0000 0.0000 1.719 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2398 2390 -1179 Optimised parameters: est. lower upper cyan_0 100.69709 NA NA log_k_JCZ38 -3.46669 NA NA log_k_J9Z38 -5.05076 NA NA log_k_JSE76 -5.55558 NA NA f_cyan_ilr_1 0.66045 NA NA f_cyan_ilr_2 0.84275 NA NA f_JCZ38_qlogis 64.22404 NA NA log_k1 -2.17715 NA NA log_k2 -4.55002 NA NA g_qlogis -0.55920 NA NA a.1 2.95785 NA NA b.1 0.04456 NA NA SD.log_k_JCZ38 1.39881 NA NA SD.log_k_J9Z38 0.67788 NA NA SD.log_k_JSE76 0.52603 NA NA SD.f_cyan_ilr_1 0.32490 NA NA SD.f_cyan_ilr_2 0.53923 NA NA SD.f_JCZ38_qlogis 2.75576 NA NA SD.log_k2 0.30694 NA NA SD.g_qlogis 0.83619 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.3988 NA NA SD.log_k_J9Z38 0.6779 NA NA SD.log_k_JSE76 0.5260 NA NA SD.f_cyan_ilr_1 0.3249 NA NA SD.f_cyan_ilr_2 0.5392 NA NA SD.f_JCZ38_qlogis 2.7558 NA NA SD.log_k2 0.3069 NA NA SD.g_qlogis 0.8362 NA NA Variance model: est. lower upper a.1 2.95785 NA NA b.1 0.04456 NA NA Backtransformed parameters: est. lower upper cyan_0 1.007e+02 NA NA k_JCZ38 3.122e-02 NA NA k_J9Z38 6.404e-03 NA NA k_JSE76 3.866e-03 NA NA f_cyan_to_JCZ38 6.187e-01 NA NA f_cyan_to_J9Z38 2.431e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA k1 1.134e-01 NA NA k2 1.057e-02 NA NA g 3.637e-01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.6187 cyan_J9Z38 0.2431 cyan_sink 0.1382 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 26.35 175.12 52.72 6.114 65.6 JCZ38 22.20 73.75 NA NA NA J9Z38 108.23 359.53 NA NA NA JSE76 179.30 595.62 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:36:26 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 527.639 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 102.0643 -2.8987 -2.7077 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.4717 -3.4008 -5.0024 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -5.8613 0.6855 1.2366 f_JCZ38_qlogis 13.7395 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 4.466 0.0000 0.000 log_k_cyan_free 0.000 0.6158 0.000 log_k_cyan_free_bound 0.000 0.0000 1.463 log_k_cyan_bound_free 0.000 0.0000 0.000 log_k_JCZ38 0.000 0.0000 0.000 log_k_J9Z38 0.000 0.0000 0.000 log_k_JSE76 0.000 0.0000 0.000 f_cyan_ilr_1 0.000 0.0000 0.000 f_cyan_ilr_2 0.000 0.0000 0.000 f_JCZ38_qlogis 0.000 0.0000 0.000 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.000 0.000 0.000 log_k_cyan_free 0.000 0.000 0.000 0.000 log_k_cyan_free_bound 0.000 0.000 0.000 0.000 log_k_cyan_bound_free 1.058 0.000 0.000 0.000 log_k_JCZ38 0.000 2.382 0.000 0.000 log_k_J9Z38 0.000 0.000 1.595 0.000 log_k_JSE76 0.000 0.000 0.000 1.245 f_cyan_ilr_1 0.000 0.000 0.000 0.000 f_cyan_ilr_2 0.000 0.000 0.000 0.000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis cyan_free_0 0.0000 0.00 0.00 log_k_cyan_free 0.0000 0.00 0.00 log_k_cyan_free_bound 0.0000 0.00 0.00 log_k_cyan_bound_free 0.0000 0.00 0.00 log_k_JCZ38 0.0000 0.00 0.00 log_k_J9Z38 0.0000 0.00 0.00 log_k_JSE76 0.0000 0.00 0.00 f_cyan_ilr_1 0.6852 0.00 0.00 f_cyan_ilr_2 0.0000 1.28 0.00 f_JCZ38_qlogis 0.0000 0.00 16.13 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2401 2394 -1181 Optimised parameters: est. lower upper cyan_free_0 102.8136 NA NA log_k_cyan_free -2.7935 NA NA log_k_cyan_free_bound -2.5440 NA NA log_k_cyan_bound_free -3.4303 NA NA log_k_JCZ38 -3.5010 NA NA log_k_J9Z38 -5.1226 NA NA log_k_JSE76 -5.6314 NA NA f_cyan_ilr_1 0.6609 NA NA f_cyan_ilr_2 0.5085 NA NA f_JCZ38_qlogis 44.0153 NA NA a.1 3.2318 NA NA SD.log_k_cyan_free 0.3211 NA NA SD.log_k_cyan_free_bound 0.8408 NA NA SD.log_k_cyan_bound_free 0.5724 NA NA SD.log_k_JCZ38 1.4925 NA NA SD.log_k_J9Z38 0.5816 NA NA SD.log_k_JSE76 0.6037 NA NA SD.f_cyan_ilr_1 0.3115 NA NA SD.f_cyan_ilr_2 0.3436 NA NA SD.f_JCZ38_qlogis 4.8937 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_cyan_free 0.3211 NA NA SD.log_k_cyan_free_bound 0.8408 NA NA SD.log_k_cyan_bound_free 0.5724 NA NA SD.log_k_JCZ38 1.4925 NA NA SD.log_k_J9Z38 0.5816 NA NA SD.log_k_JSE76 0.6037 NA NA SD.f_cyan_ilr_1 0.3115 NA NA SD.f_cyan_ilr_2 0.3436 NA NA SD.f_JCZ38_qlogis 4.8937 NA NA Variance model: est. lower upper a.1 3.232 NA NA Backtransformed parameters: est. lower upper cyan_free_0 1.028e+02 NA NA k_cyan_free 6.120e-02 NA NA k_cyan_free_bound 7.855e-02 NA NA k_cyan_bound_free 3.238e-02 NA NA k_JCZ38 3.017e-02 NA NA k_J9Z38 5.961e-03 NA NA k_JSE76 3.584e-03 NA NA f_cyan_free_to_JCZ38 5.784e-01 NA NA f_cyan_free_to_J9Z38 2.271e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.15973 0.01241 0.33124 Resulting formation fractions: ff cyan_free_JCZ38 0.5784 cyan_free_J9Z38 0.2271 cyan_free_sink 0.1945 cyan_free 1.0000 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 24.51 153.18 46.11 4.34 55.87 JCZ38 22.98 76.33 NA NA NA J9Z38 116.28 386.29 NA NA NA JSE76 193.42 642.53 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:38:31 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 651.67 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 101.3964 -2.9881 -2.7949 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.4376 -3.3626 -4.9792 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -5.8727 0.6814 6.7399 f_JCZ38_qlogis 13.7395 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 5.317 0.0000 0.000 log_k_cyan_free 0.000 0.7301 0.000 log_k_cyan_free_bound 0.000 0.0000 1.384 log_k_cyan_bound_free 0.000 0.0000 0.000 log_k_JCZ38 0.000 0.0000 0.000 log_k_J9Z38 0.000 0.0000 0.000 log_k_JSE76 0.000 0.0000 0.000 f_cyan_ilr_1 0.000 0.0000 0.000 f_cyan_ilr_2 0.000 0.0000 0.000 f_JCZ38_qlogis 0.000 0.0000 0.000 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.000 0.000 0.000 log_k_cyan_free 0.000 0.000 0.000 0.000 log_k_cyan_free_bound 0.000 0.000 0.000 0.000 log_k_cyan_bound_free 1.109 0.000 0.000 0.000 log_k_JCZ38 0.000 2.272 0.000 0.000 log_k_J9Z38 0.000 0.000 1.633 0.000 log_k_JSE76 0.000 0.000 0.000 1.271 f_cyan_ilr_1 0.000 0.000 0.000 0.000 f_cyan_ilr_2 0.000 0.000 0.000 0.000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis cyan_free_0 0.0000 0.00 0.00 log_k_cyan_free 0.0000 0.00 0.00 log_k_cyan_free_bound 0.0000 0.00 0.00 log_k_cyan_bound_free 0.0000 0.00 0.00 log_k_JCZ38 0.0000 0.00 0.00 log_k_J9Z38 0.0000 0.00 0.00 log_k_JSE76 0.0000 0.00 0.00 f_cyan_ilr_1 0.6838 0.00 0.00 f_cyan_ilr_2 0.0000 11.69 0.00 f_JCZ38_qlogis 0.0000 0.00 16.13 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2400 2392 -1180 Optimised parameters: est. lower upper cyan_free_0 100.56004 NA NA log_k_cyan_free -3.12657 NA NA log_k_cyan_free_bound -3.16825 NA NA log_k_cyan_bound_free -3.66003 NA NA log_k_JCZ38 -3.47278 NA NA log_k_J9Z38 -5.06823 NA NA log_k_JSE76 -5.54327 NA NA f_cyan_ilr_1 0.66631 NA NA f_cyan_ilr_2 0.82898 NA NA f_JCZ38_qlogis 38.31115 NA NA a.1 2.98352 NA NA b.1 0.04388 NA NA SD.log_k_cyan_free 0.49145 NA NA SD.log_k_cyan_bound_free 0.27347 NA NA SD.log_k_JCZ38 1.41193 NA NA SD.log_k_J9Z38 0.66073 NA NA SD.log_k_JSE76 0.55885 NA NA SD.f_cyan_ilr_1 0.33020 NA NA SD.f_cyan_ilr_2 0.51367 NA NA SD.f_JCZ38_qlogis 5.52122 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_cyan_free 0.4914 NA NA SD.log_k_cyan_bound_free 0.2735 NA NA SD.log_k_JCZ38 1.4119 NA NA SD.log_k_J9Z38 0.6607 NA NA SD.log_k_JSE76 0.5589 NA NA SD.f_cyan_ilr_1 0.3302 NA NA SD.f_cyan_ilr_2 0.5137 NA NA SD.f_JCZ38_qlogis 5.5212 NA NA Variance model: est. lower upper a.1 2.98352 NA NA b.1 0.04388 NA NA Backtransformed parameters: est. lower upper cyan_free_0 1.006e+02 NA NA k_cyan_free 4.387e-02 NA NA k_cyan_free_bound 4.208e-02 NA NA k_cyan_bound_free 2.573e-02 NA NA k_JCZ38 3.103e-02 NA NA k_J9Z38 6.294e-03 NA NA k_JSE76 3.914e-03 NA NA f_cyan_free_to_JCZ38 6.188e-01 NA NA f_cyan_free_to_J9Z38 2.412e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.10044 0.01124 0.36580 Resulting formation fractions: ff cyan_free_JCZ38 0.6188 cyan_free_J9Z38 0.2412 cyan_free_sink 0.1400 cyan_free 1.0000 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 26.05 164.4 49.48 6.901 61.67 JCZ38 22.34 74.2 NA NA NA J9Z38 110.14 365.9 NA NA NA JSE76 177.11 588.3 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:36:13 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ifelse(time <= tb, k1, k2) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ifelse(time <= tb, k1, k2) * cyan - k_JCZ38 * JCZ38 d_J9Z38/dt = + f_cyan_to_J9Z38 * ifelse(time <= tb, k1, k2) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 514.353 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.8845 -3.4495 -4.9355 -5.6040 0.6468 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 log_tb 1.2396 9.7220 -2.9079 -4.1810 1.7813 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.406 0.00 0.00 0.000 0.0000 log_k_JCZ38 0.000 2.33 0.00 0.000 0.0000 log_k_J9Z38 0.000 0.00 1.59 0.000 0.0000 log_k_JSE76 0.000 0.00 0.00 1.013 0.0000 f_cyan_ilr_1 0.000 0.00 0.00 0.000 0.6367 f_cyan_ilr_2 0.000 0.00 0.00 0.000 0.0000 f_JCZ38_qlogis 0.000 0.00 0.00 0.000 0.0000 log_k1 0.000 0.00 0.00 0.000 0.0000 log_k2 0.000 0.00 0.00 0.000 0.0000 log_tb 0.000 0.00 0.00 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis log_k1 log_k2 log_tb cyan_0 0.000 0.00 0.0000 0.0000 0.0000 log_k_JCZ38 0.000 0.00 0.0000 0.0000 0.0000 log_k_J9Z38 0.000 0.00 0.0000 0.0000 0.0000 log_k_JSE76 0.000 0.00 0.0000 0.0000 0.0000 f_cyan_ilr_1 0.000 0.00 0.0000 0.0000 0.0000 f_cyan_ilr_2 2.038 0.00 0.0000 0.0000 0.0000 f_JCZ38_qlogis 0.000 10.33 0.0000 0.0000 0.0000 log_k1 0.000 0.00 0.7006 0.0000 0.0000 log_k2 0.000 0.00 0.0000 0.8928 0.0000 log_tb 0.000 0.00 0.0000 0.0000 0.6773 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2427 2419 -1194 Optimised parameters: est. lower upper cyan_0 101.9660 1.005e+02 1.035e+02 log_k_JCZ38 -3.4698 -4.716e+00 -2.224e+00 log_k_J9Z38 -5.0947 -5.740e+00 -4.450e+00 log_k_JSE76 -5.5977 -6.321e+00 -4.875e+00 f_cyan_ilr_1 0.6595 3.734e-01 9.456e-01 f_cyan_ilr_2 0.5905 1.664e-01 1.015e+00 f_JCZ38_qlogis 25.8627 -4.224e+05 4.225e+05 log_k1 -3.0884 -3.453e+00 -2.723e+00 log_k2 -4.3877 -4.778e+00 -3.998e+00 log_tb 2.3057 1.715e+00 2.896e+00 a.1 3.3228 NA NA SD.log_k_JCZ38 1.4071 NA NA SD.log_k_J9Z38 0.5774 NA NA SD.log_k_JSE76 0.6214 NA NA SD.f_cyan_ilr_1 0.3058 NA NA SD.f_cyan_ilr_2 0.3470 NA NA SD.f_JCZ38_qlogis 0.0644 NA NA SD.log_k1 0.3994 NA NA SD.log_k2 0.4373 NA NA SD.log_tb 0.6419 NA NA Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.4071 NA NA SD.log_k_J9Z38 0.5774 NA NA SD.log_k_JSE76 0.6214 NA NA SD.f_cyan_ilr_1 0.3058 NA NA SD.f_cyan_ilr_2 0.3470 NA NA SD.f_JCZ38_qlogis 0.0644 NA NA SD.log_k1 0.3994 NA NA SD.log_k2 0.4373 NA NA SD.log_tb 0.6419 NA NA Variance model: est. lower upper a.1 3.323 NA NA Backtransformed parameters: est. lower upper cyan_0 1.020e+02 1.005e+02 1.035e+02 k_JCZ38 3.112e-02 8.951e-03 1.082e-01 k_J9Z38 6.129e-03 3.216e-03 1.168e-02 k_JSE76 3.706e-03 1.798e-03 7.639e-03 f_cyan_to_JCZ38 5.890e-01 NA NA f_cyan_to_J9Z38 2.318e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 0.000e+00 1.000e+00 k1 4.558e-02 3.164e-02 6.565e-02 k2 1.243e-02 8.417e-03 1.835e-02 tb 1.003e+01 5.557e+00 1.811e+01 Resulting formation fractions: ff cyan_JCZ38 5.890e-01 cyan_J9Z38 2.318e-01 cyan_sink 1.793e-01 JCZ38_JSE76 1.000e+00 JCZ38_sink 5.861e-12 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 29.02 158.51 47.72 15.21 55.77 JCZ38 22.27 73.98 NA NA NA J9Z38 113.09 375.69 NA NA NA JSE76 187.01 621.23 NA NA NA"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"pathway-2","dir":"Articles > Prebuilt","previous_headings":"Appendix > Hierarchical fit listings","what":"Pathway 2","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"Hierarchical FOMC path 2 fit two-component error Hierarchical DFOP path 2 fit constant variance Hierarchical DFOP path 2 fit two-component error Hierarchical SFORB path 2 fit constant variance Hierarchical SFORB path 2 fit two-component error","code":"saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:47:30 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 505.533 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.4477 -1.8631 -5.1087 -2.5114 0.6826 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_alpha log_beta 4.7944 15.9616 13.1566 -0.1564 2.9781 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 7.701 0.000 0.000 0.000 0.0000 log_k_JCZ38 0.000 1.448 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.000 1.724 0.000 0.0000 log_k_JSE76 0.000 0.000 0.000 3.659 0.0000 f_cyan_ilr_1 0.000 0.000 0.000 0.000 0.6356 f_cyan_ilr_2 0.000 0.000 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 0.0000 f_JSE76_qlogis 0.000 0.000 0.000 0.000 0.0000 log_alpha 0.000 0.000 0.000 0.000 0.0000 log_beta 0.000 0.000 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_alpha log_beta cyan_0 0.00 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 10.32 0.00 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 12.23 0.00 0.0000 0.0000 f_JSE76_qlogis 0.00 0.00 14.99 0.0000 0.0000 log_alpha 0.00 0.00 0.00 0.3924 0.0000 log_beta 0.00 0.00 0.00 0.0000 0.5639 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2249 2241 -1104 Optimised parameters: est. lower upper cyan_0 101.55265 9.920e+01 103.9059 log_k_JCZ38 -2.32302 -2.832e+00 -1.8142 log_k_J9Z38 -5.13082 -5.942e+00 -4.3199 log_k_JSE76 -3.01756 -4.262e+00 -1.7736 f_cyan_ilr_1 0.70850 3.657e-01 1.0513 f_cyan_ilr_2 0.95775 2.612e-01 1.6543 f_JCZ38_qlogis 3.86105 9.248e-01 6.7973 f_JSE76_qlogis 7.51583 -1.120e+02 127.0392 log_alpha -0.15308 -4.508e-01 0.1446 log_beta 2.99165 2.711e+00 3.2720 a.1 2.04034 1.843e+00 2.2382 b.1 0.06924 5.749e-02 0.0810 SD.log_k_JCZ38 0.50818 1.390e-01 0.8774 SD.log_k_J9Z38 0.86597 2.652e-01 1.4667 SD.log_k_JSE76 1.38092 4.864e-01 2.2754 SD.f_cyan_ilr_1 0.38204 1.354e-01 0.6286 SD.f_cyan_ilr_2 0.55129 7.198e-02 1.0306 SD.f_JCZ38_qlogis 1.88457 1.711e-02 3.7520 SD.f_JSE76_qlogis 2.64018 -2.450e+03 2454.9447 SD.log_alpha 0.31860 1.047e-01 0.5325 SD.log_beta 0.24195 1.273e-02 0.4712 Correlation: cyan_0 l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 f_JSE76 log_k_JCZ38 -0.0235 log_k_J9Z38 -0.0442 0.0047 log_k_JSE76 -0.0023 0.0966 0.0006 f_cyan_ilr_1 -0.0032 0.0070 -0.0536 -0.0001 f_cyan_ilr_2 -0.5189 0.0452 0.1152 0.0013 -0.0304 f_JCZ38_qlogis 0.1088 -0.0848 -0.0240 0.0040 -0.0384 -0.2303 f_JSE76_qlogis -0.0545 0.1315 0.0195 0.0020 0.0252 0.1737 -0.5939 log_alpha -0.0445 0.0056 0.0261 0.0019 -0.0055 0.0586 -0.0239 -0.0284 log_beta -0.2388 0.0163 0.0566 0.0040 -0.0078 0.2183 -0.0714 -0.0332 log_lph log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_alpha log_beta 0.2135 Random effects: est. lower upper SD.log_k_JCZ38 0.5082 1.390e-01 0.8774 SD.log_k_J9Z38 0.8660 2.652e-01 1.4667 SD.log_k_JSE76 1.3809 4.864e-01 2.2754 SD.f_cyan_ilr_1 0.3820 1.354e-01 0.6286 SD.f_cyan_ilr_2 0.5513 7.198e-02 1.0306 SD.f_JCZ38_qlogis 1.8846 1.711e-02 3.7520 SD.f_JSE76_qlogis 2.6402 -2.450e+03 2454.9447 SD.log_alpha 0.3186 1.047e-01 0.5325 SD.log_beta 0.2420 1.273e-02 0.4712 Variance model: est. lower upper a.1 2.04034 1.84252 2.238 b.1 0.06924 0.05749 0.081 Backtransformed parameters: est. lower upper cyan_0 1.016e+02 9.920e+01 103.9059 k_JCZ38 9.798e-02 5.890e-02 0.1630 k_J9Z38 5.912e-03 2.627e-03 0.0133 k_JSE76 4.892e-02 1.410e-02 0.1697 f_cyan_to_JCZ38 6.432e-01 NA NA f_cyan_to_J9Z38 2.362e-01 NA NA f_JCZ38_to_JSE76 9.794e-01 7.160e-01 0.9989 f_JSE76_to_JCZ38 9.995e-01 2.268e-49 1.0000 alpha 8.581e-01 6.371e-01 1.1556 beta 1.992e+01 1.505e+01 26.3646 Resulting formation fractions: ff cyan_JCZ38 0.6432301 cyan_J9Z38 0.2361657 cyan_sink 0.1206042 JCZ38_JSE76 0.9793879 JCZ38_sink 0.0206121 JSE76_JCZ38 0.9994559 JSE76_sink 0.0005441 Estimated disappearance times: DT50 DT90 DT50back cyan 24.759 271.61 81.76 JCZ38 7.075 23.50 NA J9Z38 117.249 389.49 NA JSE76 14.169 47.07 NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:48:20 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 555.413 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.4380 -2.3107 -5.3123 -3.7120 0.6757 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 1.1439 13.1194 12.3492 -1.9317 -4.4557 g_qlogis -0.5644 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 4.591 0.0000 0.000 0.0 0.0000 log_k_JCZ38 0.000 0.7966 0.000 0.0 0.0000 log_k_J9Z38 0.000 0.0000 1.561 0.0 0.0000 log_k_JSE76 0.000 0.0000 0.000 0.8 0.0000 f_cyan_ilr_1 0.000 0.0000 0.000 0.0 0.6349 f_cyan_ilr_2 0.000 0.0000 0.000 0.0 0.0000 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0 0.0000 f_JSE76_qlogis 0.000 0.0000 0.000 0.0 0.0000 log_k1 0.000 0.0000 0.000 0.0 0.0000 log_k2 0.000 0.0000 0.000 0.0 0.0000 g_qlogis 0.000 0.0000 0.000 0.0 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 cyan_0 0.000 0.00 0.00 0.000 0.0000 log_k_JCZ38 0.000 0.00 0.00 0.000 0.0000 log_k_J9Z38 0.000 0.00 0.00 0.000 0.0000 log_k_JSE76 0.000 0.00 0.00 0.000 0.0000 f_cyan_ilr_1 0.000 0.00 0.00 0.000 0.0000 f_cyan_ilr_2 1.797 0.00 0.00 0.000 0.0000 f_JCZ38_qlogis 0.000 13.86 0.00 0.000 0.0000 f_JSE76_qlogis 0.000 0.00 13.91 0.000 0.0000 log_k1 0.000 0.00 0.00 1.106 0.0000 log_k2 0.000 0.00 0.00 0.000 0.6141 g_qlogis 0.000 0.00 0.00 0.000 0.0000 g_qlogis cyan_0 0.000 log_k_JCZ38 0.000 log_k_J9Z38 0.000 log_k_JSE76 0.000 f_cyan_ilr_1 0.000 f_cyan_ilr_2 0.000 f_JCZ38_qlogis 0.000 f_JSE76_qlogis 0.000 log_k1 0.000 log_k2 0.000 g_qlogis 1.595 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2288 2280 -1122 Optimised parameters: est. lower upper cyan_0 102.7204 1.014e+02 1.040e+02 log_k_JCZ38 -2.8925 -4.044e+00 -1.741e+00 log_k_J9Z38 -5.1430 -5.828e+00 -4.457e+00 log_k_JSE76 -3.5577 -4.174e+00 -2.941e+00 f_cyan_ilr_1 0.6929 3.788e-01 1.007e+00 f_cyan_ilr_2 0.6066 5.342e-02 1.160e+00 f_JCZ38_qlogis 9.8071 -2.819e+03 2.838e+03 f_JSE76_qlogis 2.2229 5.684e-01 3.877e+00 log_k1 -1.9339 -2.609e+00 -1.258e+00 log_k2 -4.4709 -4.935e+00 -4.007e+00 g_qlogis -0.4987 -1.373e+00 3.757e-01 a.1 2.7368 2.545e+00 2.928e+00 SD.log_k_JCZ38 1.2747 4.577e-01 2.092e+00 SD.log_k_J9Z38 0.6758 1.418e-01 1.210e+00 SD.log_k_JSE76 0.5869 1.169e-01 1.057e+00 SD.f_cyan_ilr_1 0.3392 1.161e-01 5.622e-01 SD.f_cyan_ilr_2 0.4200 8.501e-02 7.550e-01 SD.f_JCZ38_qlogis 0.8511 -1.137e+06 1.137e+06 SD.f_JSE76_qlogis 0.3767 -5.238e-01 1.277e+00 SD.log_k1 0.7475 2.601e-01 1.235e+00 SD.log_k2 0.5179 1.837e-01 8.521e-01 SD.g_qlogis 0.9817 3.553e-01 1.608e+00 Correlation: cyan_0 l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 f_JSE76 log_k_JCZ38 -0.0351 log_k_J9Z38 -0.0541 0.0043 log_k_JSE76 -0.0078 0.0900 -0.0014 f_cyan_ilr_1 -0.0249 0.0268 -0.0962 0.0000 f_cyan_ilr_2 -0.3560 0.0848 0.1545 -0.0022 0.0463 f_JCZ38_qlogis 0.2005 -0.1226 -0.0347 0.0514 -0.1840 -0.5906 f_JSE76_qlogis -0.1638 0.1307 0.0266 0.0001 0.1645 0.5181 -0.9297 log_k1 0.0881 -0.0071 0.0005 -0.0070 -0.0064 -0.0346 0.0316 -0.0341 log_k2 0.0238 -0.0003 0.0082 -0.0022 -0.0017 -0.0017 -0.0002 -0.0076 g_qlogis 0.0198 -0.0002 -0.0109 0.0034 0.0017 -0.0176 0.0044 0.0051 log_k1 log_k2 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 0.0276 g_qlogis -0.0283 -0.0309 Random effects: est. lower upper SD.log_k_JCZ38 1.2747 4.577e-01 2.092e+00 SD.log_k_J9Z38 0.6758 1.418e-01 1.210e+00 SD.log_k_JSE76 0.5869 1.169e-01 1.057e+00 SD.f_cyan_ilr_1 0.3392 1.161e-01 5.622e-01 SD.f_cyan_ilr_2 0.4200 8.501e-02 7.550e-01 SD.f_JCZ38_qlogis 0.8511 -1.137e+06 1.137e+06 SD.f_JSE76_qlogis 0.3767 -5.238e-01 1.277e+00 SD.log_k1 0.7475 2.601e-01 1.235e+00 SD.log_k2 0.5179 1.837e-01 8.521e-01 SD.g_qlogis 0.9817 3.553e-01 1.608e+00 Variance model: est. lower upper a.1 2.737 2.545 2.928 Backtransformed parameters: est. lower upper cyan_0 102.72037 1.014e+02 104.00464 k_JCZ38 0.05544 1.752e-02 0.17539 k_J9Z38 0.00584 2.942e-03 0.01159 k_JSE76 0.02850 1.539e-02 0.05279 f_cyan_to_JCZ38 0.59995 NA NA f_cyan_to_J9Z38 0.22519 NA NA f_JCZ38_to_JSE76 0.99994 0.000e+00 1.00000 f_JSE76_to_JCZ38 0.90229 6.384e-01 0.97971 k1 0.14459 7.357e-02 0.28414 k2 0.01144 7.192e-03 0.01819 g 0.37784 2.021e-01 0.59284 Resulting formation fractions: ff cyan_JCZ38 5.999e-01 cyan_J9Z38 2.252e-01 cyan_sink 1.749e-01 JCZ38_JSE76 9.999e-01 JCZ38_sink 5.506e-05 JSE76_JCZ38 9.023e-01 JSE76_sink 9.771e-02 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 21.93 159.83 48.11 4.794 60.6 JCZ38 12.50 41.53 NA NA NA J9Z38 118.69 394.27 NA NA NA JSE76 24.32 80.78 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:51:02 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 717.675 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 101.7393 -1.4493 -5.0118 -2.1269 0.6720 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 7.3362 13.4423 13.2659 -2.0061 -4.5527 g_qlogis -0.5806 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.604 0.00 0.000 0.000 0.0000 log_k_JCZ38 0.000 2.77 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.00 1.662 0.000 0.0000 log_k_JSE76 0.000 0.00 0.000 5.021 0.0000 f_cyan_ilr_1 0.000 0.00 0.000 0.000 0.6519 f_cyan_ilr_2 0.000 0.00 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.00 0.000 0.000 0.0000 f_JSE76_qlogis 0.000 0.00 0.000 0.000 0.0000 log_k1 0.000 0.00 0.000 0.000 0.0000 log_k2 0.000 0.00 0.000 0.000 0.0000 g_qlogis 0.000 0.00 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 cyan_0 0.00 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 13.37 0.00 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 14.21 0.00 0.0000 0.0000 f_JSE76_qlogis 0.00 0.00 14.58 0.0000 0.0000 log_k1 0.00 0.00 0.00 0.8453 0.0000 log_k2 0.00 0.00 0.00 0.0000 0.5969 g_qlogis 0.00 0.00 0.00 0.0000 0.0000 g_qlogis cyan_0 0.00 log_k_JCZ38 0.00 log_k_J9Z38 0.00 log_k_JSE76 0.00 f_cyan_ilr_1 0.00 f_cyan_ilr_2 0.00 f_JCZ38_qlogis 0.00 f_JSE76_qlogis 0.00 log_k1 0.00 log_k2 0.00 g_qlogis 1.69 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2234 2226 -1095 Optimised parameters: est. lower upper cyan_0 101.25496 99.14662 103.36331 log_k_JCZ38 -2.55593 -3.32972 -1.78215 log_k_J9Z38 -5.07103 -5.85423 -4.28783 log_k_JSE76 -3.25468 -4.17577 -2.33360 f_cyan_ilr_1 0.70139 0.35924 1.04355 f_cyan_ilr_2 1.07712 0.17789 1.97636 f_JCZ38_qlogis 3.57483 0.05990 7.08976 f_JSE76_qlogis 4.54884 -7.25628 16.35395 log_k1 -2.38201 -2.51639 -2.24763 log_k2 -4.66741 -4.91865 -4.41617 g_qlogis -0.28446 -1.14192 0.57300 a.1 2.05925 1.86481 2.25369 b.1 0.06172 0.05062 0.07282 SD.log_k_JCZ38 0.81137 0.25296 1.36977 SD.log_k_J9Z38 0.83542 0.25395 1.41689 SD.log_k_JSE76 0.97903 0.30100 1.65707 SD.f_cyan_ilr_1 0.37878 0.13374 0.62382 SD.f_cyan_ilr_2 0.67274 0.10102 1.24446 SD.f_JCZ38_qlogis 1.35327 -0.42359 3.13012 SD.f_JSE76_qlogis 1.43956 -19.14972 22.02884 SD.log_k2 0.25329 0.07521 0.43138 SD.g_qlogis 0.95167 0.35149 1.55184 Correlation: cyan_0 l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 f_JSE76 log_k_JCZ38 -0.0265 log_k_J9Z38 -0.0392 0.0024 log_k_JSE76 0.0011 0.1220 -0.0016 f_cyan_ilr_1 -0.0161 0.0217 -0.0552 0.0034 f_cyan_ilr_2 -0.4718 0.0829 0.1102 0.0042 0.0095 f_JCZ38_qlogis 0.1609 -0.1318 -0.0277 0.0081 -0.1040 -0.4559 f_JSE76_qlogis -0.1289 0.1494 0.0219 0.0012 0.1004 0.4309 -0.8543 log_k1 0.2618 -0.0739 -0.0167 -0.0148 -0.0444 -0.2768 0.3518 -0.3818 log_k2 0.0603 -0.0217 0.0174 -0.0058 -0.0197 -0.0533 0.0923 -0.1281 g_qlogis 0.0362 0.0115 -0.0111 0.0040 0.0095 -0.0116 -0.0439 0.0651 log_k1 log_k2 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 0.3269 g_qlogis -0.1656 -0.0928 Random effects: est. lower upper SD.log_k_JCZ38 0.8114 0.25296 1.3698 SD.log_k_J9Z38 0.8354 0.25395 1.4169 SD.log_k_JSE76 0.9790 0.30100 1.6571 SD.f_cyan_ilr_1 0.3788 0.13374 0.6238 SD.f_cyan_ilr_2 0.6727 0.10102 1.2445 SD.f_JCZ38_qlogis 1.3533 -0.42359 3.1301 SD.f_JSE76_qlogis 1.4396 -19.14972 22.0288 SD.log_k2 0.2533 0.07521 0.4314 SD.g_qlogis 0.9517 0.35149 1.5518 Variance model: est. lower upper a.1 2.05925 1.86481 2.25369 b.1 0.06172 0.05062 0.07282 Backtransformed parameters: est. lower upper cyan_0 1.013e+02 9.915e+01 103.36331 k_JCZ38 7.762e-02 3.580e-02 0.16828 k_J9Z38 6.276e-03 2.868e-03 0.01373 k_JSE76 3.859e-02 1.536e-02 0.09695 f_cyan_to_JCZ38 6.520e-01 NA NA f_cyan_to_J9Z38 2.418e-01 NA NA f_JCZ38_to_JSE76 9.727e-01 5.150e-01 0.99917 f_JSE76_to_JCZ38 9.895e-01 7.052e-04 1.00000 k1 9.236e-02 8.075e-02 0.10565 k2 9.397e-03 7.309e-03 0.01208 g 4.294e-01 2.420e-01 0.63945 Resulting formation fractions: ff cyan_JCZ38 0.65203 cyan_J9Z38 0.24181 cyan_sink 0.10616 JCZ38_JSE76 0.97274 JCZ38_sink 0.02726 JSE76_JCZ38 0.98953 JSE76_sink 0.01047 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 24.26 185.34 55.79 7.504 73.77 JCZ38 8.93 29.66 NA NA NA J9Z38 110.45 366.89 NA NA NA JSE76 17.96 59.66 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:48:15 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 550.772 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 102.4395 -2.7673 -2.8942 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.6201 -2.3107 -5.3123 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -3.7120 0.6754 1.1448 f_JCZ38_qlogis f_JSE76_qlogis 14.8408 15.4734 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 4.589 0.0000 0.00 log_k_cyan_free 0.000 0.4849 0.00 log_k_cyan_free_bound 0.000 0.0000 1.62 log_k_cyan_bound_free 0.000 0.0000 0.00 log_k_JCZ38 0.000 0.0000 0.00 log_k_J9Z38 0.000 0.0000 0.00 log_k_JSE76 0.000 0.0000 0.00 f_cyan_ilr_1 0.000 0.0000 0.00 f_cyan_ilr_2 0.000 0.0000 0.00 f_JCZ38_qlogis 0.000 0.0000 0.00 f_JSE76_qlogis 0.000 0.0000 0.00 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.0000 0.000 0.0 log_k_cyan_free 0.000 0.0000 0.000 0.0 log_k_cyan_free_bound 0.000 0.0000 0.000 0.0 log_k_cyan_bound_free 1.197 0.0000 0.000 0.0 log_k_JCZ38 0.000 0.7966 0.000 0.0 log_k_J9Z38 0.000 0.0000 1.561 0.0 log_k_JSE76 0.000 0.0000 0.000 0.8 f_cyan_ilr_1 0.000 0.0000 0.000 0.0 f_cyan_ilr_2 0.000 0.0000 0.000 0.0 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0 f_JSE76_qlogis 0.000 0.0000 0.000 0.0 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis cyan_free_0 0.0000 0.000 0.0 0.00 log_k_cyan_free 0.0000 0.000 0.0 0.00 log_k_cyan_free_bound 0.0000 0.000 0.0 0.00 log_k_cyan_bound_free 0.0000 0.000 0.0 0.00 log_k_JCZ38 0.0000 0.000 0.0 0.00 log_k_J9Z38 0.0000 0.000 0.0 0.00 log_k_JSE76 0.0000 0.000 0.0 0.00 f_cyan_ilr_1 0.6349 0.000 0.0 0.00 f_cyan_ilr_2 0.0000 1.797 0.0 0.00 f_JCZ38_qlogis 0.0000 0.000 15.6 0.00 f_JSE76_qlogis 0.0000 0.000 0.0 17.52 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2283 2275 -1120 Optimised parameters: est. lower upper cyan_free_0 102.6517 101.40815 103.8952 log_k_cyan_free -2.8729 -3.18649 -2.5593 log_k_cyan_free_bound -2.7803 -3.60525 -1.9552 log_k_cyan_bound_free -3.5845 -4.16644 -3.0026 log_k_JCZ38 -2.3411 -2.89698 -1.7852 log_k_J9Z38 -5.2487 -6.01271 -4.4847 log_k_JSE76 -3.0259 -4.28274 -1.7690 f_cyan_ilr_1 0.7289 0.38214 1.0756 f_cyan_ilr_2 0.6891 0.18277 1.1954 f_JCZ38_qlogis 4.2162 0.47015 7.9622 f_JSE76_qlogis 5.8911 -20.19088 31.9730 a.1 2.7159 2.52587 2.9060 SD.log_k_cyan_free 0.3354 0.10979 0.5610 SD.log_k_cyan_free_bound 0.9061 0.30969 1.5025 SD.log_k_cyan_bound_free 0.6376 0.21229 1.0628 SD.log_k_JCZ38 0.5499 0.14533 0.9545 SD.log_k_J9Z38 0.7457 0.15106 1.3404 SD.log_k_JSE76 1.3822 0.47329 2.2912 SD.f_cyan_ilr_1 0.3820 0.13280 0.6313 SD.f_cyan_ilr_2 0.4317 0.06803 0.7953 SD.f_JCZ38_qlogis 1.8258 -0.25423 3.9059 SD.f_JSE76_qlogis 2.2348 -83.33679 87.8065 Correlation: cyn_f_0 lg_k_c_ lg_k_cyn_f_ lg_k_cyn_b_ l__JCZ3 l__J9Z3 log_k_cyan_free 0.1944 log_k_cyan_free_bound 0.0815 0.0814 log_k_cyan_bound_free 0.0106 0.0426 0.0585 log_k_JCZ38 -0.0231 -0.0106 -0.0089 -0.0051 log_k_J9Z38 -0.0457 -0.0108 0.0019 0.0129 0.0032 log_k_JSE76 -0.0054 -0.0024 -0.0017 -0.0005 0.1108 0.0009 f_cyan_ilr_1 0.0051 -0.0005 -0.0035 -0.0056 0.0131 -0.0967 f_cyan_ilr_2 -0.3182 -0.0771 -0.0309 -0.0038 0.0680 0.1643 f_JCZ38_qlogis 0.0834 0.0369 0.0302 0.0172 -0.1145 -0.0204 f_JSE76_qlogis -0.0553 -0.0365 -0.0441 -0.0414 0.1579 0.0175 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 log_k_cyan_free log_k_cyan_free_bound log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 -0.0002 f_cyan_ilr_2 0.0020 -0.0415 f_JCZ38_qlogis 0.0052 -0.0665 -0.3437 f_JSE76_qlogis 0.0066 0.0635 0.3491 -0.7487 Random effects: est. lower upper SD.log_k_cyan_free 0.3354 0.10979 0.5610 SD.log_k_cyan_free_bound 0.9061 0.30969 1.5025 SD.log_k_cyan_bound_free 0.6376 0.21229 1.0628 SD.log_k_JCZ38 0.5499 0.14533 0.9545 SD.log_k_J9Z38 0.7457 0.15106 1.3404 SD.log_k_JSE76 1.3822 0.47329 2.2912 SD.f_cyan_ilr_1 0.3820 0.13280 0.6313 SD.f_cyan_ilr_2 0.4317 0.06803 0.7953 SD.f_JCZ38_qlogis 1.8258 -0.25423 3.9059 SD.f_JSE76_qlogis 2.2348 -83.33679 87.8065 Variance model: est. lower upper a.1 2.716 2.526 2.906 Backtransformed parameters: est. lower upper cyan_free_0 1.027e+02 1.014e+02 103.89517 k_cyan_free 5.654e-02 4.132e-02 0.07736 k_cyan_free_bound 6.202e-02 2.718e-02 0.14153 k_cyan_bound_free 2.775e-02 1.551e-02 0.04966 k_JCZ38 9.622e-02 5.519e-02 0.16777 k_J9Z38 5.254e-03 2.447e-03 0.01128 k_JSE76 4.852e-02 1.380e-02 0.17051 f_cyan_free_to_JCZ38 6.197e-01 5.643e-01 0.84429 f_cyan_free_to_J9Z38 2.211e-01 5.643e-01 0.84429 f_JCZ38_to_JSE76 9.855e-01 6.154e-01 0.99965 f_JSE76_to_JCZ38 9.972e-01 1.703e-09 1.00000 Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.13466 0.01165 0.36490 Resulting formation fractions: ff cyan_free_JCZ38 0.619745 cyan_free_J9Z38 0.221083 cyan_free_sink 0.159172 cyan_free 1.000000 JCZ38_JSE76 0.985460 JCZ38_sink 0.014540 JSE76_JCZ38 0.997244 JSE76_sink 0.002756 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 23.293 158.67 47.77 5.147 59.5 JCZ38 7.203 23.93 NA NA NA J9Z38 131.918 438.22 NA NA NA JSE76 14.287 47.46 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 20:50:59 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 714.48 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 101.7511 -2.8370 -3.0162 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.6600 -2.2988 -5.3129 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -3.6991 0.6722 4.8596 f_JCZ38_qlogis f_JSE76_qlogis 13.4678 14.2149 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 5.629 0.000 0.000 log_k_cyan_free 0.000 0.446 0.000 log_k_cyan_free_bound 0.000 0.000 1.449 log_k_cyan_bound_free 0.000 0.000 0.000 log_k_JCZ38 0.000 0.000 0.000 log_k_J9Z38 0.000 0.000 0.000 log_k_JSE76 0.000 0.000 0.000 f_cyan_ilr_1 0.000 0.000 0.000 f_cyan_ilr_2 0.000 0.000 0.000 f_JCZ38_qlogis 0.000 0.000 0.000 f_JSE76_qlogis 0.000 0.000 0.000 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.0000 0.000 0.0000 log_k_cyan_free 0.000 0.0000 0.000 0.0000 log_k_cyan_free_bound 0.000 0.0000 0.000 0.0000 log_k_cyan_bound_free 1.213 0.0000 0.000 0.0000 log_k_JCZ38 0.000 0.7801 0.000 0.0000 log_k_J9Z38 0.000 0.0000 1.575 0.0000 log_k_JSE76 0.000 0.0000 0.000 0.8078 f_cyan_ilr_1 0.000 0.0000 0.000 0.0000 f_cyan_ilr_2 0.000 0.0000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0000 f_JSE76_qlogis 0.000 0.0000 0.000 0.0000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis cyan_free_0 0.0000 0.000 0.00 0.00 log_k_cyan_free 0.0000 0.000 0.00 0.00 log_k_cyan_free_bound 0.0000 0.000 0.00 0.00 log_k_cyan_bound_free 0.0000 0.000 0.00 0.00 log_k_JCZ38 0.0000 0.000 0.00 0.00 log_k_J9Z38 0.0000 0.000 0.00 0.00 log_k_JSE76 0.0000 0.000 0.00 0.00 f_cyan_ilr_1 0.6518 0.000 0.00 0.00 f_cyan_ilr_2 0.0000 9.981 0.00 0.00 f_JCZ38_qlogis 0.0000 0.000 14.26 0.00 f_JSE76_qlogis 0.0000 0.000 0.00 16.17 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2240 2231 -1098 Optimised parameters: est. lower upper cyan_free_0 100.73014 9.873e+01 1.027e+02 log_k_cyan_free -3.19634 -3.641e+00 -2.752e+00 log_k_cyan_free_bound -3.43533 -3.674e+00 -3.197e+00 log_k_cyan_bound_free -3.83282 -4.163e+00 -3.503e+00 log_k_JCZ38 -2.51065 -3.225e+00 -1.796e+00 log_k_J9Z38 -5.02539 -5.825e+00 -4.226e+00 log_k_JSE76 -3.24777 -4.163e+00 -2.333e+00 f_cyan_ilr_1 0.70640 3.562e-01 1.057e+00 f_cyan_ilr_2 1.42704 3.170e-01 2.537e+00 f_JCZ38_qlogis 2.84779 1.042e+00 4.654e+00 f_JSE76_qlogis 8.63674 -6.407e+02 6.580e+02 a.1 2.07082 1.877e+00 2.265e+00 b.1 0.06227 5.098e-02 7.355e-02 SD.log_k_cyan_free 0.49674 1.865e-01 8.069e-01 SD.log_k_cyan_bound_free 0.28537 6.809e-02 5.027e-01 SD.log_k_JCZ38 0.74846 2.305e-01 1.266e+00 SD.log_k_J9Z38 0.86077 2.713e-01 1.450e+00 SD.log_k_JSE76 0.97613 3.030e-01 1.649e+00 SD.f_cyan_ilr_1 0.38994 1.382e-01 6.417e-01 SD.f_cyan_ilr_2 0.82869 3.917e-02 1.618e+00 SD.f_JCZ38_qlogis 1.05000 -2.808e-02 2.128e+00 SD.f_JSE76_qlogis 0.44681 -3.985e+05 3.985e+05 Correlation: cyn_f_0 lg_k_c_ lg_k_cyn_f_ lg_k_cyn_b_ l__JCZ3 l__J9Z3 log_k_cyan_free 0.0936 log_k_cyan_free_bound 0.1302 0.1627 log_k_cyan_bound_free 0.0029 0.0525 0.5181 log_k_JCZ38 -0.0116 -0.0077 -0.0430 -0.0236 log_k_J9Z38 -0.0192 -0.0077 -0.0048 0.0229 -0.0005 log_k_JSE76 0.0007 -0.0020 -0.0134 -0.0072 0.1225 -0.0016 f_cyan_ilr_1 -0.0118 -0.0027 -0.0132 -0.0118 0.0127 -0.0505 f_cyan_ilr_2 -0.4643 -0.0762 -0.1245 0.0137 0.0497 0.1003 f_JCZ38_qlogis 0.0710 0.0371 0.1826 0.0925 -0.0869 -0.0130 f_JSE76_qlogis -0.0367 -0.0270 -0.2274 -0.1865 0.1244 0.0098 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 log_k_cyan_free log_k_cyan_free_bound log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 0.0036 f_cyan_ilr_2 0.0050 -0.0201 f_JCZ38_qlogis 0.0142 -0.0529 -0.2698 f_JSE76_qlogis 0.0064 0.0345 0.2015 -0.7058 Random effects: est. lower upper SD.log_k_cyan_free 0.4967 1.865e-01 8.069e-01 SD.log_k_cyan_bound_free 0.2854 6.809e-02 5.027e-01 SD.log_k_JCZ38 0.7485 2.305e-01 1.266e+00 SD.log_k_J9Z38 0.8608 2.713e-01 1.450e+00 SD.log_k_JSE76 0.9761 3.030e-01 1.649e+00 SD.f_cyan_ilr_1 0.3899 1.382e-01 6.417e-01 SD.f_cyan_ilr_2 0.8287 3.917e-02 1.618e+00 SD.f_JCZ38_qlogis 1.0500 -2.808e-02 2.128e+00 SD.f_JSE76_qlogis 0.4468 -3.985e+05 3.985e+05 Variance model: est. lower upper a.1 2.07082 1.87680 2.26483 b.1 0.06227 0.05098 0.07355 Backtransformed parameters: est. lower upper cyan_free_0 1.007e+02 9.873e+01 102.72898 k_cyan_free 4.091e-02 2.623e-02 0.06382 k_cyan_free_bound 3.221e-02 2.537e-02 0.04090 k_cyan_bound_free 2.165e-02 1.557e-02 0.03011 k_JCZ38 8.122e-02 3.975e-02 0.16594 k_J9Z38 6.569e-03 2.954e-03 0.01461 k_JSE76 3.886e-02 1.556e-02 0.09703 f_cyan_free_to_JCZ38 6.785e-01 6.102e-01 0.97309 f_cyan_free_to_J9Z38 2.498e-01 6.102e-01 0.97309 f_JCZ38_to_JSE76 9.452e-01 7.392e-01 0.99056 f_JSE76_to_JCZ38 9.998e-01 5.580e-279 1.00000 Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.08426 0.01051 0.41220 Resulting formation fractions: ff cyan_free_JCZ38 0.6784541 cyan_free_J9Z38 0.2498405 cyan_free_sink 0.0717054 cyan_free 1.0000000 JCZ38_JSE76 0.9452043 JCZ38_sink 0.0547957 JSE76_JCZ38 0.9998226 JSE76_sink 0.0001774 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 25.237 168.51 50.73 8.226 65.95 JCZ38 8.535 28.35 NA NA NA J9Z38 105.517 350.52 NA NA NA JSE76 17.837 59.25 NA NA NA"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"pathway-2-refined-fits","dir":"Articles > Prebuilt","previous_headings":"Appendix > Hierarchical fit listings","what":"Pathway 2, refined fits","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"Hierarchical FOMC path 2 fit reduced random effects, two-component error Hierarchical DFOP path 2 fit reduced random effects, constant variance Hierarchical DFOP path 2 fit reduced random effects, two-component error Hierarchical SFORB path 2 fit reduced random effects, constant variance Hierarchical SFORB path 2 fit reduced random effects, two-component error","code":"saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 21:04:10 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * (alpha/beta) * 1/((time/beta) + 1) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 786.038 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.4477 -1.8631 -5.1087 -2.5114 0.6826 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_alpha log_beta 4.7944 15.9616 13.1566 -0.1564 2.9781 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 7.701 0.000 0.000 0.000 0.0000 log_k_JCZ38 0.000 1.448 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.000 1.724 0.000 0.0000 log_k_JSE76 0.000 0.000 0.000 3.659 0.0000 f_cyan_ilr_1 0.000 0.000 0.000 0.000 0.6356 f_cyan_ilr_2 0.000 0.000 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.000 0.000 0.000 0.0000 f_JSE76_qlogis 0.000 0.000 0.000 0.000 0.0000 log_alpha 0.000 0.000 0.000 0.000 0.0000 log_beta 0.000 0.000 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_alpha log_beta cyan_0 0.00 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 10.32 0.00 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 12.23 0.00 0.0000 0.0000 f_JSE76_qlogis 0.00 0.00 14.99 0.0000 0.0000 log_alpha 0.00 0.00 0.00 0.3924 0.0000 log_beta 0.00 0.00 0.00 0.0000 0.5639 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2249 2242 -1106 Optimised parameters: est. lower upper cyan_0 101.24524 NA NA log_k_JCZ38 -2.85375 NA NA log_k_J9Z38 -5.07729 NA NA log_k_JSE76 -3.53511 NA NA f_cyan_ilr_1 0.67478 NA NA f_cyan_ilr_2 0.97152 NA NA f_JCZ38_qlogis 213.48001 NA NA f_JSE76_qlogis 2.02040 NA NA log_alpha -0.11041 NA NA log_beta 3.06575 NA NA a.1 2.05279 1.85495 2.2506 b.1 0.07116 0.05912 0.0832 SD.log_k_JCZ38 1.21713 0.44160 1.9927 SD.log_k_J9Z38 0.88268 0.27541 1.4900 SD.log_k_JSE76 0.59452 0.15005 1.0390 SD.f_cyan_ilr_1 0.35370 0.12409 0.5833 SD.f_cyan_ilr_2 0.78186 0.18547 1.3782 SD.log_alpha 0.27781 0.08168 0.4739 SD.log_beta 0.32608 0.06490 0.5873 Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.2171 0.44160 1.9927 SD.log_k_J9Z38 0.8827 0.27541 1.4900 SD.log_k_JSE76 0.5945 0.15005 1.0390 SD.f_cyan_ilr_1 0.3537 0.12409 0.5833 SD.f_cyan_ilr_2 0.7819 0.18547 1.3782 SD.log_alpha 0.2778 0.08168 0.4739 SD.log_beta 0.3261 0.06490 0.5873 Variance model: est. lower upper a.1 2.05279 1.85495 2.2506 b.1 0.07116 0.05912 0.0832 Backtransformed parameters: est. lower upper cyan_0 1.012e+02 NA NA k_JCZ38 5.763e-02 NA NA k_J9Z38 6.237e-03 NA NA k_JSE76 2.916e-02 NA NA f_cyan_to_JCZ38 6.354e-01 NA NA f_cyan_to_J9Z38 2.447e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA f_JSE76_to_JCZ38 8.829e-01 NA NA alpha 8.955e-01 NA NA beta 2.145e+01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.6354 cyan_J9Z38 0.2447 cyan_sink 0.1200 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 JSE76_JCZ38 0.8829 JSE76_sink 0.1171 Estimated disappearance times: DT50 DT90 DT50back cyan 25.07 259.21 78.03 JCZ38 12.03 39.96 NA J9Z38 111.14 369.19 NA JSE76 23.77 78.98 NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 21:05:56 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 892.139 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 102.4380 -2.3107 -5.3123 -3.7120 0.6757 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 1.1439 13.1194 12.3492 -1.9317 -4.4557 g_qlogis -0.5644 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 4.591 0.0000 0.000 0.0 0.0000 log_k_JCZ38 0.000 0.7966 0.000 0.0 0.0000 log_k_J9Z38 0.000 0.0000 1.561 0.0 0.0000 log_k_JSE76 0.000 0.0000 0.000 0.8 0.0000 f_cyan_ilr_1 0.000 0.0000 0.000 0.0 0.6349 f_cyan_ilr_2 0.000 0.0000 0.000 0.0 0.0000 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0 0.0000 f_JSE76_qlogis 0.000 0.0000 0.000 0.0 0.0000 log_k1 0.000 0.0000 0.000 0.0 0.0000 log_k2 0.000 0.0000 0.000 0.0 0.0000 g_qlogis 0.000 0.0000 0.000 0.0 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 cyan_0 0.000 0.00 0.00 0.000 0.0000 log_k_JCZ38 0.000 0.00 0.00 0.000 0.0000 log_k_J9Z38 0.000 0.00 0.00 0.000 0.0000 log_k_JSE76 0.000 0.00 0.00 0.000 0.0000 f_cyan_ilr_1 0.000 0.00 0.00 0.000 0.0000 f_cyan_ilr_2 1.797 0.00 0.00 0.000 0.0000 f_JCZ38_qlogis 0.000 13.86 0.00 0.000 0.0000 f_JSE76_qlogis 0.000 0.00 13.91 0.000 0.0000 log_k1 0.000 0.00 0.00 1.106 0.0000 log_k2 0.000 0.00 0.00 0.000 0.6141 g_qlogis 0.000 0.00 0.00 0.000 0.0000 g_qlogis cyan_0 0.000 log_k_JCZ38 0.000 log_k_J9Z38 0.000 log_k_JSE76 0.000 f_cyan_ilr_1 0.000 f_cyan_ilr_2 0.000 f_JCZ38_qlogis 0.000 f_JSE76_qlogis 0.000 log_k1 0.000 log_k2 0.000 g_qlogis 1.595 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2282 2274 -1121 Optimised parameters: est. lower upper cyan_0 102.6036 NA NA log_k_JCZ38 -2.9348 NA NA log_k_J9Z38 -5.1617 NA NA log_k_JSE76 -3.6396 NA NA f_cyan_ilr_1 0.6991 NA NA f_cyan_ilr_2 0.6341 NA NA f_JCZ38_qlogis 4232.3011 NA NA f_JSE76_qlogis 1.9658 NA NA log_k1 -1.9503 NA NA log_k2 -4.4745 NA NA g_qlogis -0.4967 NA NA a.1 2.7461 2.59274 2.8994 SD.log_k_JCZ38 1.3178 0.47602 2.1596 SD.log_k_J9Z38 0.7022 0.15061 1.2538 SD.log_k_JSE76 0.6566 0.15613 1.1570 SD.f_cyan_ilr_1 0.3409 0.11666 0.5652 SD.f_cyan_ilr_2 0.4385 0.09482 0.7821 SD.log_k1 0.7381 0.25599 1.2202 SD.log_k2 0.5133 0.18152 0.8450 SD.g_qlogis 0.9866 0.35681 1.6164 Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.3178 0.47602 2.1596 SD.log_k_J9Z38 0.7022 0.15061 1.2538 SD.log_k_JSE76 0.6566 0.15613 1.1570 SD.f_cyan_ilr_1 0.3409 0.11666 0.5652 SD.f_cyan_ilr_2 0.4385 0.09482 0.7821 SD.log_k1 0.7381 0.25599 1.2202 SD.log_k2 0.5133 0.18152 0.8450 SD.g_qlogis 0.9866 0.35681 1.6164 Variance model: est. lower upper a.1 2.746 2.593 2.899 Backtransformed parameters: est. lower upper cyan_0 1.026e+02 NA NA k_JCZ38 5.314e-02 NA NA k_J9Z38 5.732e-03 NA NA k_JSE76 2.626e-02 NA NA f_cyan_to_JCZ38 6.051e-01 NA NA f_cyan_to_J9Z38 2.251e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA f_JSE76_to_JCZ38 8.772e-01 NA NA k1 1.422e-01 NA NA k2 1.140e-02 NA NA g 3.783e-01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.6051 cyan_J9Z38 0.2251 cyan_sink 0.1698 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 JSE76_JCZ38 0.8772 JSE76_sink 0.1228 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 22.05 160.35 48.27 4.873 60.83 JCZ38 13.04 43.33 NA NA NA J9Z38 120.93 401.73 NA NA NA JSE76 26.39 87.68 NA NA NA saemix version used for fitting: 3.3 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:02 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan d_JCZ38/dt = + f_cyan_to_JCZ38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_to_J9Z38 * ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * cyan - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 898.534 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 101.7393 -1.4493 -5.0118 -2.1269 0.6720 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 7.3362 13.4423 13.2659 -2.0061 -4.5527 g_qlogis -0.5806 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_0 log_k_JCZ38 log_k_J9Z38 log_k_JSE76 f_cyan_ilr_1 cyan_0 5.604 0.00 0.000 0.000 0.0000 log_k_JCZ38 0.000 2.77 0.000 0.000 0.0000 log_k_J9Z38 0.000 0.00 1.662 0.000 0.0000 log_k_JSE76 0.000 0.00 0.000 5.021 0.0000 f_cyan_ilr_1 0.000 0.00 0.000 0.000 0.6519 f_cyan_ilr_2 0.000 0.00 0.000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.00 0.000 0.000 0.0000 f_JSE76_qlogis 0.000 0.00 0.000 0.000 0.0000 log_k1 0.000 0.00 0.000 0.000 0.0000 log_k2 0.000 0.00 0.000 0.000 0.0000 g_qlogis 0.000 0.00 0.000 0.000 0.0000 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis log_k1 log_k2 cyan_0 0.00 0.00 0.00 0.0000 0.0000 log_k_JCZ38 0.00 0.00 0.00 0.0000 0.0000 log_k_J9Z38 0.00 0.00 0.00 0.0000 0.0000 log_k_JSE76 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_1 0.00 0.00 0.00 0.0000 0.0000 f_cyan_ilr_2 13.37 0.00 0.00 0.0000 0.0000 f_JCZ38_qlogis 0.00 14.21 0.00 0.0000 0.0000 f_JSE76_qlogis 0.00 0.00 14.58 0.0000 0.0000 log_k1 0.00 0.00 0.00 0.8453 0.0000 log_k2 0.00 0.00 0.00 0.0000 0.5969 g_qlogis 0.00 0.00 0.00 0.0000 0.0000 g_qlogis cyan_0 0.00 log_k_JCZ38 0.00 log_k_J9Z38 0.00 log_k_JSE76 0.00 f_cyan_ilr_1 0.00 f_cyan_ilr_2 0.00 f_JCZ38_qlogis 0.00 f_JSE76_qlogis 0.00 log_k1 0.00 log_k2 0.00 g_qlogis 1.69 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2237 2229 -1099 Optimised parameters: est. lower upper cyan_0 101.00243 NA NA log_k_JCZ38 -2.80828 NA NA log_k_J9Z38 -5.04449 NA NA log_k_JSE76 -3.66981 NA NA f_cyan_ilr_1 0.72564 NA NA f_cyan_ilr_2 1.37978 NA NA f_JCZ38_qlogis 1.98726 NA NA f_JSE76_qlogis 414.80884 NA NA log_k1 -2.38601 NA NA log_k2 -4.63632 NA NA g_qlogis -0.33920 NA NA a.1 2.10837 1.91261 2.30413 b.1 0.06223 0.05085 0.07361 SD.log_k_JCZ38 1.30902 0.48128 2.13675 SD.log_k_J9Z38 0.83882 0.25790 1.41974 SD.log_k_JSE76 0.58104 0.14201 1.02008 SD.f_cyan_ilr_1 0.35421 0.12398 0.58443 SD.f_cyan_ilr_2 0.79373 0.12007 1.46739 SD.log_k2 0.27476 0.08557 0.46394 SD.g_qlogis 0.96170 0.35463 1.56878 Correlation is not available Random effects: est. lower upper SD.log_k_JCZ38 1.3090 0.48128 2.1367 SD.log_k_J9Z38 0.8388 0.25790 1.4197 SD.log_k_JSE76 0.5810 0.14201 1.0201 SD.f_cyan_ilr_1 0.3542 0.12398 0.5844 SD.f_cyan_ilr_2 0.7937 0.12007 1.4674 SD.log_k2 0.2748 0.08557 0.4639 SD.g_qlogis 0.9617 0.35463 1.5688 Variance model: est. lower upper a.1 2.10837 1.91261 2.30413 b.1 0.06223 0.05085 0.07361 Backtransformed parameters: est. lower upper cyan_0 1.010e+02 NA NA k_JCZ38 6.031e-02 NA NA k_J9Z38 6.445e-03 NA NA k_JSE76 2.548e-02 NA NA f_cyan_to_JCZ38 6.808e-01 NA NA f_cyan_to_J9Z38 2.440e-01 NA NA f_JCZ38_to_JSE76 8.795e-01 NA NA f_JSE76_to_JCZ38 1.000e+00 NA NA k1 9.200e-02 NA NA k2 9.693e-03 NA NA g 4.160e-01 NA NA Resulting formation fractions: ff cyan_JCZ38 0.68081 cyan_J9Z38 0.24398 cyan_sink 0.07521 JCZ38_JSE76 0.87945 JCZ38_sink 0.12055 JSE76_JCZ38 1.00000 JSE76_sink 0.00000 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 cyan 25.00 182.05 54.8 7.535 71.51 JCZ38 11.49 38.18 NA NA NA J9Z38 107.55 357.28 NA NA NA JSE76 27.20 90.36 NA NA NA saemix version used for fitting: 3.3 mkin version used for pre-fitting: 1.2.10 R version used for fitting: 4.5.0 Date of fit: Mon May 12 21:05:52 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 888.333 s Using 300, 100 iterations and 10 chains Variance model: Constant variance Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 102.4395 -2.7673 -2.8942 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.6201 -2.3107 -5.3123 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -3.7120 0.6754 1.1448 f_JCZ38_qlogis f_JSE76_qlogis 14.8408 15.4734 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 4.589 0.0000 0.00 log_k_cyan_free 0.000 0.4849 0.00 log_k_cyan_free_bound 0.000 0.0000 1.62 log_k_cyan_bound_free 0.000 0.0000 0.00 log_k_JCZ38 0.000 0.0000 0.00 log_k_J9Z38 0.000 0.0000 0.00 log_k_JSE76 0.000 0.0000 0.00 f_cyan_ilr_1 0.000 0.0000 0.00 f_cyan_ilr_2 0.000 0.0000 0.00 f_JCZ38_qlogis 0.000 0.0000 0.00 f_JSE76_qlogis 0.000 0.0000 0.00 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.0000 0.000 0.0 log_k_cyan_free 0.000 0.0000 0.000 0.0 log_k_cyan_free_bound 0.000 0.0000 0.000 0.0 log_k_cyan_bound_free 1.197 0.0000 0.000 0.0 log_k_JCZ38 0.000 0.7966 0.000 0.0 log_k_J9Z38 0.000 0.0000 1.561 0.0 log_k_JSE76 0.000 0.0000 0.000 0.8 f_cyan_ilr_1 0.000 0.0000 0.000 0.0 f_cyan_ilr_2 0.000 0.0000 0.000 0.0 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0 f_JSE76_qlogis 0.000 0.0000 0.000 0.0 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis cyan_free_0 0.0000 0.000 0.0 0.00 log_k_cyan_free 0.0000 0.000 0.0 0.00 log_k_cyan_free_bound 0.0000 0.000 0.0 0.00 log_k_cyan_bound_free 0.0000 0.000 0.0 0.00 log_k_JCZ38 0.0000 0.000 0.0 0.00 log_k_J9Z38 0.0000 0.000 0.0 0.00 log_k_JSE76 0.0000 0.000 0.0 0.00 f_cyan_ilr_1 0.6349 0.000 0.0 0.00 f_cyan_ilr_2 0.0000 1.797 0.0 0.00 f_JCZ38_qlogis 0.0000 0.000 15.6 0.00 f_JSE76_qlogis 0.0000 0.000 0.0 17.52 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2280 2272 -1120 Optimised parameters: est. lower upper cyan_free_0 102.6532 NA NA log_k_cyan_free -2.8547 NA NA log_k_cyan_free_bound -2.7004 NA NA log_k_cyan_bound_free -3.5078 NA NA log_k_JCZ38 -2.9255 NA NA log_k_J9Z38 -5.1089 NA NA log_k_JSE76 -3.6263 NA NA f_cyan_ilr_1 0.6873 NA NA f_cyan_ilr_2 0.6498 NA NA f_JCZ38_qlogis 3624.2149 NA NA f_JSE76_qlogis 1.9991 NA NA a.1 2.7472 2.55559 2.9388 SD.log_k_cyan_free 0.3227 0.10296 0.5423 SD.log_k_cyan_free_bound 0.8757 0.29525 1.4562 SD.log_k_cyan_bound_free 0.6128 0.20220 1.0233 SD.log_k_JCZ38 1.3431 0.48474 2.2014 SD.log_k_J9Z38 0.6881 0.14714 1.2291 SD.log_k_JSE76 0.6461 0.15321 1.1390 SD.f_cyan_ilr_1 0.3361 0.11376 0.5585 SD.f_cyan_ilr_2 0.4286 0.08419 0.7730 Correlation is not available Random effects: est. lower upper SD.log_k_cyan_free 0.3227 0.10296 0.5423 SD.log_k_cyan_free_bound 0.8757 0.29525 1.4562 SD.log_k_cyan_bound_free 0.6128 0.20220 1.0233 SD.log_k_JCZ38 1.3431 0.48474 2.2014 SD.log_k_J9Z38 0.6881 0.14714 1.2291 SD.log_k_JSE76 0.6461 0.15321 1.1390 SD.f_cyan_ilr_1 0.3361 0.11376 0.5585 SD.f_cyan_ilr_2 0.4286 0.08419 0.7730 Variance model: est. lower upper a.1 2.747 2.556 2.939 Backtransformed parameters: est. lower upper cyan_free_0 1.027e+02 NA NA k_cyan_free 5.758e-02 NA NA k_cyan_free_bound 6.718e-02 NA NA k_cyan_bound_free 2.996e-02 NA NA k_JCZ38 5.364e-02 NA NA k_J9Z38 6.042e-03 NA NA k_JSE76 2.662e-02 NA NA f_cyan_free_to_JCZ38 6.039e-01 NA NA f_cyan_free_to_J9Z38 2.285e-01 NA NA f_JCZ38_to_JSE76 1.000e+00 NA NA f_JSE76_to_JCZ38 8.807e-01 NA NA Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.1426 0.0121 0.3484 Resulting formation fractions: ff cyan_free_JCZ38 0.6039 cyan_free_J9Z38 0.2285 cyan_free_sink 0.1676 cyan_free 1.0000 JCZ38_JSE76 1.0000 JCZ38_sink 0.0000 JSE76_JCZ38 0.8807 JSE76_sink 0.1193 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 23.84 154.95 46.65 4.86 57.31 JCZ38 12.92 42.93 NA NA NA J9Z38 114.71 381.07 NA NA NA JSE76 26.04 86.51 NA NA NA saemix version used for fitting: 3.3 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:21 2025 Date of summary: Mon May 12 21:06:22 2025 Equations: d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound * cyan_free + k_cyan_bound_free * cyan_bound d_cyan_bound/dt = + k_cyan_free_bound * cyan_free - k_cyan_bound_free * cyan_bound d_JCZ38/dt = + f_cyan_free_to_JCZ38 * k_cyan_free * cyan_free - k_JCZ38 * JCZ38 + f_JSE76_to_JCZ38 * k_JSE76 * JSE76 d_J9Z38/dt = + f_cyan_free_to_J9Z38 * k_cyan_free * cyan_free - k_J9Z38 * J9Z38 d_JSE76/dt = + f_JCZ38_to_JSE76 * k_JCZ38 * JCZ38 - k_JSE76 * JSE76 Data: 433 observations of 4 variable(s) grouped in 5 datasets Model predictions using solution type deSolve Fitted in 916.619 s Using 300, 100 iterations and 10 chains Variance model: Two-component variance function Starting values for degradation parameters: cyan_free_0 log_k_cyan_free log_k_cyan_free_bound 101.7511 -2.8370 -3.0162 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 -3.6600 -2.2988 -5.3129 log_k_JSE76 f_cyan_ilr_1 f_cyan_ilr_2 -3.6991 0.6722 4.8596 f_JCZ38_qlogis f_JSE76_qlogis 13.4678 14.2149 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): cyan_free_0 log_k_cyan_free log_k_cyan_free_bound cyan_free_0 5.629 0.000 0.000 log_k_cyan_free 0.000 0.446 0.000 log_k_cyan_free_bound 0.000 0.000 1.449 log_k_cyan_bound_free 0.000 0.000 0.000 log_k_JCZ38 0.000 0.000 0.000 log_k_J9Z38 0.000 0.000 0.000 log_k_JSE76 0.000 0.000 0.000 f_cyan_ilr_1 0.000 0.000 0.000 f_cyan_ilr_2 0.000 0.000 0.000 f_JCZ38_qlogis 0.000 0.000 0.000 f_JSE76_qlogis 0.000 0.000 0.000 log_k_cyan_bound_free log_k_JCZ38 log_k_J9Z38 log_k_JSE76 cyan_free_0 0.000 0.0000 0.000 0.0000 log_k_cyan_free 0.000 0.0000 0.000 0.0000 log_k_cyan_free_bound 0.000 0.0000 0.000 0.0000 log_k_cyan_bound_free 1.213 0.0000 0.000 0.0000 log_k_JCZ38 0.000 0.7801 0.000 0.0000 log_k_J9Z38 0.000 0.0000 1.575 0.0000 log_k_JSE76 0.000 0.0000 0.000 0.8078 f_cyan_ilr_1 0.000 0.0000 0.000 0.0000 f_cyan_ilr_2 0.000 0.0000 0.000 0.0000 f_JCZ38_qlogis 0.000 0.0000 0.000 0.0000 f_JSE76_qlogis 0.000 0.0000 0.000 0.0000 f_cyan_ilr_1 f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis cyan_free_0 0.0000 0.000 0.00 0.00 log_k_cyan_free 0.0000 0.000 0.00 0.00 log_k_cyan_free_bound 0.0000 0.000 0.00 0.00 log_k_cyan_bound_free 0.0000 0.000 0.00 0.00 log_k_JCZ38 0.0000 0.000 0.00 0.00 log_k_J9Z38 0.0000 0.000 0.00 0.00 log_k_JSE76 0.0000 0.000 0.00 0.00 f_cyan_ilr_1 0.6518 0.000 0.00 0.00 f_cyan_ilr_2 0.0000 9.981 0.00 0.00 f_JCZ38_qlogis 0.0000 0.000 14.26 0.00 f_JSE76_qlogis 0.0000 0.000 0.00 16.17 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 2241 2233 -1101 Optimised parameters: est. lower upper cyan_free_0 100.95469 NA NA log_k_cyan_free -3.18706 NA NA log_k_cyan_free_bound -3.38455 NA NA log_k_cyan_bound_free -3.75788 NA NA log_k_JCZ38 -2.77024 NA NA log_k_J9Z38 -5.03665 NA NA log_k_JSE76 -3.60289 NA NA f_cyan_ilr_1 0.72263 NA NA f_cyan_ilr_2 1.45352 NA NA f_JCZ38_qlogis 2.00778 NA NA f_JSE76_qlogis 941.58570 NA NA a.1 2.11130 1.91479 2.30780 b.1 0.06299 0.05152 0.07445 SD.log_k_cyan_free 0.50098 0.18805 0.81390 SD.log_k_cyan_bound_free 0.31671 0.08467 0.54875 SD.log_k_JCZ38 1.25865 0.45932 2.05798 SD.log_k_J9Z38 0.86833 0.27222 1.46444 SD.log_k_JSE76 0.59325 0.14711 1.03940 SD.f_cyan_ilr_1 0.35705 0.12521 0.58890 SD.f_cyan_ilr_2 0.88541 0.13797 1.63286 Correlation is not available Random effects: est. lower upper SD.log_k_cyan_free 0.5010 0.18805 0.8139 SD.log_k_cyan_bound_free 0.3167 0.08467 0.5487 SD.log_k_JCZ38 1.2587 0.45932 2.0580 SD.log_k_J9Z38 0.8683 0.27222 1.4644 SD.log_k_JSE76 0.5933 0.14711 1.0394 SD.f_cyan_ilr_1 0.3571 0.12521 0.5889 SD.f_cyan_ilr_2 0.8854 0.13797 1.6329 Variance model: est. lower upper a.1 2.11130 1.91479 2.30780 b.1 0.06299 0.05152 0.07445 Backtransformed parameters: est. lower upper cyan_free_0 1.010e+02 NA NA k_cyan_free 4.129e-02 NA NA k_cyan_free_bound 3.389e-02 NA NA k_cyan_bound_free 2.333e-02 NA NA k_JCZ38 6.265e-02 NA NA k_J9Z38 6.495e-03 NA NA k_JSE76 2.724e-02 NA NA f_cyan_free_to_JCZ38 6.844e-01 NA NA f_cyan_free_to_J9Z38 2.463e-01 NA NA f_JCZ38_to_JSE76 8.816e-01 NA NA f_JSE76_to_JCZ38 1.000e+00 NA NA Estimated Eigenvalues of SFORB model(s): cyan_b1 cyan_b2 cyan_g 0.08751 0.01101 0.39586 Resulting formation fractions: ff cyan_free_JCZ38 0.68444 cyan_free_J9Z38 0.24633 cyan_free_sink 0.06923 cyan_free 1.00000 JCZ38_JSE76 0.88161 JCZ38_sink 0.11839 JSE76_JCZ38 1.00000 JSE76_sink 0.00000 Estimated disappearance times: DT50 DT90 DT50back DT50_cyan_b1 DT50_cyan_b2 cyan 25.36 163.36 49.18 7.921 62.95 JCZ38 11.06 36.75 NA NA NA J9Z38 106.71 354.49 NA NA NA JSE76 25.44 84.51 NA NA NA"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"session-info","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Session info","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"","code":"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 version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C [3] LC_TIME=de_DE.UTF-8 LC_COLLATE=de_DE.UTF-8 [5] LC_MONETARY=de_DE.UTF-8 LC_MESSAGES=de_DE.UTF-8 [7] LC_PAPER=de_DE.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=de_DE.UTF-8 LC_IDENTIFICATION=C time zone: Europe/Berlin tzcode source: system (glibc) attached base packages: [1] parallel stats graphics grDevices utils datasets methods [8] base other attached packages: [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] sass_0.4.9 generics_0.1.3 lattice_0.22-6 digest_0.6.37 [5] magrittr_2.0.3 evaluate_1.0.3 grid_4.5.0 fastmap_1.2.0 [9] cellranger_1.1.0 jsonlite_1.9.0 processx_3.8.6 pkgbuild_1.4.6 [13] deSolve_1.40 mclust_6.1.1 ps_1.9.0 gridExtra_2.3 [17] scales_1.3.0 codetools_0.2-20 textshaping_1.0.0 jquerylib_0.1.4 [21] cli_3.6.4 rlang_1.1.5 munsell_0.5.1 cachem_1.1.0 [25] yaml_2.3.10 inline_0.3.21 tools_4.5.0 dplyr_1.1.4 [29] colorspace_2.1-1 ggplot2_3.5.1 vctrs_0.6.5 R6_2.6.1 [33] zoo_1.8-13 lifecycle_1.0.4 fs_1.6.5 htmlwidgets_1.6.4 [37] MASS_7.3-65 ragg_1.3.3 pkgconfig_2.0.3 desc_1.4.3 [41] callr_3.7.6 pkgdown_2.1.1 pillar_1.10.1 bslib_0.9.0 [45] gtable_0.3.6 glue_1.8.0 systemfonts_1.2.1 xfun_0.51 [49] tibble_3.2.1 lmtest_0.9-40 tidyselect_1.2.1 htmltools_0.5.8.1 [53] nlme_3.1-168 compiler_4.5.0 readxl_1.4.4"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_cyan_pathway.html","id":"hardware-info","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Hardware info","title":"Testing hierarchical pathway kinetics with residue data on cyantraniliprole","text":"","code":"CPU model: AMD Ryzen 9 7950X 16-Core Processor MemTotal: 64927780 kB"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"introduction","dir":"Articles > Prebuilt","previous_headings":"","what":"Introduction","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"purpose document demonstrate nonlinear hierarchical models (NLHM) based parent degradation models SFO, FOMC, DFOP HS can fitted mkin package. assembled course work package 1.1 Project Number 173340 (Application nonlinear hierarchical models kinetic evaluation chemical degradation data) German Environment Agency carried 2022 2023. mkin package used version 1.2.10. contains test data functions used evaluations. saemix package used backend fitting NLHM, also loaded make convergence plot function available. document processed knitr package, also provides kable function used improve display tabular data R markdown documents. parallel processing, parallel package used.","code":"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) }"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"data","dir":"Articles > Prebuilt","previous_headings":"","what":"Data","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"test data available mkin package object class mkindsg (mkin dataset group) identifier dimethenamid_2018. following preprocessing steps still necessary: data available enantiomer dimethenamid-P (DMTAP) renamed substance name data racemic mixture dimethenamid (DMTA). reason difference degradation behaviour identified EU risk assessment. data transformation products unnecessary columns discarded observation times dataset multiplied corresponding normalisation factor also available dataset, order make possible describe datasets single set parameters independent temperature Finally, datasets observed soil (Elliot 1 Elliot 2) combined, resulting dimethenamid (DMTA) data six soils. following commented R code performs preprocessing. following tables show 6 datasets. Dataset Calke Dataset Borstel Dataset Flaach Dataset BBA 2.2 Dataset BBA 2.3 Dataset Elliot","code":"# 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 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\") }"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"separate-evaluations","dir":"Articles > Prebuilt","previous_headings":"","what":"Separate evaluations","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"order obtain suitable starting parameters NLHM fits, separate fits four models data soil generated using mmkin function mkin package. first step, constant variance assumed. Convergence checked status function. table , OK indicates convergence, C indicates failure converge. separate fits constant variance converged, sole exception HS fit BBA 2.2 data. prepare fitting NLHM using two-component error model, separate fits updated assuming two-component error. Using two-component error model, one fit converge constant variance converge, non-SFO fits failed converge.","code":"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() f_sep_tc <- update(f_sep_const, error_model = \"tc\") status(f_sep_tc) |> kable()"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"hierarchichal-model-fits","dir":"Articles > Prebuilt","previous_headings":"","what":"Hierarchichal model fits","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"following code fits eight versions hierarchical models data, using SFO, FOMC, DFOP HS parent compound, using either constant variance two-component error error model. default parameter distribution model mkin allows variation degradation parameters across assumed population soils. words, degradation parameter associated random effect first step. mhmkin function makes possible fit eight versions parallel (given sufficient number computing cores available) save execution time. Convergence plots summaries fits shown appendix. output status function shows fits terminated successfully. AIC BIC values show biphasic models DFOP HS give best fits. DFOP model preferred , better mechanistic basis batch experiments constant incubation conditions. Also, shows lowest AIC BIC values first set fits combined two-component error model. Therefore, DFOP model selected refinements fits aim make model fully identifiable.","code":"f_saem <- mhmkin(list(f_sep_const, f_sep_tc), transformations = \"saemix\") status(f_saem) |> kable() anova(f_saem) |> kable(digits = 1)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"parameter-identifiability-based-on-the-fisher-information-matrix","dir":"Articles > Prebuilt","previous_headings":"Hierarchichal model fits","what":"Parameter identifiability based on the Fisher Information Matrix","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"Using illparms function, ill-defined statistical model parameters standard deviations degradation parameters population error model parameters can found. According illparms function, fitted standard deviation second kinetic rate constant k2 ill-defined DFOP fits. suggests different values obtained standard deviation using different starting values. thus identified overparameterisation addressed removing random effect k2 parameter model. resulting fit, checked whether still ill-defined parameters, case. , refined model compared previous best model. model without random effect k2 reduced version previous model. Therefore, models nested can compared using likelihood ratio test. achieved argument test = TRUE anova function. AIC BIC criteria lower removal ill-defined random effect k2. p value likelihood ratio test much greater 0.05, indicating model higher likelihood (model random effects degradation parameters f_saem[[\"DFOP\", \"tc\"]]) fit significantly better model lower likelihood (reduced model f_saem_dfop_tc_no_ranef_k2). Therefore, AIC, BIC likelihood ratio test suggest use reduced model. convergence fit checked visually. Convergence plot NLHM DFOP fit two-component error without random effect ‘k2’ parameters appear converged satisfactory degree. final fit plotted using plot method mkin package. Plot final NLHM DFOP fit Finally, summary report fit produced.","code":"illparms(f_saem) |> kable() f_saem_dfop_tc_no_ranef_k2 <- update(f_saem[[\"DFOP\", \"tc\"]], no_random_effect = \"k2\") illparms(f_saem_dfop_tc_no_ranef_k2) anova(f_saem[[\"DFOP\", \"tc\"]], f_saem_dfop_tc_no_ranef_k2, test = TRUE) |> kable(format.args = list(digits = 4)) plot(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.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 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * DMTA Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 3.953 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: DMTA_0 k1 k2 g 98.71186 0.08675 0.01374 0.93491 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 g DMTA_0 98.71 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 g 0.00 0 0 1 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 663.7 661.8 -322.9 Optimised parameters: est. lower upper DMTA_0 98.256267 96.286112 100.22642 k1 0.064037 0.033281 0.09479 k2 0.008469 0.006002 0.01094 g 0.954167 0.914460 0.99387 a.1 1.061795 0.878608 1.24498 b.1 0.029550 0.022593 0.03651 SD.DMTA_0 2.068581 0.427178 3.70998 SD.k1 0.598285 0.258235 0.93833 SD.g 1.016689 0.360061 1.67332 Correlation: DMTA_0 k1 k2 k1 0.0213 k2 0.0541 0.0344 g -0.0521 -0.0286 -0.2744 Random effects: est. lower upper SD.DMTA_0 2.0686 0.4272 3.7100 SD.k1 0.5983 0.2582 0.9383 SD.g 1.0167 0.3601 1.6733 Variance model: est. lower upper a.1 1.06180 0.87861 1.24498 b.1 0.02955 0.02259 0.03651 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 DMTA 11.45 41.32 12.44 10.82 81.85"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"alternative-check-of-parameter-identifiability","dir":"Articles > Prebuilt","previous_headings":"Hierarchichal model fits","what":"Alternative check of parameter identifiability","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"parameter check used illparms function based quadratic approximation likelihood surface near optimum, calculated using Fisher Information Matrix (FIM). alternative way check parameter identifiability (Duchesne et al. 2021) based multistart approach recently implemented mkin. graph shows boxplots parameters obtained 50 runs saem algorithm different parameter combinations, sampled range parameters obtained individual datasets fitted separately using nonlinear regression. Scaled parameters multistart runs, full model graph clearly confirms lack identifiability variance k2 full model. overparameterisation model also indicates lack identifiability variance parameter g. parameter boxplots multistart runs reduced model shown indicate runs give similar results, regardless starting parameters. Scaled parameters multistart runs, reduced model parameters top 25% fits shown (based feature introduced mkin 1.2.2 currently development), scatter even less shown . Scaled parameters multistart runs, reduced model, fits top 25% likelihood values","code":"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_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\") 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\")"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"conclusions","dir":"Articles > Prebuilt","previous_headings":"","what":"Conclusions","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"Fitting four parent degradation models SFO, FOMC, DFOP HS part hierarchical model fits two different error models normal distributions transformed degradation parameters works without technical problems. biphasic models DFOP HS gave best fit data, default parameter distribution model fully identifiable. Removing random effect second kinetic rate constant DFOP model resulted reduced model fully identifiable showed lowest values model selection criteria AIC BIC. reliability identification model parameters confirmed using multiple starting values.","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"acknowledgements","dir":"Articles > Prebuilt","previous_headings":"","what":"Acknowledgements","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"helpful comments Janina Wöltjen German Environment Agency gratefully acknowledged.","code":""},{"path":[]},{"path":[]},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"hierarchical-model-fit-listings","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Hierarchical model fit listings","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"Hierarchical mkin fit SFO model error model const Hierarchical mkin fit SFO model error model tc Hierarchical mkin fit FOMC model error model const Hierarchical mkin fit FOMC model error model tc Hierarchical mkin fit DFOP model error model const Hierarchical mkin fit DFOP model error model tc Hierarchical mkin fit HS model error model const Hierarchical mkin fit HS model error model tc","code":"saemix version used for fitting: 3.3 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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 0.793 s Using 300, 100 iterations and 9 chains Variance model: Constant variance Starting values for degradation parameters: DMTA_0 k_DMTA 97.2953 0.0566 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k_DMTA DMTA_0 97.3 0 k_DMTA 0.0 1 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 796.3 795.3 -393.2 Optimised parameters: est. lower upper DMTA_0 97.28130 95.71113 98.8515 k_DMTA 0.05665 0.02909 0.0842 a.1 2.66442 2.35579 2.9731 SD.DMTA_0 1.54776 0.15447 2.9411 SD.k_DMTA 0.60690 0.26248 0.9513 Correlation: DMTA_0 k_DMTA 0.0168 Random effects: est. lower upper SD.DMTA_0 1.5478 0.1545 2.9411 SD.k_DMTA 0.6069 0.2625 0.9513 Variance model: est. lower upper a.1 2.664 2.356 2.973 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.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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 2.462 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: DMTA_0 k_DMTA 96.99175 0.05603 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k_DMTA DMTA_0 96.99 0 k_DMTA 0.00 1 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 798.3 797.1 -393.2 Optimised parameters: est. lower upper DMTA_0 97.271822 95.70316 98.84049 k_DMTA 0.056638 0.02911 0.08417 a.1 2.660081 2.27492 3.04525 b.1 0.001665 -0.14451 0.14784 SD.DMTA_0 1.545520 0.14301 2.94803 SD.k_DMTA 0.606422 0.26227 0.95057 Correlation: DMTA_0 k_DMTA 0.0169 Random effects: est. lower upper SD.DMTA_0 1.5455 0.1430 2.9480 SD.k_DMTA 0.6064 0.2623 0.9506 Variance model: est. lower upper a.1 2.660081 2.2749 3.0452 b.1 0.001665 -0.1445 0.1478 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.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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 1.408 s Using 300, 100 iterations and 9 chains Variance model: Constant variance Starting values for degradation parameters: DMTA_0 alpha beta 98.292 9.909 156.341 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 alpha beta DMTA_0 98.29 0 0 alpha 0.00 1 0 beta 0.00 0 1 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 734.2 732.7 -360.1 Optimised parameters: est. lower upper DMTA_0 98.3435 96.9033 99.784 alpha 7.2007 2.5889 11.812 beta 112.8745 34.8816 190.867 a.1 2.0459 1.8054 2.286 SD.DMTA_0 1.4795 0.2717 2.687 SD.alpha 0.6396 0.1509 1.128 SD.beta 0.6874 0.1587 1.216 Correlation: DMTA_0 alpha alpha -0.1125 beta -0.1227 0.3632 Random effects: est. lower upper SD.DMTA_0 1.4795 0.2717 2.687 SD.alpha 0.6396 0.1509 1.128 SD.beta 0.6874 0.1587 1.216 Variance model: est. lower upper a.1 2.046 1.805 2.286 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.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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 2.735 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: DMTA_0 alpha beta 98.772 4.663 92.597 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 alpha beta DMTA_0 98.77 0 0 alpha 0.00 1 0 beta 0.00 0 1 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 720.7 719.1 -352.4 Optimised parameters: est. lower upper DMTA_0 99.10577 97.33296 100.87859 alpha 5.46260 2.52199 8.40321 beta 81.66080 30.46664 132.85497 a.1 1.50219 1.25801 1.74636 b.1 0.02893 0.02048 0.03739 SD.DMTA_0 1.61887 -0.03843 3.27618 SD.alpha 0.58145 0.17364 0.98925 SD.beta 0.68205 0.21108 1.15302 Correlation: DMTA_0 alpha alpha -0.1321 beta -0.1430 0.2467 Random effects: est. lower upper SD.DMTA_0 1.6189 -0.03843 3.2762 SD.alpha 0.5814 0.17364 0.9892 SD.beta 0.6821 0.21108 1.1530 Variance model: est. lower upper a.1 1.50219 1.25801 1.74636 b.1 0.02893 0.02048 0.03739 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.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 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * DMTA Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 1.769 s Using 300, 100 iterations and 9 chains Variance model: Constant variance Starting values for degradation parameters: DMTA_0 k1 k2 g 98.64383 0.09211 0.02999 0.76814 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 g DMTA_0 98.64 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 g 0.00 0 0 1 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 711.8 710 -346.9 Optimised parameters: est. lower upper DMTA_0 98.092481 96.573899 99.61106 k1 0.062499 0.030336 0.09466 k2 0.009065 -0.005133 0.02326 g 0.948967 0.862080 1.03586 a.1 1.821671 1.604774 2.03857 SD.DMTA_0 1.677785 0.472066 2.88350 SD.k1 0.634962 0.270788 0.99914 SD.k2 1.033498 -0.205994 2.27299 SD.g 1.710046 0.428642 2.99145 Correlation: DMTA_0 k1 k2 k1 0.0246 k2 0.0491 0.0953 g -0.0552 -0.0889 -0.4795 Random effects: est. lower upper SD.DMTA_0 1.678 0.4721 2.8835 SD.k1 0.635 0.2708 0.9991 SD.k2 1.033 -0.2060 2.2730 SD.g 1.710 0.4286 2.9914 Variance model: est. lower upper a.1 1.822 1.605 2.039 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.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 * time)) / (g * exp(-k1 * time) + (1 - g) * exp(-k2 * time))) * DMTA Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 3.069 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: DMTA_0 k1 k2 g 98.71186 0.08675 0.01374 0.93491 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 g DMTA_0 98.71 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 g 0.00 0 0 1 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 665.7 663.6 -322.9 Optimised parameters: est. lower upper DMTA_0 98.347470 96.380815 100.31413 k1 0.064524 0.034279 0.09477 k2 0.008304 0.005843 0.01076 g 0.952128 0.909578 0.99468 a.1 1.068907 0.883665 1.25415 b.1 0.029265 0.022318 0.03621 SD.DMTA_0 2.065796 0.427951 3.70364 SD.k1 0.583703 0.251796 0.91561 SD.k2 0.004167 -7.831228 7.83956 SD.g 1.064450 0.397479 1.73142 Correlation: DMTA_0 k1 k2 k1 0.0223 k2 0.0568 0.0394 g -0.0464 -0.0269 -0.2713 Random effects: est. lower upper SD.DMTA_0 2.065796 0.4280 3.7036 SD.k1 0.583703 0.2518 0.9156 SD.k2 0.004167 -7.8312 7.8396 SD.g 1.064450 0.3975 1.7314 Variance model: est. lower upper a.1 1.06891 0.88367 1.25415 b.1 0.02927 0.02232 0.03621 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.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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 1.996 s Using 300, 100 iterations and 9 chains Variance model: Constant variance Starting values for degradation parameters: DMTA_0 k1 k2 tb 97.82176 0.06931 0.02997 11.13945 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 tb DMTA_0 97.82 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 tb 0.00 0 0 1 Starting values for error model parameters: a.1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 714 712.1 -348 Optimised parameters: est. lower upper DMTA_0 98.16102 96.47747 99.84456 k1 0.07876 0.05261 0.10491 k2 0.02227 0.01706 0.02747 tb 13.99089 -7.40049 35.38228 a.1 1.82305 1.60700 2.03910 SD.DMTA_0 1.88413 0.56204 3.20622 SD.k1 0.34292 0.10482 0.58102 SD.k2 0.19851 0.01718 0.37985 SD.tb 1.68168 0.58064 2.78272 Correlation: DMTA_0 k1 k2 k1 0.0142 k2 0.0001 -0.0025 tb 0.0165 -0.1256 -0.0301 Random effects: est. lower upper SD.DMTA_0 1.8841 0.56204 3.2062 SD.k1 0.3429 0.10482 0.5810 SD.k2 0.1985 0.01718 0.3798 SD.tb 1.6817 0.58064 2.7827 Variance model: est. lower upper a.1 1.823 1.607 2.039 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.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 Data: 155 observations of 1 variable(s) grouped in 6 datasets Model predictions using solution type analytical Fitted in 3.389 s Using 300, 100 iterations and 9 chains Variance model: Two-component variance function Starting values for degradation parameters: DMTA_0 k1 k2 tb 98.45190 0.07525 0.02576 19.19375 Fixed degradation parameter values: None Starting values for random effects (square root of initial entries in omega): DMTA_0 k1 k2 tb DMTA_0 98.45 0 0 0 k1 0.00 1 0 0 k2 0.00 0 1 0 tb 0.00 0 0 1 Starting values for error model parameters: a.1 b.1 1 1 Results: Likelihood computed by importance sampling AIC BIC logLik 667.1 665 -323.6 Optimised parameters: est. lower upper DMTA_0 97.76571 95.81350 99.71791 k1 0.05855 0.03080 0.08630 k2 0.02337 0.01664 0.03010 tb 31.09638 29.38289 32.80987 a.1 1.08835 0.90059 1.27611 b.1 0.02964 0.02261 0.03667 SD.DMTA_0 2.04877 0.42553 3.67200 SD.k1 0.59166 0.25621 0.92711 SD.k2 0.30698 0.09561 0.51835 SD.tb 0.01274 -0.10915 0.13464 Correlation: DMTA_0 k1 k2 k1 0.0160 k2 -0.0070 -0.0024 tb -0.0668 -0.0103 -0.2013 Random effects: est. lower upper SD.DMTA_0 2.04877 0.42553 3.6720 SD.k1 0.59166 0.25621 0.9271 SD.k2 0.30698 0.09561 0.5183 SD.tb 0.01274 -0.10915 0.1346 Variance model: est. lower upper a.1 1.08835 0.90059 1.27611 b.1 0.02964 0.02261 0.03667 Estimated disappearance times: DT50 DT90 DT50back DT50_k1 DT50_k2 DMTA 11.84 51.71 15.57 11.84 29.66"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"hierarchical-model-convergence-plots","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Hierarchical model convergence plots","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"Convergence plot NLHM SFO fit constant variance Convergence plot NLHM SFO fit two-component error Convergence plot NLHM FOMC fit constant variance Convergence plot NLHM FOMC fit two-component error Convergence plot NLHM DFOP fit constant variance Convergence plot NLHM DFOP fit two-component error Convergence plot NLHM HS fit constant variance Convergence plot NLHM HS fit two-component error","code":""},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"session-info","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Session info","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"","code":"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 version 3.11.0 locale: [1] LC_CTYPE=de_DE.UTF-8 LC_NUMERIC=C [3] LC_TIME=de_DE.UTF-8 LC_COLLATE=de_DE.UTF-8 [5] LC_MONETARY=de_DE.UTF-8 LC_MESSAGES=de_DE.UTF-8 [7] LC_PAPER=de_DE.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=de_DE.UTF-8 LC_IDENTIFICATION=C time zone: Europe/Berlin tzcode source: system (glibc) attached base packages: [1] parallel stats graphics grDevices utils datasets methods [8] base other attached packages: [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] 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"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_parent.html","id":"hardware-info","dir":"Articles > Prebuilt","previous_headings":"Appendix","what":"Hardware info","title":"Testing hierarchical parent degradation kinetics with residue data on dimethenamid and dimethenamid-P","text":"","code":"CPU model: AMD Ryzen 9 7950X 16-Core Processor MemTotal: 64927780 kB"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"introduction","dir":"Articles > Prebuilt","previous_headings":"","what":"Introduction","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"purpose document test demonstrate nonlinear hierarchical models (NLHM) based parent degradation models SFO, FOMC, DFOP HS, parallel formation two metabolites can fitted mkin package. assembled course work package 1.2 Project Number 173340 (Application nonlinear hierarchical models kinetic evaluation chemical degradation data) German Environment Agency carried 2022 2023. mkin package used version 1.2.10, currently development. contains test data, functions used evaluations. saemix package used backend fitting NLHM, also loaded make convergence plot function available. document processed knitr package, also provides kable function used improve display tabular data R markdown documents. parallel processing, parallel package used.","code":"library(mkin) library(knitr) library(saemix) library(parallel) n_cores <- detectCores() # We need to start a new cluster after defining a compiled model that is # saved as a DLL to the user directory, therefore we define a function # This is used again after defining the pathway model start_cluster <- function(n_cores) { if (Sys.info()[\"sysname\"] == \"Windows\") { ret <- makePSOCKcluster(n_cores) } else { ret <- makeForkCluster(n_cores) } return(ret) }"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"data","dir":"Articles > Prebuilt","previous_headings":"","what":"Data","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"test data available mkin package object class mkindsg (mkin dataset group) identifier dimethenamid_2018. following preprocessing steps done document. data available enantiomer dimethenamid-P (DMTAP) renamed substance name data racemic mixture dimethenamid (DMTA). reason difference degradation behaviour identified EU risk assessment. Unnecessary columns discarded observation times dataset multiplied corresponding normalisation factor also available dataset, order make possible describe datasets single set parameters independent temperature Finally, datasets observed soil (Elliot 1 Elliot 2) combined, resulting dimethenamid (DMTA) data six soils. following commented R code performs preprocessing. following tables show 6 datasets. Dataset Calke Dataset Borstel Dataset Flaach Dataset BBA 2.2 Dataset BBA 2.3 Dataset Elliot","code":"# 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, select = 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 for (ds_name in names(dmta_ds)) { print( kable(mkin_long_to_wide(dmta_ds[[ds_name]]), caption = paste(\"Dataset\", ds_name), booktabs = TRUE, row.names = FALSE)) cat(\"\\n\\\\clearpage\\n\") }"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"separate-evaluations","dir":"Articles > Prebuilt","previous_headings":"","what":"Separate evaluations","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"first step obtain suitable starting parameters NLHM fits, separate fits several variants pathway model used previously (Ranke et al. 2021), varying kinetic model parent compound. SFORB model often provides faster convergence DFOP model, can sometimes fitted DFOP model results errors, included set parent models tested . separate pathway fits SFO FOMC parent constant variance converged (status OK). fits DFOP SFORB parent converged well. fits HS parent converge default settings. two-component error model, set fits convergence problems slightly different, convergence problems appearing different data sets applying DFOP SFORB model additional convergence problems using FOMC model parent.","code":"if (!dir.exists(\"dmta_dlls\")) dir.create(\"dmta_dlls\") m_sfo_path_1 <- mkinmod( DMTA = mkinsub(\"SFO\", c(\"M23\", \"M27\", \"M31\")), M23 = mkinsub(\"SFO\"), M27 = mkinsub(\"SFO\"), M31 = mkinsub(\"SFO\", \"M27\", sink = FALSE), name = \"m_sfo_path\", dll_dir = \"dmta_dlls\", unload = TRUE, overwrite = TRUE, quiet = TRUE ) m_fomc_path_1 <- mkinmod( DMTA = mkinsub(\"FOMC\", c(\"M23\", \"M27\", \"M31\")), M23 = mkinsub(\"SFO\"), M27 = mkinsub(\"SFO\"), M31 = mkinsub(\"SFO\", \"M27\", sink = FALSE), name = \"m_fomc_path\", dll_dir = \"dmta_dlls\", unload = TRUE, overwrite = TRUE, quiet = TRUE ) m_dfop_path_1 <- mkinmod( DMTA = mkinsub(\"DFOP\", c(\"M23\", \"M27\", \"M31\")), M23 = mkinsub(\"SFO\"), M27 = mkinsub(\"SFO\"), M31 = mkinsub(\"SFO\", \"M27\", sink = FALSE), name = \"m_dfop_path\", dll_dir = \"dmta_dlls\", unload = TRUE, overwrite = TRUE, quiet = TRUE ) m_sforb_path_1 <- mkinmod( DMTA = mkinsub(\"SFORB\", c(\"M23\", \"M27\", \"M31\")), M23 = mkinsub(\"SFO\"), M27 = mkinsub(\"SFO\"), M31 = mkinsub(\"SFO\", \"M27\", sink = FALSE), name = \"m_sforb_path\", dll_dir = \"dmta_dlls\", unload = TRUE, overwrite = TRUE, quiet = TRUE ) m_hs_path_1 <- mkinmod( DMTA = mkinsub(\"HS\", c(\"M23\", \"M27\", \"M31\")), M23 = mkinsub(\"SFO\"), M27 = mkinsub(\"SFO\"), M31 = mkinsub(\"SFO\", \"M27\", sink = FALSE), name = \"m_hs_path\", dll_dir = \"dmta_dlls\", unload = TRUE, overwrite = TRUE, quiet = TRUE ) cl <- start_cluster(n_cores) deg_mods_1 <- list( sfo_path_1 = m_sfo_path_1, fomc_path_1 = m_fomc_path_1, dfop_path_1 = m_dfop_path_1, sforb_path_1 = m_sforb_path_1, hs_path_1 = m_hs_path_1) sep_1_const <- mmkin( deg_mods_1, dmta_ds, error_model = \"const\", quiet = TRUE) status(sep_1_const) |> kable() sep_1_tc <- update(sep_1_const, error_model = \"tc\") status(sep_1_tc) |> kable()"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"hierarchichal-model-fits","dir":"Articles > Prebuilt","previous_headings":"","what":"Hierarchichal model fits","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"following code fits two sets corresponding hierarchical models data, one assuming constant variance, one assuming two-component error. run time fits around two hours five year old hardware. recent hardware upgrade fits complete less twenty minutes. According status function, fits terminated successfully. goodness--fit models compared, warning obtained, indicating likelihood pathway fit SFORB parent compound constant variance calculated importance sampling (method ‘’). default method AIC BIC comparisons based, variant included model comparison table. Comparing goodness--fit remaining models, HS model model two-component error provides best fit. However, batch experiments performed constant conditions experiments evaluated , reason assume discontinuity, SFORB model preferable mechanistic viewpoint. addition, information criteria AIC BIC similar HS SFORB. Therefore, SFORB model selected refinements.","code":"saem_1 <- mhmkin(list(sep_1_const, sep_1_tc)) status(saem_1) |> kable() anova(saem_1) |> kable(digits = 1)"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"parameter-identifiability-based-on-the-fisher-information-matrix","dir":"Articles > Prebuilt","previous_headings":"Hierarchichal model fits","what":"Parameter identifiability based on the Fisher Information Matrix","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"Using illparms function, ill-defined statistical model parameters standard deviations degradation parameters population error model parameters can found. using constant variance, ill-defined variance parameters identified illparms function degradation models. using two-component error model, one ill-defined variance parameter variants except variant using DFOP parent compound. selected combination SFORB pathway model two-component error, random effect rate constant reversibly bound DMTA free DMTA (k_DMTA_bound_free) well-defined. Therefore, fit updated without assuming random effect parameter. expected, ill-defined parameters remain. model comparison shows reduced model preferable. convergence plot refined fit shown . parameters, example f_DMTA_ilr_1 f_DMTA_ilr_2, .e. two parameters determining formation fractions parallel formation three metabolites, movement parameters still visible second phase algorithm. However, amplitude movement range amplitude towards end first phase. Therefore, likely increase iterations improve parameter estimates much, proposed fit acceptable. numeric convergence criterion implemented saemix.","code":"illparms(saem_1) |> kable() saem_sforb_path_1_tc_reduced <- update(saem_1[[\"sforb_path_1\", \"tc\"]], no_random_effect = \"log_k_DMTA_bound_free\") illparms(saem_sforb_path_1_tc_reduced) anova(saem_1[[\"sforb_path_1\", \"tc\"]], saem_sforb_path_1_tc_reduced) |> kable(digits = 1) plot(saem_sforb_path_1_tc_reduced$so, plot.type = \"convergence\")"},{"path":"https://pkgdown.jrwb.de/mkin/articles/prebuilt/2022_dmta_pathway.html","id":"alternative-check-of-parameter-identifiability","dir":"Articles > Prebuilt","previous_headings":"Hierarchichal model fits","what":"Alternative check of parameter identifiability","title":"Testing hierarchical pathway kinetics with residue data on dimethenamid and dimethenamid-P","text":"alternative check parameter identifiability (Duchesne et al. 2021), multistart runs performed basis refined fit shown . 32 fits initiated, 17 terminated without error. reason wide variation starting parameters combination parameter variation used SAEM algorithm leads parameter combinations degradation model numerical integration routine cope . variation initial parameters, model fits take two times time original fit. Parameter boxplots multistart runs succeeded However, visual analysis boxplot parameters obtained successful fits confirms results sufficiently independent starting parameters, remaining ill-defined parameters.","code":"saem_sforb_path_1_tc_reduced_multi <- multistart(saem_sforb_path_1_tc_reduced, n = 32, cores = 10) (subscript) logical subscript too long print(saem_sforb_path_1_tc_reduced_multi)