diff options
| -rw-r--r-- | NEWS.md | 2 | ||||
| -rw-r--r-- | R/parplot.R | 16 | ||||
| -rw-r--r-- | log/test.log | 41 | ||||
| -rw-r--r-- | man/parplot.Rd | 7 | ||||
| -rw-r--r-- | tests/testthat/_snaps/multistart/parplot-for-dfop-sfo-fit.svg | 169 | ||||
| -rw-r--r-- | tests/testthat/_snaps/multistart/parplot-for-sfo-fit.svg | 4 | ||||
| -rw-r--r-- | tests/testthat/test_multistart.R | 2 | ||||
| -rw-r--r-- | vignettes/FOCUS_D.html | 37 | ||||
| -rw-r--r-- | vignettes/FOCUS_L.html | 83 |
9 files changed, 183 insertions, 178 deletions
@@ -4,6 +4,8 @@ - 'R/summary.saem.mmkin.R': List all initial parameter values in the summary, including random effects and error model parameters. Avoid redundant warnings that occurred in the calculation of correlations of the fixed effects in the case that the Fisher information matrix could not be inverted. +- 'R/parplot.R': Possibility to select the top 'llquant' fraction of the fits for the parameter plots, and improved legend text. + # mkin 1.2.1 (2022-11-19) - '{data,R}/ds_mixed.rda': Include the test data in the package instead of generating it in 'tests/testthat/setup_script.R'. Refactor the generating code to make it consistent and update tests. diff --git a/R/parplot.R b/R/parplot.R index 98579779..63306ac2 100644 --- a/R/parplot.R +++ b/R/parplot.R @@ -10,7 +10,10 @@ #' #' @param object The [multistart] object #' @param llmin The minimum likelihood of objects to be shown -#' @param scale By default, scale parameters using the best available fit. +#' @param llquant Fractional value for selecting only the fits with higher +#' likelihoods. Overrides 'llmin'. +#' @param scale By default, scale parameters using the best +#' available fit. #' If 'median', parameters are scaled using the median parameters from all fits. #' @param main Title of the plot #' @param lpos Positioning of the legend. @@ -28,7 +31,8 @@ parplot <- function(object, ...) { #' @rdname parplot #' @export -parplot.multistart.saem.mmkin <- function(object, llmin = -Inf, scale = c("best", "median"), +parplot.multistart.saem.mmkin <- function(object, llmin = -Inf, llquant = NA, + scale = c("best", "median"), lpos = "bottomleft", main = "", ...) { oldpar <- par(no.readonly = TRUE) @@ -48,6 +52,10 @@ parplot.multistart.saem.mmkin <- function(object, llmin = -Inf, scale = c("best" stop("parplot is only implemented for multistart.saem.mmkin objects") } ll <- sapply(object, llfunc) + if (!is.na(llquant[1])) { + if (llmin != -Inf) warning("Overriding 'llmin' because 'llquant' was specified") + llmin <- quantile(ll, 1 - llquant) + } selected <- which(ll > llmin) selected_parms <- all_parms[selected, ] @@ -110,7 +118,7 @@ parplot.multistart.saem.mmkin <- function(object, llmin = -Inf, scale = c("best" legend(lpos, inset = c(0.05, 0.05), bty = "n", pch = 1, col = 3:1, lty = c(NA, NA, 1), legend = c( - "Starting parameters", - "Original run", + "Original start", + "Original results", "Multistart runs")) } diff --git a/log/test.log b/log/test.log index 5764f209..7614b136 100644 --- a/log/test.log +++ b/log/test.log @@ -1,22 +1,22 @@ ℹ Testing mkin ✔ | F W S OK | Context ✔ | 5 | AIC calculation -✔ | 5 | Analytical solutions for coupled models [3.3s] +✔ | 5 | Analytical solutions for coupled models [3.2s] ✔ | 5 | Calculation of Akaike weights ✔ | 3 | Export dataset for reading into CAKE ✔ | 12 | Confidence intervals and p-values [1.1s] -✔ | 1 12 | Dimethenamid data from 2018 [31.9s] +✔ | 1 12 | Dimethenamid data from 2018 [32.0s] ──────────────────────────────────────────────────────────────────────────────── Skip ('test_dmta.R:98'): Different backends get consistent results for SFO-SFO3+, dimethenamid data Reason: Fitting this ODE model with saemix takes about 15 minutes on my system ──────────────────────────────────────────────────────────────────────────────── -✔ | 14 | Error model fitting [4.9s] +✔ | 14 | Error model fitting [4.8s] ✔ | 5 | Time step normalisation ✔ | 4 | Calculation of FOCUS chi2 error levels [0.6s] ✔ | 14 | Results for FOCUS D established in expertise for UBA (Ranke 2014) [0.8s] ✔ | 4 | Test fitting the decline of metabolites from their maximum [0.3s] ✔ | 1 | Fitting the logistic model [0.2s] -✔ | 10 | Batch fitting and diagnosing hierarchical kinetic models [41.4s] +✔ | 10 | Batch fitting and diagnosing hierarchical kinetic models [42.2s] ✔ | 1 11 | Nonlinear mixed-effects models [13.3s] ──────────────────────────────────────────────────────────────────────────────── Skip ('test_mixed.R:78'): saemix results are reproducible for biphasic fits @@ -26,32 +26,41 @@ Reason: Fitting with saemix takes around 10 minutes when using deSolve ✔ | 10 | Special cases of mkinfit calls [0.6s] ✔ | 3 | mkinfit features [0.7s] ✔ | 8 | mkinmod model generation and printing [0.2s] -✔ | 3 | Model predictions with mkinpredict [0.3s] -✔ | 12 | Multistart method for saem.mmkin models [47.5s] -✔ | 16 | Evaluations according to 2015 NAFTA guidance [2.3s] -✔ | 9 | Nonlinear mixed-effects models with nlme [9.5s] -✔ | 15 | Plotting [10.2s] +✔ | 3 | Model predictions with mkinpredict [0.4s] +✖ | 1 11 | Multistart method for saem.mmkin models [46.7s] +──────────────────────────────────────────────────────────────────────────────── +Failure ('test_multistart.R:56'): multistart works for saem.mmkin models +Snapshot of `testcase` to 'multistart/parplot-for-dfop-sfo-fit.svg' has changed +Run `testthat::snapshot_review('multistart/')` to review changes +Backtrace: + 1. vdiffr::expect_doppelganger("parplot for dfop sfo fit", parplot_dfop_sfo) + at test_multistart.R:56:2 + 3. testthat::expect_snapshot_file(...) +──────────────────────────────────────────────────────────────────────────────── +✔ | 16 | Evaluations according to 2015 NAFTA guidance [2.2s] +✔ | 9 | Nonlinear mixed-effects models with nlme [9.4s] +✔ | 15 | Plotting [10.1s] ✔ | 4 | Residuals extracted from mkinfit models -✔ | 1 36 | saemix parent models [73.0s] +✔ | 1 36 | saemix parent models [73.6s] ──────────────────────────────────────────────────────────────────────────────── Skip ('test_saemix_parent.R:143'): We can also use mkin solution methods for saem Reason: This still takes almost 2.5 minutes although we do not solve ODEs ──────────────────────────────────────────────────────────────────────────────── -✔ | 2 | Complex test case from Schaefer et al. (2007) Piacenza paper [1.5s] +✔ | 2 | Complex test case from Schaefer et al. (2007) Piacenza paper [1.4s] ✔ | 11 | Processing of residue series -✔ | 10 | Fitting the SFORB model [3.9s] +✔ | 10 | Fitting the SFORB model [3.7s] ✔ | 1 | Summaries of old mkinfit objects ✔ | 5 | Summary [0.2s] -✔ | 4 | Results for synthetic data established in expertise for UBA (Ranke 2014) [2.2s] -✔ | 9 | Hypothesis tests [8.4s] +✔ | 4 | Results for synthetic data established in expertise for UBA (Ranke 2014) [2.3s] +✔ | 9 | Hypothesis tests [8.1s] ✔ | 4 | Calculation of maximum time weighted average concentrations (TWAs) [2.2s] ══ Results ═════════════════════════════════════════════════════════════════════ -Duration: 261.1 s +Duration: 260.9 s ── Skipped tests ────────────────────────────────────────────────────────────── • Fitting this ODE model with saemix takes about 15 minutes on my system (1) • Fitting with saemix takes around 10 minutes when using deSolve (1) • This still takes almost 2.5 minutes although we do not solve ODEs (1) -[ FAIL 0 | WARN 0 | SKIP 3 | PASS 270 ] +[ FAIL 1 | WARN 0 | SKIP 3 | PASS 269 ] diff --git a/man/parplot.Rd b/man/parplot.Rd index ac9e02cf..67bf0cc1 100644 --- a/man/parplot.Rd +++ b/man/parplot.Rd @@ -10,6 +10,7 @@ parplot(object, ...) \method{parplot}{multistart.saem.mmkin}( object, llmin = -Inf, + llquant = NA, scale = c("best", "median"), lpos = "bottomleft", main = "", @@ -23,7 +24,11 @@ parplot(object, ...) \item{llmin}{The minimum likelihood of objects to be shown} -\item{scale}{By default, scale parameters using the best available fit. +\item{llquant}{Fractional value for selecting only the fits with higher +likelihoods. Overrides 'llmin'.} + +\item{scale}{By default, scale parameters using the best +available fit. If 'median', parameters are scaled using the median parameters from all fits.} \item{lpos}{Positioning of the legend.} diff --git a/tests/testthat/_snaps/multistart/parplot-for-dfop-sfo-fit.svg b/tests/testthat/_snaps/multistart/parplot-for-dfop-sfo-fit.svg index 7017908e..b01dac74 100644 --- a/tests/testthat/_snaps/multistart/parplot-for-dfop-sfo-fit.svg +++ b/tests/testthat/_snaps/multistart/parplot-for-dfop-sfo-fit.svg @@ -25,109 +25,104 @@ </clipPath> </defs> <g clip-path='url(#cpNTkuMDR8Njg5Ljc2fDU5LjA0fDUwMi41Ng==)'> -<polygon points='86.57,280.95 119.94,280.95 119.94,280.68 86.57,280.68 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='86.57' y1='280.78' x2='119.94' y2='280.78' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='103.26' y1='281.01' x2='103.26' y2='280.95' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='103.26' y1='280.63' x2='103.26' y2='280.68' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='86.57,280.91 119.94,280.91 119.94,280.26 86.57,280.26 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='86.57' y1='280.71' x2='119.94' y2='280.71' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='103.26' y1='281.01' x2='103.26' y2='280.91' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='103.26' y1='279.89' x2='103.26' y2='280.26' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='94.91' y1='281.01' x2='111.60' y2='281.01' style='stroke-width: 0.75;' /> -<line x1='94.91' y1='280.63' x2='111.60' y2='280.63' style='stroke-width: 0.75;' /> -<polygon points='86.57,280.95 119.94,280.95 119.94,280.68 86.57,280.68 ' style='stroke-width: 0.75; fill: none;' /> -<circle cx='103.26' cy='281.66' r='2.70' style='stroke-width: 0.75;' /> -<circle cx='103.26' cy='279.89' r='2.70' style='stroke-width: 0.75;' /> -<polygon points='128.29,282.46 161.66,282.46 161.66,280.33 128.29,280.33 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='128.29' y1='280.88' x2='161.66' y2='280.88' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='144.97' y1='282.96' x2='144.97' y2='282.46' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='144.97' y1='279.72' x2='144.97' y2='280.33' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='94.91' y1='279.89' x2='111.60' y2='279.89' style='stroke-width: 0.75;' /> +<polygon points='86.57,280.91 119.94,280.91 119.94,280.26 86.57,280.26 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='128.29,281.88 161.66,281.88 161.66,280.33 128.29,280.33 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='128.29' y1='280.61' x2='161.66' y2='280.61' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='144.97' y1='282.96' x2='144.97' y2='281.88' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='144.97' y1='280.24' x2='144.97' y2='280.33' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='136.63' y1='282.96' x2='153.31' y2='282.96' style='stroke-width: 0.75;' /> -<line x1='136.63' y1='279.72' x2='153.31' y2='279.72' style='stroke-width: 0.75;' /> -<polygon points='128.29,282.46 161.66,282.46 161.66,280.33 128.29,280.33 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='170.00,281.20 203.37,281.20 203.37,280.53 170.00,280.53 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='170.00' y1='280.84' x2='203.37' y2='280.84' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='186.69' y1='281.48' x2='186.69' y2='281.20' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='136.63' y1='280.24' x2='153.31' y2='280.24' style='stroke-width: 0.75;' /> +<polygon points='128.29,281.88 161.66,281.88 161.66,280.33 128.29,280.33 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='170.00,281.69 203.37,281.69 203.37,280.53 170.00,280.53 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='170.00' y1='280.75' x2='203.37' y2='280.75' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='186.69' y1='282.58' x2='186.69' y2='281.69' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='186.69' y1='280.35' x2='186.69' y2='280.53' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='178.34' y1='281.48' x2='195.03' y2='281.48' style='stroke-width: 0.75;' /> +<line x1='178.34' y1='282.58' x2='195.03' y2='282.58' style='stroke-width: 0.75;' /> <line x1='178.34' y1='280.35' x2='195.03' y2='280.35' style='stroke-width: 0.75;' /> -<polygon points='170.00,281.20 203.37,281.20 203.37,280.53 170.00,280.53 ' style='stroke-width: 0.75; fill: none;' /> -<circle cx='186.69' cy='279.30' r='2.70' style='stroke-width: 0.75;' /> -<circle cx='186.69' cy='282.58' r='2.70' style='stroke-width: 0.75;' /> -<polygon points='211.71,281.98 245.09,281.98 245.09,280.03 211.71,280.03 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='211.71' y1='281.46' x2='245.09' y2='281.46' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='228.40' y1='282.40' x2='228.40' y2='281.98' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='228.40' y1='279.21' x2='228.40' y2='280.03' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='220.06' y1='282.40' x2='236.74' y2='282.40' style='stroke-width: 0.75;' /> -<line x1='220.06' y1='279.21' x2='236.74' y2='279.21' style='stroke-width: 0.75;' /> -<polygon points='211.71,281.98 245.09,281.98 245.09,280.03 211.71,280.03 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='253.43,281.94 286.80,281.94 286.80,280.63 253.43,280.63 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='253.43' y1='281.42' x2='286.80' y2='281.42' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='270.11' y1='282.62' x2='270.11' y2='281.94' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='270.11' y1='280.34' x2='270.11' y2='280.63' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='170.00,281.69 203.37,281.69 203.37,280.53 170.00,280.53 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='211.71,281.83 245.09,281.83 245.09,280.03 211.71,280.03 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='211.71' y1='281.14' x2='245.09' y2='281.14' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='228.40' y1='282.17' x2='228.40' y2='281.83' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='228.40' y1='279.26' x2='228.40' y2='280.03' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='220.06' y1='282.17' x2='236.74' y2='282.17' style='stroke-width: 0.75;' /> +<line x1='220.06' y1='279.26' x2='236.74' y2='279.26' style='stroke-width: 0.75;' /> +<polygon points='211.71,281.83 245.09,281.83 245.09,280.03 211.71,280.03 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='253.43,282.21 286.80,282.21 286.80,280.57 253.43,280.57 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='253.43' y1='281.30' x2='286.80' y2='281.30' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='270.11' y1='282.62' x2='270.11' y2='282.21' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='270.11' y1='280.34' x2='270.11' y2='280.57' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='261.77' y1='282.62' x2='278.46' y2='282.62' style='stroke-width: 0.75;' /> <line x1='261.77' y1='280.34' x2='278.46' y2='280.34' style='stroke-width: 0.75;' /> -<polygon points='253.43,281.94 286.80,281.94 286.80,280.63 253.43,280.63 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='295.14,281.17 328.51,281.17 328.51,278.99 295.14,278.99 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='295.14' y1='280.18' x2='328.51' y2='280.18' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='311.83' y1='282.20' x2='311.83' y2='281.17' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='311.83' y1='277.75' x2='311.83' y2='278.99' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='253.43,282.21 286.80,282.21 286.80,280.57 253.43,280.57 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='295.14,281.50 328.51,281.50 328.51,278.14 295.14,278.14 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='295.14' y1='279.67' x2='328.51' y2='279.67' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='311.83' y1='282.20' x2='311.83' y2='281.50' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='311.83' y1='277.75' x2='311.83' y2='278.14' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='303.49' y1='282.20' x2='320.17' y2='282.20' style='stroke-width: 0.75;' /> <line x1='303.49' y1='277.75' x2='320.17' y2='277.75' style='stroke-width: 0.75;' /> -<polygon points='295.14,281.17 328.51,281.17 328.51,278.99 295.14,278.99 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='336.86,281.95 370.23,281.95 370.23,281.00 336.86,281.00 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='336.86' y1='281.48' x2='370.23' y2='281.48' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='353.54' y1='282.69' x2='353.54' y2='281.95' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='353.54' y1='280.80' x2='353.54' y2='281.00' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='345.20' y1='282.69' x2='361.89' y2='282.69' style='stroke-width: 0.75;' /> +<polygon points='295.14,281.50 328.51,281.50 328.51,278.14 295.14,278.14 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='336.86,281.52 370.23,281.52 370.23,280.88 336.86,280.88 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='336.86' y1='281.06' x2='370.23' y2='281.06' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='353.54' y1='281.89' x2='353.54' y2='281.52' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='353.54' y1='280.80' x2='353.54' y2='280.88' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='345.20' y1='281.89' x2='361.89' y2='281.89' style='stroke-width: 0.75;' /> <line x1='345.20' y1='280.80' x2='361.89' y2='280.80' style='stroke-width: 0.75;' /> -<polygon points='336.86,281.95 370.23,281.95 370.23,281.00 336.86,281.00 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='378.57,280.93 411.94,280.93 411.94,280.50 378.57,280.50 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='378.57' y1='280.72' x2='411.94' y2='280.72' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='395.26' y1='281.09' x2='395.26' y2='280.93' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='395.26' y1='280.44' x2='395.26' y2='280.50' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='336.86,281.52 370.23,281.52 370.23,280.88 336.86,280.88 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='378.57,281.06 411.94,281.06 411.94,280.72 378.57,280.72 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='378.57' y1='280.91' x2='411.94' y2='280.91' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='395.26' y1='281.09' x2='395.26' y2='281.06' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='395.26' y1='280.64' x2='395.26' y2='280.72' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='386.91' y1='281.09' x2='403.60' y2='281.09' style='stroke-width: 0.75;' /> -<line x1='386.91' y1='280.44' x2='403.60' y2='280.44' style='stroke-width: 0.75;' /> -<polygon points='378.57,280.93 411.94,280.93 411.94,280.50 378.57,280.50 ' style='stroke-width: 0.75; fill: none;' /> -<circle cx='395.26' cy='279.72' r='2.70' style='stroke-width: 0.75;' /> -<polygon points='420.29,409.00 453.66,409.00 453.66,-49.57 420.29,-49.57 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='420.29' y1='94.54' x2='453.66' y2='94.54' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='436.97' y1='656.81' x2='436.97' y2='409.00' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='436.97' y1='-136.13' x2='436.97' y2='-49.57' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='428.63' y1='656.81' x2='445.31' y2='656.81' style='stroke-width: 0.75;' /> +<line x1='386.91' y1='280.64' x2='403.60' y2='280.64' style='stroke-width: 0.75;' /> +<polygon points='378.57,281.06 411.94,281.06 411.94,280.72 378.57,280.72 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='420.29,106.59 453.66,106.59 453.66,-114.09 420.29,-114.09 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='420.29' y1='-49.57' x2='453.66' y2='-49.57' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='436.97' y1='280.80' x2='436.97' y2='106.59' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='436.97' y1='-136.13' x2='436.97' y2='-114.09' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='428.63' y1='280.80' x2='445.31' y2='280.80' style='stroke-width: 0.75;' /> <line x1='428.63' y1='-136.13' x2='445.31' y2='-136.13' style='stroke-width: 0.75;' /> -<polygon points='420.29,409.00 453.66,409.00 453.66,-49.57 420.29,-49.57 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='462.00,280.54 495.37,280.54 495.37,275.74 462.00,275.74 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='462.00' y1='278.14' x2='495.37' y2='278.14' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='478.69' y1='281.75' x2='478.69' y2='280.54' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='478.69' y1='272.26' x2='478.69' y2='275.74' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='420.29,106.59 453.66,106.59 453.66,-114.09 420.29,-114.09 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='462.00,281.28 495.37,281.28 495.37,276.14 462.00,276.14 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='462.00' y1='278.79' x2='495.37' y2='278.79' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='478.69' y1='281.75' x2='478.69' y2='281.28' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='478.69' y1='275.48' x2='478.69' y2='276.14' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='470.34' y1='281.75' x2='487.03' y2='281.75' style='stroke-width: 0.75;' /> -<line x1='470.34' y1='272.26' x2='487.03' y2='272.26' style='stroke-width: 0.75;' /> -<polygon points='462.00,280.54 495.37,280.54 495.37,275.74 462.00,275.74 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='503.71,282.63 537.09,282.63 537.09,280.69 503.71,280.69 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='503.71' y1='282.10' x2='537.09' y2='282.10' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='520.40' y1='283.13' x2='520.40' y2='282.63' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='520.40' y1='280.04' x2='520.40' y2='280.69' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='470.34' y1='275.48' x2='487.03' y2='275.48' style='stroke-width: 0.75;' /> +<polygon points='462.00,281.28 495.37,281.28 495.37,276.14 462.00,276.14 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='503.71,282.81 537.09,282.81 537.09,281.44 503.71,281.44 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='503.71' y1='282.28' x2='537.09' y2='282.28' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='520.40' y1='283.13' x2='520.40' y2='282.81' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='520.40' y1='280.80' x2='520.40' y2='281.44' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='512.06' y1='283.13' x2='528.74' y2='283.13' style='stroke-width: 0.75;' /> -<line x1='512.06' y1='280.04' x2='528.74' y2='280.04' style='stroke-width: 0.75;' /> -<polygon points='503.71,282.63 537.09,282.63 537.09,280.69 503.71,280.69 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='545.43,283.04 578.80,283.04 578.80,278.99 545.43,278.99 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='545.43' y1='281.15' x2='578.80' y2='281.15' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='562.11' y1='283.53' x2='562.11' y2='283.04' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='562.11' y1='275.97' x2='562.11' y2='278.99' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='512.06' y1='280.80' x2='528.74' y2='280.80' style='stroke-width: 0.75;' /> +<polygon points='503.71,282.81 537.09,282.81 537.09,281.44 503.71,281.44 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='545.43,282.16 578.80,282.16 578.80,276.63 545.43,276.63 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='545.43' y1='279.04' x2='578.80' y2='279.04' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='562.11' y1='283.53' x2='562.11' y2='282.16' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='562.11' y1='275.97' x2='562.11' y2='276.63' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='553.77' y1='283.53' x2='570.46' y2='283.53' style='stroke-width: 0.75;' /> <line x1='553.77' y1='275.97' x2='570.46' y2='275.97' style='stroke-width: 0.75;' /> -<polygon points='545.43,283.04 578.80,283.04 578.80,278.99 545.43,278.99 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='587.14,284.59 620.51,284.59 620.51,281.02 587.14,281.02 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='587.14' y1='282.38' x2='620.51' y2='282.38' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='603.83' y1='288.19' x2='603.83' y2='284.59' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='603.83' y1='279.95' x2='603.83' y2='281.02' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='595.49' y1='288.19' x2='612.17' y2='288.19' style='stroke-width: 0.75;' /> +<polygon points='545.43,282.16 578.80,282.16 578.80,276.63 545.43,276.63 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='587.14,283.70 620.51,283.70 620.51,280.38 587.14,280.38 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='587.14' y1='281.36' x2='620.51' y2='281.36' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='603.83' y1='285.50' x2='603.83' y2='283.70' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='603.83' y1='279.95' x2='603.83' y2='280.38' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='595.49' y1='285.50' x2='612.17' y2='285.50' style='stroke-width: 0.75;' /> <line x1='595.49' y1='279.95' x2='612.17' y2='279.95' style='stroke-width: 0.75;' /> -<polygon points='587.14,284.59 620.51,284.59 620.51,281.02 587.14,281.02 ' style='stroke-width: 0.75; fill: none;' /> -<polygon points='628.86,292.70 662.23,292.70 662.23,264.27 628.86,264.27 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> -<line x1='628.86' y1='282.03' x2='662.23' y2='282.03' style='stroke-width: 2.25; stroke-linecap: butt;' /> -<line x1='645.54' y1='297.26' x2='645.54' y2='292.70' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> -<line x1='645.54' y1='258.13' x2='645.54' y2='264.27' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<polygon points='587.14,283.70 620.51,283.70 620.51,280.38 587.14,280.38 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='628.86,288.92 662.23,288.92 662.23,261.71 628.86,261.71 ' style='stroke-width: 0.75; stroke: none; fill: #D3D3D3;' /> +<line x1='628.86' y1='272.96' x2='662.23' y2='272.96' style='stroke-width: 2.25; stroke-linecap: butt;' /> +<line x1='645.54' y1='297.26' x2='645.54' y2='288.92' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> +<line x1='645.54' y1='258.13' x2='645.54' y2='261.71' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' /> <line x1='637.20' y1='297.26' x2='653.89' y2='297.26' style='stroke-width: 0.75;' /> <line x1='637.20' y1='258.13' x2='653.89' y2='258.13' style='stroke-width: 0.75;' /> -<polygon points='628.86,292.70 662.23,292.70 662.23,264.27 628.86,264.27 ' style='stroke-width: 0.75; fill: none;' /> +<polygon points='628.86,288.92 662.23,288.92 662.23,261.71 628.86,261.71 ' style='stroke-width: 0.75; fill: none;' /> </g> <g clip-path='url(#cpMC4wMHw3MjAuMDB8MC4wMHw1NzYuMDA=)'> <line x1='103.26' y1='502.56' x2='645.54' y2='502.56' style='stroke-width: 0.75;' /> @@ -194,8 +189,8 @@ <circle cx='104.62' cy='437.18' r='2.70' style='stroke-width: 0.75; stroke: #61D04F;' /> <circle cx='104.62' cy='451.58' r='2.70' style='stroke-width: 0.75; stroke: #DF536B;' /> <circle cx='104.62' cy='465.98' r='2.70' style='stroke-width: 0.75;' /> -<text x='126.22' y='441.31' style='font-size: 12.00px; font-family: sans;' textLength='105.39px' lengthAdjust='spacingAndGlyphs'>Starting parameters</text> -<text x='126.22' y='455.71' style='font-size: 12.00px; font-family: sans;' textLength='62.03px' lengthAdjust='spacingAndGlyphs'>Original run</text> +<text x='126.22' y='441.31' style='font-size: 12.00px; font-family: sans;' textLength='68.03px' lengthAdjust='spacingAndGlyphs'>Original start</text> +<text x='126.22' y='455.71' style='font-size: 12.00px; font-family: sans;' textLength='80.04px' lengthAdjust='spacingAndGlyphs'>Original results</text> <text x='126.22' y='470.11' style='font-size: 12.00px; font-family: sans;' textLength='75.36px' lengthAdjust='spacingAndGlyphs'>Multistart runs</text> </g> </svg> diff --git a/tests/testthat/_snaps/multistart/parplot-for-sfo-fit.svg b/tests/testthat/_snaps/multistart/parplot-for-sfo-fit.svg index a47a585a..c8d4970b 100644 --- a/tests/testthat/_snaps/multistart/parplot-for-sfo-fit.svg +++ b/tests/testthat/_snaps/multistart/parplot-for-sfo-fit.svg @@ -103,8 +103,8 @@ <circle cx='104.62' cy='437.18' r='2.70' style='stroke-width: 0.75; stroke: #61D04F;' /> <circle cx='104.62' cy='451.58' r='2.70' style='stroke-width: 0.75; stroke: #DF536B;' /> <circle cx='104.62' cy='465.98' r='2.70' style='stroke-width: 0.75;' /> -<text x='126.22' y='441.31' style='font-size: 12.00px; font-family: sans;' textLength='105.39px' lengthAdjust='spacingAndGlyphs'>Starting parameters</text> -<text x='126.22' y='455.71' style='font-size: 12.00px; font-family: sans;' textLength='62.03px' lengthAdjust='spacingAndGlyphs'>Original run</text> +<text x='126.22' y='441.31' style='font-size: 12.00px; font-family: sans;' textLength='68.03px' lengthAdjust='spacingAndGlyphs'>Original start</text> +<text x='126.22' y='455.71' style='font-size: 12.00px; font-family: sans;' textLength='80.04px' lengthAdjust='spacingAndGlyphs'>Original results</text> <text x='126.22' y='470.11' style='font-size: 12.00px; font-family: sans;' textLength='75.36px' lengthAdjust='spacingAndGlyphs'>Multistart runs</text> </g> </svg> diff --git a/tests/testthat/test_multistart.R b/tests/testthat/test_multistart.R index 91ef71f0..dda0ea23 100644 --- a/tests/testthat/test_multistart.R +++ b/tests/testthat/test_multistart.R @@ -50,7 +50,7 @@ test_that("multistart works for saem.mmkin models", { llhist_dfop_sfo <- function() llhist(saem_dfop_sfo_m_multi) parplot_dfop_sfo <- function() parplot(saem_dfop_sfo_m_multi, - ylim = c(0.5, 2)) + ylim = c(0.5, 2), llquant = 0.5) vdiffr::expect_doppelganger("llhist for dfop sfo fit", llhist_dfop_sfo) vdiffr::expect_doppelganger("parplot for dfop sfo fit", parplot_dfop_sfo) diff --git a/vignettes/FOCUS_D.html b/vignettes/FOCUS_D.html index 9f768b06..b8a63a7b 100644 --- a/vignettes/FOCUS_D.html +++ b/vignettes/FOCUS_D.html @@ -299,8 +299,8 @@ pre code { border-radius: 4px; } -.tabset-dropdown > .nav-tabs > li.active:before { - content: ""; +.tabset-dropdown > .nav-tabs > li.active:before, .tabset-dropdown > .nav-tabs.nav-tabs-open:before { + content: "\e259"; font-family: 'Glyphicons Halflings'; display: inline-block; padding: 10px; @@ -308,16 +308,9 @@ pre code { } .tabset-dropdown > .nav-tabs.nav-tabs-open > li.active:before { - content: ""; - border: none; -} - -.tabset-dropdown > .nav-tabs.nav-tabs-open:before { - content: ""; + content: "\e258"; font-family: 'Glyphicons Halflings'; - display: inline-block; - padding: 10px; - border-right: 1px solid #ddd; + border: none; } .tabset-dropdown > .nav-tabs > li.active { @@ -364,7 +357,7 @@ pre code { <h1 class="title toc-ignore">Example evaluation of FOCUS Example Dataset D</h1> <h4 class="author">Johannes Ranke</h4> -<h4 class="date">Last change 31 January 2019 (rebuilt 2022-07-08)</h4> +<h4 class="date">Last change 31 January 2019 (rebuilt 2022-12-06)</h4> </div> @@ -438,10 +431,10 @@ print(FOCUS_2006_D)</code></pre> <p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAABgAAAAMACAMAAADVCFz+AAAC/VBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWnp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////JahlgAAAACXBIWXMAAB2HAAAdhwGP5fFlAAAgAElEQVR4nO3deZwU9Z3/8fdwDPeNB+AJUQE1AVRQAolrdFFRVzEkmpgs8QCDoIggrEFDjGv8SYzRRGIIuooa4iZeBInHz+gaE+9FXaP5GYIKIWxAjggoINOPX3VX9UzPTB9TUF3f7vq8nv90dVX1zPc7PdRr6KuUAgCYJNcDAAC4EXkA1qwAAJTBmqiP11EH4NkaAQDKoObZiA/YUQfgHnXuDwCIXGfdE/EBO/oAnBfxVwQAeM4jAABgEwEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIQPmsmD6wY+/Rt213PQ4AyIsA7Ik3L774DwU3LuigjMP/EuOIAKDFCMCe+I10V6Ft99VIn5sxaV+p/8Y4xwQALUQA9kSRAGzordaLvcttY6XL4hwTALQQAdgTRQJwrTQjs7BpX7VZG9+QAFSLLS9E+uV2PLcr7E3MBWDnlroW7LV1feFNHzUsFwnA4Wqzxl/6tjS/xcMDYMYkPRPll7tei8LexFAAvKP18g0XdFPrvqfcWx+BTd8/qU9t7yHjHw+uH61zUlvnHFjzZPrK7yce0qnX8Okrg22jdFrqpdGt1fEzF2ZWTfOf5NW7eb7Ze9LIYPFVaWw0MwCQJGfrl1F+uem6KexNbAVg2YHBMfuz6/x1y3oEK3SC/2pNLwBbP+9d9QLw8YRgU+tr/J29ADzazl9V+0CqeACWSpcHi5900EHRzABAkpQKwKoh/Rs5tfhDPASgGC8A+0mf+vqkIa2kIZ+kV73jHc/7Xzz3ilPaSJdkdvICMEEDLpm3JpU6Q2o7ZsbUUd4h/qLMtlEa0lODr1h4tfdluns7/HX5bdK1y5fvyPPN5kk/yi4fppqt0UwBQIKUCsCDauLvRXcnAMV4AVCrG9JLy7pLd6YXvL/xp2Wa+lo39cvsdLT6aEbmUf5HpKPeTC882U9aml5Ip2By+j8KW4bK/7EVfg5gthru2uOlVdFMAUCClArAA/rnFTl6EoA9kA7ANf7i76QDPvYuh6pz8D7d86QP0pdHS8dlVmz/lHqu9rc9VaNh6UsvAMf4Tx4sla4IvmSBAEyRlmWXT5XeimYKABKkdADOyr26FwHYA97Rulv2LVneMflB7+KV598IVkySMq/V9ALwaGbFw9Lc7C0/p5r0cwajgv8JpFKrpCnBlywQgPOlp7LLZ0kvRTMFAAlCAEqLMgBnZ5fvk67M2bRzWe+GAPiv3pwpPZ3dPE16LJUJwP/6K9aXDMBk/yYZY6XXo5kCgAQhAKVFGYD6g/7z0rmZhY+fumXK2CNq08+vBAFo5z/K8+XGT77cm0oHoHNw89IBuNL/L0bGP0kro5kCgAQhAKVFGYCbs8trpFO8i+3XdfNf6HncsPoA+E8Gp05uHIA7UukA9A5uXjoA86SfZJcHSZujmQKABPmcjhtfzEj1y73aTmcU3f1QzQk7AlsBmJVdfkX6YipVl36l5xn/vuTNj9OP8gQB2M/fY7z02PIG6XcGhwnAIw3frK6T+kYzAwBJ0r/pyzz31FfCjsBWAL6YXV4sXZpKLZGODF6g2SwA3ooXm3yBMAFYKZ0QLL4pjYlkAgAS5QxdcX8xMzQ892pXLSy6+2maF3YEtgLQI/tQzDjp56nU5VL2IyAmNA3AL6TvZW9505QZ6ScGwgQgNVDtg5cc3STdEs0MACQJzwGUFun7AL7jL74g7eW1YKL0ir9iXb+mAdjSR3sH/zt4uUbnpC9DBeBq6cbMwo5D1eb9aGYAIEkIQGmRBqD1D9JLj/eSfuxdzg+O46kXhnjbVqSX6gOQulPaL/Pqzdf6Zj4aqEAAbi/w3dZ3UeeXvcu6qdLEaCYAIFEIQGlRBmB/adCFU45pJR2e/gCfVT2kE364aO7pNTpMGvfoR7kB2DVCanfmNd8a6+19VWZN8wA8I31m0cL8L/HxAtJx6t23jpQOXhfNBAAkCgEoLcoAzP9m8GT5CP8ULQ8Fp+3tfveL6YvVuQFIbRoT7Fw7039rQPMAbMy8ijTfp4F65rXxbz6Qz4EAkEfpABx+Q45OBGAPZB6vefLsPrX7jFmU/VTVVRcP69Jt2Nz1qdSNfdoesq5RAFKpR796UIe9R01ZGVxtHoDUf322S+1Bawp8w+UTB7TvNeJmPgkUQD6lArC0yas8azYV3Z0AFFPkAXsAiF2pAOyYNyvX7F8U/3IEoBgCAKCSnK1fRfnlpusHYW9CAADAiYVHrIzyyy0d1PTNqyURAAAwigDsuSafG5d2XVm+EQBEiQDsuZsvbObXZflGABAlQwHY8PTThV6wCQAGGQoAACAXAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhVDQE45eWwfrjoHovuuNn1CNy4dYHrEbjxf2z+mi/8oesRuHHLb0IfCEs7peIDcJ8AAGVxX8QH7KgDsPqEo8LaV206WtRK7V0PwYUOqnE9BCfaqq3rIThh9te8TegjYQucsjriA3bUAdgNN2m66yE4MULPux6CCyt1kOshOHGDZrkeghPD9IrrIbjw/3SI6yG0CAFwhgCYQgBMIQAtRgBMIQC2EICKRgCcIQCmEABTCECLEQBTCIAtBKCiEQBnCIApBMAUAtBiBMAUAmALAahoBMAZAmAKATCFALQYATCFANhCACoaAXCGAJhCAEwhAC1GAEwhALYQgIpGAJwhAKYQAFMIQIsRAFMIgC0EoKIRAGcIgCkEwBQC0GIEwBQCYAsBqGgEwBkCYAoBMIUAtBgBMIUA2EIAKhoBcIYAmEIATCEALfbMPve7HoITMwescz0EFz46YoLrITjxf/d50PUQnJh2yAbXQ3Bhy+ALXQ+hRSogAAAAFwgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAAOyhNe/ucj2E3eI8ACumD+zYe/Rt212PIx5PfOOwrj2O/ddnG9YYmv+rtXqy/oqNee964IyDOhx57h8a1piY97MTRnQ94OQ5HzSsSfq0d3bX2pyrzadbqT8A1wFY0EEZh//F8UDisHGcAl9bH6wyNP9tg9QQABvzXvV5f5Y1k7J/HlqY94Zzg9/yXndkVyV+2kuUG4Dm063YH4DjANxXI31uxqR9pf4b3Y4kBjtP9A4FJ8y6bLT3izCmLrPK0vwnqyEANua9boDU96IrxniT/Y6/xsS8T5fanDr9gsHeb/tSf03ip71lcG4Amk+3cn8AbgOwobdaL/Yut42VLnM6kjgskGoz/yKe8H4RfpResDT/pWoIgJF5e7ObtMO7fK6zWr2WXmFi3j+RBvzRu9z1A2mvTek1iZ/2b9N/09UHoPl0K/gH4DYA10ozMgub9lWbtSV2rnpHSL/0lx6rUc9PUqbm//d91Lo+ADbm/ZR0kv8fvdulq9OXJuY9VHrGXzpbWpa+TPa0bz2tb+bRnfqJNZ9uBf8A3AbgcLVZ4y99W5rvdCjl9yfp09nlY6QXU6bmf5r6nV8fABvzHie94i99OOSIL6cvTcy7nToFz3jcIX03fZnsaR8fPONRf2BvPt0K/gE4DcB70shg8VVprMuhxOBh6ZLs8kTpIVPzn6+ax2dkA2Bj3ts6alDjNSbmXddW3YPFxdK0VOKnvWC25+CGADSfbiX/AJwGYKl0ebD4SQcd5HIoMVjUvv1N2eWvSb+3NP+3O+rSVH0AbMz7OWlS4zU25j1YesNfmuz/uWth2ic3BKD5dCv5B+A0APOC50LTDlPNVpdjidWO/VSz3tD8dwzTwG0NAbAxb/8BkJcX3rjw9WCNjXn/VBq+Ib2wrLV6pV/zYmHaOQFoPt1K/gE4DcDs+mdFM4+krXI5llhdL52dMjT/WWr7cqohADbmPUf62e8Ozzw8fOxbmTU25p2aJvW6bP4Np0s9nk9ftzDtnAA0n24l/wCcBmBK8CqBtFOlt1yOJU4/rlG79MsCrcz/6Va6LpUTABvz9o6Dk9sGTxC2vyu9xsa8U6mFwawHrc5ctTDtnAA0n24l/wCcBuB86ans8lnSSy7HEp8V3u9Aq/9MLxmZ/8YDdGz6Ra/1AbAx74mZv/1/vXr9b707vOOKlJV5p5YdmH3D+1czb3i3MO2cADSfbiX/AJwGYLL0WHZ5rPR6sX2TYsucdlLPX2WWjcz/XHV6J31ZHwAb877AOwJe5L8g8jLppJSVeS9urfZzXtr8/tJR0qfTBbAw7ZwANJ9uJf8AnAbgSunB7PI/SSsdDiUuD+zvHRbO/Jt/xcb875EWZBbqA2Bj3lOl/bf4ix/um3nk18S8369Vh+AQ53XvSykb084JQPPpVvIPwPWrgH6SXR4kbXY5llhs/pJ3+B9W/99BE/Pf3DX7yufcVwElf96pf5MuyC6Pz7zvw8S8vXv5e8Hi9v5qtdLGtBu/CqjJdCv5B+A0AI9Is4LFuk7q63Iosdh4pNTzrrr66ybmv1qNTTAy7/RHP83LLl8l3Wlk3if6b3LP+Ib0sI1p5wSg+XQr+QfgNAArpROCxTelMS6HEoePRnuT/FvOChPzzxMAE/NOPeu/DzbjosyDACbmfbT0bnZ5lrTIxrRzAtB8upX8A3D7WUAD1T74dNSbpFucDiUGM6Wv7Gy0xsL8N8/OGi6dM3t2+vVPFuad2t5dx2aXh0rp58EtzPuszF/9Pu+w+HTKxLRzApBnuhX8A3AbgKulGzMLOw5Vm/edDqX8tvbQgG2NV5maf8NzAEbmfYF0v790vzQ0fWlh3guk0cGHwb3SVt3Sp8AyMO3cADSfbgX/ANwGYH0XdX7Zu6ybKk10OpIY3JnzmHDA1PxzAmBj3u+1V/dMAR7sqVaZ98RamPf2IdLXMyeDfOZg6fb0goFp5wag+XQr+Afg+Ixg3kGx49S7bx0pHbzO7UjK75vSkcc3+HN6naX55wTAyLznSzpq8iXDlf04eBPzfquj1PPMWROP86Y9zl+V/GnnBiDPdCv3B+D6nMDz2vhPDQ6sqPdHl8WYxk+GLs+sNDT/3AAYmfcNwUdBdJ6ffemXhXm/OjL4FW93VfZzzxI/7UYByDPdiv0BuA5AavnEAe17jbi5oj4hrzwOyRcAQ/NvFAAj835j6mGdOw+aurJhjYV51917et/aroMnv9OwKunTbhyAPNOt1B+A8wAAANwgAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEZFHYC1404EAJTBuLWlj8GhRB2A+wQAKIv7Ij5gRx2Ae/SFJ0KapC+GvQkAmPMF3RPxATv6AJwX9iY3aXrEgwCA5DmPAACATYkMwA8IAACUlMgAPH/gkogHAQDJk8gAAABKIwAAYBQBAACjCAAAGEUAAMCoZAZgW8RjAIAESmQAftPmPyIeBAAkTyIDwDuBAaA0AgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMSGQDOBwAApSUyAJwPAABKS2QAAAClEQAAMIoAAIBRBAAAjCIAAGBUMgPA+QAAoKREBoDzAQBAaYkMAO8EBoDSCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAIxKZAA4HwAAlJbIAHA+AAAoLZEBAACURgAAwCgCAABGEQAAMIoAAIBRyQwA5wMAgJISGQDOBwAApSUyALwTGABKIwAAYBQBAACjCAAAGEUAAMAoAgAARhkKwNvS3AjGAwAJkcgA5D8fgJMAvHnxxX+I/ZsCQEskMgD5zwfgJAC/ke6K/ZsCQEskMgD5EQAAyEUAyosAAFXhT5f+1fUQ9tSfLl0T9iYEoJmdW+pasNfW9YU3fdSwTACAqnClfuh6CHtqpm4JexOLAfhgiHRlvj28o/XyDRd0U+u+p9xbH4FN3z+pT23vIeMfD64frXNSW+ccWPNk+srvJx7Sqdfw6SuDbaN0Wuql0a3V8TMXZlZNk+/daGYAoFyu0PfL8WXfeKKR98vxPbKm66awNzEYgA3DpBl59/ACsOzA4Jj92XX+umU9ghU6YXtmhReArZ/3rnoB+HhCsKn1Nf7OXgAebeevqn0gRQCAqlGeALxYo0Y6bynDN8kiAIG85wMIArDxaBV6n5gXgP2kT3190pBW0pBP0qve8Y7n/S+ee8UpbaRLMjt5AZigAZfMW5NKnSG1HTNj6ijvfr0os22UhvTU4CsWXu19me7eDn9dfpt07fLlO0LPAECsyhOAh7XXiTlq9bcyfJMsAuDLfz4APwCbh0vTCt3OO5K3uiG9tKy7dGd6wfsbf9qu9MJr3dQvs9PR6qMZmUf5H5GOejO98GQ/aWl6IZ2Cyen/KGwZKv8Hy3MAQFUoVwDOyL26LwEIK9p3Av/jWOmyQrdLB+Aaf/F30gEfe5dD1dl/5Mf7SemD9KX3H4jjMiu2f0o9V/vbnqrRsPSlF4Bj/CcPlkpXBF+SAACVjwBEo6ID8OFI6ZsFb+cdrbttDJZPlR70Ll55/o1gxSRpbfrSC8CjmRUP57yq6HOqST9nMCr4n0AqtUqaEnxJAgBUPgIQjUoOwIzR6cdoCt7OO1qfnV2+r/ErhXYu690QAP/FtjOlp7Obp0mPpTIB+F9/xXoCAFQTAhCNSg5AW9XUqtULhW73m5yD/vPSuZmFj5+6ZcrYI2rTz+AHAWjnP8rz5cZP79+bSgegc3BzAgBUlSv03Q3Ru0+n5F7dW2+X4ZtkTSYAGQUDoJoF10tDPylwO+9ofXN2eY10inex/bpu/gs9jxtWHwD/yeDUyY0DcEcqHYDewc0JAFBVTlUCXBx21rYCULMgtWOwCr7jzztaz8ouvyJ9MZWqS7/S84x/X/Lmx+lHeYIA7OfvMV56bHmD9DuDCQBQpf5FHXpEr7Pa5l5tpW5l+CZZ7f2DThi2ApD+8TwjdSnwoR+/yRz0fYulS1OpJdKRq/wVzQLgrXixyRcgAECV4jmAaFRAAIqfD2CCND7/7byjdY/NwfI46eep1OVS9iMgJjQNwC+k72VvedOUGeknBggAUKUIQDQqIADFzwewrmf2hZxNpd8H8B1/8QVpL68FE6VX/BXr+jUNwJY+2jv438HLNTonfUkAgCpFAKJRAQHIr/7D4H4m9c/7YRHpALT+QXrp8V7Sj73L+cFxPPXCEG/bivRSfQBSd0r7vZ5eeK1v5qOBCgTg9miGD6CMCEA0Kj8AdSOlb+Xbwzta7y8NunDKMa2kw9Mf4LOqh3TCDxfNPb1Gh0njHv0oNwC7RkjtzrzmW2O9va/KrGkegGekzyxauLn59wJQSQhANCo/AKnX26j2j3n28AIw/5vB66lGZB7vST3Uwb/a/e4X0xercwOQ2jQm2Ll2pv/WgOYB2Jh5FSmfBgpUuHIFoGP/HK0JQFjlOCHMTOn4PHtkHq958uw+tfuMWbQrWLfq4mFdug2buz6VurFP20PWNQpAKvXoVw/qsPeoKSuDq80DkPqvz3apPSj0eXoAxKs8Afjv1o1fqN8j74PPESEAe4YH7AGrynRCmLUrGtlY+ha7jwAEdrOyBACwarZudT2EPTVLPwp7k0QGIP/5AFpwOwIAGPXetwuf5rtK7MYUEhmA/O8ELo0AALDEbACafJRb2nUEAIAlZgNw84XN/JoAALDEbADy2fD007xgE4AZBAAAjCIAAGAUAQAAoxIZgPznAwAA5EpkAPKfDwAAkCuRAQAAlEYAAMAoAgAARhEAADCKAACAUckMQDnPugMACZHIAOzu+QAAwJJEBoB3AgNAaQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGJTIAnA8AAEpLZAA4HwAAlJbIAAAASiMAAGAUAQAAowgAABhFAADAqGQGgPMBAEBJiQwA5wMAgNISGQDeCQwApREAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYlcgAcD4AACgtkQHgfAAAUFoiAwAAKI0AAIBRBAAAjCIAAGAUAQAAo5IZAM4HAAAlJTIAnA8AAEpLZAB4JzAAlEYAAMAoAgAARhEAADCKAACAUQQAAIwiAABgVCIDwPkAAKC0RAaA8wEAQGmJDAAAoDQCAABGEQAAMIoAAIBRBAAAjEpmADgfAACUlMgAcD4AACgtkQHgncAAUBoBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUYkMAOcDAIDSEhkAzgcAAKUlMgAAgNIIAAAYRQAAwCgCAABGEQAAMCqZAeB8AABQUiIDwPkAAKC0RAaAdwIDQGkEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAAVWPNu7si/GqJDADnAwBQcd6W5u7p19jZXWujGEsgkQHgfAAAKk4UAVgiAgAAVSeCAGwZTAAAoPrseQB+O1oEIHLvPltqj38s+aTcgwCQbHsYgFtP6ysRgOiNrHmvxB6zdW+5BwEg2VocgHXb8609XiIA5XCE3iixx2TdVu5BAEi2+gB8MES6Mu8uo3Raav5+Us9R6T85l52+T7tDx/9PsG3BbM/BBKC0kOcDCALwwaE9Ghm0qX4PAgBgD2UDsGGYNCP/Ll4ApgV/6F9eN8VfaJP7osaTCUBJYc8HEATgBTXx3/V7EAAAeygIwMajVfCdSqPUTTVfvPX7x3nHny+oy+QF1+wr9fu4YQcCUFrYdwLXB2DYhhxHEgAA0fEDsHm4NK3QLqOkDk94l7vO8ApwyJ+9pXV7S8sbdiAApe12AIbnrh1CAABEJxOAfxwrXVZwFy8AN2YWnvMC8Exm6dvS4oYdCEBpBABAxUkH4MOR0jcL7zIqe3xfJx3qr/q5dFfDDgSgNAIAoOJ4AZiRfifX5MK7jFIPf2GLNM5feogAhEQAAFQcLwBtVVOrVi8U3GWUevsLXgCCAx8BCHuTsAHoq2NP9AxX1xNzdNGI+uX9dG3YQQBALi8AqllwvTS04AcLEICmYghA56av/8xjUthBAECut9PH/9SOwdIPC+1CAJqK4XwAg/WzJzw/0sAncgzQ7fXLZ/AQEIA94wVginfxjNTlrwV2IQBNxXA+AJ4DAFB22XcCT5DGF9iFADQV32cBEQAA5ZMNwLqe0qP5dyEATREAAElQ/2FwP5P65//AMgLQFAEAkAT1AagbKX0r7y4EoKk4A9BnVo69CQCA6DScD+D1Nqr9Y75dCEBT8QXgf5q+9PP/1e9BAADsoZwTwsyUjs+3CwFoKr7zAaTuv6GRXzXsQQAAJE8iAxD2fABH6n9K7DFZ88MOAgAqXCIDEPadwDee9lGJPR7/7FthBwEAFY4AAIBRBAAA4ndy808cuy72QRAAAIjfzRc28+vYB0EAAMAoAgAARhEAADAqkQEIez4AALAokQEIez4AALAokQEAAJRGAADAKAIAAEYRAAAwigAAgFHJDEDI8wEAgEWJDEDY8wEAgEWJDADvBAaA0ggAABhFAADAKAIAAEYRAAAwigAAgFEEAACMSmQAOB8AAJSWyABwPgAAKC2RAQAAlEYAAMAoAgAARhEAADCKAACAUckMAOcDAICSEhkAzgcAAKUlMgC8ExgASiMAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwKpEB4HwAAFBaIgPA+QAAoLREBgAAUBoBAACjCAAAGEUAAMAoAgAARiUzAJwPAABKSmQAOB8AAJSWyADwTmAAKI0AAIBRBAAAjCIAAGAUAQAAowgAABhVDQH40oaQrtMFf7HjT++4HkGc/uh6AHEyNdl3TP0erwh7UCuTL1VBAAAAZVHpAXj7wB5htZNq7GCyiWVqtsYm2yH0Ya08Dnw74gN21AHYDaaeA/hf7e16CDH6Zz3megjxuVtfdz2E+OxSK9dDiNFXdJ/rIZQLAYgXAUgsApBYBKCcCEBiEYCkIgAJQQDiRQASiwAkFgEoJwKQWAQgqQhAQhCAeBGAxCIAiUUAyokAJBYBSCoCkBAEIF4EILEIQGIRgHIiAIlFAJKKACQEAYgXAUgsApBYBKCcCEBiEYCkIgAJQQDiRQASiwAkFgEoJwKQWAQgqQhAQhCAeBGAxCIAiUUAyokAJBYBSCoCkBAEIF4EILEIQGIRgHIiAIlFAJKKACQEAYgXAUgsApBYBKCcHq6Z73oI8fl475GuhxCjybVRn8Gugv2u1bWuhxCjAYe6HkGMrmn1B9dDKJcKCEDq764HEKd/fOR6BDHatd71COJk6vd4y1bXI4hTcu/aSggAAMABAgAARhEAADCKAACAUQQAAIwiAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAKC0Ne/ucj2E6MUXgBXTB3bsPfq27S3aVGTnqhBqsk9847CuPY7912djG13UQjo2S5EAABtFSURBVM0249VaPRnHyMog1GR3PXDGQR2OPPcPsY0uYqEm++yEEV0POHnOB7GNLmrFjzo7u2tti3euGrEFYEEHZRz+lxZsKrJzVQgz2Y3jFPja+nhHGZVQd23GtkGq1gCEmuyqz/traiZV59+OYSa74dzg17jXHfEOMjIljjpLlBuAaj9EZcUVgPtqpM/NmLSv1H9jyU1Fdq4KYSa780Tv+HDCrMtGe79MY+riH+ueC3XX+iarWgMQarLrBkh9L7pijLfhO7GPNAKhJnu61ObU6RcM9n6dl8Y+0iiUOOpsGZwbgGo/RNWLKQAbeqv1Yu9y21jpslKbiuxcFUJNdoFUm/kH84T3y/SjmEcahVCz9S1VtQYg3GS9pUk7vMvnOqvVazGPNAKhJvsTacAfvctdP5D22hT3UCNQ4qjz2/RfaGtbuHMViSkA10ozMgub9lWbtSU2Fdm5KoSa7BHSL/1tj9Wo5ycxDjMioWab8fd91LpKAxBqsk9JJ/n/p7tdujrOYUYj1GSHSs/4286WlsU4yqgUO+rcelrfzOM9a1uyc3WJKQCHq80af+nb0vwSm4rsXBXCTPZP0qezG4+RXoxrjNEJdddmnKZ+51dpAEJNdpz0ir/mwyFHfDm2MUYm1GTbqVPwPMcd0ndjG2N0ih11jg+e3ljbkp2rSzwBeE8aGSy+Ko0tvqnIzlUh1GQfli7Jbp0oPRTTGKMTarYZ81Xz+IzqDECoyW7rqEGxji5ioSZb11bdgzWLpWlxjTE6RY86C2Z7Dm4IQLUfonLEE4Cl0uXB4icddFDxTUV2rgqhJruoffubslu/Jv0+pjFGJ9Rs097uqEtTVRqAUJN9TpoU6+giFu6eHSy94a+ZXJV/E5c+6pzcEIBqP0TliCcA83Ke3zxMNVuLbiqyc1UINdkcO/ZTTfW9EDTsbHcM08Bt1RqAUJP1Hwp5eeGNC1+PdZBRCXfP/lQaviF9fVlr9arCF8aUPurkBKDaD1E54gnA7PpnOjOPp60quqnIzlUh1GRzXC+dXf7RRS3sbGep7cupag1AqMnOkX72u8Mzjx0f+1asw4xGyHt2mtTrsvk3nC71eD7GUUal9FEnJwDVfojKEU8ApuS8MOBU6a2im4rsXBVCTbbBj2vUrgpfKxhytk+30nWpqg1AqMl6R8TJbYNnD9vfFes4IxH293hhMNdBq+MbY3RKH3VyAlDth6gc8QTgfOmp7PJZ0ktFNxXZuSqEmmzWCu/3qNV/xjK+aIWb7cYDdGz6pa5VGoBQk52Y+dv/16vX/9a7bzuuiHWgUQj5e7zswOw72r9afQ9ktuSokxOAaj9E5YgnAJOlx7LLY6XXi24qsnNVCDVZ35Y57aSev4ppgJEKN9tz1emd9NUqDUCoyV7gHQsv8l8aeZl0UnyjjEi4e3Zxa7Wf89Lm95eOkj5dhQUofdTJCUC1H6JyxBOAK6UHs8v/JK0suqnIzlUh1GQzHtjfO1ac+bd4hhexULO9R1qQuVqlAQg12anS/lv8FR/uW4UPFIea7Pu16hAcB73afSmmIUao9FEnJwDVfojKEdurgH6SXR4kbS66qcjOVSHUZD2bv+Qd/oc9lapOYWa7uWv2NdNVGoBQd+2/SRdk14yvwrd4hJqsd4d+L1ixvb9arYxpjNEpfdRp/Cqgqj5E5YgnAI9Is4LFuk7qW3xTkZ2rQqjJplIbj5R63lWVHwOXFma2q9XYhHiHuudC3bULpHnZrVdJd8YzxOiEmuyJOe9i/4b0cExjjE7po05OAKr9EJUjngCslE4IFt+UxhTfVGTnqhBqsqmPRnsL1fnoT0aY2VZ9AELdtc/mvCP2opzHDKpFqMkeLb2b3TpLWhTTGKNT+qiTE4BqP0TliOmzgAaqffDmkJukW0psKrJzVQg12ZnSV3bGOryIhZjt5tlZw6VzZs+uvlc9hblrt3fXsdmNQ6V34hpjZMJM9qycv/q9I+XTcY0xOiWPOjkBqPpDVIOYAnC1dGNmYcehavN+iU1Fdq4KYSa7tYcGbIt7gJEKdddmVelzAOEme4F0v7/tfmlojKOMSJjJLpBGBx8G90pbdavC82SVPOrkBqDaD1ENYgrA+i7q/LJ3WTdVmphZs+Ghhx76KP+m5muqS5jJ3pnzQHF1CnXXZlVrAEJN9r326p4pwIM91aoK3x0bZrLbh0hfz5wM8pmDpdsdjXhPFJmtLzcA1X6IahDXGcG8A13HqXffOlI6eF1mxfOSVufflGdNdQkx2W9KRx7f4M/uBr3bQt21gWoNQLjJzve2HTX5kuHKfnp8lQkz2bc6Sj3PnDXxOG+Xcc5GvCeKzDYjNwBVf4iqF9s5gee18Z/4Gxi8bzrnp9t0U7411aXlkx3T+GnR5Y4GvEdC3bW+qg1AuMneEHwUROf51fkqrzCTfXVk8Dvc7qoq/XC0IrNNaxSAqj9EZcUWgNTyiQPa9xpxc/aXI/en22RT3jXVpcWTPSQBAQh312ZUbwDCTfaNqYd17jxo6sqYxxiZMJOtu/f0vrVdB0+uvme7s4rMNtU0AFV/iArEFwAAQEUhAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEZFHYB/zJwIACiDOf+I+IAddQAWCwBQFosjPmBHHYB7NOKnIY3XiWFvAgDV7USND3uTEbon4gN29AE4L+xNbtL0iAcBABVuum4Ke5PzCAAAJAAB8BEAAOYQAB8BAGAOAfARAADmEADfQs2NeBAAUOHmamHYmyQyADuWbI14EABQ4bYu2RH2JokMAACgNAIAAEYRAAAwigAAgFEEAACMSmQA/jr17YgHAQAV7u2pfw17k0QG4BbNjHgQAFDhZuqWsDdJZAB4JzAAc3gnsI8AADCHAPgIAABzCICPAAAwhwD4CAAAcwiAjwAAMIcA+AgAAHMIgI/zAQAwh/MB+DgfAABzOB8AAKClCAAAGEUAAMAoAgAARhEAADAqkQHgfAAAzOF8AD7OBwDAHM4H4OOdwADM4Z3APgIAwBwC4CMAAMwhAD4CAMAcAuAjAADMIQA+AgDAHALgIwAAzCEAPs4HAMAczgfg43wAAMzhfAAAgJYiAABgFAEAAKMIAAAYRQAAwKhEBoDzAQAwh/MB+DgfAABzOB+Aj3cCAzCHdwL7CAAAcwiAjwAAMIcA+AgAAHMIgM9WANa8u8v1EAC4RwB8hQKw62fDOu/94h6PqLLs7K61rscAwD0C4CsUgLnyPL3HI2qRNy+++A/hN+2GJSIAAAhAVoHzAexsp46X3xPT0fI30l3hN4W3ZTABAJDifABZBc4H8LY0cc/H00IxBeC3o0UAAKQ4H0AJy6UbIvpSpcURgFtP65t+UIsAxGHjhMfi/pafXHZH3N8SxhCABju31LXgq2xdX3jTRw3L0QZg3fZ8a4+XCEBcluisuL/lX3Ro3N8SxtgJwA+Do2W+J4G9Q/LyDRd0U+u+p9xbH4FN3z+pT23vIeMfD64frXNSW+ccWPNk+srvJx7Sqdfw6SuDbaN0Wuql0a3V8TMXZlZNC77Zu82/V6NNn9xzxn61PYdeubL5fjlfef5+Us9R93rXlp2+T7tDx/9PsG3BbM/BBCAWD+lfwt5k3S/vz7VkS8jbv6NPhf2WQCgEIM0LwLIDg+2fXeevW9YjWKET/D+/vQBs/bx31QvAxxOCTa2v8Xf2DtOPtvNX1T6QanEAVn4muNLuB4WG7X3l7C0ur5viL7RZkrPDyQQgFrsRgNPUWNh3pxAAlJudAKxb/gvp0uXL8/0Z5gXA+yP7U1+fNKSVNOST9Kp3vON5/4vnXnFKG+mSzE5eACZowCXz1qRSZ0htx8yYOsr7R31RZtsoDempwVcsvNr7Mt29Hf66/Dbp2uXL8zwnk7Nps7d355Nmff1Q7+vcVmDYo9RNNV+89fvHeft8QV0mL7hmX6nfxw07EIB47EYARuqE8Q2O0YSQtycAKLdEBqDA+QAKPwfgBUCtMtuWdZfuTC94f+NPy7zD9rVu6pfZ6Wj10YzMo/yPSEe9mV54sp+0NL2QTsHk9H8UtgyV/xNtyXMA/+bV5i/eZd13W6nXpvw7e1+5wxPe5S4vOjrkz97Sur2l5Q07EIB47FYAnsu59h8EAOXF+QB8Bc4HUDwA1/iLv5MOSP99PVSdgydez5M+SF8eLR2XWbH9U+q52t/2VI2GpS+9w/Qx/pMHS6Urgi9ZKgDvt1f3jf6qGdLs/Dt7X/nGzMJz3hifySx9W1rcsAMBiAcBQKXjfAC+Au8ELhqAbsGhOHWq9KB38crzbwQrJgVHWC8Aj2ZWPKyGN5p9TjXp5wxGBf8TSKVWSVOCL1kqADc0fJ21HbR3/p1HZY/v65R9TcjPG31pAhAPAoBKxzuBfbsTgLOzy/dJV+Zs2rmsd0MA1mRWzcx5JnmalH55uHeY/l9/xfqWB+BM6bXsui9If8m78yj18Be2SOP8pYcIgAMEAJWOAPh2JwD1B/3npXMzCx8/dcuUsUfU1r/Q/mi18x/l+XLj13akX545Sp2Dm4cIwDHSP7LrLpR+l3fnUertL3gBCH4QBMCFh3TiipCG6T9zrt2os0Pe/ikCgDAIgG93AnBzdnmNdIp3sf26bv4LPY8bVh8A/8ng9CE3V/rdmvWH6TABGKhO9eu+Lf06784EoELcKgc6lR4XkEUAfLsTgFnZ5VekL6ZSdelXep7x70ve/Dj9KE8QgP38PcZLjy1vkH5n8G4FwPsfwIfZdZMKfUwpAagQd6hT/5Daq2/Otb3UJeTt988+/Ae0BAHw7U4AvphdXixdmvmQ5SNX+SuaBcBb0fSkArsVAK8xr2fXnSj9Ke/OBKBC8BwAKh0B8O1OAHpsDpbHST9PpS6Xsh8BMaFpAH4hfa/+W02ZkX5iYLcCcL10XfZWndUz/wcREYAKQQBQ6QiAr8D5AIq/D+A7/uIL0l5eCyZKr/gr1vVrGoAtfbR38L+Dl2t0TvpytwKwsl19dq6ULs+/MwGoEAQAlY7zAfgKnA+geABaZz6O5/Fe0o+9y/nBcTz1whBv24r0Un0AUndK+2UevHmtb+ajgQoE4PYCw6vfNFM66j3vsu57rdRrY/6dCUCFIACodJwPoKiiAdhfGnThlGNaSYenf4irekgn/HDR3NNrdJg07tGPcgOwa4TU7sxrvjXW2/uqzJrmAXhG+syihZubf6+cTZv6SF3HzrlwsNeYnxYYNAGoEAQACUQA0rwAzP9m8Mq7Ef7h9KEO/tXud7+YvlidG4DUpjHBzrUz/UfumwdgY+ZVpHk+DTR304ojgq/T/tZCgyYAFWK3AtClR4NOBAAVhwCkZR6UefLsPrX7jFm0K1i36uJhXboNm7s+lbqxT9tD1jUKQCr16FcP6rD3qCkrg6vNA5D6r892qT1oTd5v17Dpk7tP71fbfcis9woOmgBUiN0IwEVNXtU/L+TtCQDKzVIACivygD2QsRsBqNvQSIHPey2MAKDcCEAaAUApv9GX4/6Wq3RE3N8SxiQyAAXOB1AYAUApn/zkz7F/z5+/EPu3RBXjfAC+AucDKIwAAKh2nA/AV+CdwL4mH+WWdl3ZApDvm0WzMwA0wjuBfUUDcPOFzfy6bAHI982i2RkAGiEAvqIByGfD00/nf8EmAFQJAuALHQAAqHYEwEcAAJhDAHwEAIA5BMBHAACYQwB8Bc4HAADJxfkAfAXOBwAAycX5AAAALUUAAMAoAgAARhEAADCKAACAUYkMQOjzAQBAteN8AL7Q5wMAgGrH+QB8vBMYgDm8E9hHAACYQwB8BACAOQTARwAAmEMAfAQAgDkEwEcAAJhDAHwEAIA5BMDH+QAAmMP5AHycDwCAOZwPAADQUgQAAIwiAABgFAEAAKMIAAAYlcgAcD4AAOZwPgAf5wMAYA7nA/DxTmAA5vBOYB8BAGAOAfARAADmEAAfAQBgDgHwEQAA5hAAHwEAYA4B8BEAAOYQAB/nAwBgDucD8HE+AADmcD4AAEBLEQAAMIoAAIBRBAAAjCIAAGBUIgPA+QAAmMP5AHycDwCAOZwPwMc7gQGYwzuBfQQAgDkEwEcAAJhDAHwEAIA5BMBHAACYQwB8BACAOQkNwIHjQzpUnfYzrFdX1yNwqV+XvVwPwaW9uvRzPQSXuvZyPQKXOunQsMfKAys+AI8IAFAWj0R8wI46AHVL7w/rPB11uWHddL7rITh0pg52PQSXDtaZrofg0Pnq5noILh2l80IfLJfWRXzAjjoAu8H4cwD9tcL1EBxaqlNdD8GlU7XU9RAcWqH+rofg0m48BxA9AuAaATCMABhGAHwEwPUQHCIArofgEAFwPQQC4B4BMIwAGEYAfATA9RAcIgCuh+AQAXA9BALgHgEwjAAYRgB8BMD1EBwiAK6H4BABcD0EAuAeATCMABhGAHwEwPUQHCIArofgEAFwPQQC4B4BMIwAGEYAfATA9RAcIgCuh+AQAXA9BALgHgEwjAAYRgB8BMD1EBwiAK6H4BABcD0EAuAeATCMABhGAHwEwPUQHCIArofgEAFwPQQC4B4BMIwAGEYAfATA9RAcIgCuh+AQAXA9hIoIwI80y/UQXDpMq1wPwaEn9C+uh+DSv+gJ10NwaJUOcz0El2bpR66HUBEB+GDu+66H4NLj812PwKUd1/+36yG49N/X73A9BJfmP+56BC69P/cD10OoiAAAAFwgAABgFAEAAKMIAAAYRQAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAwigAAgFEEAACMIgAAYBQBAACjCAAAGEUAAMAoAgAARhEAADCKAACAUQQAAOKx5t1drofQmIMArJg+sGPv0bdtD7cpMYrM8YlvHNa1x7H/+mzsY4pPqXv41Vo9Ged4YlZk+rseOOOgDkee+4fYxxSfIrN/dsKIrgecPOeD2McUr53dtTbfencHvvgDsKCDMg7/S5hNiVF4jhvHKfC19S5GFodS9/C2QUpyAIpMf9Xn/U01kyrsT8ToFJ79hnOD3/xed7gYWHyWKG8AHB74Yg/AfTXS52ZM2lfqv7HlmxKj8Bx3nuj96z9h1mWjvV+EMXVuRlduJe/hyUpyAIpMf90Aqe9FV4zx9viOk7GVX5HZny61OXX6BYO9fwFLnYwtJlsG5w2AywNf3AHY0FutF3uX28ZKl7V4U2IUmeMCqTbz2/+E94vwIwdjK7+S9/BSJTkAxabvrZq0w7t8rrNaveZgbOVXZPY/kQb80bvc9QNpr00uBheP36b/uGseAKcHvrgDcK00I7OwaV+1WdvSTYlRZI5HSL/0lx6rUc9P4h5ZHErdw3/fR60THIAi039KOsn/X9/t0tWxjywORWY/VHrGXzpbWhb3wGJy62l9Mw/zNP+9d3rgizsAh6vNGn/p29L8lm5KjMJz/JP06ezyMdKL8Y4rHqXu4dPU7/wEB6DI9MdJr/hLHw454ssxjyseRWbfTp2CJz7ukL4b87jicnzwNEfzI7zTA1/MAXhPGhksviqNbeGmxCgyx4elS7LLE6WHYh1XPErdw/NV8/iM5AagyPS3ddQgByOKU5HZ17VV92BxsTQt3nHFZsFsz8F5AuD2wBdzAJZKlweLn3TQQS3clBhF5rioffubsstfk34f67jiUeIefrujLk0lOABFpv+cNMnBiOJU7M4fLL3hL01O7P/9fSfnCYDbA1/MAZiX8/zmYarZ2rJNidGyOe7YTzVJfCFo8dnvGKaB25IcgCLT9x/5eHnhjQtfdzCwWBS7838qDd+QXljWWr0S+wLAtHwBcHvgizkAs+uf6cw8KLaqZZsSo2VzvF46O64Rxan47Gep7cupJAegyPTnSD/73eGZh4iPfcvB0GJQ9M6fJvW6bP4Np0s9no99ZHHKFwC3B76YAzAl51n+U6W3WrYpMVo0xx/XqF0iXwlYdPZPt9J1qUQHoMj0vQPg5LbBk4Tt73IwtvIr/qu/MJj8oNVxjyte+QLg9sAXcwDOl57KLp8lvdSyTYnRgjmu8H4HWv1njGOKT7HZbzxAx6Zf+prgABSZ/sTM3/6/Xr3+t96933GFg8GVXdFf/WUHZt8E/9UkPvbZIF8A3B74Yg7AZOmx7PJY6fWWbUqMknPcMqed1PNXsQ4qNsVmf646vZO+THAAikz/Au/Qd5H/SsjLpJPiHlkcit35i1ur/ZyXNr+/dJT06UQXIF8A3B74Yg7AldKD2eV/kla2bFNilJrjA/t7R4Iz/xbrmOJTZPb3SAsyCwkOQJHpT5X23+IvfrhvMp//KjL792vVITjsefn7Urzjile+ALg98MX/KqCfZJcHSZtbtikxis9x85e8w/+wp5reKDEKz35z1+zrnxMcgCJ3/r9JF2SXxyfzTSBFZu/d598LFrf3V6uVsY4rXgVeBeTwwBdzAB6RZgWLdZ3Ut4WbEqPoHDceKfW8K6EfA5dWePar1dgEF8MrtyJ3/gJpXnb5KunOOIcVkyKzPzHnje/fkB6OdVzxyhcAtwe+mAOwUjohWHxTGtPCTYlRbI4fjfbWJPXRn4zCszcRgCJ3/rM5b4C9KOcRgQQpMvujpXezy7OkRbGOK175AuD2wBf3ZwENVPvgnR43Sbe0dFNiFJnjTOkrOx0MKUYFZ795dtZw6ZzZs5P5KqjCd/727jo2uzxUeifeccWj8OzPyvmr3ztCPh3vuGKVLwBuD3xxB+Bq6cbMwo5D1eb9lm5KjMJz3NpDA7Y5GVN8WnAPJ/g5gGLTv0C631+6Xxoa98BiUXj2C6TRwYfBvdJW3ZJ8PsC8AXB64Is7AOu7qPPL3mXdVGliZs2Ghx566KP8m5Kn8PTvzHkYOKmK3PlZSQ5Akem/117dMwV4sKdaJfPNsIVnv32I9PXMySCfOVi63eEYyy43AJVx4Iv9jGDega7j1LtvHSkdvC6z4nlJq/NvSqCC0/+mdOTxDf7sdphlUuTODyQ5AMWmP99bOmryJcOV/Wz45Ck8+7c6Sj3PnDXxOG/FOKdjLLfcAFTGgS/+cwLPa+M/0TcweM9zzjGg6aYkKjT9MY2fBl3udJBlU+TO9yU6AMWmf0PwURCd5yf2dWCFZ//qyODXvt1VifwQyHr5A+DywBd/AFLLJw5o32vEzdl7OvcY0GRTIhWY/iEmAlDszs9IdgCKTf+NqYd17jxo6kpHI4tD4dnX3Xt639qugycn8unvHAUC4PDA5yAAAIBKQAAAwCgCAABGEQAAMIoAAIBRBAAAjCIAAGAUAQAAowgAABhFAADAKAIAAEYRAAAw6v8D9326WQ+aolwAAAAASUVORK5CYII=" width="768" /></p> <p>A comprehensive report of the results is obtained using the <code>summary</code> method for <code>mkinfit</code> objects.</p> <pre class="r"><code>summary(fit)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.0 -## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:37 2022 -## Date of summary: Fri Jul 8 15:44:38 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:42 2022 +## Date of summary: Tue Dec 6 09:39:42 2022 ## ## Equations: ## d_parent/dt = - k_parent * parent @@ -449,7 +442,7 @@ print(FOCUS_2006_D)</code></pre> ## ## Model predictions using solution type analytical ## -## Fitted using 401 model solutions performed in 0.13 s +## Fitted using 401 model solutions performed in 0.158 s ## ## Error model: Constant variance ## @@ -492,11 +485,11 @@ print(FOCUS_2006_D)</code></pre> ## ## Parameter correlation: ## parent_0 log_k_parent log_k_m1 f_parent_qlogis sigma -## parent_0 1.000e+00 5.174e-01 -1.688e-01 -5.471e-01 -1.174e-06 -## log_k_parent 5.174e-01 1.000e+00 -3.263e-01 -5.426e-01 -8.492e-07 -## log_k_m1 -1.688e-01 -3.263e-01 1.000e+00 7.478e-01 8.220e-07 -## f_parent_qlogis -5.471e-01 -5.426e-01 7.478e-01 1.000e+00 1.307e-06 -## sigma -1.174e-06 -8.492e-07 8.220e-07 1.307e-06 1.000e+00 +## parent_0 1.000e+00 5.174e-01 -1.688e-01 -5.471e-01 -1.172e-06 +## log_k_parent 5.174e-01 1.000e+00 -3.263e-01 -5.426e-01 -8.483e-07 +## log_k_m1 -1.688e-01 -3.263e-01 1.000e+00 7.478e-01 8.205e-07 +## f_parent_qlogis -5.471e-01 -5.426e-01 7.478e-01 1.000e+00 1.305e-06 +## sigma -1.172e-06 -8.483e-07 8.205e-07 1.305e-06 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index da6c11fe..b8c9ba0c 100644 --- a/vignettes/FOCUS_L.html +++ b/vignettes/FOCUS_L.html @@ -1373,8 +1373,8 @@ pre code { border-radius: 4px; } -.tabset-dropdown > .nav-tabs > li.active:before { - content: ""; +.tabset-dropdown > .nav-tabs > li.active:before, .tabset-dropdown > .nav-tabs.nav-tabs-open:before { + content: "\e259"; font-family: 'Glyphicons Halflings'; display: inline-block; padding: 10px; @@ -1382,16 +1382,9 @@ pre code { } .tabset-dropdown > .nav-tabs.nav-tabs-open > li.active:before { - content: ""; - border: none; -} - -.tabset-dropdown > .nav-tabs.nav-tabs-open:before { - content: ""; + content: "\e258"; font-family: 'Glyphicons Halflings'; - display: inline-block; - padding: 10px; - border-right: 1px solid #ddd; + border: none; } .tabset-dropdown > .nav-tabs > li.active { @@ -1517,7 +1510,7 @@ div.tocify { <h1 class="title toc-ignore">Example evaluation of FOCUS Laboratory Data L1 to L3</h1> <h4 class="author">Johannes Ranke</h4> -<h4 class="date">Last change 18 May 2022 (rebuilt 2022-09-14)</h4> +<h4 class="date">Last change 18 May 2022 (rebuilt 2022-12-06)</h4> </div> @@ -1536,17 +1529,17 @@ FOCUS_2006_L1_mkin <- mkin_wide_to_long(FOCUS_2006_L1)</code></pre> <p>Since mkin version 0.9-32 (July 2014), we can use shorthand notation like <code>"SFO"</code> for parent only degradation models. The following two lines fit the model and produce the summary report of the model fit. This covers the numerical analysis given in the FOCUS report.</p> <pre class="r"><code>m.L1.SFO <- mkinfit("SFO", FOCUS_2006_L1_mkin, quiet = TRUE) summary(m.L1.SFO)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:35 2022 -## Date of summary: Wed Sep 14 22:28:35 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:45 2022 +## Date of summary: Tue Dec 6 09:39:45 2022 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## -## Fitted using 133 model solutions performed in 0.032 s +## Fitted using 133 model solutions performed in 0.033 s ## ## Error model: Constant variance ## @@ -1637,10 +1630,10 @@ summary(m.L1.SFO)</code></pre> <pre><code>## Warning in sqrt(1/diag(V)): NaNs produced</code></pre> <pre><code>## Warning in cov2cor(ans$covar): diag(.) had 0 or NA entries; non-finite result is ## doubtful</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:35 2022 -## Date of summary: Wed Sep 14 22:28:35 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:45 2022 +## Date of summary: Tue Dec 6 09:39:45 2022 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent @@ -1742,17 +1735,17 @@ plot(m.L2.FOMC, show_residuals = TRUE, main = "FOCUS L2 - FOMC")</code></pre> <p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAqAAAAJACAMAAABlpiR1AAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////isF19AAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO2deWAM5//Hn0hC5BARouJoE3GUVBH31ZZQZ6mjLdXqgRLfOuqIuHpQTVF1taSt4yvqiKPqiFYULV+qXz+Uto4Kpa2vEqGiiBzPb2Z2E8lmZzOzM0/2mdn364+Z2dlnPp+PzcvMzuzM8xAKAMcQVxcAgCMgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoLaMozksVd4tePFR/0aDFyba3kvN6FNsE/Y4B/F5XGEbBDnGwgZLczSpzUP8Y987of8OPGELMl/ceypGgFN3r2rKl20sK6PuDBWWHjANr27AEFtKWjMP8Oti10yxLcudbK88vyQ2gp6rqK15Yy8OAUFTfGS3mtyT0266DwvW1kXCqV3FyCoLYIxMUskLtEXCfHvM6mvPyGDxLfaElJhwNS+noR8aSvoY4T0XfHv18uSUoescQoKWpv4LtrVh5D31KQTBSXnKc30sQpaKL27AEFtEYzZZl1MIaT2aWF+uhYhKZSuI6TeBeHlVkIa2gh615u0F1+sIeRd68YFBL1EyARKb5cjHdWkEwR9gHxO6SFS0VMUtHB6dwGC2lLAmCaE7JYWdgqHZ0rbEPK19PLZho3vFRb0CiHNsoUXd1etOmjduICg3xFRNFqf1FKTThC0NxlB6TzylJcoaOH07gIEtUUw5vXPRH6+5kFqWlc+RDzSaDCpWKBd4UN8OCER49efzb3/fgFB76WlZVJ6059Eq0knCPoBaST4SOIlQQundxcgqC35Zy2LjgpnK9aVwunJkTRhL1mgXWFBD1SVtinfd3fe+4XO4gVyXiAkSUU6UdCDxPMmrUG+EwW1Se8uQFBb7htzhJCu1pVdCflBOIy3KtDO5jLTrdWDG3iLmy20vm8j6I1ehIyyLm/9l8BHxaQTBT1bh+y6RLxvi4LapHcXIKgt978Upt0/5tYk5C8aREItr+7dvZtrK6hI5n9eFc7DcywvCgv6f+HEc1beiymikJ2LSycK+gp5Z6Ow45QO8YXTuwsQ1JYCZy1RlqvnlO4l4pfBloRYriE1It7/0OmELBVfLCXkLXoyPt5yciQcnP+wbFxI0JU+pHL+wV9GUNt0oqBLSedxgv+SoIXTuwsQ1JYCxnxNSN1UYf5rXUJ2UPpvQlpeF15u8yCtKd1ISA+xUQ9C1tPThDwmnsXnNiN+2ZaNCwq6xYPU+V1tOlHQUySwpfDNVRK0cHp3AYLaUsAYOpCQcv3f7B9AyEAq2UceGjPneU9SSvDnRnVBmKlThP1ateviXo20XLQxoT0hz1i3FQQdKJ2ef/ZPTg1CxohLG9WkEwWlFYkH+d0iaOH07gIEtaWgMTc7Ws9gOt4UX15oa3nl/YH46ttgy6uK3wovfqtibRl5w7ptfP5vmH/sz1uKVJNOEvQpQmpQi6A26d0ECGpLQWNo7uoBj/g9MmC19aQkZ2GvhwIav/qr5dX12JbBwS1jr0svbs1uU61MeIflWXmbFhB0hUJBC6eTBJ1FyLN5gtqkdw8gKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGt0FvSNKAAU0OayawRtsfwwAMVTS+mY4noLerD4NgDQBhAU8AxHgh6bPXVDlp2mwJ3hRtCc16uPe6dD5Fl90wCjw42gCa0zhOnCpvqmAUaHG0Fbp4jT3PBf9M0DDA43gj50Xpp1+lrfPMDgcCNoyz3SrNYJffMAg8ONoAufuCNMlzZ0nwHOjUHOtynsSZPPz42g2S/XfGtez5r4CsoZxwOimRMxST4/N4JSun/K4OV39c0CNHO0Ifsc702Uf48fQTPnNik/DJdBeQOCWrnRsO+SsJkh2+y2Bi4DglqZOJjeCKAHQ+/pmwdoBIJaefQwpQF/08aH9M0DNAJBrYSdo/Thn2jHnfrmARqBoFa6rKW001d3Q37XNw/QCAS18tVDP9FXPn6ln75pgFYgaB6JVdrV9h90U980QCsQNJ+Mb0f01zcJ0A4ELcD2LvomAdqBoAX4MVLfJEA7ELQA6eX1TQK0w1DQo/+1LhhGUOqPcyTeYCjowF7WBeMIWuekvlmAZooImrtp3Ji1OQ63yc0usiqz0E/Ylod3DShohxR9swDN2Ap6vW2r2XMfb/KX7AZnyY4aHnXfF5ZujorwCZ+eS+kD68YF/UZXRPnWXyasrvZRNAnsd422IITckDYxjqCDluubBWjGVtAhMeIzD7HPym5wlpQdsnmcx5uU9g6as2UUWS0I2qTn+jsLvKclj/b4WBA06Pm9s7xG0Gt9Ol+2PD9hHEEnT9c3C9CMJOjd9DyuBZwRZ3/4/y9/lc1N5mdJH2E6PuAm7bNEWKgzQRC0QS69FSz+bYeECIKKj/X0aG3IQ/zHr+mbBWhGErRrUB7lPSzzUoH5q1oV3uAs2SBMTxDxtrSMo595jxMEjaX0EDmQlpaWSC7SanHCO+NaGFLQrd30zQI0Y3uIr/CHOL1eTvbhnLNkv9hA0HR/A4/QLqGioHMpXUcsHKfVZlMuBM28mu7gXfuCHn1UbRbAGFtBx/QXzsezh74qu4FlD/oL+U+6d4xwKtVCFHQepbvJFWuDanOo6wW9OCnMg5AyERPPyTSwL+jVYFVZAHtsBb3du25s3COd/5bd4Cx5RphO8r2+k5yl9E4Vq6BXyiwVVk/twIegR3yrD5+fuHLBv8KDjtlvYV/Q3LL/qEkD2FP0Qv2++Jm7HWxwlvjEbI8rNYn+5t1v35fN/bqmSYLS2NIzksd7zC8g6Mv1D1sumJa8oI93uW1ZyOofbb+FTP+gEWfUpAHsUf1L0lnyRbfytWcKZ+pr6/g12/rvoGkWQXPnRPrWE8/r8wXdExZg2RGXvKDlNuQt7ZP5eV1G0Mcd/d8ELsAJQQ+rzVHygjaJyVt6q7n9FjKCDlypJg1gjzkF3eDRZemBk6e+X9nbc4P9FjKCTnxXTRrAHnMKSrc9IV308mifLNNARtBFMXZXA5ehWtCs3zLV5nDJddD0n1JSjl+zWZmR359ZxHd2t9r8lMo0gDGmvmE551ebfkJ+yO/QzOsDu1scbuxEGsAQkwq6rW/35TQhmPjMkOnv0/9Tu6svh6hKA5hjTkGTSNOu3jH+03dM8V5qv4WMoLk+d9TkAcwxp6CNhlO6mLwnLMU1st9CRlAalqomD2COOQX13U7pVemulh1+9lvICdp2r5o8gDnmFLTW+5R+T1YIS/Pr2m8hJ+iAVWryAOaYU9D3fEa/XbVp6M60L4Mn228hJ+iEeDV5AHPMKWjWlNBKw+4NJoR0k7k9SU7Q+a+ryQOYY05BrZxcKdsdrZygG592Ig9gx0+lw5lTYap8flc98iEn6CEM1skZ51PZ4+AuYN4E/bOKvnmAweFN0Jwyqm81AGaGN0Hpg+f1TQSMDXeCtt6nbyJgbLgT9Nk1+iYCxoY7QcfO1jcRMDbcCTp3tL6JgLHhTtCkvvomAsaGO0EPyjwGCtwT7gT9s7K+iYCx4U5QWuGK3DvADeFP0LZ79M0EDA1/gg5foG8mYGj4E3QRelkG9+FP0L2t9c0EDA1/gl4LlHmUHrgj/AlKq1zUNxUwMhwK2lGuwzHghnAo6JhZ+qYCRoZDQZcO0jcVMDIcCnooSt9UwMhwKOgtv6KD5QJ3xVUDeTkQlIZjqA+Qh6sG8nIk6FObVOUCZsZVA3k5EnTSO2pyAVPjqoG8HAn6+TNqcgFT45SgF08U07j4gbwcCfpjPYU1AfPjhKBfPUAI7eDwprjiB/JyJGhmWdmxnoG7oV7QRK+hiYRO9Uhw0Lj4gbwcCUrrKS0KmB71gtYbQ9OEF+MjHbUudiAvh4I+87nCooDpUS+ob7IkaLJM5/N52B3I68TQPLw/dLDt25MUFgVMj3pBG78pCfpOg+K3uX4mp/CKq58mWPFZ7GC7jRhvDlhRL+hy7+kHyJWlpec6ap3Y8ji93J2QgHkyDRwe4k+HKywKmB4nzuIXBAvfLstMdHTf+3zS9jLtVHXxjvFeMuNrOxQ02++mwqqA2XHmOuitQ+t2X3XYOCyO0stEHC92osqBvCxEyfZtD9wMNr8kBW8UTojILWEpOcB+C8eCDpIZQBG4HeoF7Z6Hg8a9umbSrIBdwlKcExfqKZ31hsKqgNlRL+hLIr1CSv3LQeOTIZGz9swIXbYnzlumP1rHgiZ3UlgVMDvOHuJvdXI4XsyZ4RWkK/X118k0cCzoxVCFVQGz4/R30L3EcSdf984d3H5A7m7Q4gSlwZcUlgVMjtOCriirqX+FYgTttVZLcGAe1Au6SuKtKo9ryluMoB8O0xQdmAb1gvpI+LY4qSlvMYIekxnFG7gbHD7VKZJb8U998wGDwqmg9GkMlwREVApasQCa8hYn6Dz0EgpEVAq6ZMmSef5VR8+dXKf6Nk15ixP0WB1N4YFZUH+IH/aY+MRQdmdtu7jiBMWXUCChXtCqlqeMNmv7sac4QWnv1ZriA5PghKCWe5AX1NCUt1hB5w/VFB+YBPWCDgncKky3BbI9xNMfa2uKD0yCekEzHicVIiuQ9rc05S1W0NxKf2hKAMyBM9dBd8ePnv2dxrzFCkr74NljwO+FeuFL7hB9MwJDolLQz1LoZ3loylu8oMdraUoAzIFKQf2fo+Xz0JS3eEFzK+NLKOD4EE/7rtI3JTAiTgp6aodtrzYqUSDowsHaUgAzoF7Qix3H0Y2epKJM18kKUSDoiQhNGYApUC9or9DNtOGT5zo6euy4eBQImhuCMRGBekErfEIvkRS6upKmvAoEpYMwcjxQL2j5tXSZzx261V9TXiWC7mipKQUwA+oF7dRq1yM96a1uTTTlVSJoVmX5x5aBm6Be0GOVSbkjtHbprcVs4PxAXvkMe09RacDEOHGZ6fbha5QmOR4NTtNAXvnsDF+yE+MpuDdshqHRNpBXHnvDfZ55LOI/CgsEpoTNMDTaBvKycinkm7FT6fbKjrsiBeaGzTA02gbysjJ7GD1cM5cO+khhhcCMsBmGRttAXlaGf0xp7f/SuWMUVgjMCJthaDQO5GVh8jRKp46lE2YorBCYEUbD0GgbyMvCf2vMahLsPTu0uDMyYGYYDUMjM5DXbpLP+8VmzK3tM2pVRW8F4zEB88JmGJo8Nst2cqtgD7ot6kBMt1a963+jqD5gTtgMQ5O/0S65dxQIOlbcyV6oNG2yolTAnKgVNOPrdReEo++fP+/t5qDxpoEWSIeBA+23UCDoiIXitNWrYxVWCMyISkF/qS4c3of/WEM8AXLQeLsviWohQB5u0cJ+CwWCftJbnC6viCc/3BmVgvYITjy+plLVlmt3HnLYccOpqGanqcZD/D8RM+/Ru1O9DyqsEJgRlYJWEr8Xvk9+K7b5vVj/xRoFpRd6BjUq33fqAIUVAjOiUlAijnu0UdETnnurd/mfNkEpTT9ynd6seFZJU2BO1Aoq/jC0WdkjyNefrahVUIm41xU3BaaDoaCUrhl9Wu4tFYJeroD7mdwXpoI6QIWgdMhb2vMBg6JWUG8fHx9vIg2VpCmvGkFPV9LW1SMwMCoFnVIATXnVCEp7L9KUCxgYjvtmus8PYVn6ZgeGwRCC0nYY1stdMYagKTXv6JseGAVjCEp7v6NvemAUDCLoxYroZMQ9MYigdPrT+uYHBsEogmbWlnu6CZgaowhKv45AJzjuiGEEpT3Qk5g7YhxBUysUfxcqMB3GEZROw3mSG2IgQe82Ri9N7oeBBKVnQ/5P3yIA/xhJUJoUcUPfKgD3GEpQGtNP3yoA9xhL0LuNP9a3DMA7xhIUX0PdDoMJStfVvKxvIYBvjCYofesRR6PbALPBTlAdxkmyyxstCz1Cl7tp7MhV2c4GA7zDSFB9xkmyS+7L3Qs8ofT3Yy3nzO/46CVnowHOYSOoPuMkyZDd+/mc/BcxQ8WOdN/s5XQ0wDdsBNVlnCRZbrcZlb8c/Kc4vVMuw/lwgGfYCKrLOEnyXI96zfqt8563ZR6GJ0JMChtBdRknyQG3ujz1j2XpAcnMjHL/aAkH+IWNoLqMk+SIrMHNLOMzxPb5bvZ7KTEvaIoG+IXRWbwe4yQ5JPfNmr+K8xs1vJu2LBuMe5nNCrProHbHSbryaYIVn8Uq4xVhYVWxb/DYvt/NmrlrBPagZoWVoNnnLH2B3Cn0y+SJoXl4f6gunh22hMzOxXdQs8NG0KypZUnZCeKZ9gqZ7bQe4kUutu7+P5zFmxw2gs7yGrthjNfLlK2gNOvNGoG4Dmpu2AhaO06YrCJfMhaU0q2+jcXH5fE8nWlhI6jfdnE6MOwOa0HpqYplRyyIbojf4s0KG0GbviFO/woZwVxQmjupXN0E3M1kWtgIuoCM3Ckcerd7vhjLWlBK/x4RukrJ0MvAiDC6zDS9HBGH39oaStgLSunB5o1TdAwHOILVddDMVKmvr6xvlth/X1dBKd1SKxoPK5kSwz3yIcO9j6oMOK5vSMADZhGU0oz3q3Tfr3dQ4GrMIyild5dEtNmaU3w7YCDMJCil2WubhsdfYREZuAhzCSpweHCF/t/iqpNpMJ2glN5YGBk+9RSz8KBEMaGgAkfHhjad/wfLDKCEMKegwrfRlJeCm79/hm0SwB6zCipwL2V4aGTctxiH1tCYWFCBnINTmgT1+eRCSeQCTDC3oCJ/JQ6oHP5qIr6RGhPzCyry86K+lWoNSjiBq/iGwz0EFcj96ZNBtQM7T/0S9zYbCrcRVOLKlmldKlXt+eYXeMbOKLiXoBK/rZ/cvVpguxGLv7N9bB/whxsKKpG2a/6QVoGhHWIW7ryAH0Y5xl0FtXBx54KYDtXKPvL0hE92n8eDTTzi3oJayDi6/r1XHqtepmb0kBmJ+y7gyj5PQNB8pjxc5dHXn29drUy1Vs++MW/9/t8yXV0RgKD5ZFQp+/TgCM9ESu9d2L9mzqg+rWqUDnmky6C4+Wv3/pLm6urcFwhqpVuIuMMcV7rgAf7yj8krZo589rGHg0uHNuj4/KgZCZu+/fkyvgKUJBDUip9lrO8ySXbfzfzj6I7EuXGv9mr7cIhnYHjTzv1HTJm7fPO3P164WZJFuiEQ1IrXPmkWPLP4ptd/PZS8auHbo198qu0j1f09K4Q3bv/0SyOnzEpYm7zvWGr6PScr2D9t1LK7Tm5rXiColfLviNMsr51qN8xKO3t418Zl894ZN/SZLq0bhAV5lQkOf7RVdL+XYmJnfJCw5ouU74+l/nW9uDDZL0W8/WGvmj87U7uZgaBWXvM9Q2lOdKAOoe5cTT26f2fSskXxk8YMfbZndLMG4SGBxCfowYio1tG9+g0ZFjslfl7CqqTklO8Pn069YtnlLmgv7j2XP4pfDQoDQfNoW6pOU7+AH5jFv51+/szh/Smbkj75OP6d2JFDn+/XObp5VK3wSkHexDeoik/1B4OCH+ri32PosNjYGfEfJCSsTtqYkvKfw4dTUy+mp7vtsR+C5rP7he7vu+h+vFvpl6pGNp8587nAmhMSPo6Pnxz7xtChz/XrHR3dKioqPLxaUFAZQvyDgmqEhzeMimofLeyH+704dGhMbGzszHhhb5zwSVJS0paUFGGXfPhIamrq+fT09Buu+ac4z/G5b20q+mue8QaTNSnV24vTXV4nZFtkpKdfSE09evjwNykpXyQl/TshYVF8fHxcrLA3Hjq4X79+PaKjo5tFRTUKDw9/KCgoKFAcZkWYBz0orIiMiopqLrwf3VNo2E8cJECUO3ZSvIg0qsUqQfGkDSkiwl5b4EyqyOV0kRLovzp3uH9wueCIVNv1xhtM1qRUbyBasMD3v3oGzRXt+k3w7Lhg3PeifJtFD0UhPxLVfFe0NFYa1eJ50dw+osPiXlugVrhIZVHxID9pUCESKL2oIL0RXjvKQntpm+ge/Sy8Yh0lY1SshcnxVj7MG+BF3Ntb2JGSx77Do0v3Szq6qny47T/BgIPJmpOwQdXHzegQ2Ub1VYQS5Ia0O70m7VpTTx+28I1FsS1W6ZZaNZxn9XKGVdTY0XkDvAzpl0fn6DzaRHkHRUVNp7942vpixMFkTUnHTcdmTV5/u8p5VxfiIrx3SLPAuTbrS3Yw2R/y/8t4faAmnhuwubbwbejeqO6ursNV+KyWZr62O66SHUw2I/9LR8R3auK5A4tCug8M6+W2N/k3eFi8FyLB87LtehcNJtvioJp4bsHVLYnyp/Cm5z9lq727qJvvs7brXTWYLAQFhUmpUb6Sb2yRe3BLdjDZ+0BQYEPuuSN3iq511S9JEBQoAoICroGggGtcJWj7gKAi+BEPZjAMzRKDll2q6B/XWbyU9vGqs6B30ouyusM5ZlT/llnoKa8wC33Gi1nocwPfYRZ6V4Sdv66TKL4pS2dB7bGlB7vYYex6Y5o7hlnoLC9moWnMR8xCn6zLLLQ8EFQOCGoLBFUPBLUFgqoGgtoCQZUDQeWAoLZAUPVAUFsgqGogqC0QVDkQVA4IaotZBf11CbvYb7J7MHHfF8xC545jFpqu/55Z6BvTmYWWpwQEBcB5ICjgGggKuAaCAq6BoIBrICjgGggKuAaCAq6BoIBrICjgGggKuAaCAq5hL+iGpoFPHGUSeZPUUdRgBpFTvpRmLEq3hGZR+sLm/nVmi+Pl6V92Xmh2n7gszAXd5jF8fWe/iyxCzwlZIrBX/8A5zaTbjViUbg3NoPTp5I1tE72msig7PzSzT1we5oI+0ZnS29UnsQgd04FFVPr7R+2IZJH+peeH1r/0zHIjhenYstn6l30/NKNP3BGsBU0nS4XpsDAWsTsPZRGVJrdp4yNaxKD0vNAMSk8lYjf4G8g5/cvOD83qE3cEa0F/IgeE6XwPFoO0136ysd+jCQwC0wjRIjalS6EZlH73rDgo2Jiyd/QvOz80w09cFtaC7iInhWkiuap/6JzSwfM3DyZz9I9ssYhN6VJoVqWv8hrP6hMXQzP8xGVhLWgKOSVMVxJHo385SeZacXioF8sxGEZOsohN6VJoNqVfeYG8lMWmbEtohp+4LKwFPU7EZ2QWlGGWYBM5q39QySI2pUfcfx5J59K3h4RtpmzKtoa2wOQTl4W1oNc8VgrT12syCP3XYXEI4S3EdgAJHZAsYlO6FJpF6ds9Y6SuthmUnRea4ScuC/PLTO17U5oVHssgcgr5XJi+VoNBaMtujknp1m8PupeeFfqCdUn3svNDM/zEZWEuaLLn2/sHBLF4fD27ecj07SNLrWcQ2iIok9Kl0AxK/4ZMWCFyR/+y80Mz/MRlYf9T5/pmgR3Y/NR5e3TdgFY7WES2flFkUboltP6lL7EMDCsefvUu+35odp+4LLhZBHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVco7Ogm+MBUMDsf1wjaIsXYgEonko/MhT0+hnZLnZbHHQiHnA/GrARNLHlcXq5OyEB82QaQFCgCDaCzidtL9NOVRfvGO+10n4LCAoUwUbQsDhKL5PvhKWJjey3gKBAEWwEDd5I6QlyS1hKDii4/tvwPDyXqIkH3BY2gvbqmkmzAnYJS3HNC67PSc3D71M18YDbwkbQkyGRs/bMCF22J857jf0W/hAUKIHRWfyZ4RWk3srrr5NpAEGBIhgJSum9cwe3H5AfQQKCAkUwE7QY9BN0z4w3t+fqFQzwhtEFvd2z/uS3mrZN0yca4A6jCzp2QBaluWOe0yca4A6jCxpyQZzeKqf0nhdgMAwu6D1vyzyMxRCKgAMMLiit+Ls4vV0uQ59wgDeMLuioQeIdfRP76RMNcIfRBc3o1GjmrHZN/9InGuAOowtKaXLc+A2y90UDo2N8QYGpgaCAayAo4BoICrgGggKuKVlBTw0basVb7nFPAApSsoL+mZBHmUV6xAOmB4d4wDUQFHANBAVcA0EB10BQwDUQFHANBAVco1LQRQXQlBeCAkWoFLRiATTlhaBAETjEA65xWtC0aZry2hM057amkMCMqBc09/MJ4wS6VtCUt6igp3v4+dZeiU5sQCHUC/o2qeMd2qKKn0zf3gopIuiZyotu537fcIamqMB0qBc0bDT9tCu93WybprxFBH1xlji9VP5vtZG2x41PwkNzpkW9oGW20iNVKd3c3FFrmn3ujjS/c9n++0UErX1SmrXbq7AeKxmdGuOxYzOjXtCaH9CbHlfod/4OGmdNLUvKTsgWllbInP0XFfSUNGu3R2E9Vka+LO494/qo2woYBvWCjgxaQhvEnO1X10HjWV5jN4zxepmqEPSF2eL0f0E3FNZjpeIf4vR2uVvqNgNGQb2gf/fuQ/d5E+8kB41rxwmTVeRLW0EzUvIo+4nNJqcrL75Lf2g8XWE5VtB5mNlx8jpo+s7fHDX22y5OB4bdsRH0h+g8POfYbvNLV9+AiOVqLzOFSHWg+0XTwuaXpKZviNO/QkYoP8QLZN1Uk8PCuOfEDmxHD1C/JTAE6gXtnoeDxgvIyJ13Kd3u+WKsCkGd4XavepPfbtLumj7RAHeoF/QlkV4hpf7lqPX0cuSsMNsaShgLSuneGW8m4/cn0+LsIf5Wp6YOm2em3hVnWd/IDHmIm0WAIpz+DrqXXNGSF4ICRTgt6Iqymo6rEBQoQr2gqyTeqvK4prwQFChCvaA+Er4tTmrKC0GBInBHPeAaCAq4Bg/NAa5RKeiSJUvm+VcdPXdyneo63xEoeD8AAAsSSURBVLAMgD3UH+KHPSZegs/u/JqmvBAUKEK9oFU3SLPNoZryQlCgCCcEtXTevaCGprwQFChCvaBDArcK022BJXCIT+pav/sXmtIAo6Ne0IzHSYXICqS9tocslAj6UpPNP214dLimPMDgOHMddHf86NnfacyrQNBv6okPht6KOKAxlXHIdnUBHFKyF+r3BOXhUTYoKPgipYOCZOc+vtLcx6eYdiaZB5b2KNN4o+vr4G2uTtDPUuhneThj6PV0K/7z09PFPhoy02Xnr34gzWfEFNPOHPOjld/6LXNXvVmuroO3eaQ6Qf2fo+XzULihfRQc4hc9L816L9OUyCi8InX6czEow9WFuIyvxr6yoOi/nuPf4m9U+zSX5iwMc49H3uv8Is3afOviOlzFvb6Rcz4bWP2I7XonBT21Q+NjakrO4n9uHdb5wcdPa0tkFCIs/87Hd7u4Dlcxp0uWMF33sO1t8OoFvdhxHN3oSSoe01SPouuguad3/Kopi4HoP1+cXg1Kd3UhLqLZ3uXPdY37X11bH9UL2it0M2345LmOjh47Lh78kmTDTyErc+gvLSe7ug5X8WCTDmu2jw9p8rXNevWCVviEXiIpdHUlTfVAUFuOtPN/4IFFbtuR5IPSI0QpXj/brFcvaPm1dJnPHbrVUe92xQNBi/KPTFeVbsGDdcUvN+/5HbZZr17QTq12PdKT3urWRFM9EBQU4qHhVYbFtWjaVvsh/lhlUu4IrV16q6Z6ICgoRPSm04tn7rhb5bzNeicuM90+fI3SpDPa6oGgoBBf1D5P6b1RRU69nboOevGE5nogKCjMopAeL4T3TLNd7YSgXz1ACO2wwHFz1X3Ug5Mfv7vdbU/iBa58ufJ40bXqBU30GppI6FSPBAeNneij3t3Jja0yfFKrJg77BXZH1AtabwxNE16Mj3TQ2Ik+6ik9OmtKUpbCaszHsqbXhen7rVxdB2+oF9Q3WRI02c9BY7k+6u9TRNCcf9UYPz26vtv8tGlLu6/Eae6Dp1xdCGeoF7Txm5Kg7zRw0Fiuj/r7FBF0SWvxrqVF2q6uGpiHzkuzTrbXAd0d9YIu955+gFxZWnqug8ZO9FHfepc4zQ3/RWE9+vN34hzb39lKkGaWh2iK3Czh7jhxFr8gmBBSZqKj7kHl+qg//kw/K162frt8DzLO0zeoVL2rrko/90nx+/ea+ujNvDDOXAe9dWjd7mL+kDJ91Kcl5eGz2GaLlnukWS3tl1idY45XIqV/PBTmovT03jP1Z302oLq2mxhNiFpBM75ed0E4Ev/5895uDpur7qN+4RPihdOlDV21BwmRHm++VqrILd0lxs5xry5yj6cH1KBS0F+qC4f34T/WEKYexW8zUb4b+yKCZr9c8615PWu67Cuol2UU2wrvu6oAYBeVgvYITjy+plLVlmt3HlLwn53IXzOxcx10/7RRy+8WWpPyfOuBexTWpxWf1dLMDz8g8IVKQSuJO5j3icLfO9QJWoTRdZfv+zQiTlkurTSVfnhILOW+T1XyiUpByTphslHpE57aBN0fIcpyvcb/KcymjXNlwhasf7pUbIkkA4pRK6jY9+JmpYLukh/iVYGgE9+RZpPeVphNI9eiA31qauuVF+gPU0EdoEDQ4R9Js7ljtGcDhoVjQee9Is36O7ptCpgdtYJ6+/j4eBNpqCRNeRUIevUB8X/Dqqo3NCUCxkaloFMKoCmvkrP4/z7SqH+DRvhtxa3huG8mSrMOrf7Bne8xB5wLCgAEBVwDQQHXQFDANRAUcA0EBVxTsoJek7+jHgB7lKygx/rJPpMEgD1wiAdcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDVMBb1+RvaJIggKFMFI0MSWx+nl7oQEzJNpAEGBItgIOp+0vUw7VV28Y7zXSvstIChQBBtBw+IovUzEXtcnNrLfAoICRbARNHgjpSeI2INocoD9FhAUKIKNoL26ZtKsAHHcjrjm9ltAUKAINoKeDImctWdG6LI9cd5r7LeAoEARjM7izwyvQETqryu0eo8HseKxSFU84K4wuw5679zB7QfOyb7d4qDKeMA9YXqh3sEoHxAUKIKpoA76qIegQBEQFHCNywSdmFCEd58cyIz2zzML3as7s9ADn2AXultvZqH79yr6x3WWqiwFdTDKx6dDi/K4/8PM8I5gFrpyBWah63owC/1w0APMQoeXtfPXdZKRfzMUVCVberCLHSZ/QUErDIchyfJiFprGfMQs9Mm6zELLA0HlgKC2QFD1QFBbIKhqIKgtEFQ5EFQOCGoLBFUPBLUFgqoGgtoCQZUDQeWAoLaYVdAdT7OLXesis9ALxzMLne3HLDQd+Qmz0L9GMgstTwkImn2dXew0dqHv3GIXm2HZGXfZxWZYtiwlICgAzgNBAddAUMA1EBRwDQQFXANBAddAUMA1EBRwDQQFXANBAddAUMA1EBRwDXtBNzQNfOIok8ibpP7KBjOInPKlNGNRuiU0i9IXNvevMzuLsig7LzS7T1wW5oJu8xi+vrMfk7vi5oQsEdirf+CcZuPEGYvSraEZlD6dvLFtotdUFmXnh2b2icvDXNAnOlN6u/okFqFjOrCISn//qB2RLNK/9PzQ+peeWW6kMB1bNlv/su+HZvSJO4K1oOlkqTAdFsYiduehLKLS5DZtfESLGJSeF5pB6alkpzDdQM7pX3Z+aFafuCNYC/oTOSBM53tkMohd+8nGfo8mMAhMI0SL2JQuhWZQ+t2z4p3KY8re0b/s/NAMP3FZWAu6i5wUponkqv6hc0oHz988mMzRP7LFIjalS6FZlb7KazyrT1wMzfATl4W1oClSV3grSbr+oTPXpgrTF8vJDnvnPJJFbEqXQrMp/coL5KUsNmVbQjP8xGVhLehx8r0wXVCGWYJN5Kz+QSWL2JRuOcRL6Fz69pCwzZRN2dbQFph84rKwFvSahzgm3es1GYT+63CuMN1CLusfWrKITelSaBalb/eMuSPOGZSdF5rhJy4L88tM7XtTmhUeyyByCvlcmL5Wg0Foy26OSenWbw+6l54V+oJ1Sfey80Mz/MRlYS5osufb+wcEsehfIbt5yPTtI0utZxDaIiiT0qXQDEr/hkxYIXJH/7LzQzP8xGVh/1Pn+maBHdj81Hl7dN2AVjtYRLZ+UWRRuiW0/qUvsY5UdVn/su+HZveJy4KbRQDXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRBUb/paO4ohM6j/Z64uxvhAUL05uGHDhuqthclJ+lyKq4sxPhCUBZHPuboC0wBBWWAVtLxwiK/2YUvf8I8vdS//oNi35ooo3/rLXFubwYCgLCgoqPe0b7p7VPswuZlPBl3gPS15tMfHLq7OUEBQFhQUtB+lp8goSreRE7eCpwtrh4S4tjhjAUFZUFDQ9ynNJqtFS48dIgfS0tISCZOBIU0KBGVBQUHniIJukARdZ70AddzF5RkJCMoCGUF3kysuLsx4QFAWyAh6pYw4jObUkh+R1cBAUBbICEpjS89IHu8x38XVGQoIygI5QXPnRPrWW+Li4owFBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVc8/9HUPmXE6PYSgAAAABJRU5ErkJggg==" /><!-- --></p> <pre class="r"><code>summary(m.L2.FOMC, data = FALSE)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:35 2022 -## Date of summary: Wed Sep 14 22:28:35 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:45 2022 +## Date of summary: Tue Dec 6 09:39:45 2022 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted using 239 model solutions performed in 0.049 s +## Fitted using 239 model solutions performed in 0.048 s ## ## Error model: Constant variance ## @@ -1820,10 +1813,10 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, main = "FOCUS L2 - DFOP")</code></pre> <p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAqAAAAJACAMAAABlpiR1AAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////isF19AAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO2dB3zU5BvH39IWuhilUCi0SEtZUgFZgsJf9kZQQREBFzKKIrtMQYoKgshSi7JkKKMgu0IREJChCAoiqBRZIkIpQqGDjvef5O5Ke3e5Xi55e2/uft/PxySXe/M8T8PXyyWXvC+hAHAMcXYBANgCggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayAo4BoICrgGggKugaCAayCoBYOJiX3Cq4T+9fzr9l2Ta3gvd1HzIJ/wAb+Iy6MJiRfn8YQMF2Ypbz8WHBDV+4e8ODMIict78fNTVUo2ejfDSrpxUirf8Nev204+WlpbLPzpkyz+aG6BoBbkd+TeEONip1TxravtDa88P6Lmgp4vZ2w53RQnv6CJXtJ7je5bphtnylbhgM3ko01veB5gvAO4AoJaIDgSHSdxlfYnJODZCT0DCHlJfKsFIWX7TO7pSchmc0GfJKTn8i/e9CXFjhrj5Be0BvFbuPtZQt63TCcI+krcpxPrEVI721ZyId3AuLj3GxLSgPUe4AkIaoHgyDbjYiIhNX4X5r9XJySR0rWEPHxReLmVkPpmgmZ4k9bii68Iede4cT5BrxIyltK0UqSdZTpB0FXCLEP4cF5qK7kx3d1ypJi1bwquCgS1IJ8jjQjZIy3sEg7PlDYnZKf08vn6De4XFPQ6IU2ED0CasWrVYePG+QTdT8hqYVaHVLdMZxSUHiKko63kpnSPE3JZ07+XbyCoBYIjby4WOX3Tg1QzrqxKPJJpECmXr13BQ3wEIZFj1p/LffB+PkHvJydnUnongLS1TGcSlPqIyeSTG9PdK098crT8czkHglqQd56y8IRwfmJcKRyAjycLn5L52hUU9FBlaZsyPfeY3i9wFi+Q04+QdZbp8gSNIJ735ZOL6QYvXjy7CSFdtPtb+QeCWvDAkeOEdDau7EzID8Jh/PF87cwuM939ckBdb3GzBcb3zQT9rwchbxmXt74h8LFhOU/QagUFNUv+4Cw+4Iz2fzO/QFALHnwNTH5wlK1GyL80kFQyvLqfkZFrLqhI5vevCQIZj8AFBf1J+ID8wPRikihaR8NynqC+xkO8THKjoJWfOaflH8s9ENSCfOcpDQ3XyyndR8ijlDYjxHAN6VHifY/GErJEfLGEkKn0zIwZhpMj4XB8xbBxAUFX+JAKeQd/q4IeNT9JMktu+v/BzYCgFuRzZCchtZKE+Z+1CEmg9AtCmt0SXm7zIE9QuoGQbmKjboSsp78T8qR4Fp/bhPhnGzbOL+gWD1JT5tzbKGhmJ0IW20oOQYGBfI7QvoSUemHKCyUJ6Usl+0jVEbNf9CTFBGP+CxN8nTxJ+FgNvSV+qJJmCzcsak3Ic8ZtBUH7Sifki+/lVCFkhLi0wTKdIOhriz+f2kDQMdtWcggKDOR35E4745lJuzviy4stDK+8PxRffRdkeFXuO+HFhRBjy6j/jNvOyPvV8spB01KUZbq8nzqDD9hMDkGBgfyO0Nwv+zzi/0ifL43XN3MW9KhassFrfxpe3YppFhTULOaW9OLurOahJSLaLMsybZpP0OWFCupTdYDpZhGZ5BAUAP6AoIBrICjgGggKuAaCAq6BoIBrICjgGggKuAaCOkBS9/5xhbcCWgBBHWDMcdrC2TW4CxDUAdJyLz9dSJN9PexZJUfKX8oKcmEgqCOs6/5PIS3sEfSzmhX6pokLNTy9vDbQsaEho+iFpjXnURp9VqNC9Q8EdYBtI3MLrrhv/O/BKzsEPVb1n9udpwoLudKjeDsbpd2q+c2gb7Ib3P49WuOCdQwEVcTI6kfp+lqvdOnZV3o5PzxyZM7BPr0nif/R6RHhY6RXVLJRerPD15Q2+k5aNBd04yxKl4hxLj8ivvw+gdJnN46d89fDaf3/LuI/i2MgqDI2dEh7KMH0Yl+963d7Tz9Y+iwV//um5u2M5ivEJfGtHoY3l/WjSRF7pUVJ0E4PSRj7V7rZWPCX7g9rU+EV4Vi/oXPXjNuD2u76cbwz/jJOgaDKyK3+Ure8F5OrtmxZp8vBJykV/xs7kdKPXxaXqCio4c3/Ktx/7x3DouVRf1Wk9DTSLx/eT+4gHOuvrq8nud/z1pYOo92p9xBbQFCFzCz+54Pl9ylNv32wqyCo8N8Y4dAe109coqKghjdp9531zxsWJUHbV5YQu8fL6feU4VQrK5vSFU9vPETph4OFlzs/SG5wcfKHRf+ncQkEVUjvoAefbSdqXc/ovMMkaELt1MwWy/IENbxJv2rZytjO4jtoG/FUK/ds1twO6be7fbig1e3bbYQz+Nyn0i+2yV03oQj/Jp6BoMo4GNJt1YNXcdUqv0lNgtLYyIhROXmCGt6kd/2WG9uZCzqieEBAwEsZ5EJWTETYiKyctx4KHZpF6ZrFlE6p3eJa0f1NXANBFZHTcNnehlmFtwNaAUEVsaRRLu0bet7ZZbgREBRwDQQFXANBAddAUMA1EBRwDQQFXANBAddAUMA1EBRwDQQFXANBAddoLOjIhgDYQXN779bSWNCmy44BUDjVf3GSoIcLbwMArQtBAc9wJOjPsybH4xZfUBBuBM15M2z0tDZR7jV8HygUbgRd9ESqMF3QWNs0QO9wI+gTieI0N+I3bfMAncONoFX/kmbtd2qbB+gcbgRttleaVT+lbR6gc7gRdEGrdGG6pH6utdbAaeR8l8ieZPn83Aia/Uq1qXO7V8NXUM44WbItcyJt9I7CjaCUHnyrwzL0hMUbJ+qzz/H+OPn3+BE0c06j0oNxGZQ3IKiR/+r3nFfrveBtVlsDpwFBjYwbQP+oTg9Xum+1OXAWENRIvWP03/KUNjiqbR6gEghqJPw8zShOabtd2uYBKoGgRjqtobREekbwZW3zAJVAUCPfVP2VVjj/ai9t0wC1QFATK0P+51/upTvapgFqgaB5pH5XfbO2SYB6IGg+2uIMiTsgaD56rtM2CVAPBM3Ha59rmwSoh6GgJ340LuhG0JGztU0C1MNQ0L6m4XZ0I+g7k7VNAtRjIWjuxtEj1uTY3CY322JVZoGfsA0P7+pP0LnDtE0C1GMu6K0Wj8+a07LRv7IbnCMJVTxqzRSW7rwV6RMRm0tpxbWjAy/Q5Q396iwVVod+3JaU7nWTNiWE/CdtohtBl/fXNglQj7mgr0eLzzzEPC+7wTni+/qm0R5TKH0mcPaWt8iXgqCNuq9Pn+/99o7hHp8Igga+uO8Dr6H05rMdrxmen9CNoF931zYJUI8kaEaKiZsl/xBnVwL+yVtldpP5OfKsMB1T8g59Nk5YqDlWELRuLr0bFCu8ej1YEFR8rKfbE3o8xO9pqW0SoB5J0M6BJsp4GObFSueterzgBudIvDA9RcTb0lJPLPYeLQgaQ+lRcig5OXkluURDxVHuRzfVo6A/PaptEqAe80N82Svi9FYp2YdzzpGDYgNB04N1PSp1qiQKOofStcTASRo6i3IhaOaNFBvvWhU0KVxpEsAac0FHvCCcj2cPfE12A8Mn6G/k+xTvaOFUqqko6Fzh6EiuGxuEihcTnS3opQnhHoSUiBwnN5yqVUGTyypKAooAc0HTnqkVM/6RjrdlNzhHnhOmE/xu7SLnKE0PMQp6vcQSYfXkNnwIetwvbMi8lSvmvxER+LP1FlYFzfLCU/G8YXmh/sCM9/bY2OAc8YnePr7YBHrBu9eBzY/5d06WBKUxxafvGOMxL5+gr9Q5ZrhgWvSCtuyUZljIeqGt9RbW+wf1v6skCygCFP+SdI583aVMjfeEj5o1Nf2bbP0i8G2DoLmzo/weFs/r8wTdG17S8EFc9IKWijctHShjvYV1QStdUZIFFAEOCHpMaY6iF7RRtGlp6mPWW1gXtPZpJVlAEeCagsZ7dFpy6MzZIyue8Yy33sK6oE0PKckCigDXFJRuayVd9PJovUOmgXVBOyQoygLYo1jQrAuZSnM45Tpoyq+JiSdvmq1MzevPLHK/tY2e/0phFsAal75hOedPs35Cfsjr0MzrQ2sbDIxzIAtgiYsKuq1n12V0URDxmS5zZTPA6s3zY2YqygLY45qCriONO3tHB8QmTPJeYr2FdUGn2+goEjgF1xT00SGUfkreF5bGy9z/YV3QBUOVZAFFgGsK6red0hvSXS0J/tZbWBd05YtKsoAiwDUFrS58lzxClgtL82pZb2Fd0C1dlWQBRYBrCvq+z/B3KjeutCt5c9BE6y2sC/pdCyVZQBHgmoJmTapUfvD9AYSQLvest7Au6M91lWQBRYBrCmrkzArZ7mitC3qhigNZAEt+LR7BnLI2Hjd31iMf1gW9VVrbLEA9fyWxR+YwK8KXoDmetnsEAG4HX4LSUv9pmwboHc4EDbuobRqgdzgTNOqktmmA3uFM0CcOaJsG6B3OBO2CkeZAATgTtM8qbdMAvcOZoEM+1jYN0DucCTruPW3TAL3DmaDvx2ibBugdzgT9ZLC2aYDe4UzQL3trmwboHc4E3d5J2zRA73Am6MHHra4Gbgtngp6qo20aoHc4E/RyZW3TAL3DmaB3ArRNA/QOZ4LmemVpmwfoHM4EpYHm/Y0B94Y3QavKjboA3BPeBK13Qts8QOfwJuj/9mmbB+gcZw3kJSfoU5uU5gEujbMG8pITtN8XivIAV8dZA3nJCfrGfCV5gMvjrIG85ASdGKskD3B5HBL00qlCGhc+kJecoB+MtrMe4B44IOg3FQmhbWweigsfyEtO0M9et7Me4B4oF3Sl18CVhE72WGSjceEDeckJuraXnfUA90C5oA+PoMnCizFRtloXOpCXnKDftLezHuAeKBfUb4ck6A6ZzudNWB3I69RAE94fWd/ssMxXAuCmKBe0wRRJ0Gl2dIZ86w+z3hRvfL7IiM+n1jf5TaZPe+CmKBd0mXfsIXJ9SfE5tlqvbHaSXutKSMm5Mg3kDvFXK9pZD3APHDiLnx8kfLssMU5mCDmJeaTFNdq+8qcJY7xWWG8hJ+g9XzvrAe6BI9dB7x5du+eGzcbh4ym9RsTxYscpGshLwFvxWLnAlWHzS1LQBuGEiNwVlnaUtN5CVtBy15UkAq6OckG7mrDRuEfnTJpVcrewNF7hhXpa7U87CwJugXJBXxbpEVzsDRuNzwRHfbB3eqWle8d7y4wALytog2N2FgTcAkcP8XfbN7bV+o8hZaUr9XXWyjSQFbTVt3YWBNwCh7+D7iO2vyzeP394+yH5B4xkBX16g50FAbfAYUGX+9q6zlQosoK+vFRNWOBqKBd0lcTUkJaq8soK+pbMb6DAPVEuqI+EX9MzqvLKCvr2VFVxgYvB21Od9MMR2iYC+oY7QZe8om0ioG8UClouH6ryygoa/4yquMDFUChoXFzc3IDKw+dMrBmmbsgtWUET26iKC1wM5Yf4wU9mCNPsjoNU5ZUV9IdGquICF0O5oJUNTxltqqQqr6ygf1ZRdX0VuBgOCGq4B3l+FVV5ZQXNfVTmMTvgligX9PXSW4XpttKMDvF0Z437qiIDl0K5oKktSdmosqT1XVV55QWlbWUeVwLuiCPXQffMGD5rv8q8NgT9ueIdlcGB68DdhXqBF6ZpmwvoGIWCLk6ki02oymtL0PNB11TFBi6EQkEDetMyJlTltSUofXOYqtjAheDxEE9vlDunbTagWxwU9GyCyuFibApK3+mjLjpwGZQLeqndaLrBk5ST6TrZTmwLmlr5sKrowGVQLmiPSpto/Q7n29l67LhwbAtKVzXIVhUeuArKBS37Gb1KEumX5VXlLURQ2mqhqvDAVVAuaJk1dKlPOt2qbtjXwgQ9Xe6qqvjARVAuaPvHdz/Snd7tou62uMIEpaNfVhUfuAjKBf25Ail1nNYovrWQDRwbyCuPO6EYcw44dJkp7dhNStf9YbO1wwN5PWBNVFbaN3G7MuysD7gmbIahcXwgr3x0fDOi1aAnI7+3s0DgkrAZhsbxgbzysb/YOmG6vYLtrkiBa8NmGBrHB/LKx6yG3cXZSx/bWSFwRdgMQ+P4QF75GDL/kdXCbA56cnBn2AxD4/hAXvmY+Pbxiv9QOna6nRUCV4TRMDQOD+SVjx+rfBBSovHsSoWdkQFXhtEwNDIDee0hecwsNGNuDZ83QtuXsGM8JuC6sBmGxsQm2U5u7fgE3dbwUHRzvx110OWyO8NmGJq8jXbLvWOHoKPED9mY7u9MtCsVcE2UCpq6c+1F4ej79+l9XWw03tjXAGnTt6/1FnYIOnSBMEmv1XeUnRUCV0ShoL+FCYf3Ib9UEU+AbDTe7kcaNhUgtZs2td7CDkE/k/q5+6k4ulx2ZxQK2i1o5cmvyldutmbXUZsdN5xt2OR3qvIQfy/yvfs0Y0qFxuhpxI1RKGh58XvhTHKh0Ob3YwI+VSkovdg98NEyPS93m2RnicAFUSgoEcc92mDXE577wjr9o05QSlOO36L0eqhsFODyKBVU/GFok32PIN96vpxaQQ3sq4SOHNwWhoJS+tXw3+XeUiIoHdcJfYa6K0wFtYEiQe83LeRnK+CyKBXU28fHx5tIQyWpyqtIUHohZI+qbEC3KBR0Uj5U5VUmKD1QMUlVOqBXuOybyQof1bunbQFAH+hFUPoqumtyS3QjaHrj2dpWAHSBbgSlF0MStS0B6AH9CEr3VTitbQ1AB+hIUPpV2CVtiwD8oydB6XtRt7StAnCPrgSlb7bK1LYMwDv6EjT76T74Vd690Jeg9F6TydrWAThHZ4LS6zVxOdSt0Jug9Eq1WZoWAvhGd4LSSxEYbNaN0J+g9ELVzzQsBPCNDgWlf4Su1q4QwDd6FJT+GrJKs0IA3+hSUHqmCnpzcBP0KSi9UD1Go0IA3+hUUHqt/jD8puQOsBNU5ThJhZHSbKA0nGfuxlHDVmFgT5eFkaAajJNUGHfbP5VK6e0nm82e164exk10VdgIqsk4SYWRFV33Io0eKB7pp/RQHQ3wCRtBNRknqXDmhv4Y9Le4kF4qVYNwgEPYCKrJOEl2sKOip2EhXO6rBNA5bATVZJwkezjmOVY8xKeWwlPzLgobQTUZJ8kuooOaT3s/MbqfNtEAdzA6i9dinCS7+K9KseKP+AYV3qMu0CfMroNaHSfp+ueLjPhodM9cTM/9LwWMHIpPUFeFlaDZ59OleXqBvmdPDTThrdGP6RWFs6Pfavcvie+gLgobQbMm+xLfseLvO8tlttPoEH/fW5zefdV7gybhAHewEfQDr1HxI7xeocwFpcbroH7BMRgLxCVhI2iN8cJkFdnMXtChA8TLTG8//Xf7pue0iQi4go2g/tvFad/wdOaC3m7ZdPb8tvWv0tx55T7EPSOuBxtBG48Up/8GD2UuKM39etSw1ZKZSe3q/6hRUMANbASdT4btyqB0u2f/GNaC5mddxWE2x78D+oPRZabYUkT8Sri1EilKQem156tvYxEXOA1W10EzkzLEWda3cdbfZyMopQm1Op1lExk4Bb0+8iHL/bnlhv3HKDYoelxOUOHcbEDIQvTS6Cq4oKCUnugcviKHYXxQdLikoJTufyJqM9MEoIhwUUEp3Vqv4UZ8iuoflxWU0sSmdb7IYp4FsMWFBRUUbVntE1y51zcuLSil3z9bbuzFIskE2ODiglL616ig5w4UUS6gPS4vKKWpC2o//JH5wydAJ7iBoAL7+5V5cS/O6fWIewhKacq8eg9NwFif+sNdBBU4FRPWYA7OmHSGGwlKac6e18o9Ngu95OgJtxJUIGvXwPINpx1H57d6wd0EFcjeM7J66KCteJJeF7ihoCK/z24V0Or9Yzix5x43FVQgdduw2uV6LTyFoz3XuK+gIldWvFa93DNzj0idPkypX607uhLnDfcWVOTv1dH1ApqP+TLY9+kBkZ4rnV0OKAgEFUn9Nja4WHCXKVsHFccNenwBQY34f3z564kdK3o8MvKLn9KcXQzIA4Ia8TLc8hT4yowX6vlG9piw8oeifjb04NtvLc0o4pz8A0GNlJkmTrO8donTs/HTejfwD2k5cGb8z0U0fkj2y5HvfNSjGm4XMAOCGhnk9welOW1L51t1KfHTUT2ifCs07T1u0TenGd+aP7+1+Om5rB6uehUEgppoUaxmY/+SP1i+cfX71dMHtKvlW7Zu18HvLE04dYNJ+mbfLnmuY8zf1U8xia5fIGgee/p1nWnrp6UbJ7Z8MvnlDnXKFQ9t0u3ViXNX7/7lb+2+Mz7UqN2ahHHBjXZqFlFvnJwzdaNlB5q6HUzWiWRePrLls9hhL7SMCvEuGdGkY583p8z7YsuBU1fuqAga1lqc7vZy10/Q3CEBQaWCIpPM1+t3MFk+uP3nkR0r5015s1/XJ+pUCigWGPHok117DxwbO2fR2m17jp29kGJvz+RhdcWzsfl+7trF6ewS0b9eXlUmwny9jgeT5ZDsm0k/7d365aIZk0YM7NW5ZcOaVQK9SWCFiLoNW7fv1Wfg0JhJM2YtWrRuXULid8eOJSVdSknJ60Qq/KWw0dPbRDXf5cz6nUjIc+L0N09zX3Q9mKwuyE35J+mXY7t3rlu1aMGM2JhRAwf26tWxbYuGDSMiQgMDixMSEBgYFhERUO3hsPC6rX069+rVf+DAgW/ECLw3QyBOHFZq2TqRjYkSB44ZOJ1k5GpKHmq+ZTgV7wRpVnqO2fqiHUz2h7YmvD5UEs+lSU1JuZSUNL/KlqOJCU83Fjz8QjBygejmeNHSQeKwUi/1EnnasO+aNzRQO8JISGAeASQ/3oEWhETIU6uhOh5r6zjFotq2FZzwM//gKtrBZFMTTUTuVxLPHVgY3LVveA9tH4++n2LB1SR5zhxTx5FEx4mokpB4mS7yvGb2JzhrMNmmh5XEcwtubFnprqfwAt/7hr67sIvf8+brnTWYLAQFBUmsUqa8X4xFz8NFO5jsAyAoMCP3/PF0y7XO+iUJggK7gKCAayAo4BpnCdq6pOUlOn/iwQyGoVmi07KLWf7jOorXH84RNN3yCl3Kl23OMyPsO2ahJ73KLPQfXsxCn+87jVno3ZFW/nUdxO7nGjQW1BpburGLHc6uR6Y5I5iFzvJiFppGf8ws9JlazELLA0HlgKDmQFDlQFBzIKhiIKg5ENR+IKgcENQcCKocCGoOBFUMBDUHgtoPBJUDgprjqoL+Gccu9hR2HYQc+JpZ6NzRzELT9UeYhf4vllloeYpAUAAcB4ICroGggGsgKOAaCAq4BoICroGggGsgKOAaCAq4BoICroGggGsgKOAa9oLGNy7d6gSTyBuljqIGMIicuFmasSjdEJpF6QseC6g5SxxRT/uyTaHZ7XFZmAu6zWPI+o7+l1iEnh0cJ7BP+8A5TaTbjViUbgzNoPRYMnLbOK/JLMrOC81sj8vDXNBWHSlNC5vAInR0GxZR6eWP/0cki7QvPS+09qVnlhomTEf5Zmtf9oPQjPa4LVgLmkKWCNPB4SxidxzIIird0by5j2gRg9JNoRmUnkTEbvDjyXnty84LzWqP24K1oL+SQ8J0nodFB5EaUKNDA/96ixgEppGiRWxKl0IzKD3jnDi40wjfdO3LzgvNcI/LwlrQ3eSMMF1JGIzjllM8aN6mAWS29pENFrEpXQrNqvRVXmNY7XExNMM9LgtrQRPJWWG6gtga/ctBMteIw0P1L2VrGDkHkSxiU7oUmk3p1/uRl7PYlG0IzXCPy8Ja0JNEfEZmfglmCTaSc9oHlSxiU3rkg+eRNC59e3D4JsqmbGNoA0z2uCysBb3psUKYvlmNQeh/j4lDCG8h5gNIaIBkEZvSpdAsSt/uGS11tc2gbFNohntcFuaXmVo/Q2lWRAyDyIlktTAdVIVBaMPHHJPSjd8eNC89q1I/45LmZeeFZrjHZWEu6A7Pdw72CWTx+Hr2Y8Gx24cVW88gtEFQJqVLoRmU/i0Zu1wkXfuy80Iz3OOysP+pc32T0m3Y/NSZNrxWyccTWEQ2flFkUbohtPalxxlHobumfdkPQrPb47LgZhHANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcI3Ggm6aAYAdzLrnHEGb9osBoHDK/+IkQQ9rGw+4KHUhKOAZCAq4BoICroGggGsgKOAaCAq4BoICroGggGs4EvT4zIlrs7TNAnQPN4LmDK0ydnq7h//UNo2euLN3M4tRonQON4LGPXFXmC5spG0aHbG8YsunQvrednYZvMFK0Ozz0qCNNF1muEYLQZ/YLU5zI04ry+My7Ij4jdK01591dh28wUbQrMm+xHdstrC0XGY7C0Gr/iXN2u9UkseF6LhWnGZWuOTsQjiDjaAfeI2KH+H1ClUgaLO90qz6KSV5XIhww/fPdrucXAdvsBG0xnhhsopsViDogpbid4LF9XOV5HEh6v0ozR496uQ6eIONoP7bxWnf8HQzQX/Puw81MN5sk+xXq035qHvkb0rSSPz99TpXOPmdIB5v6MFQXGgrCBtBG48Up/8GDzUT9MpM0538JRZabPT91OHLM5RkEcmdGNyjZ8XBmUq3447bDbpvOxAbXLQjBOsANoLOJ8N2Ca5t9+wfI7NdwOdK4skzu8V14R+3x0htojmT+/M6Nx/6l7Or4A5Gl5liS5FzwmxrJcJY0Konxen1Mvr/CAVWYXUdNDNJOlpnfRtn/X2NBL3vbTipCv9Lk3CAO5z1S5JWn6ABKeI0q/QtbcIB3tC7oH2miNPPWmoTDXCH3gX9u3qfTTsGVVZ+dQroA70LStNmde807T+NggHu0L2gwLWBoIBrICjgGggKuIZvQVNP39U2LdAbPAt6pVfJhwP6/KNtYqAvOBY0rcY7GTRtfBR+ZndnFAq6MB+q8toh6K8d898AAAxySURBVOLu0qz9alWJgL5RKGi5fKjKa4egb8yXZrNGq0oE9A3Hh/iRH0izaRO0zQx0hcOCJr+tKq8dgm5uLD4Wer9uoqpEQN8oFzR39djRAp3Lqsprh6C5nTr/mHqk7dOq8gCdo1zQd0hN70pNQ/xXqMprz2Wm+3PqBtRfgKfI3BrlgoYPp593pmlNttlovHn0A6y3wC9JwC6UC1piKz1emdJNj9lonNiE+EQayb/+SFsTnrMdqRa4HcoFrfYhveNxne4PsNU6u0Vza6vv7U404vuZgiKB+6Jc0GGBcbRu9LletWw2X2hV0AfgEA/sQrmgt595lh7wJt7rbDa/vMN2OAhqwT/rF//o7Br4w8HroCm7LqjLC0HN+TC452s1Ol53dhm8wfEvSe7F+oevUJo1toOz63AiGTesrFQuaFcTqoqBoGa02SxOsyonObsQZ/FLa79yleNyzFcrF/RlkR7Bxd5QVQ4ENcPdO/A9XWFZFj312CTz9Y4e4u+2b6yqHghqRkNDh6mP/OTkOpxFn4/E6fUy5l3EOPwddB9R9X0egpoR21PsZSohwuIY5yZU/12atdxrtt5hQZf7quoLGYKakda6+fJNb1U44Ow6nEVNQ+cwzfeZrVcu6CqJqSHqukOCoObkrn6p+xRr57HuwSvvidPLgXfM1isX1EfCr+kZVfVAUFCApIozD323pc4s8/W4Dgr4YJGfZzHPzmnmqyEo4ILvQg/TzJTnXjJfz/FDc8Cd6PaFOL0XdNVsvUJB4+Li5gZUHj5nYs0wWzcsFw4EBQWIMPyE1tb8CTTlh/jBT4q9z2d3HKSqHggKChB1Qpo1PmS2XrmglQ1DcG2qpKoeCAoKMHKYOP25gnk/Mg4IOleaza9SyAaZN1JsvAtBQQGSq7925OynIRZ3GSsX9PXSW4XpttI2D/GXJoR7EFIicpzcKIUQFBQk9e0mtXr/bLFauaCpLUnZqLKkta2OEY/7hQ2Zt3LF/DciAi1TSkBQYBeOXAfdM2P4rP02G7fsZLzgmvVCW+stICiwCzYX6kvljWV8oEz+9XsDTXj4BgYGXaL0pUDMMbc1Vybo4kS62ISNxo2iTUtTCz4+fyvFSMC8lJTbworMFMwxtzGPUiZoQG9axoSNxvEenZYcOnP2yIpnPM3HhTeCQzywC0a/xW9rRUQ8Wss9fQxBgV04KOjZhJuFtE/5NTHxpHwjCArsQrmgl9qNphs8STmZ60d2AkGBXSgXtEelTbR+h/Pt8NgxKAKUC1r2M3qVJNIvy6vKC0GBXSgXtMwautQnnW612btdoUBQYBfKBW3/+O5HutO7XRqpygtBgV0oF/TnCqTUcVqj+FZVeSEosAsHLjOlHbtJ6bo/1OWFoMAuHLoOeumU6rw8CnrpsLt268ExDgj6TUVCaJv56vLyJ+jSAA8vj24YGJQzlAu60mvgSkIneyxSlZc7QdcVG0HpvjJNnF0HKIhyQR8eQZOFF2OiVOXlTtCw58XpOQ+V362BxigX1G+HJOgOf1V5uRPUe5c0C5zj5DpAQZQL2mCKJOi0uqrycido8U3SrJS6UcaB1igXdJl37CFyfUlxdR813AlaW+qt7zuPf51dCCiAA2fx84MIISXGqeoelD9BD3u2OnFjqncvZ9cBCuLIddC7R9fuUduRJXeC0mPVipGSU51dBTBDqaCpO9depDT379P7uqjKy5+ggEsUCvpbmHB4H/JLFfFxDgeyZSWZ8IegwB4UCtotaOXJr8pXbrZm11FbHTfIcSDCRLGZDmxujbPPBpftdEyjYIA7FApaXvRqJlE5DiLV7hD/Y/CCf28ur5igTTTAHQoFJWuFyQYNelvWStDWK8XpnpraRAPcoVRQ8TH3TfwImlMiXZoH/6NJOMAdOhc0q3i2NA+9pEk4wB06F5TWl35CPx2COzldFKWCevv4+HgTaagkVXntEjTp3dffL+R8bHPV/ZSeiIpTVQzgF4WCTsqHqrz2CPpphdGfjQhearvR5siQhyotV1UL4BiOx0k6XfGiMD0X/Gch7a5e1KQiwCUcCzp1nDR7S6tL+kCPcCzokE+k2UfDtc0MdAXHgr47QpoNxj3u7gzHgp4vL9b2Y/kr2mYGuoJjQWl8hb7TXqi4RdvEQF/wLChNXjx5ma3RwIDrw7WgALASNPu84S6O9GvW34egwC7YCJo12Zf4jhXv41gusx0EBXbBRtAPvEbFj/B6hUJQoBI2gtYYL0xWkc0QFKiEjaD+28Vp3/B0CArUwUbQxiPF6b/BQyEoUAcbQeeTYbsyKN3u2T+mwHbJ60z4fGqx0ZlP3t2O+45BQRhdZootRc4Js62VSIHtfullImC12Ra5MSFDJjzeSP3zosClYHUdNDMpQ5xlfStzr3vTw2Yrlja+JUxnPq4sDXB1mP6SNO667FsWgv7vG3Ga+9BZB/IA14WpoETeNgtBq/4lzdrvdCAPcF24EbTJfmlWy956gHvAjaBzOmQJ06/qqOt2FLgaTAXdfU/2LQtB7z9X54PFfcLUDfINXA5n3W5nISilu0a/ttCRLvOAK8ORoABYAkEB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVc4zRBxy2y4N0OfZnR+kVmoXt0ZRa6byt2obs8wyz0Cz0s/3EdpbKTBP18oCUtA2ozwzuSWegKZZmFruXBLHTtwIrMQkf4WvnXdZBht50jqDW2dGMXO/w8s9BzRjALneXFLDSN/phZ6DO1mIWWB4LKAUHNgaDKgaDmQFDFQFBzIKj9QFA5IKg5EFQ5ENQcCKoYCGoOBLUfCCoHBDXHVQVNeJpd7OrshkNeMIZZ6Gx/ZqHpsM+Yhf4zilloeYpA0Oxb7GInswudzvBJP4Zlp2awi82wbFmKQFAAHAeCAq6BoIBrICjgGggKuAaCAq6BoIBrICjgGggKuAaCAq6BoIBrICjgGvaCxjcu3eoEk8gbicgABpETN0szFqUbQrMofcFjATVniSMCaF+2KTS7PS4Lc0G3eQxZ39GfyV1xs4PjBPZpHzinyWhxxqJ0Y2gGpceSkdvGeU1mUXZeaGZ7XB7mgrbqSGla2AQWoaPbsIhKL3/8PyJZpH3peaG1Lz2z1DBhOso3W/uyH4RmtMdtwVrQFLJEmA4OZxG740AWUemO5s19RIsYlG4KzaD0JLJLmMaT89qXnRea1R63BWtBfyWHhOk8j0wGsWt0aOBfbxGDwDRStIhN6VJoBqVnnBPvVB7hm6592XmhGe5xWVgLupucEaYryQ3tQ+cUD5q3aQCZrX1kg0VsSpdCsyp9ldcYVntcDM1wj8vCWtBEaVi6FSRF+9CZa5KEaf9SDAajlyxiU7oUmk3p1/uRl7PYlG0IzXCPy8Ja0JPkiDCdX4JZgo3SCPYaI1nEpnTDIV5C49K3B4dvomzKNoY2wGSPy8Ja0JseK4Tpm9UYhP73mDgA6BZyTfvQkkVsSpdCsyh9u2d0ujhnULYpNMM9Lgvzy0ytn6E0KyKGQeREslqYDqrCILThY45J6cZvD5qXnlWpn3FJ87LzQjPc47IwF3SH5zsH+wSy6F8h+7Hg2O3Diq1nENogKJPSpdAMSv+WjF0ukq592XmhGe5xWdj/1Lm+Sek2bH7qTBteq+TjCSwiG78osijdEFr70uOIgWval/0gNLs9LgtuFgFcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQwDUQFHANBAVcA0EB10BQrelp7CiGTKcBi51djP6BoFpzOD4+PuwJYXKG9k50djH6B4KyIKq3sytwGSAoC4yClhEO8aEfNfOL+ORq1zIPiX1rLm/oV2epc2vTGRCUBfkF9X77264eoR/taOKTSud7v71juMcnTq5OV0BQFuQXtBelZ8lblG4jp+4GxQprXw92bnH6AoKyIL+gMynNJl+Klv58lBxKTk5eSZgMDOmiQFAW5Bd0tihovCToWuMFqJNOLk9PQFAWyAi6h1x3cmH6A4KyQEbQ6yXEYTQnF/2IrDoGgrJARlAaU3z6jjEe85xcna6AoCyQEzR3dpTfw3FOLk5fQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANRAUcA0EBVwDQQHXQFDANf8Hpqlws4GXNvMAAAAASUVORK5CYII=" /><!-- --></p> <pre class="r"><code>summary(m.L2.DFOP, data = FALSE)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:36 2022 -## Date of summary: Wed Sep 14 22:28:36 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:46 2022 +## Date of summary: Tue Dec 6 09:39:46 2022 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1832,7 +1825,7 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, ## ## Model predictions using solution type analytical ## -## Fitted using 581 model solutions performed in 0.135 s +## Fitted using 581 model solutions performed in 0.131 s ## ## Error model: Constant variance ## @@ -1920,10 +1913,10 @@ plot(mm.L3)</code></pre> <p>The objects returned by mmkin are arranged like a matrix, with models as a row index and datasets as a column index.</p> <p>We can extract the summary and plot for <em>e.g.</em> the DFOP fit, using square brackets for indexing which will result in the use of the summary and plot functions working on mkinfit objects.</p> <pre class="r"><code>summary(mm.L3[["DFOP", 1]])</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:36 2022 -## Date of summary: Wed Sep 14 22:28:36 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:46 2022 +## Date of summary: Tue Dec 6 09:39:46 2022 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1932,7 +1925,7 @@ plot(mm.L3)</code></pre> ## ## Model predictions using solution type analytical ## -## Fitted using 376 model solutions performed in 0.081 s +## Fitted using 376 model solutions performed in 0.078 s ## ## Error model: Constant variance ## @@ -2028,17 +2021,17 @@ plot(mm.L4)</code></pre> <p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAqAAAAJACAMAAABlpiR1AAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5fX19gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7////isF19AAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO2dB3wURdvAjxdFhSBCOgkkhJJQBAwBQYogKK8BQpDQA6hUKVKkRIp0FLBRpAmiEBRBglggCJ8EAQVFlKLwKr2XBEivl+fbvcuFsHd7O3szdzt7ef6/Xy53s09mZ5/89+Z2b3fGAAjCMQatG4Ag9kBBEa5BQRGuQUERrkFBEa5BQRGuQUERrkFBEa5BQRGuQUERrkFBEa5BQRGuQUERrkFBEa5BQRGuQUERrkFB1ZFrMPEQwKHgh8tU7C+WbfMr+x+f94UnBkMqwAqDB8DOwLKGcq1umP+mjOGfor/e16+fUZt26xYUVB2CoP4CdWFTGUPZ8gZDTYAFgq/lDIZXSwh6rqzhMe8yhkDz39wXtLLBkKtZ0/UJCqoOQVDT7/wKhnYASQ8bZmSXM0QCTDWUyb0vaJzBD+CEwWB+Cy0WdHSZMiioSlBQdVgE3WB4ROysRxn8NhoeFp8N7nnjvqATDQ99kg9nTphttAh6o2y3h1BQlaCg6ij6DDpjkKGO+HK/4eHRhuqWhcWC3vUwGMpUbLHfXGwRtNljuSioWlBQdRR9Bl3/qlnQQ4ayowzBloX3D5Ly32nhWcZQxmxokaDfG1YDCqoWFFQdli5+vbmLH2vw3WB4VCxo2fAfQdAL4jHT43Ds6zMAd58wdDPFFgka8p+mTcsYmuRr1HKdgoKqwyJobnnDCwD7HjZMzXzY8AoIh0VljPCIYSBAPeHIvrfpXfVJw1BTbJGg/uZPB/gWqgoUVB0WQWFjGcNDwifNGgBzDIbHHjWdZhpoMJQrazCsh1NlDQ8FPmYoe8oUWsbgGxgY2EJ8il28WlBQdRQLCgerP1TGo5/4bKPPf8pUmiE+m1yxzH+8PxWe/BjykOGhmj+aI8uY3jk9xKcoqFpQUIRrUFCEa1BQhGtQUIRrUFCEa1BQhGt4FPRs1wErtW4Dwgk8CjrxKLTWug0IJ/AoaFbh5W5atwHhBC0E3dfMr+EiYdWeXl5e74FxXgPvlt8/ELC563WFKpKiSYrkuHOeOBTRGA0EzfT9wXim3RowZJhevtLpkvFgyKYSAd+NL1Sqg0TQSYH+b5ierA71jc2CedWD34cLzUMXA4w47XDjERejgaCn/AX/fvqoSNA/q2aB+MV2wf2AQZ1iYh/4i7yin/uvCATdFZF1NzRReHIk+Hpq5MxvGly72ujAsMSC8NT/jWC0JTpBsceysCdGfFwTV7JMRcfkFDQQ1Fit2x7xkgmzoEuHmQqrnxUfx9c+DFvCil1cUqPWeOOBvr2niT8wN6TGRNMrMKXNtLDjNoCIfaan0kwe3AnQPUF4kiD8c9bGThb+bP64Se+fr5c14KprNpQTFHusYtLPi4+lXlC4PaF2heF3TXt0CEycYyp7xnzlz9aOWUE7LXFJjW5l9J57oNJpEH8SQ1NzWq0Xn4mLos0L1/WHsyF7TU/NmXwxyMRxU2WRnXPMNaU03bah6Z3kp3unDuvww29vunJjtcduj5VnHY+CCpz4b5+ihC0ebioIvmD6VVh7YJfioOnBbdvW73TgWQDxZ9JUgI9eFp+BmDbzwnu+efNnmZ9aZ/LalkZm2eNrxYPxdb+6/QeJr2LuftNxQo7zto035HssU3dk6n7SXq3T4htIiil4vdpTg+MSewPErYHJQT69cpOizcu0QgNBV/cVHo5XK0rYkYBs4fFQ1aKFC8r9Wxy44G2A7NQDnYVMCj8ThT56ZX/xGYiCmhdC112Nz5mfmgV9IcDEMaFn/xngPdF+Y/+o6wAF+QBvzRZe7lqYHH5x+nsu3GCtke2xxO7I3BOtiCn8ZZAg6MctslKCiwS93CE7v9m+pGjzMq3QQNB/nthbmDch1vKZKLbrtcLDNROKFvb2vP/e9kfYrZzIHRZBd9ZNz229rlhQ80L4om27ojjpO+jSdqmp7RcXns5PaC+eFPi7XsbNEOFtozAq+2L7ws1TXLChHGG7xxK7I3P3czJg3K5CQdA+KwFet7yDXto43TMxKdq8TCu06OIPPONXrX+KRdCCWfW8G8dbFvl3ib8fuLJmwGiwCApzaoW8YSwW1LwQMsp/WhQnFdQ4JihwZH6O4cK4ch4eHgNhVtVGXwvFm9YAzKjb+obTt5Ib5HssMZlFPVHK5wOaCoL2Wy30MyZBh685WHfFb90FQc3LtIKvb5KMTdbtxdseGSPfY4mCmruf2RPzb1e8mxSztmVOamhcUv38e/5r3h4F5wK3J0Wbl2nVer4EXRtRCLGB57Ruhpsh22OZuiNT93OllV+rNaaDpKDwKXE5A2oEdllzKbJBlzdDkqLNy7RCUdCTMwYPnHbUFU1BEGuUBF0cNuuTT2bX/9AljUEQKUqChqSJj1kNXNAUBLFGSdA6F8XHm41c0BQEsUZJ0ITA2KlTB1Td4pLGIIgUxYOklPXz5nx62xVNQRBr8Cge4RoHj+KvryrFrE5j/m/AfMqh9ig+Z/pkkdY+Q0svgbtYOPkAH4dpvVEaYjefao/i8z94RyQ8iMF/Ra90TGRe5eohzKvUD3bz6eBRfP8wykbpGRSULVSCyhzFo6BsQUHlcPBiEYug54a1G3MN4MaxLMfq0SUoKFucKOjNwPk/xtVOHeLd0E+7mwJcDgrKFhpB7zUMNiMpLxJ0uXgrQNTw5zLhiK9mVwy6HBSULVTvoKlhCVdEJMVFgs6aUVj/w+HNxKsLO+yhaKK+QEHZQtfFzztiq7RI0L1hKaefffj5RQCF9f8AuPLB23861kZdQSXozuz06dELpVMpoKByUB4kzfBq7jmo6qMLfxrcxghHfUZMDljvWIXO48MOLFhaokYaQV95PqdXt/hO0rFNUFA5KAWFW4fuQtaQR4In3gOIFOQ85a3dDYC26TV1Nz1v9itRI42gfjngmwFp1YoLjh0RmdzT8Sp1j1MFNXP+pZrbAGqdEZ568XblUy+bo7yoJJ6VoE+fgPBzcKa25XVGyyYinv4UjdM7LhAUYE+DDiej1gAc93OsQufBl6BHg3tG+Q8I+FJSjF98yMFKUMhf6tPXe+AoX+4ubeZLUMjesnDGisvSUhRUDmaCAqSM9Or94T+O1edEOBPUNiioHAwFBfirY13bA09qCQrKO64TFOD7sI4nHavRaaCgvONKQSFvic+wm47V6SRQUN5xqaAAd8Z5vS2OTwWXNv9fgWyU60BBecfFggL82z0ovhA+9e3RLIL93TuqQUF5x+WCAux/OmKH5xmA4RwMw1lC0LyFLZ79WN03XVmTTb9QUCeigaBQuCngib8AdnV0rHqWlBB0fIeDu59aovQHDwzbfs98oSEK6kS0EBTgsofPkKuLhjpWPUsEQU9vNvFlhbWbNy+obn6x+VJxxMretQL6ZZiHZL8/bPuhTs+FdL7Z8xFx+FcUlIb/vTlmp73l2ggKrzX972MVbF6r51oEQTf0MNH9oZgePTqXN7/osaM4YuXjlwq6zzAPyX5/2PZDZU4WRn+E76DUHPaZ8WG9eXYCNBK08Itho3v6fqD5dBoluvjnFkHB669aRax8CWBHa/OQ7PeHbT/UBGDuuygoNVEbAG48bmO+GwsaCWriROfg9UbHVsGKEoKeb1rN54U7VhErYwB2NzcPyX5/2PZD/0VBmdDwmPBQ/ZJ8gJaCAvz0TMPvHFsHIx44zXTR1uQJKytdNfacYh6S/f6w7UWCBppCUFDHeVno3n/3s3PyRFtBAbY3aLXfsbUwQfk86MrolkG9081Dst8ftt0sqLGZafIQFNRxrtR8vrePvZt+tRYUjOtrdPrDsfUwgEDQOKUIFJSG7G+/uGZvueaCAuQuq9pLqymwUVDe4UBQgMy3vV8+79i6KFEWNNnu/m0CBXUiXAgKcO8tz9ekt9e7AvwunnfoBLU9wrJDCb090XOc6y/FQ0F5h0pQmRGWHUzo9TGek1191ycKyjtUgsrMk+RwQi+P8JyS4uDfOgYKyjtUgsrMk0SR0IvDvKa5UlEUlHe4G2H5wlCv6a5TdIBPCD0+Jb/BR0HZwnSE5eyRpnHvawdSNenCEM8pyVQ1kJNxlgWZJWqkFjTSqgQFlUNJ0MKt2yA+Zna25bVxvWnmkOYhlI26MMxz8i3KOrSCQtDbsSIVYmMl5SioHEqCjolo1a3VZ12kg68xSOilkVXGK58h5xEKQfMneS3ZtMl7U/Hn4ox6po8QFX2ZtEyfUAlaNSOv0mlIDZAUM9njr47zHHmRQT2uhqqL39v0Wyg5XvUl0yeI6Dq0jdIxVIIGpmQ89hvcqi4pZtQl3YzzHPQvk5pcCd1n0Lt9hlgPZYddvBxKgr7rX31C2JTG0yTFzBJ6Z5ZP3+OM6nIVtAdJ8dFWRSioHIpH8X8eh/1xG6WlDBOavtC/y8/ManMFeJqJLZxcLCJlS+/+RXM0Zi+v0Zb9/9x5oKBs4VPQBQ03rg2JL3qRv+HJ8C95GCeHCBSULXwK6nMB4OfGxS8Lv21Za2W2nXiOQEHZwqWgueWFN8ybPiWLDnTxm2d9yyWHoKBs4VJQaCYceC3s9mDZX694jtXBiVEUlC18Cnos+JnwelY2Xp3k1c/q6mjeQEHZwqegkL3/cL6N4tR3q7XfwdtcSw+CgrKFU0FlydvQuP4azUfMsQMKyha9CSqwJ9JvFr/XOqGgbNGhoAB/D608hLfpGCygoGzRpaAAt+f4v8Dnh1EUlC06FRQgd32T0I/SXbAilaCgbNGtoAL7YzzfOOuaVZGDgrJFz4ICXJzkFbWHr54eBWWLvgUFyFz1ZL3lPPX0KChb9C6owN7unmP+59I12gMFZYsbCApwaYrvC9s5uSAPBWWLWwgKkBPfImg+F9OAoqBscRNBBY4OrtxXy9HEi3CVoP8s/+we8zVxiPsICnB3cd0nl2n9X3ORoF/4De9ZVX/3vKrHnQQFKPyxZ+XBv2m1dhOuEbTQ+2+AJT2Zr4o/3EtQgRvzazRZpeE8yjSCZs2M/Vr4NVBSbCOfV8XBMv6u6/iqdIPbCQpgTHyp8lDN3kZpBI2NXBsRD+AhKbaRT+MTFwA+6er4qnSDGwoqcG1eyFPLtfk0SiOoTxok17p+X9C8998RCQ+yDl0R9NZon2OOr0o3uKegwlvMDz0rD9TioJ5G0DrXATZEFRYLmvvWZJH6D4wtlH/adDHs4VmLtJh1wuVwM4kCe269Vzdska3ZDZ0KjaCraiwS+vn25STFD+TzYHBolX52Jl91NziaRMEJHHy1crdvbd3c5DyojuKPJAjv/VvHSkpL5jO/2k7I7voOxTp0BleTKLDhgaub0te29J90yoVrd/JpphNPCg8/PM98HdzC2SQK9HxVu2yjpAdKTk3yf2a1y46YnCzolaoFAPGl4QRoEdxNokDJn1UPG3f6XX2wMP+77k/02+WaqemdfaK+c7+fN1f7P+br4BamkyhkRnUQ8a/KomWOMVccrLT/eqvy5KURgXGu6OqdLWj6lBadfmC+Cn6hE3Rndvr06IW5lpeFe3eLdKjJpGkOMf9N4aGv1ZClIicnVX16mdOnD1ESNNoCeZX8HHRqAImgxrmdUnfYWv7K8zm9usV3GiEp1jKhJ/1/zNzsL3NyqSCx3xPRW5077IOSoEkWyKtEQeUoEnRc7/rZkaNtLPfLAd8MSKsmKdY0oTsaV2h1SH5x2rrnPIftd+JtTARdfO6+xMTtKqbqQUHlKBK0ZmZbyJZO5SHy9AkIPwdnakuKOU/o5QVPBk/9y1m1Ewjas5ln98Ap5FVynk/nQiJo8O22cKuWjeVHg3tG+Q8I+FJSzH9Cj00MbLzoslOqJhDUO3/akXNR5FXyn08nQiLo4kbB88NW2QrI3rJwxgqrf7QeEmr8cVCVtqudMCkogaDVc+JXgPTbDTvoIZ9Og+gofu/0BX+oqFMnCc3Z1rNSl88zGNdKIOjkDhfrjupEXqVO8ukcSARdY4K8Tv0kNG19ZKVe25ge1pOcB70MB99W8QFDMZ9fhAf2cc4HFu0hETQuLm64v7t+Zkpe1a7KwB1Flwfl/32etj4CQVnv8Ikh+y6/9RQnt12zhvRE/V3puSQ76EpQgetLWnoN/iFfOHaqXTfg+VS6yggEZb3DD/hUeGjgplcvkwqa797n7S6939x76O4Gn0PB8FF0NRF+1clyh++5WXhodpi8Qj1BImhvgfoDyevUoaACF95t8p8hu/JO1qerhlBQljv8+qbX4fNgnsdFp4BE0ESBpFw7cRL0KShAboV5zT07NqD7TxMIynyHn1nRown30584iKKgcRbI69SroDDihR9XeoVVjt2W5XgdBIKy3+GNrE+WOUoK85tRFAUVjjcHB78+vc5U8jp1K2jeu+06b4dry9tX6vG5o8dKSoK68w7/+5NVPCYxvs6BpIt/KgUgNZy8Tt0kVJbktZ0fj1zt0GBkSoK68Q6fF7wZklutZlspiaCBQo+XHURep14Sape0L3s90fo99SOME3Txju7wf3St14ufgVCt+eMp4eFrFRe6kkAi6OhnVq5uOYa8TrcQVCBnx1C/htOPEPRZ+9pW7/q3+SmBoA7u8Jf81v71YXUnXD7Aiv+Jl7Vt6Mu2UhJBCz7u02+tii8q3EVQAePPk0OrjfhB4YjmlF/CxWXV75qeEwgqv8PbHWdg8UjhoefnitVrhrHFhDN7ajC+5UVR0JjzEGOCvE43ElTk9IKWlXtssPfONVe8uLNrguk5gaCyO7z9cQZMt1sNW0nSZI248WrN5tI7KGlRFHRPOuAtCrfWvVSpzULZe+6mzBce+plvhFIS1N4ObzXOwNWzItEhwkMWHAref3aXz3Hhneqc6XXp+N2B9Lt4+2l/APcTVCD7++HVao7dbbOz3x9yFn72Mw+VpCSovR1eOs5ARliISMWHhYf3Acb856GyPsLv/wsxvy4Vv8MIBP1mkLFpRbe83E4lf8x9+omYT2zckbfC0yuo6LZC0u/ibRUqjTOQd8lNL1myA8lBUvNfEgfcbkVep9sKKnDz056VI2Yclg4CYSweHIBAUPkdXjLOQBHunE9FSAStUzj2i3x/8jrdPKH5eyfW9+n/hcwd9gSC4g6vAhJBX+kRkDKhNXmdpSChF1Z2rdR8ltUbKRAJiju8CkgETVtxEmarGC21dCQ0d8+EBl5910s/kRIIiju8CqhGFpGh9CT08pqYKk/F7S15aE8gKO7wKqAaWUSGUpXQgoNvNa/UeUnxl/YEguIOrwKqkUVkKHUJvbNlSPGInQSC4g6vAqqRRWTAhNoHd3gV0I0swvUkCtpAICju8CqgGlmE+0kUNIBAUPkd3jaYTzlMgp7/JQuMN88ss7Gcz0kUtIXkq073HErIOSgK+kGZ6t7f+z5e7Vkby3mcREFrFAW1s8PLgPmUQxS06q+wz7De9gQEHE6ioDlKgtrb4WXAfMohClpF+KksN0GG5OKGjGeaiHiq+BrP7VAS1N4OLwMKKocoqJfw42evhsgSz48dEYmUDrlcmlAS1O4ObxsUVA7TO2hycrKP8GPj4p3bsSIVYmMl5ZhQOyju8NZgPuUQBX3UgvXi/EleSzZt8t60SVKOCbWDnR1eDsynHIrzJO1t+i0EW5ViQu1gZ4eXA/Mph/J88Xf7DLE+IsKEsgXzKYeyoADx1iNJYELZgvmUg0RQG2BC2YL5lAMFVQ8KyhYUlDE0gmbNjP1a+DVQUoz5lAMFVQ+NoLGRayPiATwkxZhPOVBQ9dAI6pMGybWuo6AlQUEZQyNonesAG6IKUdASoKCMoRF0VY1FQj/fvpzldWZkBxF/vPhGBhRUPVRH8UcSAIxbx1peFv60W6SDittD3A4UlDG0p5nOWBdhPuVAQdVDK2iEdRHmUw4UVD0oKFtQUMbQCmpjFhfMpxwoqHrwq062oKCMQUHZgoIyBgVlCwrKGBSULSgoY1BQtqCgjEFB2YKCMgYFZQsKyhgUlC0oKGNQULagoIxBQdmCgjIGBWULCsoYFJQtKChjUFC20AmKkyhYgYKyhUpQnETBGhSULVSC4iQK1qCgbKESFCdRsAYFZQuVoDiJgjUoKFvoDpKkkyjUDRGp6MuiZTqFWtBIqxIUVA4lQa0Gu7p8ViS6DoOG6RUKQXHMfxtQCYqDXVlDISiO+W8DKkFxsCtrqLp4HPPfCrqjeBzsygq6z6APjvlvjF8l0jyEtlE6hkpQ6WBXRaCgFJQc8z9rxFCR2oF0VeoauqN4yWBXRaCgFODYTA+CF4swhrehb5LfaNP/BMXfaw0KyhjOBM0Nf33fEj8bb8t6AQVlDGdjM+1rITxMfcvxCrQGBWUMZ191bu0uPCwfwawtLgcFZQxngl7xPQWpzTaza4yrQUEZw5mg8LlPuNc4Zk1xPSgoY3gTFNKP3mLVEC1AQRnDnaA6BwVlDArKFhSUMSgoW1BQxqCgbEFBGYOCsgUFZYx7CJqSeNjo8pXaBAVljFsI+p1vxwbNU129VpugoIxxB0HzfX+DwhETXLxW26CgjHEHQU/VFR4OtHLxWm2DgjLGHQS990Q6wNpeLl6rbVBQxriDoDCi9cZFPr+6eq02QUEZ4xaCGtf1HcnJZfgoKGPcQlCOQEEZg4KyBQVlDArKFhSUMVSC7sxOnx69MFdSivmUAwVVD42grzyf06tbfCfpHUSYTzlQUPXQCOqXA74ZkFZNUoz5lAMFVQ+NoE+fgPBzcKa2pBjzKQcKqh4aQY8G94zyHxDwpeV1hofBhB+LhukUFJQxVAdJ2VsWzlhxWVq6eghFlXoH50lijBMGD0NB5cB5ktTjhLGZUFA51M6TlOFfWeQR/MzkOCjog7CdJyn9jsiHmFDHsTF4GAoqh4PzJGFC2YL5lEPtPElFYELZgvmUw8HTTJhQtmA+5UBB1YOCsgUFZQwKyhZnCLrOq4mZ8IfKE1FW/2Hlija5yeP/51jSWOWz3COKIY8RVPPwo8rVPKxczUOPKYY8WqmJPezm00FBc48UceCRDSR8/CjLsNXlicJWViAKW+FBFLa8kmWbj7EfkUNVPtu/rBjyrq9yNeHjFUOmhilXEzJbMWRE+yP2sJtPBwUtJrM8UVi6dKo626RVJAq7V4ko7E5lorBkT6KwW95EYZSQ5PO15Yoh/9ZSriZqu2JI0rPK1UT8phiysa9yNXKgoCioPCioBBQUBZWAgqKg8qCgElBQFFQCraB5bYnCctuRhT1HFJbTnigsqwNRWObzRGEZHYnCKCHJ50Jls252U65m/GHFkL+GKlfT/5xiyI/TlauRg1ZQBHEqKCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNTwJmm9165Ozucc0jBtclUiCvFCnjlLQbaEhi5ViVof6xmaRRE6KJajw6NP+gwjC3q0WMLZQKSy1IViqshdpCiPeCioU6ydqhnIiCbKonEGC9NGnjk7QlMCLqaEn7cccCb6eGjmTIDLJK1a5wuzg4znP71QM+9M/OSNiq0LYvLBgyybYizSFEW8FFYr1EzVDOZEEWVTOIEH6GKSOTtD1wq46c4b9mIRFAGtjlSPvNl0Sq1zhlp4AOdmKYcerZ+a12qYQdiQh2LIJ9iJNYaRbQYdi/STNIEgkQRaVM0iQPgapoxN03lSATwYrhqU03aYc2eebTbHKFS7sGu7bJ0O5tgmPVYpWrO1KsGUT7EaKYaRbQQdJ/YrNIEgkSRaVM0iQPvrU0Qk6R1zrQKWo+FrxypHxQ0DIq2LYDP+z6d3eUgw7EPbrqZbxSmFi+swxdiPNWSbaCkoI6ldsBkkiCbJIkEGC9NGnjk7QdcIaZ0+zH2PsH3WdIDLaP9i7QqRi2AohID5aMSxuCsBnMUphYvrMMXYjxTDCraBEsX6CZpAkkiCLBBkkSB996ugEvRVwK7v+7/ZjEtoXEkYKO75i2OWAs+mdFyqGbWx8Ky1qtlKYmD5zjN1IMYx8K2hQrJ+sGYqJJMgiQQYJ0kefOsrTTJsb1f9AIWRcOQ8Pj4EkkUJelcO+rBX4Wq5iWOG0alWHZCuFmTogc4y9SDGMfCuoUKqfrBnKiVTOIkEGCdJHnzqeTtQjiBUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjXcCnq264CVWrcB0R5uBZ14FFpr3QZEe7gVNKvwMsFsaYi742JBDZ5eXl7vgXFeA++W3wOcNIhjmt76zzTIfbNWlVZHSoZu7npdobKkaJIiOe6cJw5FtMPVgmaYfr3S6ZLxYMgmOFlBnJxzXYVpMKJ7csH6gMz7kd+NL1SqjERQ8+i+YBkDeFKg/xtwoXmo8HTEacc3A3EZmgj6Z1XRmYPVC042Cr0D8NJL0y55pQslr/x8P3JQp5jYB/40r+jn/isCQc2j+4JlJOBdEVl3QxOHJRaEp/5vBJst4hsVPZaFPTHi45q4kmUqOibWaCLo0mGmF9XPnmw8Lh5yas+ftq3L/ZjxtQ/DlrBiF5fUqDXeeKBv72niD8wNqTHR9ApMaTMt7LgNIGKf6ak0k+bRfcEyEvDBnQDdEya9f75e1oCrTt9YDiDvsYpJPy8+llpBhT06BCbOMb145seTjZN6wM6h70xbXHLs3a0ds4J2Wl4kNbqV0XvugUqnQfxJDE3NabVefCYuijYvXNcfzobsNT01Z/LFIBPHTRWkNBX8LR4qeGtk55zUYR1++O1NV22yppD1WHnWf1hqBTUlbPFw04vgCycbFwTnjvjunWlfRYkFv5r7nMLaA++/oU4Pbtu2fqcDzwKIP5MEzz56WXwGYtrMC+/55s2fZX5qnUlxdF+4P3bxtS2NTO7H3P2m44QcZ20mNxD0WKbuyNT9pL1ap8U3kBRT8Hq1pwbHJfYGiFsDk4N8euUmRZuXabEFLl6dKWFHArKFx0NV4WRjGLAjNOudaWerpAol7T4zR+RDVNMAABftSURBVC0o92/xXyx4GyA79UBnIZPCz0Sha1/ZX3wGoqDmhdB1V+Nz5qdmQV8IMHHMMrovWAb4TRDeMt4T945dC5PDL05/z2XbrRUEPZbYHZl7ohUxhb8MEgT9uEVWSnCRoJc7ZOc325cUbV6mxRa4eHXmz0SxXa8VHq6ZIAq6tUU0vDMNBnW9nbfa7445qrfn/fe2P8Ju5UTusAi6s256but1xYKaF8IXbdsVxVl9BjWN7lt4Ot88wO/Sdqmp7YXjhMKo7IvtCzdPccUmawpBjyV2R+bu52TAuF2FgqB9VgK8bnkHvbRxumdiUrR5mRZb4OLVmQUtmFXPu7HQ9QqCZjy6VhQ0Y0x17zZFI0Qf8O8Sf/9PVtYMGA0WQWFOrZA3jMWCmhdCRvlPi+KkgppH980xXDAP8GscExQ4Mh9g0xqAGXVb33DBFmsLQY8lJrOoJ0r5fEBTQdB+qwHeMgk6fM3Buit+6y4Ial6mxRZosVL7GJus29skX+tWuAcEPZYoqLn7mT0x/3bFu0kxa1vmpIbGJdXPv+e/5u1RcC5we1K0eZkWW6DBOhVYG1EIsYHntG6GW0DQY5m6I1P3c6WVX6s1poOkoPApcTkDagR2WXMpskGXN0OSos3LtNgCLVaKIKSgoAjXoKAI16CgCNegoAjXKAp6csbggdOOuqIpCGKNkqCLw2Z98sns+h+6pDEIIkVJ0JA08TGrgaQ4dXcpZg/7i0wwn3IoCVrnovh4s5GkeE1Ah9KL5w9UMtoC8ymHkqAJgbFTpw6oukVSvHoI5X9Ez3RMZF4l5lMOxYOklPXz5nx6W1qKCWUL5lMOB4/iMaFswXzK4eBRPCaULZhPOdQexWd2Mn2urVGPRct0CgrKFipBpUfxhftMZwY61LIUXKNrnB5BQdlCJajMUXz/sKInhdU1uUpQU1BQtjjjKL5Y0Jyfaix2tGF6BQVlC52gtrEIuuSJGpV9Ft2YNWK90bGK9AiVoIkJOTuHL5He0YKCykEn6C81LsNp3xCPcR8928+xivQIjaBvhHbv0ui99qMkxSioHAqC3msYbEZSXiToHHEEmuHP+Y005gSVnrG4aAT1vwcvfgX3qhUXHDQddI6OYdAuvUL1DpoalnBFRFJcJOgScQSuHs980rZXrmkthTYGUXE/aAQNzIKvrsG94h0+87+m03b+/iwaplPouvh5NodAKxL0sv+qE+8Gzeqb9VIb7+uQO9KjfFQpOO1EI+jEiJ0AP7cbJykuPugsjTjzIOlY9wax/2a3iuhTPuAqTI+6kxP3gmMV6gmqg6RdBwESlkoPKVFQOSgFNWPcu/Hs/OC/Gh4HKKiU7liNOsIJp5lQUDmYCGoi3re28Gkg7/Esx2rUESgoW1wkKOyrWOf0tUGl4HAUBWWLqwSFU16V/IZqMX6Pi0FB2eIyQSGlfedUx+rTFSgoW1wnKOSPrPuPYxXqCRSULS4UFOBj3x2O1agjUFC2uFRQOBg4u2gk3uw/LjlWOe+goGxxraBwrVVn00jeuwMa+/V0y68+UVC2uFhQyBsbckR4//TfD7nRCx2rnm9QULa4WlCALd7L4fcmwpMdkY5VzzcoKFtcLyj881SPk4H5AKsHOFY932gh6JUZIz/XZJYNF6CBoJAzKqR1rx/X+v/iWPV8o4Gg5/wmLntmsP0Y3aKFoADbfDq063XQsdo5RwNBX58LkOV/kfl6uUAbQeFyu7alMaGOoZTPqO3CQ5t9zNfLBRoJCsYFPhsdq51zNBB05iCAi14pzNfLBVoJCnC0fs9kx+rnGg0ETW/WrK/vKuar5QPtBIXsNwI0mSLXuWhxFF+wZ+N55mvlBA0FBdhfa8Adx1bBL3gelC2aCgqZY6omOLYObkFB2eKM8UHVJPRgWIx73erpKkFzjpxy13PzD+CM8UFV7fHZU31Wu1OmXSTor0FNarQpBTcosB0ftAiVXdLx5i1PqvoDrqERNGtm7NfCr4GSYlv5rLMdjKNHOL4q3cB0fNAi1H5mMi73npyh7k/4hUbQ2Mi1EfEAHpJiG/m84Sc8nJC+M7gjTh0flJgbsUFb1f4Np9AI6pMGybWu3xfUuH6VSPMQq8icikLn9f1zjq9KNzh1fFAVJDV4wT0GGKMRtM51gA1RhcWCZo8cKlI70Dp0dLvvNgRtd3xVuoFK0MKt2yA+Zna25XWGh8GEtwMNyf/Qe/w9B/6ON2gEXVVjkdDPty8nKba1w+cve7HHbsfXpB+oBB0T0apbq8+6SIevdPC83c0hfqsLHPpLnqA6ij+SIPTrW8dKSvE8qBxKglbNyKt0GlIDJMUOJ/Ro24a6f1fAE/VsoRI0MCXjsd/gVnVJMUVCE2q/eMLhP+YCFJQtVIK+6199QtiUxtMkxTQJzVvs8+plx/9ce1BQttAdxf95HPbHWV3XSZfQe296TtbxtY0oKFu0vVjENleHeU/5dItOD+lRULbwKCjAxkcfbeCvz+8/UVC28Clo/R+Od6sclkNbjRagoCpJvmD3WiEuBc19VGjzj+UCl2Qrx/IGCqqKrO5Vqte311dyKSgE/gXwQ7MjUQHvZ9JW5WpQUFVMjs2FdU/aCeBT0HVBC+b47QT4M8Z3ns4OllBQVTT5XXjwuy4fwKegcGBC3HHTk78HeMbZaT5/oKCqeG4XQE4lOxdbcipoCc6PqjJMRyMzo6Cq+CJ0528v9bcTwL+gALdm+Lykm3FyUFB1fNEuYqa9yYn0IChA5rJazTdLZ6nmExSULfoQFMCY0CZooR5uokdB2aIXQQV+7195KP+XOqGgbNGRoAA3ZlVt9xXnPT0KyhZdCQqQt6lNwEzp/PRcgYKyRWeCCpwYUaVbonTGan5AQdmiP0EB0ldHBM/l9W0UBWWLHgUV+H14lS7buJxnCQVli04FBcj8tLXvGxxeMoqCskW3ggr8MyWw6TLeRmlGQdmiZ0EBCnb1eyI6gasLm1FQtuhbUIHUde08hyTxc1RPIei9hsFmLAW5MyaL1Jfe112a0L2gApcWNKw20WoYXY2geQdNDUu4ImJ5nffeOyLhQQzapVfcQVCBk1NCQmf87fLV2oCqi593xFYpdvFy6EdQgcPjqz05R/sR8vAzKFu0HqOeJYUHxgQ8OVPjU08oKFs0H6OeLcYDY6uHTflNw0HvUVC2cDBGPWMKf40Lrf76/2l10RMKyhYexqhnSvr7I1fm/jW/qWf/LWlarB8FZQsfY9SzI6Puy8uini0AuLLiv4+/sPS8yxugJGi0BfIqUVA5igQ1zu2UusNmgGSM+sKk3SIdajFsoEpW9Raa0cK8Uelfvezz5OSfXDtqs5KgSRbIq0RB5SgSdFzv+tmRo20F7MxOnx69MNfyMrNzBxH/qkybqIo33xYehq+wvDT+Mi28Sq/PbriuAQRdfO6+xMTt1jN3yIKCylEkaM3MtpAtHedb5JXnc3p1i+8knU9Ky4Rubl0AmbUPlyy6tjamcpMpP7noqIlA0J7NPLsHTiGvEgWVo0jQ4Ntt4ZatbtsvB3wzIK2apFjLhBq71x9QY5y0NH//1Ignopf/64IGEAjqnT/tyLko8ipRUDmKBF3cKHh+2Coby58+AeHn4ExtSbG2CT342Z82y29//nJAjSGbnX15HoGg1XPiV4CKWeJQUDksR/F7py/4w9byo8E9o/wHBHwpKeY3oSc/7Px4+MREZ867SCDo5A4X647qRF4lv/l0ASSCrjFhKyB7y8IZK6zmPOA6ofkHZrXxaPXWXmcNPUpyHvQyHHxbxUwRXOfT2ZAIGhcXN9zfnT4zZf7wZnOPdjOTnCEpgaDyO7xtuM+nMyE9UX9XeiRkB10kNH3H5OYebd7azbq7JxDU7XZ4Z0IqaL47nrdL3zWttUfTcdtuWQqMZ6nPmBJ+1en4Dn/py92cD67CFBJBewvUH0hep24EFcn+ad6LlcJe/US8kvSv+jW8uqTT1UcoqMM7fLxvzxaN7qprk54hETRRICnXTpwEXQkqUvDnsr5B3lELaq+EvFelU7mqhEBQmh0+y/N/AK9PUN8uvaIoaJwF8jp1J6iJK5uHP1S+6esL6tBVQyAozQ7/x1PCw9526tulVxQFFY43Bwe/Pr3OVPI69SkoQE7FlAPvtnnEp8uc3Y7P3KAkKOUOn1IlDWDxK461TY+QdPFPpQCkhpPXqVdBYUjng9+Fxl/ZOvHZimEDlh5y6CyUkqC0O/yYpium+HBxe6BrIBE0MEs4lggir1O3gmbPa9HB/LVYwfG1w8PLhw/7+KjaAaAIuniqHb5w8/CpF1S2Sc+QCDr6mZWrW44hr1O3gkrIPrTs5QYVmg3/+HcVHxgJBC01OzwLSAQt+LhPv7UqLvt1q4RmHlg8oMFj4YM/+tneKf3cS5b8EAgqv8PzdZesau6MDX+B9R0vioLGnIcYE+R16iahxGT9snxIk/J1e7+zy/aJ/HefqOa31fyUQFDZHZ6/u2RVUdhu2O9bA/eyrVRR0D3pgLcomMg79tn456r4d5y88YTkq5wf61yFP/xMdxAqCmpvh5feJZsV20MkqILwsBng9x58/94utvO1/mzrbUT6Xbz9tD+Amwpq5tK3c3uEPhY+8N0frhWXme4zid1geq4kqL0dXnqXrPHrzSJt/ISHfwCSN/P9+8cQ4ffKKLb1tiYQ9JtBxqYV8eqbEmT++vGY57y8nhu96udU4eX8ScJDp+2mRaTfxdsq5PAuWTXkBO6E9A7L2FZKcpDU/JfEAbdbkdepl4TScn33+4OaVgh6ceICz9V/vVPTfBs+gaDyO7zkLtkidJPPAyE1Kr3GeCRMEkHrFI79It+fvE7dJJQFhWe/ead/WNkK/c+ZXxMI6r47vPFcKusqSQR9pUdAyoTW5HXqJ6HMMBafgyIQFHd4FZAImrbiJMxWMe0LJtQ+uMOrgG5kEdtgQu2DO7wK6EYWsQ0m1D64w6uAamQRGTCh9sEdXgVUI4vIgAm1D+7wKqAaWUQGTKh9cIdXAdXIInq/+sYZEAiKO7wKlAU9/0sWGG+esfUNls6vvnEKJF91yu7wtsF8yiEK+kGZ6t7f+z5e7Vkby7kco15jFAW1s8PLgPmUQxS06q+wz7De9hesPI5RrzVKgtrb4WXAfMohClpF+KkscwGAzq++cQpKgtrb4WXAfMohCuol/PjJBUiuvsmZap78VMWwLm6HkqD2dngZUFA5TO+gycnJPsKP7MCvkfef5i81T34aLBdbClAS1O4ObxsUVA5R0EctWC++HStSITZWUo4JtYPiDm8N5lMOhXmS8id5Ldm0yXvTJkk5JtQOdnZ4OTCfcihOJru36bdg3Z9jQtmC+ZRDeTruu32GWF95iwllC+ZTDpL54uOt5/TDhLIF8ykHiaA2wISyBfMpBwqqHt4Ezf6g35vXlMO4BQVlDIWg9xoGm5GU0+TT2D4mflLAdccr0BoUlDE076CpYQlXRCTFNPn8ubERYOxcxyvQGhSUMVRd/LwjD7zMaNFExFPlN08l2dJDeFg5jKJNGoOCMobpZ9DjR0QipZOhquCs/3XI77COWYtcDgrKGFpBz1gXUeXzA9+omj0YD0fjSlBQxtAKGmFdRJfPS9uP0fy51qCgjOFOUJ2DgjKGVtDV1kWYTzlQUPXwdqJe76CgjEFB2YKCMgYFZQsKyhgUlC0oKGNQULagoIxBQdmCgjIGBWULCsoYFJQtKChjUFC2oKCMQUHZgoIyBgVlCwrKGBSULSgoY1BQtqCgjEFB2UInKI5RbwUKyhYqQXGMemtQULZQCYpj1FuDgrKFSlAco94aFJQtVILiGPXW0AiaNTP2a+HXQEkx5lMOxYMkyRj1GfVCRCr6smiZTqERNDZybUQ8gIekGAWVQ0lQqz3+0lmR6DoMGqZXaAT1SYPkWtdR0JJQCYp7vDU0gta5DrAhqhDzWQIqQXGPt4ZG0FU1Fgl7fftyltcZ/pVFynmzaJhOoTuKxz3eCqqj+CMJAMatY4tfp90R6YX5lEFJUOkeXwQKSgHrsZn0Dt1RvHSPN4MJpQCHvnkQvFiEMSgoW1BQxuDYTGxBQRmDX3WyBQVlDArKFhSUMSgoW1BQxqCgbEFBGYOCsgUFZQwKyhYUlDHOFPTUt2eZV847KChjnCjoa9U7+09nXjvnoKCMcZ6g34dnw72avzKvnm9QUMY4T9AZs4SHUUuZV2+LS5/vzHPJihRBQRnjPEFX9xUent/GvHobfOHbp02928pxLgAFZYzzBE2tOebzV5/KYV69NbmepwAmjXbBmpRBQRnjxIOk5Lf6zE9jXrsN/q4nPBx8xhWrUgQFZYw7nAdNr3QH4KN+Ll6rbVBQxriDoBDXeOkUnxOuXqtNUFDGuIWgsH3UtPMuX6lNUFDGuIeg/ICCMoZK0MSEnJ3Dl+RLSjGfcqCg6qER9I3Q7l0avdd+lKQY8ykHCqoeGkH978GLX8G9apJizKccKKh6aAQNzIKvrsG9YEkx5lMOFFQ9NIJOjNgJ8HO7ccUFqeaRRUIZtEuvoKCMoTpI2nUQIGGp0fIyI8A0NtMj0j6/NIGCMsYJQ9+sHkJXpa7BWT4Y44SRRVBQOXCWD/WgoGxhO8vHZRxh2QlD36Cgcqid5SMjzDRG/ePVWbRMpzjhq04UVA4HZ/nAhLIF8ymH2lk+isCEsgXzKYeDp5kwoWzBfMqBgqoHBWULCsoYFJQtThF0wJ0irp4rLdywbPJzThDUvfOZcscedvPpoKDbKlswlCHCDcLKFG8z+6E/VOeTwTaRrIlRNWUr28VePh0UtJjM8kRh6dKJlmyTVpEo7F4lorA7lYnCkj2Jwm65ZK4tkny+tlwx5N9aytVEbVcMSXpWuZqI3xRDNvZVrkYOFBQFlQcFlYCCoqASUFAUVB4UVAIKioJKQEFRUHncQFDjRKKwgkksw/InE4XlxRGF5b7JMowSknxuPqQYkjZTuZqV/yiGXP5AuZoFNxVDjn+qXI0ctIIiiFNBQRGuQUERrkFBEa5BQRGuQUERrkFBEa5BQRGu4UnQfJdP23OPaRg3uCqRBHmhTh2loNtCQxYrxawO9Y3NIomcFEtQ4dGn/QcRhL1bLWBsoVJYakOwVGUv0hRGvBVUKNZP1AzlRBJkUTmDBOmjTx2doCmBF1NDT9qPORJ8PTVyJkFkklescoXZwcdznt+pGPanf3JGxFaFsHlhwZZNsBdpCiPeCioU6ydqhnIiCbKonEGC9DFIHZ2g64VddeYM+zEJiwDWxipH3m26JFa5wi09AXKyFcOOV8/Ma7VNIexIQrBlE+xFmsJIt4IOxfpJmkGQSIIsKmeQIH0MUkcn6LypAJ8MVgxLabpNObLPN5tilStc2DXct0+Gcm0THqsUrVjblWDLJtiNvBJMvhV0kNSv2AyCRJJkUTmDBOmjTx2doHPEtQ5UioqvFa8cGT8EhLwqhs3wP5ve7S3FsANhv55qGa8UJqbPHGM30pxloq2ghKB+xWaQJJIgiwQZJEgfferoBF0nrHH2NPsxxv5R1wkio/2DvStEKoatEALioxXD4qYAfBajFCamzxxjN1IMI9wKShTrJ2gGSSIJskiQQYL00aeOTtBbAbey6/9uPyahfSFhpLDjK4ZdDjib3nmhYtjGxrfSomYrhYnpM8fYjRTDyLeCBsX6yZqhmEiCLBJkkCB99KmjPM20uVF9pWtax5Xz8PAYSBIp5FU57Mtaga/lKoYVTqtWdUi2UpipAzLH2IsUw8i3ggql+smaoZxI5SwSZJAgffSp4+lEPYJYgYIiXIOCIlyDgiJcg4IiXIOCIlyDgiJcg4IiXIOCIlyDgiJcg4IiXIOCIlyDgiJcg4IiXIOCIlyDgiJcg4IiXIOCIlyjH0E7eVUxeHkN3xOjdUPcBJ3kUz+CAiSLU4Wkn9e6GW6DLvKpO0ERZugin7oTNCnmUOtG/sNn1ml5GZbUqDXeqHWr9Isu8qlHQctfSKs4A4YsSmp0K6P3XK1bpV90kU89Cvo8QON/YFnc9OC2bet30rpV+kUX+dSjoP8VEnpBSOiCtwGyU7VulX7RRT71LOgfYbdyIndo3Sr9oot86llQWFkzYLTWjdIxusinngRFSiEoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1KCjCNSgowjUoKMI1/w8EAGl9XZFhIAAAAABJRU5ErkJggg==" /><!-- --></p> <p>The <span class="math inline"><em>χ</em><sup>2</sup></span> error level of 3.3% as well as the plot suggest that the SFO model fits very well. The error level at which the <span class="math inline"><em>χ</em><sup>2</sup></span> test passes is slightly lower for the FOMC model. However, the difference appears negligible.</p> <pre class="r"><code>summary(mm.L4[["SFO", 1]], data = FALSE)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:36 2022 -## Date of summary: Wed Sep 14 22:28:37 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:47 2022 +## Date of summary: Tue Dec 6 09:39:47 2022 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## -## Fitted using 142 model solutions performed in 0.034 s +## Fitted using 142 model solutions performed in 0.03 s ## ## Error model: Constant variance ## @@ -2092,10 +2085,10 @@ plot(mm.L4)</code></pre> ## DT50 DT90 ## parent 106 352</code></pre> <pre class="r"><code>summary(mm.L4[["FOMC", 1]], data = FALSE)</code></pre> -<pre><code>## mkin version used for fitting: 1.1.2 -## R version used for fitting: 4.2.1 -## Date of fit: Wed Sep 14 22:28:37 2022 -## Date of summary: Wed Sep 14 22:28:37 2022 +<pre><code>## mkin version used for fitting: 1.2.2 +## R version used for fitting: 4.2.2 +## Date of fit: Tue Dec 6 09:39:47 2022 +## Date of summary: Tue Dec 6 09:39:47 2022 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent |
