From fe411af3dbed474d9e5b3cc1d509f66bd76f7319 Mon Sep 17 00:00:00 2001
From: Johannes Ranke
summary(f_saem_dfop_sfo, data = TRUE)
#> saemix version used for fitting: 3.4
-#> mkin version used for pre-fitting: 1.2.11
-#> R version used for fitting: 4.5.1
-#> Date of fit: Fri Sep 12 20:34:55 2025
-#> Date of summary: Fri Sep 12 20:34:55 2025
+#> mkin version used for pre-fitting: 1.2.10
+#> R version used for fitting: 4.5.2
+#> Date of fit: Fri Nov 28 09:18:42 2025
+#> Date of summary: Fri Nov 28 09:18:42 2025
#>
#> Equations:
#> d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -420,7 +420,7 @@ using mmkin.
#>
#> Model predictions using solution type analytical
#>
-#> Fitted in 4.017 s
+#> Fitted in 3.833 s
#> Using 300, 100 iterations and 10 chains
#>
#> Variance model: Constant variance
@@ -682,28 +682,11 @@ using mmkin.
#> Dataset 10 A1 120 14.3 12.79238 1.50762 1.882 0.80111
#> Dataset 10 A1 120 12.1 12.79238 -0.69238 1.882 -0.36791
-# The following takes about 6 minutes
-f_saem_dfop_sfo_deSolve <- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve",
- nbiter.saemix = c(200, 80))
-#> DINTDY- T (=R1) illegal
-#> In above message, R1 = 70
-#>
-#> T not in interval TCUR - HU (= R1) to TCUR (=R2)
-#> In above message, R1 = 53.1042, R2 = 56.6326
-#>
-#> DINTDY- T (=R1) illegal
-#> In above message, R1 = 91
-#>
-#> T not in interval TCUR - HU (= R1) to TCUR (=R2)
-#> In above message, R1 = 53.1042, R2 = 56.6326
-#>
-#> DLSODA- Trouble in DINTDY. ITASK = I1, TOUT = R1
-#> In above message, I1 = 1
-#>
-#> In above message, R1 = 91
-#>
-#> Error in deSolve::lsoda(y = odeini, times = outtimes_deSolve, func = lsoda_func, :
-#> illegal input detected before taking any integration steps - see written message
+# The following took about 6 minutes and currently fails with
+# "Error in deSolve::lsoda(y = odeini, times = outtimes_deSolve, func = lsoda_func, :
+# illegal input detected before taking any integration steps - see written message"
+#f_saem_dfop_sfo_deSolve <- saem(f_mmkin["DFOP-SFO", ], solution_type = "deSolve",
+# nbiter.saemix = c(200, 80))
#anova(
# f_saem_dfop_sfo,
--
cgit v1.2.3