From cf54ccca37d27480dbf8d59eb027300518f7ad75 Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 19 May 2023 17:08:43 +0200 Subject: Prepare release of v1.2.4 - Update DESCRIPTION - Update Makefile to document how to use R-patched - Remove markup from two URLs to avoid CRAN NOTE - Switch two vignettes from html_document to html_vignette to save space in the docs directory, also avoiding a CRAN NOTE - Complete rebuild of pkgdown docs for release --- docs/articles/prebuilt/2022_cyan_pathway.html | 3262 +++++++++----------- .../figure-html/unnamed-chunk-13-1.png | Bin 363943 -> 363589 bytes .../figure-html/unnamed-chunk-14-1.png | Bin 365867 -> 364310 bytes .../figure-html/unnamed-chunk-15-1.png | Bin 363662 -> 365199 bytes .../figure-html/unnamed-chunk-20-1.png | Bin 378667 -> 378245 bytes .../figure-html/unnamed-chunk-21-1.png | Bin 372548 -> 371558 bytes .../figure-html/unnamed-chunk-22-1.png | Bin 373913 -> 373392 bytes .../figure-html/unnamed-chunk-7-1.png | Bin 373737 -> 373137 bytes .../figure-html/unnamed-chunk-8-1.png | Bin 373913 -> 373392 bytes docs/articles/prebuilt/2022_dmta_parent.html | 6 +- docs/articles/prebuilt/2022_dmta_pathway.html | 168 +- .../saem-sforb-path-1-tc-reduced-convergence-1.png | Bin 160921 -> 160917 bytes .../figure-html/unnamed-chunk-2-1.png | Bin 108560 -> 114625 bytes .../figure-html/unnamed-chunk-3-1.png | Bin 403164 -> 403548 bytes .../figure-html/unnamed-chunk-4-1.png | Bin 403636 -> 401796 bytes .../figure-html/unnamed-chunk-6-1.png | Bin 393985 -> 394501 bytes .../figure-html/unnamed-chunk-7-1.png | Bin 400319 -> 399532 bytes .../figure-html/unnamed-chunk-8-1.png | Bin 403164 -> 403548 bytes 18 files changed, 1574 insertions(+), 1862 deletions(-) (limited to 'docs/articles/prebuilt') diff --git a/docs/articles/prebuilt/2022_cyan_pathway.html b/docs/articles/prebuilt/2022_cyan_pathway.html index cd63fa3c..c22c6735 100644 --- a/docs/articles/prebuilt/2022_cyan_pathway.html +++ b/docs/articles/prebuilt/2022_cyan_pathway.html @@ -33,7 +33,7 @@ mkin - 1.2.3.1 + 1.2.4 @@ -135,7 +135,7 @@ residue data on cyantraniliprole Ranke

Last change on 20 April 2023, -last compiled on 20 April 2023

+last compiled on 19 Mai 2023 Source: vignettes/prebuilt/2022_cyan_pathway.rmd @@ -155,7 +155,7 @@ be fitted with the mkin package.

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

-

The mkin package is used in version 1.2.3 which is currently under +

The mkin package is used in version 1.2.4 which is currently under development. The newly introduced functionality that is used here is a simplification of excluding random effects for a set of fits based on a related set of fits with a reduced model, and the documentation of the @@ -1356,9 +1356,9 @@ when the two-component error model is used.

SFORB tc 10 -700.1 -696.2 --340.1 +700.0 +696.1 +-340.0 HS tc @@ -1447,7 +1447,7 @@ pathway fits are performed for all datasets.

OK OK OK -OK +C OK @@ -1511,7 +1511,7 @@ pathway fits are performed for all datasets.

OK OK OK -C +OK dfop_path_1 @@ -1524,7 +1524,7 @@ pathway fits are performed for all datasets.

sforb_path_1 OK -C +OK OK OK OK @@ -1535,7 +1535,7 @@ pathway fits are performed for all datasets.

OK C OK -OK +C @@ -1564,7 +1564,7 @@ which is currently under development.

sfo_path_1 -Fth, FO +FO Fth, FO @@ -1584,8 +1584,8 @@ which is currently under development.

hs_path_1 -Fth, FO -Fth, FO +FO +E @@ -1634,7 +1634,7 @@ sd(f_JCZ38_qlogis) hs_path_1 NA -NA +E @@ -1654,30 +1654,30 @@ or SFORB for the parent compound provide the best fit.

sfo_path_1 const 16 -2692.8 -2686.6 --1330.4 +2693.0 +2686.8 +-1330.5 sfo_path_1 tc 17 -2657.7 -2651.1 --1311.9 +2657.6 +2651.0 +-1311.8 fomc_path_1 const 18 -2427.8 -2420.8 --1195.9 +2427.9 +2420.9 +-1196.0 fomc_path_1 tc 19 -2423.4 -2416.0 --1192.7 +2423.6 +2416.2 +-1192.8 dfop_path_1 const @@ -1696,30 +1696,23 @@ or SFORB for the parent compound provide the best fit.

hs_path_1 const 20 -2427.3 -2419.5 --1193.7 +2427.2 +2419.4 +-1193.6 dfop_path_1 tc 20 2398.0 -2390.2 +2390.1 -1179.0 sforb_path_1 tc 20 -2399.8 -2392.0 --1179.9 - - -hs_path_1 tc -21 -2422.3 -2414.1 --1190.2 +2399.9 +2392.1 +-1180.0 @@ -1827,7 +1820,7 @@ alternative pathway fits.

OK OK OK -C +OK OK @@ -1850,7 +1843,7 @@ the fits to the Sassafras soil data.

fomc_path_2 OK -C +OK OK C OK @@ -1858,7 +1851,7 @@ the fits to the Sassafras soil data.

dfop_path_2 OK -OK +C OK C OK @@ -1868,7 +1861,7 @@ the fits to the Sassafras soil data.

OK OK OK -OK +C OK @@ -1891,8 +1884,8 @@ parent and the Sassafras dataset.

fomc_path_2 +E OK -FO dfop_path_2 @@ -1924,8 +1917,8 @@ successfully.

fomc_path_2 -sd(f_JCZ38_qlogis), sd(f_JSE76_qlogis) -NA +E +sd(f_JSE76_qlogis) dfop_path_2 @@ -1954,46 +1947,39 @@ JCZ38 are ill-defined.

-fomc_path_2 const -20 -2308.3 -2300.5 --1134.2 - - fomc_path_2 tc 21 -2248.3 -2240.1 --1103.2 +2249.0 +2240.8 +-1103.5 - + dfop_path_2 const 22 -2289.6 -2281.0 --1122.8 +2288.4 +2279.8 +-1122.2 - + sforb_path_2 const 22 -2284.1 -2275.5 --1120.0 +2283.3 +2274.7 +-1119.7 - + dfop_path_2 tc 22 2234.4 2225.8 -1095.2 - + sforb_path_2 tc 22 -2240.4 -2231.8 --1098.2 +2239.7 +2231.1 +-1097.9 @@ -2120,37 +2106,37 @@ confidence intervals for the optimised parameters are available.

fomc_path_2 tc 19 -2250.9 -2243.5 --1106.5 +2249.1 +2241.6 +-1105.5 dfop_path_2 const 20 -2281.7 -2273.9 --1120.8 +2282.2 +2274.4 +-1121.1 sforb_path_2 const 20 -2279.5 -2271.7 --1119.7 +2279.7 +2271.9 +-1119.9 dfop_path_2 tc 20 -2231.5 -2223.7 --1095.8 +2237.3 +2229.5 +-1098.6 sforb_path_2 tc 20 -2235.7 -2227.9 --1097.9 +2241.3 +2233.5 +-1100.7 @@ -2219,10 +2205,10 @@ Hierarchical SFO path 1 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:44:55 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:27:54 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - k_cyan * cyan
@@ -2235,7 +2221,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 431.793 s
+Fitted in 422.743 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Constant variance 
@@ -2244,7 +2230,7 @@ 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        22.5404        10.4289 
+        0.8158        23.5335        11.8774 
 
 Fixed degradation parameter values:
 None
@@ -2266,8 +2252,8 @@ 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        8.684           0.00
-f_JCZ38_qlogis       0.0000        0.000          13.48
+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 
@@ -2277,68 +2263,68 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2693 2687  -1330
+  2693 2687  -1331
 
 Optimised parameters:
-                     est. lower upper
-cyan_0            95.0946    NA    NA
-log_k_cyan        -3.8544    NA    NA
-log_k_JCZ38       -3.0402    NA    NA
-log_k_J9Z38       -5.0109    NA    NA
-log_k_JSE76       -5.2857    NA    NA
-f_cyan_ilr_1       0.8069    NA    NA
-f_cyan_ilr_2      16.6623    NA    NA
-f_JCZ38_qlogis     1.3602    NA    NA
-a.1                4.8326    NA    NA
-SD.log_k_cyan      0.5842    NA    NA
-SD.log_k_JCZ38     1.2680    NA    NA
-SD.log_k_J9Z38     0.3626    NA    NA
-SD.log_k_JSE76     0.5244    NA    NA
-SD.f_cyan_ilr_1    0.2752    NA    NA
-SD.f_cyan_ilr_2    2.3556    NA    NA
-SD.f_JCZ38_qlogis  0.2400    NA    NA
+                     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.5842    NA    NA
-SD.log_k_JCZ38    1.2680    NA    NA
-SD.log_k_J9Z38    0.3626    NA    NA
-SD.log_k_JSE76    0.5244    NA    NA
-SD.f_cyan_ilr_1   0.2752    NA    NA
-SD.f_cyan_ilr_2   2.3556    NA    NA
-SD.f_JCZ38_qlogis 0.2400    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
 
 Variance model:
      est. lower upper
-a.1 4.833    NA    NA
+a.1 4.839    NA    NA
 
 Backtransformed parameters:
-                      est. lower upper
-cyan_0           95.094581    NA    NA
-k_cyan            0.021186    NA    NA
-k_JCZ38           0.047825    NA    NA
-k_J9Z38           0.006665    NA    NA
-k_JSE76           0.005063    NA    NA
-f_cyan_to_JCZ38   0.757885    NA    NA
-f_cyan_to_J9Z38   0.242115    NA    NA
-f_JCZ38_to_JSE76  0.795792    NA    NA
+                      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.579e-01
-cyan_J9Z38  2.421e-01
-cyan_sink   5.877e-10
-JCZ38_JSE76 7.958e-01
-JCZ38_sink  2.042e-01
+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.72 108.68
-JCZ38  14.49  48.15
-J9Z38 103.99 345.46
-JSE76 136.90 454.76
+cyan   32.66 108.50
+JCZ38  14.46  48.05
+J9Z38 103.86 345.00
+JSE76 143.91 478.04

@@ -2347,10 +2333,10 @@ Hierarchical SFO path 1 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:44:53 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:27:49 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - k_cyan * cyan
@@ -2363,7 +2349,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 429.526 s
+Fitted in 417.436 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Two-component variance function 
@@ -2372,7 +2358,7 @@ 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        20.0030        15.1336 
+        0.7937        22.3422        17.8932 
 
 Fixed degradation parameter values:
 None
@@ -2394,8 +2380,8 @@ 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        7.334           0.00
-f_JCZ38_qlogis       0.0000        0.000          11.78
+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 
@@ -2409,66 +2395,66 @@ Likelihood computed by importance sampling
 
 Optimised parameters:
                       est. lower upper
-cyan_0            94.72923    NA    NA
-log_k_cyan        -3.91670    NA    NA
-log_k_JCZ38       -3.12917    NA    NA
-log_k_J9Z38       -5.06070    NA    NA
-log_k_JSE76       -5.09254    NA    NA
-f_cyan_ilr_1       0.81116    NA    NA
-f_cyan_ilr_2      39.97850    NA    NA
-f_JCZ38_qlogis     3.09728    NA    NA
-a.1                3.95044    NA    NA
-b.1                0.07998    NA    NA
-SD.log_k_cyan      0.58855    NA    NA
-SD.log_k_JCZ38     1.29753    NA    NA
-SD.log_k_J9Z38     0.62851    NA    NA
-SD.log_k_JSE76     0.37235    NA    NA
-SD.f_cyan_ilr_1    0.37346    NA    NA
-SD.f_cyan_ilr_2    1.41667    NA    NA
-SD.f_JCZ38_qlogis  1.81467    NA    NA
+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.5886    NA    NA
-SD.log_k_JCZ38    1.2975    NA    NA
-SD.log_k_J9Z38    0.6285    NA    NA
-SD.log_k_JSE76    0.3724    NA    NA
-SD.f_cyan_ilr_1   0.3735    NA    NA
-SD.f_cyan_ilr_2   1.4167    NA    NA
-SD.f_JCZ38_qlogis 1.8147    NA    NA
+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.95044    NA    NA
-b.1 0.07998    NA    NA
+a.1 3.93879    NA    NA
+b.1 0.08057    NA    NA
 
 Backtransformed parameters:
-                      est. lower upper
-cyan_0           94.729229    NA    NA
-k_cyan            0.019907    NA    NA
-k_JCZ38           0.043754    NA    NA
-k_J9Z38           0.006341    NA    NA
-k_JSE76           0.006142    NA    NA
-f_cyan_to_JCZ38   0.758991    NA    NA
-f_cyan_to_J9Z38   0.241009    NA    NA
-f_JCZ38_to_JSE76  0.956781    NA    NA
+                     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.75899
-cyan_J9Z38  0.24101
+cyan_JCZ38  0.75807
+cyan_J9Z38  0.24193
 cyan_sink   0.00000
-JCZ38_JSE76 0.95678
-JCZ38_sink  0.04322
+JCZ38_JSE76 0.95409
+JCZ38_sink  0.04591
 
 Estimated disappearance times:
         DT50   DT90
-cyan   34.82 115.67
-JCZ38  15.84  52.63
-J9Z38 109.31 363.12
-JSE76 112.85 374.87
+cyan   34.78 115.54
+JCZ38  15.81  52.52
+J9Z38 107.97 358.68
+JSE76 114.20 379.35

@@ -2477,10 +2463,10 @@ Hierarchical FOMC path 1 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:45:50 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:28:29 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan
@@ -2495,7 +2481,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 477.996 s
+Fitted in 457.122 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Constant variance 
@@ -2504,7 +2490,7 @@ 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.3870        15.7604        -0.1791         2.9811 
+        7.0229        14.9234        -0.1791         2.9811 
 
 Fixed degradation parameter values:
 None
@@ -2526,8 +2512,8 @@ 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.33           0.00    0.0000   0.0000
-f_JCZ38_qlogis         0.00          20.42    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
 
@@ -2542,78 +2528,78 @@ Likelihood computed by importance sampling
   2428 2421  -1196
 
 Optimised parameters:
-                      est.     lower    upper
-cyan_0            101.0225 98.306270 103.7387
-log_k_JCZ38        -3.3786 -4.770657  -1.9866
-log_k_J9Z38        -5.2603 -5.902085  -4.6186
-log_k_JSE76        -6.1427 -7.318336  -4.9671
-f_cyan_ilr_1        0.7437  0.421215   1.0663
-f_cyan_ilr_2        0.9108  0.267977   1.5537
-f_JCZ38_qlogis      2.0487  0.524897   3.5724
-log_alpha          -0.2268 -0.618049   0.1644
-log_beta            2.8986  2.700701   3.0964
-a.1                 3.4058  3.169913   3.6416
-SD.cyan_0           2.5279  0.454190   4.6016
-SD.log_k_JCZ38      1.5636  0.572824   2.5543
-SD.log_k_J9Z38      0.5316 -0.004405   1.0677
-SD.log_k_JSE76      0.9903  0.106325   1.8742
-SD.f_cyan_ilr_1     0.3464  0.112066   0.5807
-SD.f_cyan_ilr_2     0.2804 -0.393900   0.9546
-SD.f_JCZ38_qlogis   0.9416 -0.152986   2.0362
-SD.log_alpha        0.4273  0.161044   0.6936
+                      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.0156                                                        
-log_k_J9Z38    -0.0493  0.0073                                                
-log_k_JSE76    -0.0329  0.0018  0.0069                                        
-f_cyan_ilr_1   -0.0086  0.0180 -0.1406  0.0012                                
-f_cyan_ilr_2   -0.2629  0.0779  0.2826  0.0274  0.0099                        
-f_JCZ38_qlogis  0.0713 -0.0747 -0.0505  0.1169 -0.1022 -0.4893                
-log_alpha      -0.0556  0.0120  0.0336  0.0193  0.0036  0.0840 -0.0489        
-log_beta       -0.2898  0.0460  0.1305  0.0768  0.0190  0.4071 -0.1981  0.2772
+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.5279  0.454190 4.6016
-SD.log_k_JCZ38    1.5636  0.572824 2.5543
-SD.log_k_J9Z38    0.5316 -0.004405 1.0677
-SD.log_k_JSE76    0.9903  0.106325 1.8742
-SD.f_cyan_ilr_1   0.3464  0.112066 0.5807
-SD.f_cyan_ilr_2   0.2804 -0.393900 0.9546
-SD.f_JCZ38_qlogis 0.9416 -0.152986 2.0362
-SD.log_alpha      0.4273  0.161044 0.6936
+                    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.406  3.17 3.642
+a.1 3.416  3.18 3.652
 
 Backtransformed parameters:
                       est.     lower     upper
-cyan_0           1.010e+02 9.831e+01 1.037e+02
-k_JCZ38          3.409e-02 8.475e-03 1.372e-01
-k_J9Z38          5.194e-03 2.734e-03 9.867e-03
-k_JSE76          2.149e-03 6.633e-04 6.963e-03
-f_cyan_to_JCZ38  6.481e-01        NA        NA
-f_cyan_to_J9Z38  2.264e-01        NA        NA
-f_JCZ38_to_JSE76 8.858e-01 6.283e-01 9.727e-01
-alpha            7.971e-01 5.390e-01 1.179e+00
-beta             1.815e+01 1.489e+01 2.212e+01
+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.6481
-cyan_J9Z38  0.2264
-cyan_sink   0.1255
-JCZ38_JSE76 0.8858
-JCZ38_sink  0.1142
+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.15  308.01    92.72
-JCZ38  20.33   67.54       NA
-J9Z38 133.46  443.35       NA
-JSE76 322.53 1071.42       NA
+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

@@ -2622,10 +2608,10 @@ Hierarchical FOMC path 1 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:45:45 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:28:21 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan
@@ -2640,41 +2626,41 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 480.648 s
+Fitted in 449.531 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.13827       -3.32493       -5.08921       -5.93478        0.71330 
+     101.13294       -3.32499       -5.09097       -5.93566        0.71359 
   f_cyan_ilr_2 f_JCZ38_qlogis      log_alpha       log_beta 
-      10.05989       12.79248       -0.09621        3.10646 
+      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.643       0.000       0.000        0.00       0.0000
-log_k_JCZ38     0.000       2.319       0.000        0.00       0.0000
-log_k_J9Z38     0.000       0.000       1.731        0.00       0.0000
-log_k_JSE76     0.000       0.000       0.000        1.86       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000        0.00       0.7186
-f_cyan_ilr_2    0.000       0.000       0.000        0.00       0.0000
-f_JCZ38_qlogis  0.000       0.000       0.000        0.00       0.0000
-log_alpha       0.000       0.000       0.000        0.00       0.0000
-log_beta        0.000       0.000       0.000        0.00       0.0000
+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.49           0.00    0.0000   0.0000
-f_JCZ38_qlogis         0.00          20.19    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.7331
+log_beta               0.00           0.00    0.0000   0.7333
 
 Starting values for error model parameters:
 a.1 b.1 
@@ -2684,74 +2670,74 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2423 2416  -1193
+  2424 2416  -1193
 
 Optimised parameters:
                        est. lower upper
-cyan_0            100.57649    NA    NA
-log_k_JCZ38        -3.46250    NA    NA
-log_k_J9Z38        -5.24442    NA    NA
-log_k_JSE76        -5.75229    NA    NA
-f_cyan_ilr_1        0.68480    NA    NA
-f_cyan_ilr_2        0.61670    NA    NA
-f_JCZ38_qlogis     87.97407    NA    NA
-log_alpha          -0.15699    NA    NA
-log_beta            3.01540    NA    NA
-a.1                 3.11518    NA    NA
-b.1                 0.04445    NA    NA
-SD.log_k_JCZ38      1.40732    NA    NA
-SD.log_k_J9Z38      0.56510    NA    NA
-SD.log_k_JSE76      0.72067    NA    NA
-SD.f_cyan_ilr_1     0.31199    NA    NA
-SD.f_cyan_ilr_2     0.36894    NA    NA
-SD.f_JCZ38_qlogis   6.92892    NA    NA
-SD.log_alpha        0.25662    NA    NA
-SD.log_beta         0.35845    NA    NA
+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.4073    NA    NA
-SD.log_k_J9Z38    0.5651    NA    NA
-SD.log_k_JSE76    0.7207    NA    NA
-SD.f_cyan_ilr_1   0.3120    NA    NA
-SD.f_cyan_ilr_2   0.3689    NA    NA
-SD.f_JCZ38_qlogis 6.9289    NA    NA
-SD.log_alpha      0.2566    NA    NA
-SD.log_beta       0.3585    NA    NA
+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.11518    NA    NA
-b.1 0.04445    NA    NA
+a.1 3.11051    NA    NA
+b.1 0.04508    NA    NA
 
 Backtransformed parameters:
                       est. lower upper
-cyan_0           1.006e+02    NA    NA
-k_JCZ38          3.135e-02    NA    NA
-k_J9Z38          5.277e-03    NA    NA
-k_JSE76          3.175e-03    NA    NA
-f_cyan_to_JCZ38  5.991e-01    NA    NA
-f_cyan_to_J9Z38  2.275e-01    NA    NA
+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.547e-01    NA    NA
-beta             2.040e+01    NA    NA
+alpha            8.652e-01    NA    NA
+beta             2.074e+01    NA    NA
 
 Resulting formation fractions:
                 ff
-cyan_JCZ38  0.5991
-cyan_J9Z38  0.2275
-cyan_sink   0.1734
+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.50 281.29    84.68
-JCZ38  22.11  73.44       NA
-J9Z38 131.36 436.35       NA
-JSE76 218.28 725.11       NA
+        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

@@ -2760,10 +2746,10 @@ Hierarchical DFOP path 1 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:46:41 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:29:15 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -2782,16 +2768,16 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 528.713 s
+Fitted in 503.737 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.0644        -3.4008        -5.0024        -5.8613         0.6855 
+      102.0643        -3.4008        -5.0024        -5.8612         0.6855 
   f_cyan_ilr_2 f_JCZ38_qlogis         log_k1         log_k2       g_qlogis 
-        1.2365        13.7245        -1.8641        -4.5063        -0.6468 
+        1.2366        13.6901        -1.8641        -4.5063        -0.6468 
 
 Fixed degradation parameter values:
 None
@@ -2815,7 +2801,7 @@ 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.11 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
@@ -2832,72 +2818,72 @@ Likelihood computed by importance sampling
 
 Optimised parameters:
                       est. lower upper
-cyan_0            102.6079    NA    NA
-log_k_JCZ38        -3.4855    NA    NA
-log_k_J9Z38        -5.1686    NA    NA
-log_k_JSE76        -5.6697    NA    NA
-f_cyan_ilr_1        0.6714    NA    NA
-f_cyan_ilr_2        0.4986    NA    NA
-f_JCZ38_qlogis     55.4760    NA    NA
-log_k1             -1.8409    NA    NA
-log_k2             -4.4915    NA    NA
-g_qlogis           -0.6403    NA    NA
-a.1                 3.2387    NA    NA
-SD.log_k_JCZ38      1.4524    NA    NA
-SD.log_k_J9Z38      0.5151    NA    NA
-SD.log_k_JSE76      0.6514    NA    NA
-SD.f_cyan_ilr_1     0.3023    NA    NA
-SD.f_cyan_ilr_2     0.2959    NA    NA
-SD.f_JCZ38_qlogis   1.9984    NA    NA
-SD.log_k1           0.5188    NA    NA
-SD.log_k2           0.3894    NA    NA
-SD.g_qlogis         0.8579    NA    NA
+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.4524    NA    NA
-SD.log_k_J9Z38    0.5151    NA    NA
-SD.log_k_JSE76    0.6514    NA    NA
-SD.f_cyan_ilr_1   0.3023    NA    NA
-SD.f_cyan_ilr_2   0.2959    NA    NA
-SD.f_JCZ38_qlogis 1.9984    NA    NA
-SD.log_k1         0.5188    NA    NA
-SD.log_k2         0.3894    NA    NA
-SD.g_qlogis       0.8579    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
 
 Variance model:
-     est. lower upper
-a.1 3.239    NA    NA
+    est. lower upper
+a.1 3.24    NA    NA
 
 Backtransformed parameters:
                       est. lower upper
 cyan_0           1.026e+02    NA    NA
-k_JCZ38          3.064e-02    NA    NA
-k_J9Z38          5.692e-03    NA    NA
-k_JSE76          3.449e-03    NA    NA
-f_cyan_to_JCZ38  5.798e-01    NA    NA
-f_cyan_to_J9Z38  2.243e-01    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.587e-01    NA    NA
-k2               1.120e-02    NA    NA
-g                3.452e-01    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.5798
-cyan_J9Z38  0.2243
-cyan_sink   0.1958
+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.21 167.73    50.49   4.368   61.87
-JCZ38  22.62  75.15       NA      NA      NA
-J9Z38 121.77 404.50       NA      NA      NA
-JSE76 200.98 667.64       NA      NA      NA
+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

@@ -2906,10 +2892,10 @@ Hierarchical DFOP path 1 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:49:05 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:31:24 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -2928,7 +2914,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 673.139 s
+Fitted in 632.55 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Two-component variance function 
@@ -2937,7 +2923,7 @@ 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.7799        13.7245        -1.9222        -4.5035        -0.7172 
+        6.8713        13.6901        -1.9222        -4.5035        -0.7172 
 
 Fixed degradation parameter values:
 None
@@ -2948,7 +2934,7 @@ 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.6838
+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
@@ -2960,8 +2946,8 @@ 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.77           0.00 0.0000 0.0000    0.000
-f_JCZ38_qlogis         0.00          16.11 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
@@ -2977,73 +2963,73 @@ Likelihood computed by importance sampling
   2398 2390  -1179
 
 Optimised parameters:
-                      est. lower upper
-cyan_0            100.8076    NA    NA
-log_k_JCZ38        -3.4684    NA    NA
-log_k_J9Z38        -5.0844    NA    NA
-log_k_JSE76        -5.5743    NA    NA
-f_cyan_ilr_1        0.6669    NA    NA
-f_cyan_ilr_2        0.7912    NA    NA
-f_JCZ38_qlogis     84.1825    NA    NA
-log_k1             -2.1671    NA    NA
-log_k2             -4.5447    NA    NA
-g_qlogis           -0.5631    NA    NA
-a.1                 2.9627    NA    NA
-b.1                 0.0444    NA    NA
-SD.log_k_JCZ38      1.4044    NA    NA
-SD.log_k_J9Z38      0.6410    NA    NA
-SD.log_k_JSE76      0.5391    NA    NA
-SD.f_cyan_ilr_1     0.3203    NA    NA
-SD.f_cyan_ilr_2     0.5038    NA    NA
-SD.f_JCZ38_qlogis   3.5865    NA    NA
-SD.log_k2           0.3119    NA    NA
-SD.g_qlogis         0.8276    NA    NA
+                       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.4044    NA    NA
-SD.log_k_J9Z38    0.6410    NA    NA
-SD.log_k_JSE76    0.5391    NA    NA
-SD.f_cyan_ilr_1   0.3203    NA    NA
-SD.f_cyan_ilr_2   0.5038    NA    NA
-SD.f_JCZ38_qlogis 3.5865    NA    NA
-SD.log_k2         0.3119    NA    NA
-SD.g_qlogis       0.8276    NA    NA
+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.9627    NA    NA
-b.1 0.0444    NA    NA
+       est. lower upper
+a.1 2.95785    NA    NA
+b.1 0.04456    NA    NA
 
 Backtransformed parameters:
                       est. lower upper
-cyan_0           1.008e+02    NA    NA
-k_JCZ38          3.117e-02    NA    NA
-k_J9Z38          6.193e-03    NA    NA
-k_JSE76          3.794e-03    NA    NA
-f_cyan_to_JCZ38  6.149e-01    NA    NA
-f_cyan_to_J9Z38  2.395e-01    NA    NA
+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.145e-01    NA    NA
-k2               1.062e-02    NA    NA
-g                3.628e-01    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.6149
-cyan_J9Z38  0.2395
-cyan_sink   0.1456
+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.26 174.32    52.47   6.053   65.25
-JCZ38  22.24  73.88       NA      NA      NA
-J9Z38 111.93 371.82       NA      NA      NA
-JSE76 182.69 606.88       NA      NA      NA
+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

@@ -3052,10 +3038,10 @@ Hierarchical SFORB path 1 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:46:35 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:29:23 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -3073,7 +3059,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 531.17 s
+Fitted in 511.715 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Constant variance 
@@ -3086,7 +3072,7 @@ log_k_cyan_bound_free           log_k_JCZ38           log_k_J9Z38
           log_k_JSE76          f_cyan_ilr_1          f_cyan_ilr_2 
               -5.8613                0.6855                1.2366 
        f_JCZ38_qlogis 
-              13.7418 
+              13.7395 
 
 Fixed degradation parameter values:
 None
@@ -3124,7 +3110,7 @@ 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.14
+f_JCZ38_qlogis              0.0000         0.00          16.13
 
 Starting values for error model parameters:
 a.1 
@@ -3138,77 +3124,77 @@ Likelihood computed by importance sampling
 
 Optimised parameters:
                              est. lower upper
-cyan_free_0              102.7803    NA    NA
-log_k_cyan_free           -2.8068    NA    NA
-log_k_cyan_free_bound     -2.5714    NA    NA
-log_k_cyan_bound_free     -3.4426    NA    NA
-log_k_JCZ38               -3.4994    NA    NA
-log_k_J9Z38               -5.1148    NA    NA
-log_k_JSE76               -5.6335    NA    NA
-f_cyan_ilr_1               0.6597    NA    NA
-f_cyan_ilr_2               0.5132    NA    NA
-f_JCZ38_qlogis            37.2090    NA    NA
-a.1                        3.2367    NA    NA
-SD.log_k_cyan_free         0.3161    NA    NA
-SD.log_k_cyan_free_bound   0.8103    NA    NA
-SD.log_k_cyan_bound_free   0.5554    NA    NA
-SD.log_k_JCZ38             1.4858    NA    NA
-SD.log_k_J9Z38             0.5859    NA    NA
-SD.log_k_JSE76             0.6195    NA    NA
-SD.f_cyan_ilr_1            0.3118    NA    NA
-SD.f_cyan_ilr_2            0.3344    NA    NA
-SD.f_JCZ38_qlogis          0.5518    NA    NA
+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.3161    NA    NA
-SD.log_k_cyan_free_bound 0.8103    NA    NA
-SD.log_k_cyan_bound_free 0.5554    NA    NA
-SD.log_k_JCZ38           1.4858    NA    NA
-SD.log_k_J9Z38           0.5859    NA    NA
-SD.log_k_JSE76           0.6195    NA    NA
-SD.f_cyan_ilr_1          0.3118    NA    NA
-SD.f_cyan_ilr_2          0.3344    NA    NA
-SD.f_JCZ38_qlogis        0.5518    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
 
 Variance model:
      est. lower upper
-a.1 3.237    NA    NA
+a.1 3.232    NA    NA
 
 Backtransformed parameters:
                           est. lower upper
 cyan_free_0          1.028e+02    NA    NA
-k_cyan_free          6.040e-02    NA    NA
-k_cyan_free_bound    7.643e-02    NA    NA
-k_cyan_bound_free    3.198e-02    NA    NA
-k_JCZ38              3.022e-02    NA    NA
-k_J9Z38              6.007e-03    NA    NA
-k_JSE76              3.576e-03    NA    NA
-f_cyan_free_to_JCZ38 5.787e-01    NA    NA
-f_cyan_free_to_J9Z38 2.277e-01    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.15646 0.01235 0.33341 
+0.15973 0.01241 0.33124 
 
 Resulting formation fractions:
                     ff
-cyan_free_JCZ38 0.5787
-cyan_free_J9Z38 0.2277
-cyan_free_sink  0.1936
+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.48 153.7    46.26         4.43        56.15
-JCZ38  22.94  76.2       NA           NA           NA
-J9Z38 115.39 383.3       NA           NA           NA
-JSE76 193.84 643.9       NA           NA           NA
+        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

@@ -3217,10 +3203,10 @@ Hierarchical SFORB path 1 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:49:08 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:31:23 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -3238,7 +3224,7 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 675.659 s
+Fitted in 630.627 s
 Using 300, 100 iterations and 10 chains
 
 Variance model: Two-component variance function 
@@ -3249,9 +3235,9 @@ Starting values for degradation parameters:
 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.8139 
+              -5.8727                0.6814                6.7399 
        f_JCZ38_qlogis 
-              13.7419 
+              13.7395 
 
 Fixed degradation parameter values:
 None
@@ -3288,8 +3274,8 @@ 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.84           0.00
-f_JCZ38_qlogis              0.0000         0.00          16.14
+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 
@@ -3303,77 +3289,77 @@ Likelihood computed by importance sampling
 
 Optimised parameters:
                               est. lower upper
-cyan_free_0              100.69983    NA    NA
-log_k_cyan_free           -3.11584    NA    NA
-log_k_cyan_free_bound     -3.15216    NA    NA
-log_k_cyan_bound_free     -3.65986    NA    NA
-log_k_JCZ38               -3.47811    NA    NA
-log_k_J9Z38               -5.08835    NA    NA
-log_k_JSE76               -5.55514    NA    NA
-f_cyan_ilr_1               0.66764    NA    NA
-f_cyan_ilr_2               0.78329    NA    NA
-f_JCZ38_qlogis            25.35245    NA    NA
-a.1                        2.99088    NA    NA
-b.1                        0.04346    NA    NA
-SD.log_k_cyan_free         0.48797    NA    NA
-SD.log_k_cyan_bound_free   0.27243    NA    NA
-SD.log_k_JCZ38             1.42450    NA    NA
-SD.log_k_J9Z38             0.63496    NA    NA
-SD.log_k_JSE76             0.55951    NA    NA
-SD.f_cyan_ilr_1            0.32687    NA    NA
-SD.f_cyan_ilr_2            0.48056    NA    NA
-SD.f_JCZ38_qlogis          0.43818    NA    NA
+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.4880    NA    NA
-SD.log_k_cyan_bound_free 0.2724    NA    NA
-SD.log_k_JCZ38           1.4245    NA    NA
-SD.log_k_J9Z38           0.6350    NA    NA
-SD.log_k_JSE76           0.5595    NA    NA
-SD.f_cyan_ilr_1          0.3269    NA    NA
-SD.f_cyan_ilr_2          0.4806    NA    NA
-SD.f_JCZ38_qlogis        0.4382    NA    NA
+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.99088    NA    NA
-b.1 0.04346    NA    NA
+a.1 2.98352    NA    NA
+b.1 0.04388    NA    NA
 
 Backtransformed parameters:
                           est. lower upper
-cyan_free_0          1.007e+02    NA    NA
-k_cyan_free          4.434e-02    NA    NA
-k_cyan_free_bound    4.276e-02    NA    NA
-k_cyan_bound_free    2.574e-02    NA    NA
-k_JCZ38              3.087e-02    NA    NA
-k_J9Z38              6.168e-03    NA    NA
-k_JSE76              3.868e-03    NA    NA
-f_cyan_free_to_JCZ38 6.143e-01    NA    NA
-f_cyan_free_to_J9Z38 2.389e-01    NA    NA
+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.10161 0.01123 0.36636 
+0.10044 0.01124 0.36580 
 
 Resulting formation fractions:
-                       ff
-cyan_free_JCZ38 6.143e-01
-cyan_free_J9Z38 2.389e-01
-cyan_free_sink  1.468e-01
-cyan_free       1.000e+00
-JCZ38_JSE76     1.000e+00
-JCZ38_sink      9.763e-12
+                    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   25.91 164.4    49.49        6.822        61.72
-JCZ38  22.46  74.6       NA           NA           NA
-J9Z38 112.37 373.3       NA           NA           NA
-JSE76 179.22 595.4       NA           NA           NA
+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

@@ -3382,10 +3368,10 @@ Hierarchical HS path 1 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:46:30 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:28:57 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ifelse(time <= tb, k1, k2) * cyan
@@ -3400,27 +3386,27 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 525.846 s
+Fitted in 485.304 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.8738        -3.4490        -4.9348        -5.5989         0.6469 
+      102.8845        -3.4495        -4.9355        -5.6040         0.6468 
   f_cyan_ilr_2 f_JCZ38_qlogis         log_k1         log_k2         log_tb 
-        1.2854         9.7193        -2.9084        -4.1810         1.7813 
+        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.409        0.00        0.00       0.000       0.0000
+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.006       0.0000
-f_cyan_ilr_1    0.000        0.00        0.00       0.000       0.6371
+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
@@ -3432,11 +3418,11 @@ 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.167           0.00 0.0000 0.0000 0.0000
-f_JCZ38_qlogis        0.000          10.22 0.0000 0.0000 0.0000
-log_k1                0.000           0.00 0.7003 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.6774
+log_tb                0.000           0.00 0.0000 0.0000 0.6773
 
 Starting values for error model parameters:
 a.1 
@@ -3446,220 +3432,76 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2427 2420  -1194
-
-Optimised parameters:
-                       est. lower upper
-cyan_0            101.84849    NA    NA
-log_k_JCZ38        -3.47365    NA    NA
-log_k_J9Z38        -5.10562    NA    NA
-log_k_JSE76        -5.60318    NA    NA
-f_cyan_ilr_1        0.66127    NA    NA
-f_cyan_ilr_2        0.60283    NA    NA
-f_JCZ38_qlogis     45.06408    NA    NA
-log_k1             -3.10124    NA    NA
-log_k2             -4.39028    NA    NA
-log_tb              2.32256    NA    NA
-a.1                 3.32683    NA    NA
-SD.log_k_JCZ38      1.41427    NA    NA
-SD.log_k_J9Z38      0.54767    NA    NA
-SD.log_k_JSE76      0.62147    NA    NA
-SD.f_cyan_ilr_1     0.30189    NA    NA
-SD.f_cyan_ilr_2     0.34960    NA    NA
-SD.f_JCZ38_qlogis   0.04644    NA    NA
-SD.log_k1           0.39534    NA    NA
-SD.log_k2           0.43468    NA    NA
-SD.log_tb           0.60781    NA    NA
-
-Correlation is not available
-
-Random effects:
-                     est. lower upper
-SD.log_k_JCZ38    1.41427    NA    NA
-SD.log_k_J9Z38    0.54767    NA    NA
-SD.log_k_JSE76    0.62147    NA    NA
-SD.f_cyan_ilr_1   0.30189    NA    NA
-SD.f_cyan_ilr_2   0.34960    NA    NA
-SD.f_JCZ38_qlogis 0.04644    NA    NA
-SD.log_k1         0.39534    NA    NA
-SD.log_k2         0.43468    NA    NA
-SD.log_tb         0.60781    NA    NA
-
-Variance model:
-     est. lower upper
-a.1 3.327    NA    NA
-
-Backtransformed parameters:
-                      est. lower upper
-cyan_0           1.018e+02    NA    NA
-k_JCZ38          3.100e-02    NA    NA
-k_J9Z38          6.063e-03    NA    NA
-k_JSE76          3.686e-03    NA    NA
-f_cyan_to_JCZ38  5.910e-01    NA    NA
-f_cyan_to_J9Z38  2.320e-01    NA    NA
-f_JCZ38_to_JSE76 1.000e+00    NA    NA
-k1               4.499e-02    NA    NA
-k2               1.240e-02    NA    NA
-tb               1.020e+01    NA    NA
-
-Resulting formation fractions:
-               ff
-cyan_JCZ38  0.591
-cyan_J9Z38  0.232
-cyan_sink   0.177
-JCZ38_JSE76 1.000
-JCZ38_sink  0.000
-
-Estimated disappearance times:
-        DT50   DT90 DT50back DT50_k1 DT50_k2
-cyan   29.09 158.91    47.84   15.41   55.91
-JCZ38  22.36  74.27       NA      NA      NA
-J9Z38 114.33 379.80       NA      NA      NA
-JSE76 188.04 624.66       NA      NA      NA
-
-
-

- -Hierarchical HS path 1 fit with two-component error - -

-saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:46:19 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
-
-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.968 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.168         -3.358         -4.941         -5.794          0.676 
-  f_cyan_ilr_2 f_JCZ38_qlogis         log_k1         log_k2         log_tb 
-         5.740         13.863         -3.147         -4.262          2.173 
-
-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.79       0.000       0.000       0.000       0.0000
-log_k_JCZ38      0.00       2.271       0.000       0.000       0.0000
-log_k_J9Z38      0.00       0.000       1.614       0.000       0.0000
-log_k_JSE76      0.00       0.000       0.000       1.264       0.0000
-f_cyan_ilr_1     0.00       0.000       0.000       0.000       0.6761
-f_cyan_ilr_2     0.00       0.000       0.000       0.000       0.0000
-f_JCZ38_qlogis   0.00       0.000       0.000       0.000       0.0000
-log_k1           0.00       0.000       0.000       0.000       0.0000
-log_k2           0.00       0.000       0.000       0.000       0.0000
-log_tb           0.00       0.000       0.000       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.000
-log_k_JCZ38           0.000           0.00 0.0000 0.0000  0.000
-log_k_J9Z38           0.000           0.00 0.0000 0.0000  0.000
-log_k_JSE76           0.000           0.00 0.0000 0.0000  0.000
-f_cyan_ilr_1          0.000           0.00 0.0000 0.0000  0.000
-f_cyan_ilr_2          9.572           0.00 0.0000 0.0000  0.000
-f_JCZ38_qlogis        0.000          19.19 0.0000 0.0000  0.000
-log_k1                0.000           0.00 0.8705 0.0000  0.000
-log_k2                0.000           0.00 0.0000 0.9288  0.000
-log_tb                0.000           0.00 0.0000 0.0000  1.065
-
-Starting values for error model parameters:
-a.1 b.1 
-  1   1 
-
-Results:
-
-Likelihood computed by importance sampling
-   AIC  BIC logLik
-  2422 2414  -1190
+  2427 2419  -1194
 
 Optimised parameters:
-                      est. lower upper
-cyan_0            100.9521    NA    NA
-log_k_JCZ38        -3.4629    NA    NA
-log_k_J9Z38        -5.0346    NA    NA
-log_k_JSE76        -5.5722    NA    NA
-f_cyan_ilr_1        0.6560    NA    NA
-f_cyan_ilr_2        0.7983    NA    NA
-f_JCZ38_qlogis     42.7949    NA    NA
-log_k1             -3.1721    NA    NA
-log_k2             -4.4039    NA    NA
-log_tb              2.3994    NA    NA
-a.1                 3.0586    NA    NA
-b.1                 0.0380    NA    NA
-SD.log_k_JCZ38      1.3754    NA    NA
-SD.log_k_J9Z38      0.6703    NA    NA
-SD.log_k_JSE76      0.5876    NA    NA
-SD.f_cyan_ilr_1     0.3272    NA    NA
-SD.f_cyan_ilr_2     0.5300    NA    NA
-SD.f_JCZ38_qlogis   6.4465    NA    NA
-SD.log_k1           0.4135    NA    NA
-SD.log_k2           0.4182    NA    NA
-SD.log_tb           0.6035    NA    NA
+                      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.3754    NA    NA
-SD.log_k_J9Z38    0.6703    NA    NA
-SD.log_k_JSE76    0.5876    NA    NA
-SD.f_cyan_ilr_1   0.3272    NA    NA
-SD.f_cyan_ilr_2   0.5300    NA    NA
-SD.f_JCZ38_qlogis 6.4465    NA    NA
-SD.log_k1         0.4135    NA    NA
-SD.log_k2         0.4182    NA    NA
-SD.log_tb         0.6035    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
 
 Variance model:
      est. lower upper
-a.1 3.059    NA    NA
-b.1 0.038    NA    NA
+a.1 3.323    NA    NA
 
 Backtransformed parameters:
-                      est. lower upper
-cyan_0           1.010e+02    NA    NA
-k_JCZ38          3.134e-02    NA    NA
-k_J9Z38          6.509e-03    NA    NA
-k_JSE76          3.802e-03    NA    NA
-f_cyan_to_JCZ38  6.127e-01    NA    NA
-f_cyan_to_J9Z38  2.423e-01    NA    NA
-f_JCZ38_to_JSE76 1.000e+00    NA    NA
-k1               4.191e-02    NA    NA
-k2               1.223e-02    NA    NA
-tb               1.102e+01    NA    NA
+                      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  0.6127
-cyan_J9Z38  0.2423
-cyan_sink   0.1449
-JCZ38_JSE76 1.0000
-JCZ38_sink  0.0000
+                   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.94 161.54    48.63   16.54   56.68
-JCZ38  22.12  73.47       NA      NA      NA
-J9Z38 106.50 353.77       NA      NA      NA
-JSE76 182.30 605.60       NA      NA      NA
+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

@@ -3668,14 +3510,14 @@ JSE76 182.30 605.60 NA NA NA

Pathway 2

-Hierarchical FOMC path 2 fit with constant variance +Hierarchical FOMC path 2 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:58:00 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:39:30 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan
@@ -3690,88 +3532,89 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 522.351 s
+Fitted in 474.942 s
 Using 300, 100 iterations and 10 chains
 
-Variance model: Constant variance 
+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.8173        -1.8998        -5.1449        -2.5415         0.6705 
+      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.4669        16.1281        13.3327        -0.2314         2.8738 
+        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          5.742       0.000       0.000        0.00       0.0000
-log_k_JCZ38     0.000       1.402       0.000        0.00       0.0000
-log_k_J9Z38     0.000       0.000       1.718        0.00       0.0000
-log_k_JSE76     0.000       0.000       0.000        3.57       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000        0.00       0.5926
-f_cyan_ilr_2    0.000       0.000       0.000        0.00       0.0000
-f_JCZ38_qlogis  0.000       0.000       0.000        0.00       0.0000
-f_JSE76_qlogis  0.000       0.000       0.000        0.00       0.0000
-log_alpha       0.000       0.000       0.000        0.00       0.0000
-log_beta        0.000       0.000       0.000        0.00       0.0000
+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.56           0.00           0.00    0.0000   0.0000
-f_JCZ38_qlogis         0.00          12.04           0.00    0.0000   0.0000
-f_JSE76_qlogis         0.00           0.00          15.26    0.0000   0.0000
-log_alpha              0.00           0.00           0.00    0.4708   0.0000
-log_beta               0.00           0.00           0.00    0.0000   0.4432
+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 
-  1 
+a.1 b.1 
+  1   1 
 
 Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2308 2301  -1134
+  2249 2241  -1104
 
 Optimised parameters:
-                      est.    lower     upper
-cyan_0            101.9586 99.22024 104.69700
-log_k_JCZ38        -2.4861 -3.17661  -1.79560
-log_k_J9Z38        -5.3926 -6.08842  -4.69684
-log_k_JSE76        -3.1193 -4.12904  -2.10962
-f_cyan_ilr_1        0.7368  0.42085   1.05276
-f_cyan_ilr_2        0.6196  0.06052   1.17861
-f_JCZ38_qlogis      4.8970 -4.68003  14.47398
-f_JSE76_qlogis      4.4066 -1.02087   9.83398
-log_alpha          -0.3021 -0.68264   0.07838
-log_beta            2.7438  2.57970   2.90786
-a.1                 2.9008  2.69920   3.10245
-SD.cyan_0           2.7081  0.64216   4.77401
-SD.log_k_JCZ38      0.7043  0.19951   1.20907
-SD.log_k_J9Z38      0.6248  0.05790   1.19180
-SD.log_k_JSE76      1.0750  0.33157   1.81839
-SD.f_cyan_ilr_1     0.3429  0.11688   0.56892
-SD.f_cyan_ilr_2     0.4774  0.09381   0.86097
-SD.f_JCZ38_qlogis   1.5565 -7.83970  10.95279
-SD.f_JSE76_qlogis   1.6871 -1.25577   4.63000
-SD.log_alpha        0.4216  0.15913   0.68405
+                       est.      lower      upper
+cyan_0            101.55265  9.920e+01  103.90593
+log_k_JCZ38        -2.32302 -2.832e+00   -1.81416
+log_k_J9Z38        -5.13082 -5.942e+00   -4.31990
+log_k_JSE76        -3.01756 -4.262e+00   -1.77360
+f_cyan_ilr_1        0.70850  3.657e-01    1.05135
+f_cyan_ilr_2        0.95775  2.612e-01    1.65432
+f_JCZ38_qlogis      3.86105  9.248e-01    6.79733
+f_JSE76_qlogis      7.51583 -1.120e+02  127.03921
+log_alpha          -0.15308 -4.508e-01    0.14462
+log_beta            2.99165  2.711e+00    3.27202
+a.1                 2.04034  1.811e+00    2.26968
+b.1                 0.06924  5.745e-02    0.08104
+SD.log_k_JCZ38      0.50818  1.390e-01    0.87736
+SD.log_k_J9Z38      0.86597  2.652e-01    1.46671
+SD.log_k_JSE76      1.38092  4.864e-01    2.27541
+SD.f_cyan_ilr_1     0.38204  1.354e-01    0.62864
+SD.f_cyan_ilr_2     0.55129  7.198e-02    1.03060
+SD.f_JCZ38_qlogis   1.88457  1.710e-02    3.75205
+SD.f_JSE76_qlogis   2.64018 -2.450e+03 2455.27887
+SD.log_alpha        0.31860  1.047e-01    0.53249
+SD.log_beta         0.24195  1.273e-02    0.47117
 
 Correlation: 
                cyan_0  l__JCZ3 l__J9Z3 l__JSE7 f_cy__1 f_cy__2 f_JCZ38 f_JSE76
-log_k_JCZ38    -0.0167                                                        
-log_k_J9Z38    -0.0307  0.0057                                                
-log_k_JSE76    -0.0032  0.1358  0.0009                                        
-f_cyan_ilr_1   -0.0087  0.0206 -0.1158 -0.0009                                
-f_cyan_ilr_2   -0.1598  0.0690  0.1770  0.0002 -0.0007                        
-f_JCZ38_qlogis  0.0966 -0.1132 -0.0440  0.0182 -0.1385 -0.4583                
-f_JSE76_qlogis -0.0647  0.1157  0.0333 -0.0026  0.1110  0.3620 -0.8586        
-log_alpha      -0.0389  0.0113  0.0209  0.0021  0.0041  0.0451 -0.0605  0.0412
-log_beta       -0.2508  0.0533  0.0977  0.0098  0.0220  0.2741 -0.2934  0.1999
+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           
@@ -3781,199 +3624,54 @@ f_cyan_ilr_2
 f_JCZ38_qlogis        
 f_JSE76_qlogis        
 log_alpha             
-log_beta        0.2281
+log_beta        0.2135
 
 Random effects:
-                    est.    lower   upper
-SD.cyan_0         2.7081  0.64216  4.7740
-SD.log_k_JCZ38    0.7043  0.19951  1.2091
-SD.log_k_J9Z38    0.6248  0.05790  1.1918
-SD.log_k_JSE76    1.0750  0.33157  1.8184
-SD.f_cyan_ilr_1   0.3429  0.11688  0.5689
-SD.f_cyan_ilr_2   0.4774  0.09381  0.8610
-SD.f_JCZ38_qlogis 1.5565 -7.83970 10.9528
-SD.f_JSE76_qlogis 1.6871 -1.25577  4.6300
-SD.log_alpha      0.4216  0.15913  0.6840
-
-Variance model:
-     est. lower upper
-a.1 2.901 2.699 3.102
-
-Backtransformed parameters:
-                      est.     lower     upper
-cyan_0           101.95862 99.220240 1.047e+02
-k_JCZ38            0.08323  0.041727 1.660e-01
-k_J9Z38            0.00455  0.002269 9.124e-03
-k_JSE76            0.04419  0.016098 1.213e-01
-f_cyan_to_JCZ38    0.61318        NA        NA
-f_cyan_to_J9Z38    0.21630        NA        NA
-f_JCZ38_to_JSE76   0.99259  0.009193 1.000e+00
-f_JSE76_to_JCZ38   0.98795  0.264857 9.999e-01
-alpha              0.73924  0.505281 1.082e+00
-beta              15.54568 13.193194 1.832e+01
-
-Resulting formation fractions:
-                  ff
-cyan_JCZ38  0.613182
-cyan_J9Z38  0.216298
-cyan_sink   0.170519
-JCZ38_JSE76 0.992586
-JCZ38_sink  0.007414
-JSE76_JCZ38 0.987950
-JSE76_sink  0.012050
-
-Estimated disappearance times:
-         DT50   DT90 DT50back
-cyan   24.157 334.68    100.7
-JCZ38   8.328  27.66       NA
-J9Z38 152.341 506.06       NA
-JSE76  15.687  52.11       NA
-
-
-

- -Hierarchical FOMC path 2 fit with two-component error - -

-saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:57:52 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
-
-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 514.301 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.9028        -1.9055        -5.0249        -2.5646         0.6807 
-  f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis      log_alpha       log_beta 
-        4.8883        16.0676         9.3923        -0.1346         3.0364 
-
-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          6.321       0.000       0.000       0.000       0.0000
-log_k_JCZ38     0.000       1.392       0.000       0.000       0.0000
-log_k_J9Z38     0.000       0.000       1.561       0.000       0.0000
-log_k_JSE76     0.000       0.000       0.000       3.614       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000       0.000       0.6339
-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.41           0.00           0.00    0.0000   0.0000
-f_JCZ38_qlogis         0.00          12.24           0.00    0.0000   0.0000
-f_JSE76_qlogis         0.00           0.00          15.13    0.0000   0.0000
-log_alpha              0.00           0.00           0.00    0.3701   0.0000
-log_beta               0.00           0.00           0.00    0.0000   0.5662
-
-Starting values for error model parameters:
-a.1 b.1 
-  1   1 
-
-Results:
-
-Likelihood computed by importance sampling
-   AIC  BIC logLik
-  2248 2240  -1103
-
-Optimised parameters:
-                       est.      lower      upper
-cyan_0            101.55545  9.920e+01  1.039e+02
-log_k_JCZ38        -2.37354 -2.928e+00 -1.819e+00
-log_k_J9Z38        -5.14736 -5.960e+00 -4.335e+00
-log_k_JSE76        -3.07802 -4.243e+00 -1.913e+00
-f_cyan_ilr_1        0.71263  3.655e-01  1.060e+00
-f_cyan_ilr_2        0.95202  2.701e-01  1.634e+00
-f_JCZ38_qlogis      3.58473  1.251e+00  5.919e+00
-f_JSE76_qlogis     19.03623 -1.037e+07  1.037e+07
-log_alpha          -0.15297 -4.490e-01  1.431e-01
-log_beta            2.99230  2.706e+00  3.278e+00
-a.1                 2.04816         NA         NA
-b.1                 0.06886         NA         NA
-SD.log_k_JCZ38      0.56174         NA         NA
-SD.log_k_J9Z38      0.86509         NA         NA
-SD.log_k_JSE76      1.28450         NA         NA
-SD.f_cyan_ilr_1     0.38705         NA         NA
-SD.f_cyan_ilr_2     0.54153         NA         NA
-SD.f_JCZ38_qlogis   1.65311         NA         NA
-SD.f_JSE76_qlogis   7.51468         NA         NA
-SD.log_alpha        0.31586         NA         NA
-SD.log_beta         0.24696         NA         NA
-
-Correlation is not available
-
-Random effects:
-                    est. lower upper
-SD.log_k_JCZ38    0.5617    NA    NA
-SD.log_k_J9Z38    0.8651    NA    NA
-SD.log_k_JSE76    1.2845    NA    NA
-SD.f_cyan_ilr_1   0.3870    NA    NA
-SD.f_cyan_ilr_2   0.5415    NA    NA
-SD.f_JCZ38_qlogis 1.6531    NA    NA
-SD.f_JSE76_qlogis 7.5147    NA    NA
-SD.log_alpha      0.3159    NA    NA
-SD.log_beta       0.2470    NA    NA
+                    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.710e-02    3.7520
+SD.f_JSE76_qlogis 2.6402 -2.450e+03 2455.2789
+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.04816    NA    NA
-b.1 0.06886    NA    NA
+       est.   lower   upper
+a.1 2.04034 1.81101 2.26968
+b.1 0.06924 0.05745 0.08104
 
 Backtransformed parameters:
-                      est.    lower    upper
-cyan_0           1.016e+02 99.20301 103.9079
-k_JCZ38          9.315e-02  0.05349   0.1622
-k_J9Z38          5.815e-03  0.00258   0.0131
-k_JSE76          4.605e-02  0.01436   0.1477
-f_cyan_to_JCZ38  6.438e-01       NA       NA
-f_cyan_to_J9Z38  2.350e-01       NA       NA
-f_JCZ38_to_JSE76 9.730e-01  0.77745   0.9973
-f_JSE76_to_JCZ38 1.000e+00  0.00000   1.0000
-alpha            8.582e-01  0.63824   1.1538
-beta             1.993e+01 14.97621  26.5262
+                      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  6.438e-01
-cyan_J9Z38  2.350e-01
-cyan_sink   1.212e-01
-JCZ38_JSE76 9.730e-01
-JCZ38_sink  2.700e-02
-JSE76_JCZ38 1.000e+00
-JSE76_sink  5.403e-09
+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.771 271.70    81.79
-JCZ38   7.441  24.72       NA
-J9Z38 119.205 395.99       NA
-JSE76  15.052  50.00       NA
+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

@@ -3982,10 +3680,10 @@ Hierarchical DFOP path 2 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:58:43 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:40:29 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -4005,25 +3703,25 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 565.562 s
+Fitted in 533.901 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.4358        -2.3107        -5.3123        -3.7120         0.6753 
+      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.1462        12.4095        12.3630        -1.9317        -4.4557 
+        1.1439        13.1194        12.3492        -1.9317        -4.4557 
       g_qlogis 
-       -0.5648 
+       -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.594      0.0000       0.000         0.0       0.0000
+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
@@ -4035,17 +3733,17 @@ 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.0  0.000 0.0000
-log_k_JCZ38           0.000           0.00            0.0  0.000 0.0000
-log_k_J9Z38           0.000           0.00            0.0  0.000 0.0000
-log_k_JSE76           0.000           0.00            0.0  0.000 0.0000
-f_cyan_ilr_1          0.000           0.00            0.0  0.000 0.0000
-f_cyan_ilr_2          1.797           0.00            0.0  0.000 0.0000
-f_JCZ38_qlogis        0.000          13.85            0.0  0.000 0.0000
-f_JSE76_qlogis        0.000           0.00           14.1  0.000 0.0000
-log_k1                0.000           0.00            0.0  1.106 0.0000
-log_k2                0.000           0.00            0.0  0.000 0.6141
-g_qlogis              0.000           0.00            0.0  0.000 0.0000
+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
@@ -4067,45 +3765,45 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2290 2281  -1123
+  2288 2280  -1122
 
 Optimised parameters:
-                      est.     lower    upper
-cyan_0            102.6903 101.44420 103.9365
-log_k_JCZ38        -2.4018  -2.98058  -1.8230
-log_k_J9Z38        -5.1865  -5.92931  -4.4437
-log_k_JSE76        -3.0784  -4.25226  -1.9045
-f_cyan_ilr_1        0.7157   0.37625   1.0551
-f_cyan_ilr_2        0.7073   0.20136   1.2132
-f_JCZ38_qlogis      4.6797   0.43240   8.9269
-f_JSE76_qlogis      5.0080  -1.01380  11.0299
-log_k1             -1.9620  -2.62909  -1.2949
-log_k2             -4.4894  -4.94958  -4.0292
-g_qlogis           -0.4658  -1.34443   0.4129
-a.1                 2.7158   2.52576   2.9059
-SD.log_k_JCZ38      0.5818   0.15679   1.0067
-SD.log_k_J9Z38      0.7421   0.16751   1.3167
-SD.log_k_JSE76      1.2841   0.43247   2.1356
-SD.f_cyan_ilr_1     0.3748   0.13040   0.6192
-SD.f_cyan_ilr_2     0.4550   0.08396   0.8261
-SD.f_JCZ38_qlogis   2.0862  -0.73390   4.9062
-SD.f_JSE76_qlogis   1.9585  -3.14773   7.0647
-SD.log_k1           0.7389   0.25761   1.2201
-SD.log_k2           0.5132   0.18143   0.8450
-SD.g_qlogis         0.9870   0.35773   1.6164
+                      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.0170                                                        
-log_k_J9Z38    -0.0457  0.0016                                                
-log_k_JSE76    -0.0046  0.1183  0.0005                                        
-f_cyan_ilr_1    0.0079  0.0072 -0.0909  0.0003                                
-f_cyan_ilr_2   -0.3114  0.0343  0.1542  0.0023 -0.0519                        
-f_JCZ38_qlogis  0.0777 -0.0601 -0.0152  0.0080 -0.0520 -0.2524                
-f_JSE76_qlogis -0.0356  0.0817  0.0073  0.0051  0.0388  0.1959 -0.6236        
-log_k1          0.0848 -0.0028  0.0010 -0.0010 -0.0014 -0.0245  0.0121 -0.0177
-log_k2          0.0274 -0.0001  0.0075  0.0000 -0.0023 -0.0060  0.0000 -0.0130
-g_qlogis        0.0159  0.0002 -0.0095  0.0002  0.0029 -0.0140 -0.0001  0.0149
+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                   
@@ -4115,56 +3813,56 @@ f_cyan_ilr_2
 f_JCZ38_qlogis                
 f_JSE76_qlogis                
 log_k1                        
-log_k2          0.0280        
-g_qlogis       -0.0278 -0.0310
+log_k2          0.0276        
+g_qlogis       -0.0283 -0.0309
 
 Random effects:
-                    est.    lower  upper
-SD.log_k_JCZ38    0.5818  0.15679 1.0067
-SD.log_k_J9Z38    0.7421  0.16751 1.3167
-SD.log_k_JSE76    1.2841  0.43247 2.1356
-SD.f_cyan_ilr_1   0.3748  0.13040 0.6192
-SD.f_cyan_ilr_2   0.4550  0.08396 0.8261
-SD.f_JCZ38_qlogis 2.0862 -0.73390 4.9062
-SD.f_JSE76_qlogis 1.9585 -3.14773 7.0647
-SD.log_k1         0.7389  0.25761 1.2201
-SD.log_k2         0.5132  0.18143 0.8450
-SD.g_qlogis       0.9870  0.35773 1.6164
+                    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.716 2.526 2.906
+a.1 2.737 2.545 2.928
 
 Backtransformed parameters:
                       est.     lower     upper
-cyan_0           1.027e+02 1.014e+02 103.93649
-k_JCZ38          9.056e-02 5.076e-02   0.16154
-k_J9Z38          5.591e-03 2.660e-03   0.01175
-k_JSE76          4.603e-02 1.423e-02   0.14890
-f_cyan_to_JCZ38  6.184e-01        NA        NA
-f_cyan_to_J9Z38  2.248e-01        NA        NA
-f_JCZ38_to_JSE76 9.908e-01 6.064e-01   0.99987
-f_JSE76_to_JCZ38 9.934e-01 2.662e-01   0.99998
-k1               1.406e-01 7.214e-02   0.27393
-k2               1.123e-02 7.086e-03   0.01779
-g                3.856e-01 2.068e-01   0.60177
+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  0.618443
-cyan_J9Z38  0.224770
-cyan_sink   0.156787
-JCZ38_JSE76 0.990803
-JCZ38_sink  0.009197
-JSE76_JCZ38 0.993360
-JSE76_sink  0.006640
+                   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.674 161.70    48.68   4.931   61.74
-JCZ38   7.654  25.43       NA      NA      NA
-J9Z38 123.966 411.81       NA      NA      NA
-JSE76  15.057  50.02       NA      NA      NA
+        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

@@ -4173,10 +3871,10 @@ Hierarchical DFOP path 2 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:01:24 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:43:04 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -4196,59 +3894,59 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 726.501 s
+Fitted in 688.913 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.7523        -1.5948        -5.0119        -2.2723         0.6719 
+      101.7393        -1.4493        -5.0118        -2.1269         0.6720 
   f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis         log_k1         log_k2 
-        5.1681        12.8238        12.4130        -2.0057        -4.5526 
+        7.3362        13.4423        13.2659        -2.0061        -4.5527 
       g_qlogis 
-       -0.5805 
+       -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.627       0.000       0.000       0.000       0.0000
-log_k_JCZ38     0.000       2.327       0.000       0.000       0.0000
-log_k_J9Z38     0.000       0.000       1.664       0.000       0.0000
-log_k_JSE76     0.000       0.000       0.000       4.566       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000       0.000       0.6519
-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_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
+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.0           0.00           0.00 0.0000 0.0000
-log_k_JCZ38             0.0           0.00           0.00 0.0000 0.0000
-log_k_J9Z38             0.0           0.00           0.00 0.0000 0.0000
-log_k_JSE76             0.0           0.00           0.00 0.0000 0.0000
-f_cyan_ilr_1            0.0           0.00           0.00 0.0000 0.0000
-f_cyan_ilr_2           10.1           0.00           0.00 0.0000 0.0000
-f_JCZ38_qlogis          0.0          13.99           0.00 0.0000 0.0000
-f_JSE76_qlogis          0.0           0.00          14.15 0.0000 0.0000
-log_k1                  0.0           0.00           0.00 0.8452 0.0000
-log_k2                  0.0           0.00           0.00 0.0000 0.5968
-g_qlogis                0.0           0.00           0.00 0.0000 0.0000
+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.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.691
+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 
@@ -4261,42 +3959,42 @@ Likelihood computed by importance sampling
   2234 2226  -1095
 
 Optimised parameters:
-                       est.      lower      upper
-cyan_0            101.10667  9.903e+01  103.18265
-log_k_JCZ38        -2.49437 -3.297e+00   -1.69221
-log_k_J9Z38        -5.08171 -5.875e+00   -4.28846
-log_k_JSE76        -3.20072 -4.180e+00   -2.22163
-f_cyan_ilr_1        0.71059  3.639e-01    1.05727
-f_cyan_ilr_2        1.15398  2.981e-01    2.00984
-f_JCZ38_qlogis      3.18027  1.056e+00    5.30452
-f_JSE76_qlogis      5.61578 -2.505e+01   36.28077
-log_k1             -2.38875 -2.517e+00   -2.26045
-log_k2             -4.67246 -4.928e+00   -4.41715
-g_qlogis           -0.28231 -1.135e+00    0.57058
-a.1                 2.08190  1.856e+00    2.30785
-b.1                 0.06114  5.015e-02    0.07214
-SD.log_k_JCZ38      0.84622  2.637e-01    1.42873
-SD.log_k_J9Z38      0.84564  2.566e-01    1.43464
-SD.log_k_JSE76      1.04385  3.242e-01    1.76351
-SD.f_cyan_ilr_1     0.38568  1.362e-01    0.63514
-SD.f_cyan_ilr_2     0.68046  7.166e-02    1.28925
-SD.f_JCZ38_qlogis   1.25244 -4.213e-02    2.54700
-SD.f_JSE76_qlogis   0.28202 -1.515e+03 1515.87968
-SD.log_k2           0.25749  7.655e-02    0.43843
-SD.g_qlogis         0.94535  3.490e-01    1.54174
+                       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.83267   2.28582
+b.1                 0.06172   0.05076   0.07268
+SD.log_k_JCZ38      0.81137   0.25296   1.36977
+SD.log_k_J9Z38      0.83542   0.25396   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.42361   3.13015
+SD.f_JSE76_qlogis   1.43956 -19.15140  22.03052
+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.0086                                                        
-log_k_J9Z38    -0.0363 -0.0007                                                
-log_k_JSE76     0.0015  0.1210 -0.0017                                        
-f_cyan_ilr_1   -0.0048  0.0095 -0.0572  0.0030                                
-f_cyan_ilr_2   -0.4788  0.0328  0.1143  0.0027 -0.0316                        
-f_JCZ38_qlogis  0.0736 -0.0664 -0.0137  0.0145 -0.0444 -0.2175                
-f_JSE76_qlogis -0.0137  0.0971  0.0035  0.0009  0.0293  0.1333 -0.6767        
-log_k1          0.2345 -0.0350 -0.0099 -0.0113 -0.0126 -0.1652  0.1756 -0.2161
-log_k2          0.0440 -0.0133  0.0199 -0.0040 -0.0097 -0.0119  0.0604 -0.1306
-g_qlogis        0.0438  0.0078 -0.0123  0.0029  0.0046 -0.0363 -0.0318  0.0736
+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                   
@@ -4306,56 +4004,56 @@ f_cyan_ilr_2
 f_JCZ38_qlogis                
 f_JSE76_qlogis                
 log_k1                        
-log_k2          0.3198        
-g_qlogis       -0.1666 -0.0954
+log_k2          0.3269        
+g_qlogis       -0.1656 -0.0928
 
 Random effects:
-                    est.      lower     upper
-SD.log_k_JCZ38    0.8462  2.637e-01    1.4287
-SD.log_k_J9Z38    0.8456  2.566e-01    1.4346
-SD.log_k_JSE76    1.0439  3.242e-01    1.7635
-SD.f_cyan_ilr_1   0.3857  1.362e-01    0.6351
-SD.f_cyan_ilr_2   0.6805  7.166e-02    1.2893
-SD.f_JCZ38_qlogis 1.2524 -4.213e-02    2.5470
-SD.f_JSE76_qlogis 0.2820 -1.515e+03 1515.8797
-SD.log_k2         0.2575  7.655e-02    0.4384
-SD.g_qlogis       0.9453  3.490e-01    1.5417
+                    est.     lower   upper
+SD.log_k_JCZ38    0.8114   0.25296  1.3698
+SD.log_k_J9Z38    0.8354   0.25396  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.42361  3.1301
+SD.f_JSE76_qlogis 1.4396 -19.15140 22.0305
+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.08190 1.85595 2.30785
-b.1 0.06114 0.05015 0.07214
+a.1 2.05925 1.83267 2.28582
+b.1 0.06172 0.05076 0.07268
 
 Backtransformed parameters:
                       est.     lower     upper
-cyan_0           1.011e+02 9.903e+01 103.18265
-k_JCZ38          8.255e-02 3.701e-02   0.18411
-k_J9Z38          6.209e-03 2.809e-03   0.01373
-k_JSE76          4.073e-02 1.530e-02   0.10843
-f_cyan_to_JCZ38  6.608e-01        NA        NA
-f_cyan_to_J9Z38  2.419e-01        NA        NA
-f_JCZ38_to_JSE76 9.601e-01 7.419e-01   0.99506
-f_JSE76_to_JCZ38 9.964e-01 1.322e-11   1.00000
-k1               9.174e-02 8.070e-02   0.10430
-k2               9.349e-03 7.243e-03   0.01207
-g                4.299e-01 2.432e-01   0.63890
+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.660808
-cyan_J9Z38  0.241904
-cyan_sink   0.097288
-JCZ38_JSE76 0.960085
-JCZ38_sink  0.039915
-JSE76_JCZ38 0.996373
-JSE76_sink  0.003627
+                 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.359 186.18    56.05   7.555   74.14
-JCZ38   8.397  27.89       NA      NA      NA
-J9Z38 111.631 370.83       NA      NA      NA
-JSE76  17.017  56.53       NA      NA      NA
+        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

@@ -4364,10 +4062,10 @@ Hierarchical SFORB path 2 fit with constant variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 07:58:46 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:40:32 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -4385,20 +4083,20 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 568.562 s
+Fitted in 536.94 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.4394               -2.7673               -2.8942 
+             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 
-              13.2672               13.3538 
+              14.8408               15.4734 
 
 Fixed degradation parameter values:
 None
@@ -4429,17 +4127,17 @@ 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.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.6349        0.000           0.00           0.00
-f_cyan_ilr_2                0.0000        1.797           0.00           0.00
-f_JCZ38_qlogis              0.0000        0.000          13.84           0.00
-f_JSE76_qlogis              0.0000        0.000           0.00          14.66
+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 
@@ -4449,45 +4147,45 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2284 2275  -1120
+  2283 2275  -1120
 
 Optimised parameters:
-                             est.      lower      upper
-cyan_free_0              102.7730  1.015e+02  1.041e+02
-log_k_cyan_free           -2.8530 -3.167e+00 -2.539e+00
-log_k_cyan_free_bound     -2.7326 -3.543e+00 -1.922e+00
-log_k_cyan_bound_free     -3.5582 -4.126e+00 -2.990e+00
-log_k_JCZ38               -2.3810 -2.921e+00 -1.841e+00
-log_k_J9Z38               -5.2301 -5.963e+00 -4.497e+00
-log_k_JSE76               -3.0286 -4.286e+00 -1.771e+00
-f_cyan_ilr_1               0.7081  3.733e-01  1.043e+00
-f_cyan_ilr_2               0.5847  7.846e-03  1.162e+00
-f_JCZ38_qlogis             9.5676 -1.323e+03  1.342e+03
-f_JSE76_qlogis             3.7042  7.254e-02  7.336e+00
-a.1                        2.7222  2.532e+00  2.913e+00
-SD.log_k_cyan_free         0.3338  1.086e-01  5.589e-01
-SD.log_k_cyan_free_bound   0.8888  3.023e-01  1.475e+00
-SD.log_k_cyan_bound_free   0.6220  2.063e-01  1.038e+00
-SD.log_k_JCZ38             0.5221  1.334e-01  9.108e-01
-SD.log_k_J9Z38             0.7104  1.371e-01  1.284e+00
-SD.log_k_JSE76             1.3837  4.753e-01  2.292e+00
-SD.f_cyan_ilr_1            0.3620  1.248e-01  5.992e-01
-SD.f_cyan_ilr_2            0.4259  8.145e-02  7.704e-01
-SD.f_JCZ38_qlogis          3.5332 -1.037e+05  1.037e+05
-SD.f_JSE76_qlogis          1.6990 -2.771e-01  3.675e+00
+                             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.2126                                                
-log_k_cyan_free_bound  0.0894  0.0871                                        
-log_k_cyan_bound_free  0.0033  0.0410  0.0583                                
-log_k_JCZ38           -0.0708 -0.0280 -0.0147      0.0019                    
-log_k_J9Z38           -0.0535 -0.0138  0.0012      0.0148      0.0085        
-log_k_JSE76           -0.0066 -0.0030 -0.0021     -0.0005      0.1090  0.0010
-f_cyan_ilr_1          -0.0364 -0.0157 -0.0095     -0.0015      0.0458 -0.0960
-f_cyan_ilr_2          -0.3814 -0.1104 -0.0423      0.0146      0.1540  0.1526
-f_JCZ38_qlogis         0.2507  0.0969  0.0482     -0.0097     -0.2282 -0.0363
-f_JSE76_qlogis        -0.1648 -0.0710 -0.0443     -0.0087      0.2002  0.0226
+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                                
@@ -4495,63 +4193,63 @@ log_k_cyan_bound_free
 log_k_JCZ38                                          
 log_k_J9Z38                                          
 log_k_JSE76                                          
-f_cyan_ilr_1           0.0001                        
-f_cyan_ilr_2           0.0031  0.0586                
-f_JCZ38_qlogis         0.0023 -0.1867 -0.6255        
-f_JSE76_qlogis         0.0082  0.1356  0.4519 -0.7951
+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.3338  1.086e-01 5.589e-01
-SD.log_k_cyan_free_bound 0.8888  3.023e-01 1.475e+00
-SD.log_k_cyan_bound_free 0.6220  2.063e-01 1.038e+00
-SD.log_k_JCZ38           0.5221  1.334e-01 9.108e-01
-SD.log_k_J9Z38           0.7104  1.371e-01 1.284e+00
-SD.log_k_JSE76           1.3837  4.753e-01 2.292e+00
-SD.f_cyan_ilr_1          0.3620  1.248e-01 5.992e-01
-SD.f_cyan_ilr_2          0.4259  8.145e-02 7.704e-01
-SD.f_JCZ38_qlogis        3.5332 -1.037e+05 1.037e+05
-SD.f_JSE76_qlogis        1.6990 -2.771e-01 3.675e+00
+                           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.722 2.532 2.913
+a.1 2.716 2.526 2.906
 
 Backtransformed parameters:
                           est.     lower     upper
-cyan_free_0          1.028e+02 1.015e+02 104.06475
-k_cyan_free          5.767e-02 4.213e-02   0.07894
-k_cyan_free_bound    6.505e-02 2.892e-02   0.14633
-k_cyan_bound_free    2.849e-02 1.614e-02   0.05028
-k_JCZ38              9.246e-02 5.390e-02   0.15859
-k_J9Z38              5.353e-03 2.572e-03   0.01114
-k_JSE76              4.838e-02 1.376e-02   0.17009
-f_cyan_free_to_JCZ38 6.011e-01 5.028e-01   0.83792
-f_cyan_free_to_J9Z38 2.208e-01 5.028e-01   0.83792
-f_JCZ38_to_JSE76     9.999e-01 0.000e+00   1.00000
-f_JSE76_to_JCZ38     9.760e-01 5.181e-01   0.99935
+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.13942 0.01178 0.35948 
+0.13466 0.01165 0.36490 
 
 Resulting formation fractions:
-                       ff
-cyan_free_JCZ38 6.011e-01
-cyan_free_J9Z38 2.208e-01
-cyan_free_sink  1.780e-01
-cyan_free       1.000e+00
-JCZ38_JSE76     9.999e-01
-JCZ38_sink      6.996e-05
-JSE76_JCZ38     9.760e-01
-JSE76_sink      2.403e-02
+                      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.390 157.60    47.44        4.971        58.82
-JCZ38   7.497  24.90       NA           NA           NA
-J9Z38 129.482 430.13       NA           NA           NA
-JSE76  14.326  47.59       NA           NA           NA
+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

@@ -4560,10 +4258,10 @@ Hierarchical SFORB path 2 fit with two-component error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:01:30 2023 
-Date of summary: Thu Apr 20 20:01:30 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:42:47 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -4581,20 +4279,20 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 732.212 s
+Fitted in 671.849 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.751                -2.837                -3.016 
+             101.7511               -2.8370               -3.0162 
 log_k_cyan_bound_free           log_k_JCZ38           log_k_J9Z38 
-               -3.660                -2.299                -5.313 
+              -3.6600               -2.2988               -5.3129 
           log_k_JSE76          f_cyan_ilr_1          f_cyan_ilr_2 
-               -3.699                 0.672                 5.873 
+              -3.6991                0.6722                4.8596 
        f_JCZ38_qlogis        f_JSE76_qlogis 
-               13.216                13.338 
+              13.4678               14.2149 
 
 Fixed degradation parameter values:
 None
@@ -4625,17 +4323,17 @@ 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.00           0.00           0.00
-log_k_cyan_free             0.0000         0.00           0.00           0.00
-log_k_cyan_free_bound       0.0000         0.00           0.00           0.00
-log_k_cyan_bound_free       0.0000         0.00           0.00           0.00
-log_k_JCZ38                 0.0000         0.00           0.00           0.00
-log_k_J9Z38                 0.0000         0.00           0.00           0.00
-log_k_JSE76                 0.0000         0.00           0.00           0.00
-f_cyan_ilr_1                0.6519         0.00           0.00           0.00
-f_cyan_ilr_2                0.0000        10.78           0.00           0.00
-f_JCZ38_qlogis              0.0000         0.00          13.96           0.00
-f_JSE76_qlogis              0.0000         0.00           0.00          14.69
+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 
@@ -4645,45 +4343,45 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2240 2232  -1098
+  2240 2231  -1098
 
 Optimised parameters:
-                              est.     lower    upper
-cyan_free_0              101.10205  98.99221 103.2119
-log_k_cyan_free           -3.16929  -3.61395  -2.7246
-log_k_cyan_free_bound     -3.38259  -3.63022  -3.1350
-log_k_cyan_bound_free     -3.81075  -4.13888  -3.4826
-log_k_JCZ38               -2.42057  -3.00756  -1.8336
-log_k_J9Z38               -5.07501  -5.85138  -4.2986
-log_k_JSE76               -3.12442  -4.21277  -2.0361
-f_cyan_ilr_1               0.70577   0.35788   1.0537
-f_cyan_ilr_2               1.14824   0.15810   2.1384
-f_JCZ38_qlogis             3.52245   0.43257   6.6123
-f_JSE76_qlogis             5.65140 -21.22295  32.5257
-a.1                        2.07062   1.84329   2.2980
-b.1                        0.06227   0.05124   0.0733
-SD.log_k_cyan_free         0.49468   0.18566   0.8037
-SD.log_k_cyan_bound_free   0.28972   0.07188   0.5076
-SD.log_k_JCZ38             0.58852   0.16800   1.0090
-SD.log_k_J9Z38             0.82500   0.24730   1.4027
-SD.log_k_JSE76             1.19201   0.40313   1.9809
-SD.f_cyan_ilr_1            0.38534   0.13640   0.6343
-SD.f_cyan_ilr_2            0.72463   0.10076   1.3485
-SD.f_JCZ38_qlogis          1.38223  -0.20997   2.9744
-SD.f_JSE76_qlogis          2.07989 -72.53027  76.6901
+                              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.846e+00  2.296e+00
+b.1                        0.06227  5.120e-02  7.334e-02
+SD.log_k_cyan_free         0.49674  1.865e-01  8.069e-01
+SD.log_k_cyan_bound_free   0.28537  6.808e-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.809e-02  2.128e+00
+SD.f_JSE76_qlogis          0.44681 -3.986e+05  3.986e+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.1117                                                
-log_k_cyan_free_bound  0.1763  0.1828                                        
-log_k_cyan_bound_free  0.0120  0.0593  0.5030                                
-log_k_JCZ38           -0.0459 -0.0230 -0.0931     -0.0337                    
-log_k_J9Z38           -0.0381 -0.0123 -0.0139      0.0237      0.0063        
-log_k_JSE76           -0.0044 -0.0038 -0.0175     -0.0072      0.1120  0.0003
-f_cyan_ilr_1          -0.0199 -0.0087 -0.0407     -0.0233      0.0268 -0.0552
-f_cyan_ilr_2          -0.4806 -0.1015 -0.2291     -0.0269      0.1156  0.1113
-f_JCZ38_qlogis         0.1805  0.0825  0.3085      0.0963     -0.1674 -0.0314
-f_JSE76_qlogis        -0.1586 -0.0810 -0.3560     -0.1563      0.2025  0.0278
+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                                
@@ -4691,63 +4389,63 @@ log_k_cyan_bound_free
 log_k_JCZ38                                          
 log_k_J9Z38                                          
 log_k_JSE76                                          
-f_cyan_ilr_1           0.0024                        
-f_cyan_ilr_2           0.0087  0.0172                
-f_JCZ38_qlogis        -0.0016 -0.1047 -0.4656        
-f_JSE76_qlogis         0.0119  0.1034  0.4584 -0.8137
+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.4947   0.18566  0.8037
-SD.log_k_cyan_bound_free 0.2897   0.07188  0.5076
-SD.log_k_JCZ38           0.5885   0.16800  1.0090
-SD.log_k_J9Z38           0.8250   0.24730  1.4027
-SD.log_k_JSE76           1.1920   0.40313  1.9809
-SD.f_cyan_ilr_1          0.3853   0.13640  0.6343
-SD.f_cyan_ilr_2          0.7246   0.10076  1.3485
-SD.f_JCZ38_qlogis        1.3822  -0.20997  2.9744
-SD.f_JSE76_qlogis        2.0799 -72.53027 76.6901
+                           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.808e-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.809e-02 2.128e+00
+SD.f_JSE76_qlogis        0.4468 -3.986e+05 3.986e+05
 
 Variance model:
-       est.   lower  upper
-a.1 2.07062 1.84329 2.2980
-b.1 0.06227 0.05124 0.0733
+       est.  lower   upper
+a.1 2.07082 1.8458 2.29588
+b.1 0.06227 0.0512 0.07334
 
 Backtransformed parameters:
-                          est.     lower     upper
-cyan_free_0          1.011e+02 9.899e+01 103.21190
-k_cyan_free          4.203e-02 2.695e-02   0.06557
-k_cyan_free_bound    3.396e-02 2.651e-02   0.04350
-k_cyan_bound_free    2.213e-02 1.594e-02   0.03073
-k_JCZ38              8.887e-02 4.941e-02   0.15984
-k_J9Z38              6.251e-03 2.876e-03   0.01359
-k_JSE76              4.396e-02 1.481e-02   0.13054
-f_cyan_free_to_JCZ38 6.590e-01 5.557e-01   0.95365
-f_cyan_free_to_J9Z38 2.429e-01 5.557e-01   0.95365
-f_JCZ38_to_JSE76     9.713e-01 6.065e-01   0.99866
-f_JSE76_to_JCZ38     9.965e-01 6.067e-10   1.00000
+                          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.08749 0.01063 0.40855 
+0.08426 0.01051 0.41220 
 
 Resulting formation fractions:
-                     ff
-cyan_free_JCZ38 0.65905
-cyan_free_J9Z38 0.24291
-cyan_free_sink  0.09805
-cyan_free       1.00000
-JCZ38_JSE76     0.97132
-JCZ38_sink      0.02868
-JSE76_JCZ38     0.99650
-JSE76_sink      0.00350
+                       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   24.91 167.16    50.32        7.922        65.19
-JCZ38   7.80  25.91       NA           NA           NA
-J9Z38 110.89 368.36       NA           NA           NA
-JSE76  15.77  52.38       NA           NA           NA
+         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

@@ -4761,10 +4459,10 @@ error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:15:01 2023 
-Date of summary: Thu Apr 20 20:01:31 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:55:35 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - (alpha/beta) * 1/((time/beta) + 1) * cyan
@@ -4779,27 +4477,27 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 808.728 s
+Fitted in 748.54 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.9028        -1.9055        -5.0249        -2.5646         0.6807 
+      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.8883        16.0676         9.3923        -0.1346         3.0364 
+        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          6.321       0.000       0.000       0.000       0.0000
-log_k_JCZ38     0.000       1.392       0.000       0.000       0.0000
-log_k_J9Z38     0.000       0.000       1.561       0.000       0.0000
-log_k_JSE76     0.000       0.000       0.000       3.614       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000       0.000       0.6339
+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
@@ -4811,11 +4509,11 @@ 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.41           0.00           0.00    0.0000   0.0000
-f_JCZ38_qlogis         0.00          12.24           0.00    0.0000   0.0000
-f_JSE76_qlogis         0.00           0.00          15.13    0.0000   0.0000
-log_alpha              0.00           0.00           0.00    0.3701   0.0000
-log_beta               0.00           0.00           0.00    0.0000   0.5662
+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 
@@ -4825,76 +4523,76 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2251 2244  -1106
+  2249 2242  -1106
 
 Optimised parameters:
-                      est.   lower   upper
-cyan_0           101.05768      NA      NA
-log_k_JCZ38       -2.73252      NA      NA
-log_k_J9Z38       -5.07399      NA      NA
-log_k_JSE76       -3.52863      NA      NA
-f_cyan_ilr_1       0.72176      NA      NA
-f_cyan_ilr_2       1.34610      NA      NA
-f_JCZ38_qlogis     2.08337      NA      NA
-f_JSE76_qlogis  1590.31880      NA      NA
-log_alpha         -0.09336      NA      NA
-log_beta           3.10191      NA      NA
-a.1                2.08557 1.85439 2.31675
-b.1                0.06998 0.05800 0.08197
-SD.log_k_JCZ38     1.20053 0.43329 1.96777
-SD.log_k_J9Z38     0.85854 0.26708 1.45000
-SD.log_k_JSE76     0.62528 0.16061 1.08995
-SD.f_cyan_ilr_1    0.35190 0.12340 0.58039
-SD.f_cyan_ilr_2    0.85385 0.15391 1.55378
-SD.log_alpha       0.28971 0.08718 0.49225
-SD.log_beta        0.31614 0.05938 0.57290
+                     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.82393 2.28166
+b.1               0.07116 0.05910 0.08322
+SD.log_k_JCZ38    1.21713 0.44160 1.99266
+SD.log_k_J9Z38    0.88268 0.27541 1.48995
+SD.log_k_JSE76    0.59452 0.15005 1.03898
+SD.f_cyan_ilr_1   0.35370 0.12409 0.58331
+SD.f_cyan_ilr_2   0.78186 0.18547 1.37824
+SD.log_alpha      0.27781 0.08168 0.47394
+SD.log_beta       0.32608 0.06490 0.58726
 
 Correlation is not available
 
 Random effects:
                   est.   lower  upper
-SD.log_k_JCZ38  1.2005 0.43329 1.9678
-SD.log_k_J9Z38  0.8585 0.26708 1.4500
-SD.log_k_JSE76  0.6253 0.16061 1.0900
-SD.f_cyan_ilr_1 0.3519 0.12340 0.5804
-SD.f_cyan_ilr_2 0.8538 0.15391 1.5538
-SD.log_alpha    0.2897 0.08718 0.4923
-SD.log_beta     0.3161 0.05938 0.5729
+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.08557 1.854 2.31675
-b.1 0.06998 0.058 0.08197
+       est.  lower   upper
+a.1 2.05279 1.8239 2.28166
+b.1 0.07116 0.0591 0.08322
 
 Backtransformed parameters:
                       est. lower upper
-cyan_0           1.011e+02    NA    NA
-k_JCZ38          6.506e-02    NA    NA
-k_J9Z38          6.257e-03    NA    NA
-k_JSE76          2.935e-02    NA    NA
-f_cyan_to_JCZ38  6.776e-01    NA    NA
-f_cyan_to_J9Z38  2.442e-01    NA    NA
-f_JCZ38_to_JSE76 8.893e-01    NA    NA
-f_JSE76_to_JCZ38 1.000e+00    NA    NA
-alpha            9.109e-01    NA    NA
-beta             2.224e+01    NA    NA
+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.67761
-cyan_J9Z38  0.24417
-cyan_sink   0.07822
-JCZ38_JSE76 0.88928
-JCZ38_sink  0.11072
-JSE76_JCZ38 1.00000
-JSE76_sink  0.00000
+                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.36 256.37    77.18
-JCZ38  10.65  35.39       NA
-J9Z38 110.77 367.98       NA
-JSE76  23.62  78.47       NA
+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

@@ -4904,10 +4602,10 @@ variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:16:32 2023 
-Date of summary: Thu Apr 20 20:01:31 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:57:10 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -4927,25 +4625,25 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 900.061 s
+Fitted in 843.793 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.4358        -2.3107        -5.3123        -3.7120         0.6753 
+      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.1462        12.4095        12.3630        -1.9317        -4.4557 
+        1.1439        13.1194        12.3492        -1.9317        -4.4557 
       g_qlogis 
-       -0.5648 
+       -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.594      0.0000       0.000         0.0       0.0000
+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
@@ -4957,17 +4655,17 @@ 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.0  0.000 0.0000
-log_k_JCZ38           0.000           0.00            0.0  0.000 0.0000
-log_k_J9Z38           0.000           0.00            0.0  0.000 0.0000
-log_k_JSE76           0.000           0.00            0.0  0.000 0.0000
-f_cyan_ilr_1          0.000           0.00            0.0  0.000 0.0000
-f_cyan_ilr_2          1.797           0.00            0.0  0.000 0.0000
-f_JCZ38_qlogis        0.000          13.85            0.0  0.000 0.0000
-f_JSE76_qlogis        0.000           0.00           14.1  0.000 0.0000
-log_k1                0.000           0.00            0.0  1.106 0.0000
-log_k2                0.000           0.00            0.0  0.000 0.6141
-g_qlogis              0.000           0.00            0.0  0.000 0.0000
+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
@@ -4993,74 +4691,74 @@ Likelihood computed by importance sampling
 
 Optimised parameters:
                      est.   lower  upper
-cyan_0           102.5254      NA     NA
-log_k_JCZ38       -2.9358      NA     NA
-log_k_J9Z38       -5.1424      NA     NA
-log_k_JSE76       -3.6458      NA     NA
-f_cyan_ilr_1       0.6957      NA     NA
-f_cyan_ilr_2       0.6635      NA     NA
-f_JCZ38_qlogis  4984.8163      NA     NA
-f_JSE76_qlogis     1.9415      NA     NA
-log_k1            -1.9456      NA     NA
-log_k2            -4.4705      NA     NA
-g_qlogis          -0.5117      NA     NA
-a.1                2.7455 2.55392 2.9370
-SD.log_k_JCZ38     1.3163 0.47635 2.1563
-SD.log_k_J9Z38     0.7162 0.16133 1.2711
-SD.log_k_JSE76     0.6457 0.15249 1.1390
-SD.f_cyan_ilr_1    0.3424 0.11714 0.5677
-SD.f_cyan_ilr_2    0.4524 0.09709 0.8077
-SD.log_k1          0.7353 0.25445 1.2161
-SD.log_k2          0.5137 0.18206 0.8453
-SD.g_qlogis        0.9857 0.35651 1.6148
+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.3163 0.47635 2.1563
-SD.log_k_J9Z38  0.7162 0.16133 1.2711
-SD.log_k_JSE76  0.6457 0.15249 1.1390
-SD.f_cyan_ilr_1 0.3424 0.11714 0.5677
-SD.f_cyan_ilr_2 0.4524 0.09709 0.8077
-SD.log_k1       0.7353 0.25445 1.2161
-SD.log_k2       0.5137 0.18206 0.8453
-SD.g_qlogis     0.9857 0.35651 1.6148
+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.745 2.554 2.937
+a.1 2.746 2.593 2.899
 
 Backtransformed parameters:
                       est. lower upper
-cyan_0           1.025e+02    NA    NA
-k_JCZ38          5.309e-02    NA    NA
-k_J9Z38          5.844e-03    NA    NA
-k_JSE76          2.610e-02    NA    NA
-f_cyan_to_JCZ38  6.079e-01    NA    NA
-f_cyan_to_J9Z38  2.272e-01    NA    NA
+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.745e-01    NA    NA
-k1               1.429e-01    NA    NA
-k2               1.144e-02    NA    NA
-g                3.748e-01    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.6079
-cyan_J9Z38  0.2272
-cyan_sink   0.1649
+cyan_JCZ38  0.6051
+cyan_J9Z38  0.2251
+cyan_sink   0.1698
 JCZ38_JSE76 1.0000
 JCZ38_sink  0.0000
-JSE76_JCZ38 0.8745
-JSE76_sink  0.1255
+JSE76_JCZ38 0.8772
+JSE76_sink  0.1228
 
 Estimated disappearance times:
         DT50   DT90 DT50back DT50_k1 DT50_k2
-cyan   22.29 160.20    48.22    4.85   60.58
-JCZ38  13.06  43.37       NA      NA      NA
-J9Z38 118.61 394.02       NA      NA      NA
-JSE76  26.56  88.22       NA      NA      NA
+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

@@ -5070,10 +4768,10 @@ error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:16:47 2023 
-Date of summary: Thu Apr 20 20:01:31 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:57:32 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 *
@@ -5093,59 +4791,59 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 914.763 s
+Fitted in 865.636 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.7523        -1.5948        -5.0119        -2.2723         0.6719 
+      101.7393        -1.4493        -5.0118        -2.1269         0.6720 
   f_cyan_ilr_2 f_JCZ38_qlogis f_JSE76_qlogis         log_k1         log_k2 
-        5.1681        12.8238        12.4130        -2.0057        -4.5526 
+        7.3362        13.4423        13.2659        -2.0061        -4.5527 
       g_qlogis 
-       -0.5805 
+       -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.627       0.000       0.000       0.000       0.0000
-log_k_JCZ38     0.000       2.327       0.000       0.000       0.0000
-log_k_J9Z38     0.000       0.000       1.664       0.000       0.0000
-log_k_JSE76     0.000       0.000       0.000       4.566       0.0000
-f_cyan_ilr_1    0.000       0.000       0.000       0.000       0.6519
-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_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
+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.0           0.00           0.00 0.0000 0.0000
-log_k_JCZ38             0.0           0.00           0.00 0.0000 0.0000
-log_k_J9Z38             0.0           0.00           0.00 0.0000 0.0000
-log_k_JSE76             0.0           0.00           0.00 0.0000 0.0000
-f_cyan_ilr_1            0.0           0.00           0.00 0.0000 0.0000
-f_cyan_ilr_2           10.1           0.00           0.00 0.0000 0.0000
-f_JCZ38_qlogis          0.0          13.99           0.00 0.0000 0.0000
-f_JSE76_qlogis          0.0           0.00          14.15 0.0000 0.0000
-log_k1                  0.0           0.00           0.00 0.8452 0.0000
-log_k2                  0.0           0.00           0.00 0.0000 0.5968
-g_qlogis                0.0           0.00           0.00 0.0000 0.0000
+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.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.691
+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 
@@ -5155,78 +4853,78 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2232 2224  -1096
+  2237 2229  -1099
 
 Optimised parameters:
                      est.   lower   upper
-cyan_0          101.20051      NA      NA
-log_k_JCZ38      -2.93542      NA      NA
-log_k_J9Z38      -5.03151      NA      NA
-log_k_JSE76      -3.67679      NA      NA
-f_cyan_ilr_1      0.67290      NA      NA
-f_cyan_ilr_2      0.99787      NA      NA
-f_JCZ38_qlogis  348.32484      NA      NA
-f_JSE76_qlogis    1.87846      NA      NA
-log_k1           -2.32738      NA      NA
-log_k2           -4.61295      NA      NA
-g_qlogis         -0.38342      NA      NA
-a.1               2.06184 1.83746 2.28622
-b.1               0.06329 0.05211 0.07447
-SD.log_k_JCZ38    1.29042 0.47468 2.10617
-SD.log_k_J9Z38    0.84235 0.25903 1.42566
-SD.log_k_JSE76    0.56930 0.13934 0.99926
-SD.f_cyan_ilr_1   0.35183 0.12298 0.58068
-SD.f_cyan_ilr_2   0.77269 0.17908 1.36631
-SD.log_k2         0.28549 0.09210 0.47888
-SD.g_qlogis       0.93830 0.34568 1.53093
+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.88051 2.33623
+b.1               0.06223 0.05108 0.07338
+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.46740
+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.2904 0.4747 2.1062
-SD.log_k_J9Z38  0.8423 0.2590 1.4257
-SD.log_k_JSE76  0.5693 0.1393 0.9993
-SD.f_cyan_ilr_1 0.3518 0.1230 0.5807
-SD.f_cyan_ilr_2 0.7727 0.1791 1.3663
-SD.log_k2       0.2855 0.0921 0.4789
-SD.g_qlogis     0.9383 0.3457 1.5309
+                  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.06184 1.83746 2.28622
-b.1 0.06329 0.05211 0.07447
+a.1 2.10837 1.88051 2.33623
+b.1 0.06223 0.05108 0.07338
 
 Backtransformed parameters:
                       est. lower upper
-cyan_0           1.012e+02    NA    NA
-k_JCZ38          5.311e-02    NA    NA
-k_J9Z38          6.529e-03    NA    NA
-k_JSE76          2.530e-02    NA    NA
-f_cyan_to_JCZ38  6.373e-01    NA    NA
-f_cyan_to_J9Z38  2.461e-01    NA    NA
-f_JCZ38_to_JSE76 1.000e+00    NA    NA
-f_JSE76_to_JCZ38 8.674e-01    NA    NA
-k1               9.755e-02    NA    NA
-k2               9.922e-03    NA    NA
-g                4.053e-01    NA    NA
+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.6373
-cyan_J9Z38  0.2461
-cyan_sink   0.1167
-JCZ38_JSE76 1.0000
-JCZ38_sink  0.0000
-JSE76_JCZ38 0.8674
-JSE76_sink  0.1326
+                 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   24.93 179.68    54.09   7.105   69.86
-JCZ38  13.05  43.36       NA      NA      NA
-J9Z38 106.16 352.67       NA      NA      NA
-JSE76  27.39  91.00       NA      NA      NA
+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

@@ -5236,10 +4934,10 @@ variance 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:16:33 2023 
-Date of summary: Thu Apr 20 20:01:31 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:57:01 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -5257,20 +4955,20 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 901.179 s
+Fitted in 834.906 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.4394               -2.7673               -2.8942 
+             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 
-              13.2672               13.3538 
+              14.8408               15.4734 
 
 Fixed degradation parameter values:
 None
@@ -5301,17 +4999,17 @@ 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.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.6349        0.000           0.00           0.00
-f_cyan_ilr_2                0.0000        1.797           0.00           0.00
-f_JCZ38_qlogis              0.0000        0.000          13.84           0.00
-f_JSE76_qlogis              0.0000        0.000           0.00          14.66
+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 
@@ -5321,83 +5019,83 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2279 2272  -1120
+  2280 2272  -1120
 
 Optimised parameters:
                               est.   lower  upper
-cyan_free_0               102.5621      NA     NA
-log_k_cyan_free            -2.8531      NA     NA
-log_k_cyan_free_bound      -2.6916      NA     NA
-log_k_cyan_bound_free      -3.5032      NA     NA
-log_k_JCZ38                -2.9436      NA     NA
-log_k_J9Z38                -5.1140      NA     NA
-log_k_JSE76                -3.6472      NA     NA
-f_cyan_ilr_1                0.6887      NA     NA
-f_cyan_ilr_2                0.6874      NA     NA
-f_JCZ38_qlogis           4063.6389      NA     NA
-f_JSE76_qlogis              1.9556      NA     NA
-a.1                         2.7460 2.55451 2.9376
-SD.log_k_cyan_free          0.3131 0.09841 0.5277
-SD.log_k_cyan_free_bound    0.8850 0.29909 1.4710
-SD.log_k_cyan_bound_free    0.6167 0.20391 1.0295
-SD.log_k_JCZ38              1.3555 0.49101 2.2200
-SD.log_k_J9Z38              0.7200 0.16166 1.2783
-SD.log_k_JSE76              0.6252 0.14619 1.1042
-SD.f_cyan_ilr_1             0.3386 0.11447 0.5627
-SD.f_cyan_ilr_2             0.4699 0.09810 0.8417
+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.3131 0.09841 0.5277
-SD.log_k_cyan_free_bound 0.8850 0.29909 1.4710
-SD.log_k_cyan_bound_free 0.6167 0.20391 1.0295
-SD.log_k_JCZ38           1.3555 0.49101 2.2200
-SD.log_k_J9Z38           0.7200 0.16166 1.2783
-SD.log_k_JSE76           0.6252 0.14619 1.1042
-SD.f_cyan_ilr_1          0.3386 0.11447 0.5627
-SD.f_cyan_ilr_2          0.4699 0.09810 0.8417
+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.746 2.555 2.938
+a.1 2.747 2.556 2.939
 
 Backtransformed parameters:
                           est. lower upper
-cyan_free_0          1.026e+02    NA    NA
-k_cyan_free          5.767e-02    NA    NA
-k_cyan_free_bound    6.777e-02    NA    NA
-k_cyan_bound_free    3.010e-02    NA    NA
-k_JCZ38              5.267e-02    NA    NA
-k_J9Z38              6.012e-03    NA    NA
-k_JSE76              2.606e-02    NA    NA
-f_cyan_free_to_JCZ38 6.089e-01    NA    NA
-f_cyan_free_to_J9Z38 2.299e-01    NA    NA
+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.761e-01    NA    NA
+f_JSE76_to_JCZ38     8.807e-01    NA    NA
 
 Estimated Eigenvalues of SFORB model(s):
 cyan_b1 cyan_b2  cyan_g 
- 0.1434  0.0121  0.3469 
+ 0.1426  0.0121  0.3484 
 
 Resulting formation fractions:
                     ff
-cyan_free_JCZ38 0.6089
-cyan_free_J9Z38 0.2299
-cyan_free_sink  0.1612
+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.8761
-JSE76_sink      0.1239
+JSE76_JCZ38     0.8807
+JSE76_sink      0.1193
 
 Estimated disappearance times:
         DT50   DT90 DT50back DT50_cyan_b1 DT50_cyan_b2
-cyan   23.94 155.06    46.68        4.832        57.28
-JCZ38  13.16  43.71       NA           NA           NA
-J9Z38 115.30 383.02       NA           NA           NA
-JSE76  26.59  88.35       NA           NA           NA
+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

@@ -5407,10 +5105,10 @@ error 

 saemix version used for fitting:      3.2 
-mkin version used for pre-fitting:  1.2.3 
-R version used for fitting:         4.2.3 
-Date of fit:     Thu Apr 20 08:16:19 2023 
-Date of summary: Thu Apr 20 20:01:31 2023 
+mkin version used for pre-fitting:  1.2.4 
+R version used for fitting:         4.3.0 
+Date of fit:     Fri May 19 09:57:17 2023 
+Date of summary: Fri May 19 09:57:33 2023 
 
 Equations:
 d_cyan_free/dt = - k_cyan_free * cyan_free - k_cyan_free_bound *
@@ -5428,20 +5126,20 @@ Data:
 
 Model predictions using solution type deSolve 
 
-Fitted in 887.343 s
+Fitted in 850.751 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.751                -2.837                -3.016 
+             101.7511               -2.8370               -3.0162 
 log_k_cyan_bound_free           log_k_JCZ38           log_k_J9Z38 
-               -3.660                -2.299                -5.313 
+              -3.6600               -2.2988               -5.3129 
           log_k_JSE76          f_cyan_ilr_1          f_cyan_ilr_2 
-               -3.699                 0.672                 5.873 
+              -3.6991                0.6722                4.8596 
        f_JCZ38_qlogis        f_JSE76_qlogis 
-               13.216                13.338 
+              13.4678               14.2149 
 
 Fixed degradation parameter values:
 None
@@ -5472,17 +5170,17 @@ 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.00           0.00           0.00
-log_k_cyan_free             0.0000         0.00           0.00           0.00
-log_k_cyan_free_bound       0.0000         0.00           0.00           0.00
-log_k_cyan_bound_free       0.0000         0.00           0.00           0.00
-log_k_JCZ38                 0.0000         0.00           0.00           0.00
-log_k_J9Z38                 0.0000         0.00           0.00           0.00
-log_k_JSE76                 0.0000         0.00           0.00           0.00
-f_cyan_ilr_1                0.6519         0.00           0.00           0.00
-f_cyan_ilr_2                0.0000        10.78           0.00           0.00
-f_JCZ38_qlogis              0.0000         0.00          13.96           0.00
-f_JSE76_qlogis              0.0000         0.00           0.00          14.69
+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 
@@ -5492,83 +5190,83 @@ Results:
 
 Likelihood computed by importance sampling
    AIC  BIC logLik
-  2236 2228  -1098
+  2241 2233  -1101
 
 Optimised parameters:
                               est.   lower   upper
-cyan_free_0              100.72760      NA      NA
-log_k_cyan_free           -3.18281      NA      NA
-log_k_cyan_free_bound     -3.37924      NA      NA
-log_k_cyan_bound_free     -3.77107      NA      NA
-log_k_JCZ38               -2.92811      NA      NA
-log_k_J9Z38               -5.02759      NA      NA
-log_k_JSE76               -3.65835      NA      NA
-f_cyan_ilr_1               0.67390      NA      NA
-f_cyan_ilr_2               1.15106      NA      NA
-f_JCZ38_qlogis           827.82299      NA      NA
-f_JSE76_qlogis             1.83064      NA      NA
-a.1                        2.06921 1.84443 2.29399
-b.1                        0.06391 0.05267 0.07515
-SD.log_k_cyan_free         0.50518 0.18962 0.82075
-SD.log_k_cyan_bound_free   0.30991 0.08170 0.53813
-SD.log_k_JCZ38             1.26661 0.46578 2.06744
-SD.log_k_J9Z38             0.88272 0.27813 1.48730
-SD.log_k_JSE76             0.53050 0.12561 0.93538
-SD.f_cyan_ilr_1            0.35547 0.12461 0.58633
-SD.f_cyan_ilr_2            0.91446 0.20131 1.62761
+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.88299 2.33960
+b.1                        0.06299 0.05176 0.07421
+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.5052 0.1896 0.8207
-SD.log_k_cyan_bound_free 0.3099 0.0817 0.5381
-SD.log_k_JCZ38           1.2666 0.4658 2.0674
-SD.log_k_J9Z38           0.8827 0.2781 1.4873
-SD.log_k_JSE76           0.5305 0.1256 0.9354
-SD.f_cyan_ilr_1          0.3555 0.1246 0.5863
-SD.f_cyan_ilr_2          0.9145 0.2013 1.6276
+                           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.5488
+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.06921 1.84443 2.29399
-b.1 0.06391 0.05267 0.07515
+a.1 2.11130 1.88299 2.33960
+b.1 0.06299 0.05176 0.07421
 
 Backtransformed parameters:
                           est. lower upper
-cyan_free_0          1.007e+02    NA    NA
-k_cyan_free          4.147e-02    NA    NA
-k_cyan_free_bound    3.407e-02    NA    NA
-k_cyan_bound_free    2.303e-02    NA    NA
-k_JCZ38              5.350e-02    NA    NA
-k_J9Z38              6.555e-03    NA    NA
-k_JSE76              2.578e-02    NA    NA
-f_cyan_free_to_JCZ38 6.505e-01    NA    NA
-f_cyan_free_to_J9Z38 2.508e-01    NA    NA
-f_JCZ38_to_JSE76     1.000e+00    NA    NA
-f_JSE76_to_JCZ38     8.618e-01    NA    NA
+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.08768 0.01089 0.39821 
+0.08751 0.01101 0.39586 
 
 Resulting formation fractions:
                      ff
-cyan_free_JCZ38 0.65053
-cyan_free_J9Z38 0.25082
-cyan_free_sink  0.09864
+cyan_free_JCZ38 0.68444
+cyan_free_J9Z38 0.24633
+cyan_free_sink  0.06923
 cyan_free       1.00000
-JCZ38_JSE76     1.00000
-JCZ38_sink      0.00000
-JSE76_JCZ38     0.86184
-JSE76_sink      0.13816
+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.32 164.79    49.61        7.906        63.64
-JCZ38  12.96  43.04       NA           NA           NA
-J9Z38 105.75 351.29       NA           NA           NA
-JSE76  26.89  89.33       NA           NA           NA
+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

@@ -5577,13 +5275,13 @@ JSE76 26.89 89.33 NA NA NA

Session info

-
R version 4.2.3 (2023-03-15)
+
R version 4.3.0 Patched (2023-05-18 r84448)
 Platform: x86_64-pc-linux-gnu (64-bit)
 Running under: Debian GNU/Linux 12 (bookworm)
 
 Matrix products: default
-BLAS:   /usr/lib/x86_64-linux-gnu/openblas-serial/libblas.so.3
-LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/libopenblas-r0.3.21.so
+BLAS:   /home/jranke/svn/R/r-patched/build/lib/libRblas.so 
+LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/liblapack.so.3;  LAPACK version 3.11.0
 
 locale:
  [1] LC_CTYPE=de_DE.UTF-8       LC_NUMERIC=C              
@@ -5593,34 +5291,40 @@ locale:
  [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] saemix_3.2 npde_3.3   knitr_1.42 mkin_1.2.3
+[1] saemix_3.2 npde_3.3   knitr_1.42 mkin_1.2.4
 
 loaded via a namespace (and not attached):
- [1] pillar_1.9.0      bslib_0.4.2       compiler_4.2.3    jquerylib_0.1.4  
- [5] tools_4.2.3       mclust_6.0.0      digest_0.6.31     tibble_3.2.1     
- [9] jsonlite_1.8.4    evaluate_0.20     memoise_2.0.1     lifecycle_1.0.3  
-[13] nlme_3.1-162      gtable_0.3.3      lattice_0.21-8    pkgconfig_2.0.3  
-[17] rlang_1.1.0       DBI_1.1.3         cli_3.6.1         yaml_2.3.7       
-[21] pkgdown_2.0.7     xfun_0.38         fastmap_1.1.1     gridExtra_2.3    
-[25] dplyr_1.1.1       stringr_1.5.0     generics_0.1.3    desc_1.4.2       
-[29] fs_1.6.1          vctrs_0.6.1       sass_0.4.5        systemfonts_1.0.4
-[33] tidyselect_1.2.0  rprojroot_2.0.3   lmtest_0.9-40     grid_4.2.3       
-[37] inline_0.3.19     glue_1.6.2        R6_2.5.1          textshaping_0.3.6
-[41] fansi_1.0.4       rmarkdown_2.21    purrr_1.0.1       ggplot2_3.4.2    
-[45] magrittr_2.0.3    scales_1.2.1      htmltools_0.5.5   colorspace_2.1-0 
-[49] ragg_1.2.5        utf8_1.2.3        stringi_1.7.12    munsell_0.5.0    
-[53] cachem_1.0.7      zoo_1.8-12       

+ [1] sass_0.4.6        utf8_1.2.3        generics_0.1.3    stringi_1.7.12   
+ [5] lattice_0.21-8    digest_0.6.31     magrittr_2.0.3    evaluate_0.21    
+ [9] grid_4.3.0        fastmap_1.1.1     cellranger_1.1.0  rprojroot_2.0.3  
+[13] jsonlite_1.8.4    processx_3.8.1    pkgbuild_1.4.0    deSolve_1.35     
+[17] DBI_1.1.3         mclust_6.0.0      ps_1.7.5          gridExtra_2.3    
+[21] purrr_1.0.1       fansi_1.0.4       scales_1.2.1      codetools_0.2-19 
+[25] textshaping_0.3.6 jquerylib_0.1.4   cli_3.6.1         crayon_1.5.2     
+[29] rlang_1.1.1       munsell_0.5.0     cachem_1.0.8      yaml_2.3.7       
+[33] inline_0.3.19     tools_4.3.0       memoise_2.0.1     dplyr_1.1.2      
+[37] colorspace_2.1-0  ggplot2_3.4.2     vctrs_0.6.2       R6_2.5.1         
+[41] zoo_1.8-12        lifecycle_1.0.3   stringr_1.5.0     fs_1.6.2         
+[45] ragg_1.2.5        callr_3.7.3       pkgconfig_2.0.3   desc_1.4.2       
+[49] pkgdown_2.0.7     bslib_0.4.2       pillar_1.9.0      gtable_0.3.3     
+[53] glue_1.6.2        systemfonts_1.0.4 highr_0.10        xfun_0.39        
+[57] tibble_3.2.1      lmtest_0.9-40     tidyselect_1.2.0  htmltools_0.5.5  
+[61] nlme_3.1-162      rmarkdown_2.21    compiler_4.3.0    prettyunits_1.1.1
+[65] readxl_1.4.2

Hardware info 

CPU model: AMD Ryzen 9 7950X 16-Core Processor
-
MemTotal:       64936316 kB
+
MemTotal:       64925476 kB
diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png index b969f2ff..b3d6066b 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-13-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png index 60393da3..341d48b9 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-14-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png index b9a410f7..92df5b6c 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-15-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png index cf921dab..1a17b3bc 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-20-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png index ff732730..15905501 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-21-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png index e30011bc..cc7d0126 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-22-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png index 4aad76df..48bfc544 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-7-1.png differ diff --git a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png index e30011bc..cc7d0126 100644 Binary files a/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png and b/docs/articles/prebuilt/2022_cyan_pathway_files/figure-html/unnamed-chunk-8-1.png differ diff --git a/docs/articles/prebuilt/2022_dmta_parent.html b/docs/articles/prebuilt/2022_dmta_parent.html index 378a7e8e..2da41981 100644 --- a/docs/articles/prebuilt/2022_dmta_parent.html +++ b/docs/articles/prebuilt/2022_dmta_parent.html @@ -33,7 +33,7 @@ mkin - 1.2.3 + 1.2.4 @@ -135,7 +135,7 @@ with residue data on dimethenamid and dimethenamid-P Ranke

Last change on 5 January -2023, last compiled on 20 April 2023

+2023, last compiled on 19 Mai 2023 Source: vignettes/prebuilt/2022_dmta_parent.rmd @@ -154,7 +154,7 @@ FOMC, DFOP and HS can be fitted with the mkin package.

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

-

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

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

diff --git a/docs/articles/prebuilt/2022_dmta_pathway.html b/docs/articles/prebuilt/2022_dmta_pathway.html index c8323add..ea4bd087 100644 --- a/docs/articles/prebuilt/2022_dmta_pathway.html +++ b/docs/articles/prebuilt/2022_dmta_pathway.html @@ -33,7 +33,7 @@ mkin - 1.2.3.1 + 1.2.4 @@ -135,7 +135,7 @@ residue data on dimethenamid and dimethenamid-P Ranke

Last change on 20 April 2023, -last compiled on 20 April 2023

+last compiled on 19 Mai 2023 Source: vignettes/prebuilt/2022_dmta_pathway.rmd @@ -155,7 +155,7 @@ can be fitted with the mkin package.

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

-

The mkin package is used in version 1.2.3, which is currently under +

The mkin package is used in version 1.2.4, which is currently under development. It contains the test data, and the functions used in the evaluations. The saemix package is used as a backend for fitting the NLHM, but is also loaded to make the convergence plot @@ -1496,7 +1496,7 @@ of parent models tested here.

dfop_path_1 OK OK -C +OK OK OK OK @@ -1514,7 +1514,7 @@ of parent models tested here.

hs_path_1 C C -C +OK C C C @@ -1552,24 +1552,24 @@ not converge with default settings.

fomc_path_1 OK OK -C OK OK -C +OK +OK dfop_path_1 OK C OK -OK +C OK OK sforb_path_1 OK -C +OK OK OK OK @@ -1582,7 +1582,7 @@ not converge with default settings.

C C C -OK +C @@ -1643,8 +1643,6 @@ than twenty minutes.

successfully.

 anova(saem_1) |> kable(digits = 1)
-
Warning in FUN(X[[i]], ...): Could not obtain log likelihood with 'is' method
-for sforb_path_1 const
@@ -1664,58 +1662,65 @@ for sforb_path_1 const - - - + + + - - - + + + - - - + + + - - - + + + - + - - - + + + + + + + + + + - - - + + + - + - - - + + + - + - - - + + +
sfo_path_1 tc 182276.32272.5-1120.12276.42272.7-1120.2
fomc_path_1 const 192099.02095.0-1030.52095.92091.9-1028.9
fomc_path_1 tc 201939.61935.5-949.81939.01934.8-949.5
dfop_path_1 const 212038.82034.4-998.42039.72035.3-998.8
hs_path_1 constsforb_path_1 const 212024.22019.8-991.12017.72013.4-987.9
hs_path_1 const212023.72019.3-990.9
dfop_path_1 tc 221879.81875.2-917.91881.71877.1-918.9
sforb_path_1 tc 221832.91828.3-894.41832.71828.1-894.3
hs_path_1 tc 221831.41826.8-893.71831.61827.0-893.8
@@ -1740,7 +1745,7 @@ Matrix
+
 illparms(saem_1) |> kable()

 
@@ -1766,7 +1771,7 @@ found.

-
+
@@ -1786,13 +1791,13 @@ two-component error, the random effect for the rate constant from
 reversibly bound DMTA to the free DMTA (k_DMTA_bound_free)
 is not well-defined. Therefore, the fit is updated without assuming a
 random effect for this parameter.

-
+
 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)

 
As expected, no ill-defined parameters remain. The model comparison
 below shows that the reduced model is preferable.

-
+
 anova(saem_1[["sforb_path_1", "tc"]], saem_sforb_path_1_tc_reduced) |> kable(digits = 1)

 

 

 

 

 sforb_path_1sd(log_k_DMTA_bound_free)
 sd(log_k_DMTA_bound_free)
 

 

@@ -1806,21 +1811,21 @@ below shows that the reduced model is preferable.

-
-
+
+
-
-
-
+
+
+

 

 

 saem_sforb_path_1_tc_reduced
 211830.31825.91830.41826.0
 -894.2
 

 

 saem_1[[“sforb_path_1”, “tc”]]
 221832.91828.3-894.41832.71828.1-894.3
 

 
 

 
The convergence plot of the refined fit is shown below.

-
+
 plot(saem_sforb_path_1_tc_reduced$so, plot.type = "convergence")

 

 
For some parameters, for example for f_DMTA_ilr_1 and
@@ -1839,14 +1844,14 @@ saemix.

 
 
As an alternative check of parameter identifiability (Duchesne et al. 2021), multistart runs were
 performed on the basis of the refined fit shown above.

-
+
 saem_sforb_path_1_tc_reduced_multi <- multistart(saem_sforb_path_1_tc_reduced,
   n = 32, cores = 10)

-
+
 print(saem_sforb_path_1_tc_reduced_multi)

 
<multistart> object with 32 fits:
  E OK 
-15 17 
+ 7 25 
 OK: Fit terminated successfully
 E: Error

 
Out of the 32 fits that were initiated, only 17 terminated without an
@@ -1856,7 +1861,7 @@ the SAEM algorithm leads to parameter combinations for the degradation
 model that the numerical integration routine cannot cope with. Because
 of this variation of initial parameters, some of the model fits take up
 to two times more time than the original fit.

-
+
 par(mar = c(12.1, 4.1, 2.1, 2.1))
 parplot(saem_sforb_path_1_tc_reduced_multi, ylim = c(0.5, 2), las = 2)

 

@@ -1875,14 +1880,14 @@ ill-defined parameters.

 
 
The SFORB pathway fits with full and reduced parameter distribution
 model are shown below.

-
+
 plot(saem_1[["sforb_path_1", "tc"]])

 

 SFORB pathway fit with two-component error

 SFORB pathway fit with two-component error
 

 

-
+
 plot(saem_sforb_path_1_tc_reduced)

 

 SFORB pathway fit with two-component error, reduced parameter model

@@ -1891,7 +1896,7 @@ SFORB pathway fit with two-component error, reduced parameter model
 

 
Plots of the remaining fits and listings for all successful fits are
 shown in the Appendix.

-
 

@@ -1931,21 +1936,21 @@ Level with Nonlinear Mixed-Effects Models.” Environments
 

 
Plots of hierarchical fits not selected for refinement
 

-
+
 plot(saem_1[["sfo_path_1", "tc"]])

 

 SFO pathway fit with two-component error

 SFO pathway fit with two-component error
 

 

-
+
 plot(saem_1[["fomc_path_1", "tc"]])

 

 FOMC pathway fit with two-component error

 FOMC pathway fit with two-component error
 

 

-
+
 plot(saem_1[["sforb_path_1", "tc"]])

 

 HS pathway fit with two-component error

@@ -1971,13 +1976,13 @@ error
 
Session info
 

-
R version 4.2.3 (2023-03-15)
+
R version 4.3.0 Patched (2023-05-18 r84448)
 Platform: x86_64-pc-linux-gnu (64-bit)
 Running under: Debian GNU/Linux 12 (bookworm)
 
 Matrix products: default
-BLAS:   /usr/lib/x86_64-linux-gnu/openblas-serial/libblas.so.3
-LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/libopenblas-r0.3.21.so
+BLAS:   /home/jranke/svn/R/r-patched/build/lib/libRblas.so 
+LAPACK: /usr/lib/x86_64-linux-gnu/openblas-serial/liblapack.so.3;  LAPACK version 3.11.0
 
 locale:
  [1] LC_CTYPE=de_DE.UTF-8       LC_NUMERIC=C              
@@ -1987,36 +1992,39 @@ locale:
  [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] saemix_3.2 npde_3.3   knitr_1.42 mkin_1.2.3
+[1] saemix_3.2 npde_3.3   knitr_1.42 mkin_1.2.4
 
 loaded via a namespace (and not attached):
- [1] deSolve_1.35      zoo_1.8-12        tidyselect_1.2.0  xfun_0.38        
- [5] bslib_0.4.2       purrr_1.0.1       lattice_0.21-8    colorspace_2.1-0 
- [9] vctrs_0.6.1       generics_0.1.3    htmltools_0.5.5   yaml_2.3.7       
-[13] utf8_1.2.3        rlang_1.1.0       pkgbuild_1.4.0    pkgdown_2.0.7    
-[17] jquerylib_0.1.4   pillar_1.9.0      glue_1.6.2        DBI_1.1.3        
-[21] lifecycle_1.0.3   stringr_1.5.0     munsell_0.5.0     gtable_0.3.3     
-[25] ragg_1.2.5        codetools_0.2-19  memoise_2.0.1     evaluate_0.20    
-[29] inline_0.3.19     callr_3.7.3       fastmap_1.1.1     ps_1.7.4         
-[33] lmtest_0.9-40     fansi_1.0.4       highr_0.10        scales_1.2.1     
-[37] cachem_1.0.7      desc_1.4.2        jsonlite_1.8.4    systemfonts_1.0.4
-[41] fs_1.6.1          textshaping_0.3.6 gridExtra_2.3     ggplot2_3.4.2    
-[45] digest_0.6.31     stringi_1.7.12    processx_3.8.0    dplyr_1.1.1      
-[49] grid_4.2.3        rprojroot_2.0.3   cli_3.6.1         tools_4.2.3      
-[53] magrittr_2.0.3    sass_0.4.5        tibble_3.2.1      crayon_1.5.2     
-[57] pkgconfig_2.0.3   prettyunits_1.1.1 rmarkdown_2.21    R6_2.5.1         
-[61] mclust_6.0.0      nlme_3.1-162      compiler_4.2.3   

+ [1] sass_0.4.6        utf8_1.2.3        generics_0.1.3    stringi_1.7.12   
+ [5] lattice_0.21-8    digest_0.6.31     magrittr_2.0.3    evaluate_0.21    
+ [9] grid_4.3.0        fastmap_1.1.1     rprojroot_2.0.3   jsonlite_1.8.4   
+[13] processx_3.8.1    pkgbuild_1.4.0    deSolve_1.35      DBI_1.1.3        
+[17] mclust_6.0.0      ps_1.7.5          gridExtra_2.3     purrr_1.0.1      
+[21] fansi_1.0.4       scales_1.2.1      codetools_0.2-19  textshaping_0.3.6
+[25] jquerylib_0.1.4   cli_3.6.1         crayon_1.5.2      rlang_1.1.1      
+[29] munsell_0.5.0     cachem_1.0.8      yaml_2.3.7        inline_0.3.19    
+[33] tools_4.3.0       memoise_2.0.1     dplyr_1.1.2       colorspace_2.1-0 
+[37] ggplot2_3.4.2     vctrs_0.6.2       R6_2.5.1          zoo_1.8-12       
+[41] lifecycle_1.0.3   stringr_1.5.0     fs_1.6.2          ragg_1.2.5       
+[45] callr_3.7.3       pkgconfig_2.0.3   desc_1.4.2        pkgdown_2.0.7    
+[49] bslib_0.4.2       pillar_1.9.0      gtable_0.3.3      glue_1.6.2       
+[53] systemfonts_1.0.4 highr_0.10        xfun_0.39         tibble_3.2.1     
+[57] lmtest_0.9-40     tidyselect_1.2.0  htmltools_0.5.5   nlme_3.1-162     
+[61] rmarkdown_2.21    compiler_4.3.0    prettyunits_1.1.1

 

 

 
Hardware info
 

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

-
MemTotal:       64936316 kB

+
MemTotal:       64925476 kB

 

 

   

diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png
index 206c424d..a7187abc 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/saem-sforb-path-1-tc-reduced-convergence-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png
index 0fe084d3..ba36beb2 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-2-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png
index 1c81601e..80376942 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-3-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png
index e0961dce..15a30cfb 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-4-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png
index 00db0c76..5445ddb5 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-6-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png
index ac5271ec..067b144d 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-7-1.png differ
diff --git a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png
index 1c81601e..80376942 100644
Binary files a/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png and b/docs/articles/prebuilt/2022_dmta_pathway_files/figure-html/unnamed-chunk-8-1.png differ
-- 
cgit v1.2.1