diff options
| -rw-r--r-- | DESCRIPTION | 5 | ||||
| -rw-r--r-- | NAMESPACE | 8 | ||||
| -rw-r--r-- | NEWS.md | 2 | ||||
| -rw-r--r-- | R/aw.R | 36 | ||||
| -rw-r--r-- | R/multistart.R | 59 | ||||
| -rw-r--r-- | R/saem.R | 8 | ||||
| -rw-r--r-- | log/test.log | 28 | ||||
| -rw-r--r-- | man/aw.Rd | 6 | ||||
| -rw-r--r-- | man/multistart.Rd | 21 | ||||
| -rw-r--r-- | man/saem.Rd | 2 | ||||
| -rw-r--r-- | tests/testthat/print_sfo_saem_1.txt | 6 | ||||
| -rw-r--r-- | tests/testthat/summary_hfit_sfo_tc.txt | 2 | ||||
| -rw-r--r-- | vignettes/FOCUS_L.html | 151 |
13 files changed, 223 insertions, 111 deletions
diff --git a/DESCRIPTION b/DESCRIPTION index 6734a43c..c653be1c 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -2,7 +2,7 @@ Package: mkin Type: Package Title: Kinetic Evaluation of Chemical Degradation Data Version: 1.1.2 -Date: 2022-09-14 +Date: 2022-09-16 Authors@R: c( person("Johannes", "Ranke", role = c("aut", "cre", "cph"), email = "johannes.ranke@jrwb.de", @@ -24,7 +24,8 @@ Description: Calculation routines based on the FOCUS Kinetics Report (2006, purpose. Depends: R (>= 2.15.1), Imports: stats, graphics, methods, parallel, deSolve, R6, inline (>= 0.3.19), - numDeriv, lmtest, pkgbuild, nlme (>= 3.1-151), saemix (>= 3.1), rlang, vctrs + numDeriv, lmtest, pkgbuild, nlme (>= 3.1-151), saemix (>= 3.1), rlang, vctrs, + KernSmooth Suggests: knitr, rbenchmark, tikzDevice, testthat, rmarkdown, covr, vdiffr, benchmarkme, tibble, stats4 License: GPL @@ -6,8 +6,10 @@ S3method(AIC,mhmkin) S3method(AIC,mmkin) S3method(BIC,mhmkin) S3method(BIC,mmkin) +S3method(aw,mixed.mmkin) S3method(aw,mkinfit) S3method(aw,mmkin) +S3method(aw,multistart) S3method(confint,mkinfit) S3method(convergence,mhmkin) S3method(convergence,mmkin) @@ -32,6 +34,7 @@ S3method(nlme,mmkin) S3method(nobs,mkinfit) S3method(parms,mkinfit) S3method(parms,mmkin) +S3method(parms,multistart) S3method(parms,saem.mmkin) S3method(plot,mixed.mmkin) S3method(plot,mkinfit) @@ -59,7 +62,6 @@ S3method(residuals,mkinfit) S3method(saem,mmkin) S3method(summary,mkinfit) S3method(summary,mmkin) -S3method(summary,multistart.saem.mmkin) S3method(summary,nlme.mmkin) S3method(summary,saem.mmkin) S3method(update,mkinfit) @@ -85,6 +87,7 @@ export(illparms) export(ilr) export(intervals) export(invilr) +export(llhist) export(loftest) export(logistic.solution) export(lrtest) @@ -115,6 +118,7 @@ export(nafta) export(nlme) export(nlme_data) export(nlme_function) +export(parhist) export(parms) export(plot_err) export(plot_res) @@ -127,6 +131,7 @@ export(transform_odeparms) import(deSolve) import(graphics) import(nlme) +importFrom(KernSmooth,bkde) importFrom(R6,R6Class) importFrom(grDevices,dev.cur) importFrom(lmtest,lrtest) @@ -147,6 +152,7 @@ importFrom(stats,dist) importFrom(stats,dnorm) importFrom(stats,lm) importFrom(stats,logLik) +importFrom(stats,median) importFrom(stats,na.fail) importFrom(stats,nlminb) importFrom(stats,nobs) @@ -1,6 +1,6 @@ # mkin 1.1.2 -- 'R/multistart.R': New method for testing multiple start parameters for hierarchical model fits. +- 'R/multistart.R': New method for testing multiple start parameters for hierarchical model fits, with diagnostic plotting functions 'llhist' and 'parhist'. - 'R/mhmkin.R': New method for performing multiple hierarchical mkin fits in one function call, optionally in parallel. @@ -30,6 +30,14 @@ #' @export aw <- function(object, ...) UseMethod("aw") +.aw <- function(all_objects) { + AIC_all <- sapply(all_objects, AIC) + delta_i <- AIC_all - min(AIC_all) + denom <- sum(exp(-delta_i/2)) + w_i <- exp(-delta_i/2) / denom + return(w_i) +} + #' @export #' @rdname aw aw.mkinfit <- function(object, ...) { @@ -43,11 +51,7 @@ aw.mkinfit <- function(object, ...) { } } all_objects <- list(object, ...) - AIC_all <- sapply(all_objects, AIC) - delta_i <- AIC_all - min(AIC_all) - denom <- sum(exp(-delta_i/2)) - w_i <- exp(-delta_i/2) / denom - return(w_i) + .aw(all_objects) } #' @export @@ -56,3 +60,25 @@ aw.mmkin <- function(object, ...) { if (ncol(object) > 1) stop("Please supply an mmkin column object") do.call(aw, object) } + +#' @export +#' @rdname aw +aw.mixed.mmkin <- function(object, ...) { + oo <- list(...) + data_object <- object$data[c("ds", "name", "time", "value")] + for (i in seq_along(oo)) { + if (!inherits(oo[[i]], "mixed.mmkin")) stop("Please supply objects inheriting from mixed.mmkin") + data_other_object <- oo[[i]]$data[c("ds", "name", "time", "value")] + if (!identical(data_object, data_other_object)) { + stop("It seems that the mixed.mmkin objects have not all been fitted to the same data") + } + } + all_objects <- list(object, ...) + .aw(all_objects) +} + +#' @export +#' @rdname aw +aw.multistart <- function(object, ...) { + do.call(aw, object) +} diff --git a/R/multistart.R b/R/multistart.R index a3afa08b..94292e82 100644 --- a/R/multistart.R +++ b/R/multistart.R @@ -14,7 +14,12 @@ #' @param n How many different combinations of starting parameters should be #' used? #' @param cores How many fits should be run in parallel? -#' @param \dots Passed to the update function. +#' @param \dots Passed to the update function, or to the basic plotting +#' function in the case of the graphical function. +#' @param x The multistart object to print +#' @param breaks Passed to [hist] +#' @param main title of the plot +#' @param lpos Positioning of the legend. #' @return A list of [saem.mmkin] objects, with class attributes #' 'multistart.saem.mmkin' and 'multistart'. #' @@ -39,6 +44,7 @@ multistart.saem.mmkin <- function(object, n = 50, cores = 1, ...) { res <- parallel::mclapply(1:n, function(x) { update(object, degparms_start = start_parms[x, ], ...) }, mc.cores = cores) + attr(res, "orig") <- object attr(res, "start_parms") <- start_parms class(res) <- c("multistart.saem.mmkin", "multistart") return(res) @@ -53,8 +59,53 @@ print.multistart <- function(x, ...) { #' @rdname multistart #' @export -summary.multistart.saem.mmkin <- function(object) { +parms.multistart <- function(object, ...) { + t(sapply(object, parms)) +} + +#' @rdname multistart +#' @importFrom stats median +#' @export +parhist <- function(object, lpos = "topleft", ...) { + orig <- attr(object, "orig") + orig_parms <- parms(orig) + start_parms <- orig$mean_dp_start + all_parms <- parms(object) + median_parms <- apply(all_parms, 2, median) + all_scaled_parms <- t(apply(all_parms, 1, function(x) x / median_parms)) + orig_scaled_parms <- orig_parms / median_parms + start_scaled_parms <- rep(NA_real_, length(orig_parms)) + names(start_scaled_parms) <- names(orig_parms) + start_scaled_parms[names(start_parms)] <- + start_parms / median_parms[names(start_parms)] + + boxplot(all_scaled_parms, log = "y", ...) + points(orig_scaled_parms, col = 2, cex = 2) + points(start_scaled_parms, col = 3, cex = 3) + legend(lpos, inset = c(0.05, 0.05), bty = "n", + pch = 1, col = 3:1, lty = c(NA, NA, 1), + legend = c( + "Starting parameters", + "Converged parameters", + "Multistart runs")) +} + +#' @rdname multistart +#' @importFrom KernSmooth bkde +#' @export +llhist <- function(object, breaks = "Sturges", main = "", lpos = "topleft", ...) { + ll <- sapply(object, logLik) + kde <- KernSmooth::bkde(ll) + h <- hist(ll, freq = TRUE, + xlim = range(kde$x), + xlab = "", main = main, + ylab = "Frequency of log likelihoods", breaks = breaks, ...) - parm_matrix <- sapply(object, parms) - parm_matrix + freq_factor <- h$counts[1] / h$density[1] + lines(kde$x, freq_factor * kde$y) + abline(v = logLik(attr(object, "orig")), col = 2) + legend(lpos, inset = c(0.05, 0.05), bty = "n", + lty = 1, col = c(2, 1), + legend = c("original log likelihood", + "kernel density estimate")) } @@ -568,7 +568,7 @@ saemix_data <- function(object, verbose = FALSE, ...) { ds_list <- lapply(object, function(x) x$data[c("time", "variable", "observed")]) names(ds_list) <- ds_names - ds_saemix_all <- purrr::map_dfr(ds_list, function(x) x, .id = "ds") + ds_saemix_all <- vctrs::vec_rbind(!!!ds_list, .names_to = "ds") ds_saemix <- data.frame(ds = ds_saemix_all$ds, name = as.character(ds_saemix_all$variable), time = ds_saemix_all$time, @@ -617,9 +617,9 @@ update.saem.mmkin <- function(object, ..., evaluate = TRUE) { #' @rdname saem #' @param ci Should a matrix with estimates and confidence interval boundaries #' be returned? If FALSE (default), a vector of estimates is returned. -parms.saem.mmkin <- function(x, ci = FALSE, ...) { - conf.int <- x$so@results@conf.int[c("estimate", "lower", "upper")] - rownames(conf.int) <- x$so@results@conf.int[["name"]] +parms.saem.mmkin <- function(object, ci = FALSE, ...) { + conf.int <- object$so@results@conf.int[c("estimate", "lower", "upper")] + rownames(conf.int) <- object$so@results@conf.int[["name"]] conf.int.var <- grepl("^Var\\.", rownames(conf.int)) conf.int <- conf.int[!conf.int.var, ] estimate <- conf.int[, "estimate"] diff --git a/log/test.log b/log/test.log index 8be4a512..cb123776 100644 --- a/log/test.log +++ b/log/test.log @@ -1,23 +1,23 @@ ℹ Testing mkin ✔ | F W S OK | Context ✔ | 5 | AIC calculation -✔ | 5 | Analytical solutions for coupled models [3.4s] +✔ | 5 | Analytical solutions for coupled models [3.3s] ✔ | 5 | Calculation of Akaike weights ✔ | 3 | Export dataset for reading into CAKE -✔ | 12 | Confidence intervals and p-values [1.0s] -✔ | 1 12 | Dimethenamid data from 2018 [32.4s] +✔ | 12 | Confidence intervals and p-values [1.1s] +✔ | 1 12 | Dimethenamid data from 2018 [32.5s] ──────────────────────────────────────────────────────────────────────────────── Skip (test_dmta.R:98:3): 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 [5.0s] +✔ | 14 | Error model fitting [5.1s] ✔ | 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.4s] ✔ | 1 | Fitting the logistic model [0.2s] -✔ | 5 | Batch fitting and diagnosing hierarchical kinetic models [14.5s] -✔ | 1 12 | Nonlinear mixed-effects models [0.2s] +✔ | 7 | Batch fitting and diagnosing hierarchical kinetic models [14.6s] +✔ | 1 12 | Nonlinear mixed-effects models [0.3s] ──────────────────────────────────────────────────────────────────────────────── Skip (test_mixed.R:68:3): saemix results are reproducible for biphasic fits Reason: Fitting with saemix takes around 10 minutes when using deSolve @@ -26,25 +26,25 @@ Reason: Fitting with saemix takes around 10 minutes when using deSolve ✔ | 10 | Special cases of mkinfit calls [0.4s] ✔ | 3 | mkinfit features [0.7s] ✔ | 8 | mkinmod model generation and printing [0.2s] -✔ | 3 | Model predictions with mkinpredict [0.4s] +✔ | 3 | Model predictions with mkinpredict [0.3s] ✔ | 16 | Evaluations according to 2015 NAFTA guidance [1.8s] -✔ | 9 | Nonlinear mixed-effects models with nlme [8.7s] -✔ | 16 | Plotting [10.1s] +✔ | 9 | Nonlinear mixed-effects models with nlme [8.5s] +✔ | 16 | Plotting [10.0s] ✔ | 4 | Residuals extracted from mkinfit models -✔ | 28 | saemix parent models [181.2s] +✔ | 28 | saemix parent models [183.3s] ✔ | 2 | Complex test case from Schaefer et al. (2007) Piacenza paper [1.4s] ✔ | 7 | 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.1s] -✔ | 9 | Hypothesis tests [7.8s] -✔ | 4 | Calculation of maximum time weighted average concentrations (TWAs) [2.1s] +✔ | 9 | Hypothesis tests [8.1s] +✔ | 4 | Calculation of maximum time weighted average concentrations (TWAs) [2.2s] ══ Results ═════════════════════════════════════════════════════════════════════ -Duration: 279.6 s +Duration: 282.2 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) -[ FAIL 0 | WARN 0 | SKIP 2 | PASS 233 ] +[ FAIL 0 | WARN 0 | SKIP 2 | PASS 235 ] @@ -4,6 +4,8 @@ \alias{aw} \alias{aw.mkinfit} \alias{aw.mmkin} +\alias{aw.mixed.mmkin} +\alias{aw.multistart} \title{Calculate Akaike weights for model averaging} \usage{ aw(object, ...) @@ -11,6 +13,10 @@ aw(object, ...) \method{aw}{mkinfit}(object, ...) \method{aw}{mmkin}(object, ...) + +\method{aw}{mixed.mmkin}(object, ...) + +\method{aw}{multistart}(object, ...) } \arguments{ \item{object}{An \link{mmkin} column object, containing two or more diff --git a/man/multistart.Rd b/man/multistart.Rd index 347b12fb..4dd95503 100644 --- a/man/multistart.Rd +++ b/man/multistart.Rd @@ -4,7 +4,9 @@ \alias{multistart} \alias{multistart.saem.mmkin} \alias{print.multistart} -\alias{summary.multistart.saem.mmkin} +\alias{parms.multistart} +\alias{parhist} +\alias{llhist} \title{Perform a hierarchical model fit with multiple starting values} \usage{ multistart(object, n = 50, cores = 1, ...) @@ -13,7 +15,11 @@ multistart(object, n = 50, cores = 1, ...) \method{print}{multistart}(x, ...) -\method{summary}{multistart.saem.mmkin}(object) +\method{parms}{multistart}(object, ...) + +parhist(object, lpos = "topleft", ...) + +llhist(object, breaks = "Sturges", main = "", lpos = "topleft", ...) } \arguments{ \item{object}{The fit object to work with} @@ -23,7 +29,16 @@ used?} \item{cores}{How many fits should be run in parallel?} -\item{\dots}{Passed to the update function.} +\item{\dots}{Passed to the update function, or to the basic plotting +function in the case of the graphical function.} + +\item{x}{The multistart object to print} + +\item{lpos}{Positioning of the legend.} + +\item{breaks}{Passed to \link{hist}} + +\item{main}{title of the plot} } \value{ A list of \link{saem.mmkin} objects, with class attributes diff --git a/man/saem.Rd b/man/saem.Rd index dfb652b6..dd787475 100644 --- a/man/saem.Rd +++ b/man/saem.Rd @@ -42,7 +42,7 @@ saemix_model( saemix_data(object, verbose = FALSE, ...) -\method{parms}{saem.mmkin}(x, ci = FALSE, ...) +\method{parms}{saem.mmkin}(object, ci = FALSE, ...) } \arguments{ \item{object}{An \link{mmkin} row object containing several fits of the same diff --git a/tests/testthat/print_sfo_saem_1.txt b/tests/testthat/print_sfo_saem_1.txt index a11e1e96..34eea058 100644 --- a/tests/testthat/print_sfo_saem_1.txt +++ b/tests/testthat/print_sfo_saem_1.txt @@ -7,13 +7,13 @@ Data: Likelihood computed by importance sampling AIC BIC logLik - 1311 1315 -649 + 1312 1316 -650 Fitted parameters: estimate lower upper -parent_0 1e+02 99.13 1e+02 +parent_0 1e+02 99.28 1e+02 k_parent 4e-02 0.03 4e-02 a.1 9e-01 0.75 1e+00 b.1 5e-02 0.04 5e-02 -SD.parent_0 7e-01 -1.09 3e+00 +SD.parent_0 2e-01 -6.33 7e+00 SD.k_parent 3e-01 0.20 4e-01 diff --git a/tests/testthat/summary_hfit_sfo_tc.txt b/tests/testthat/summary_hfit_sfo_tc.txt index ad701cee..e3e2f7e4 100644 --- a/tests/testthat/summary_hfit_sfo_tc.txt +++ b/tests/testthat/summary_hfit_sfo_tc.txt @@ -45,7 +45,7 @@ SD.log_k_parent 0.3 0.1 0.4 Variance model: est. lower upper -a.1 0.91 0.64 1.17 +a.1 0.90 0.64 1.17 b.1 0.05 0.04 0.06 Backtransformed parameters: diff --git a/vignettes/FOCUS_L.html b/vignettes/FOCUS_L.html index 13d0d9eb..da6c11fe 100644 --- a/vignettes/FOCUS_L.html +++ b/vignettes/FOCUS_L.html @@ -1517,7 +1517,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-07-08)</h4> +<h4 class="date">Last change 18 May 2022 (rebuilt 2022-09-14)</h4> </div> @@ -1536,17 +1536,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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:39 2022 -## Date of summary: Fri Jul 8 15:44:39 2022 +## Date of fit: Wed Sep 14 22:28:35 2022 +## Date of summary: Wed Sep 14 22:28:35 2022 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## -## Fitted using 133 model solutions performed in 0.029 s +## Fitted using 133 model solutions performed in 0.032 s ## ## Error model: Constant variance ## @@ -1578,9 +1578,9 @@ summary(m.L1.SFO)</code></pre> ## ## Parameter correlation: ## parent_0 log_k_parent sigma -## parent_0 1.000e+00 6.186e-01 -1.712e-09 -## log_k_parent 6.186e-01 1.000e+00 -3.237e-09 -## sigma -1.712e-09 -3.237e-09 1.000e+00 +## parent_0 1.000e+00 6.186e-01 -1.516e-09 +## log_k_parent 6.186e-01 1.000e+00 -3.124e-09 +## sigma -1.516e-09 -3.124e-09 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. @@ -1627,25 +1627,27 @@ summary(m.L1.SFO)</code></pre> <pre class="r"><code>mkinresplot(m.L1.SFO, ylab = "Observed", xlab = "Time")</code></pre> <p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAkAAAAHgCAMAAAB6sCJ3AAAC/VBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+Pj4/Pz9AQEBBQUFCQkJDQ0NERERFRUVGRkZHR0dISEhJSUlKSkpLS0tMTExNTU1OTk5PT09QUFBRUVFSUlJTU1NUVFRVVVVWVlZXV1dYWFhZWVlaWlpbW1tcXFxdXV1eXl5gYGBhYWFiYmJjY2NkZGRlZWVmZmZnZ2doaGhpaWlqampra2tsbGxtbW1ubm5vb29wcHBxcXFycnJzc3N0dHR1dXV2dnZ3d3d4eHh5eXl6enp7e3t8fHx9fX1+fn5/f3+AgICBgYGCgoKDg4OEhISFhYWGhoaHh4eIiIiJiYmKioqLi4uMjIyNjY2Ojo6Pj4+QkJCRkZGSkpKTk5OUlJSVlZWWlpaXl5eYmJiZmZmampqbm5ucnJydnZ2enp6fn5+goKChoaGioqKjo6OkpKSlpaWmpqanp6eoqKipqamqqqqrq6usrKytra2urq6vr6+wsLCxsbGysrKzs7O0tLS1tbW2tra3t7e4uLi5ubm6urq7u7u8vLy9vb2+vr6/v7/AwMDBwcHCwsLDw8PExMTFxcXGxsbHx8fIyMjJycnKysrLy8vMzMzNzc3Ozs7Pz8/Q0NDR0dHS0tLT09PU1NTV1dXW1tbX19fY2NjZ2dna2trb29vc3Nzd3d3e3t7f39/g4ODh4eHi4uLj4+Pk5OTl5eXm5ubn5+fo6Ojp6enq6urr6+vs7Ozt7e3u7u7v7+/w8PDx8fHy8vLz8/P09PT19fX29vb39/f4+Pj5+fn6+vr7+/v8/Pz9/f3+/v7///+7nL8pAAAACXBIWXMAAA7DAAAOwwHHb6hkAAAW6UlEQVR4nO3deXjU5KLH8ZTdspRSabEsSgHZKoIglMWjCAhVuCi1goKgQhWUTcoiCEeFqyA9IFU5RRaRCrIUZClFqSLniCxaEUGFoy1y8YhQ2Xe6TJ6bZFoYsNOZvL9kMk1/nz+StjNv3jzM96GdTCYjyUQAyeodoJKNARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARGEARHE4IBebE220OmoNQFFvZ9BdtDoe4sC2mHs9sgiLRgQIRgQQRgQQRgQQRgQQRgQQRgQQRgQQRgQQRgQQRgQQRgQQRiQTX0UZ7RfipyHAdlUbNw8Y7VeWeQ8DMimYot+vA3fIAOyKQZEEAZEEAZEEAZEEAZEEAZEEAZEEAZEEAZEEAZEEAZEEKMC+u6b4jfIgGzKqIAGPFz8BhmQTbk+3kcH1w2MWn/DHRx5fxlzJcf1u1xtyYBKKZfH+0jdSYcvbGqS4HJrprSpXkCTGcpXZ0c2rBQx1SHLtVbEBx+SF7cObL5I+XGdd7tKQbEn5ChJkk7fuEFXDMimXB7vERPU5W/Bp6/dmindNGRtfMDfZblPcML6kdIyJaA2vVddSiw/JW1UwFwloOD+W98s97x8IqbHUceNG3TFgGzK5fFumaGtumy+dmumFKMsx1Y9K8ckKV80HqcE1MIhnw+Zqnw3JFQJqKWSTa+O/BVWark83pF7tVX0xmu3ZkopynKftEtZnvtuQfl4JaDxsrxL2n78+PFk6bBc5yXllvgoBlRquTzeA/+hLs+E/HHt1kxpm7I8pWS0rUVAeHS4GtAsWV4hOe2V68yUGVCp5vJ4/xi6zCEf6vqCy63O/4F+kr46WX7YMeXRUAN6S5a3SNkFd6ij/snNgEox18c7o0P1iOA3cl1uzZQeU5YTA09tljJl+dItBQFlV1yo/HhyFwZE1z/epw7mX3drplRp2MaXykyUD5WP/XJdu8oPHtcCksdXmJY2NmCOS0BPNc/IK2KDVzEgmyr+SHSm9PFD1W9/XXmmtbxx5bYbPgie4gzIkRAZ2Ex9XnY1oC/qVz1T3AYZkE15CijDoA0yIJtiQARhQAQpPqDcQ1cM2iADsimeUEYQBkQQBkQQBkQQBkQQBkQQBkQQBkSQ2OkGfzhqVwZUqkw3+uOZ2xT94gcDIggDIggDIggDIggDIggDIggDIggDIggDIggDIoiZAeVnHnV7GwOyCXMCGrxVWcysJkn1Vrm5BwOyCXMCkpJkOUkauGZ9XJlNRd+DAdmEeQHdEad+9UJH15+fWFmo3hY92yO/ZV5AgevUr1Kruf58T2yhcrP1bI/8lnkBtUpUv3qtadH3qDJfz/bIb5kUUGiXuAdCsmTH0tBRRd+DAdmEOQGtSRjW4/YKy+UMqd3Fou/BgGzCxONA+VfkY+n5bm5kQDZh1ZFoBmQTDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDEjc8UMOq3fBegxI1NY7Q26t+XapT4gBCdp2S6osH4h61er9sBoDEtQjWV0eCTpv9Y5YjAEJCnVeOKL1Nxbvh9UYkKDww9oqco/F+2E1BiSon/rB6vIPYTlW74jFGJCgn8NGvfHaK7cusXo/rMaABOU/FRhYoWoL95fQKiUYkKA5fzsry46J0Vbvh9UYkKCWX6nL3JpHrN4RizEgQTVOaKtSf6U1BiSoifb83VH7oNU7YjEGJGjKo+qFR+a3s3o/rMaABF2Mbj17wWP191u9H1YrUQFdntW728Q/jd0RcetHPjPXzeWzShG/DyjtuUenFDzTyW7Wd+Wa0WFFf3YnWcPPA8rrd8fcVWNC07RvhvSPqly94ciWxu4JQfw8oPn3qq817bjlgvpNjdCUfHlH86DDxu4KIfw8oOi12ur+T9RlmQ/U5S/l+DvMj/h5QG13aasB2tlbZT5Vl2fKfG3srhDCzwPqu1Bb3bldXd7U8ndZPv94UKl/6uxP/Dygz247KMuO2S21qwVHx4T2fCz84bA8Y3eFEH4ekJwUGju8ddss7etvw2atXLaw/gdF3vG/F4zaNdLD3wOSj65I/KLwvTP7oqtUvruoz4/KnRFcu8r9ew3aN9LB7wO6Xn7RZ5AO7XZIzl8YlonsEQkpYQEV7fDN59TV/8YZuE3yji0CWtdLW317l4HbJO/oDOhmF9C8hgaU6jyxdOfdBm6TvKMzoKSkpLeq1B41a1LjuqnQvIYGlF1DO7c9foyB2yTv6P8V9ty9l5VlXo9noXmNPR/otRZfXDn6ct1S/xYJC+gPqHaKtlobDs1r8AllK1pVrPXMH4ZukrwiENBb2iqxHjRvcQHte6n/y7/o3WCpv86KRfQHNCRog7JMDTLtV9iM8Fc/nBQ2D9o8+Yr+gM7dJ9WIrCHdj13XxH1AGXWOKctDoTwqWCKIHAfaMn3UzH+D87oP6KVXtNXwBHAG8gmhA4mH98Hzug9osPOWN8fCc5APCAT0SS1JkrskYvO6D2jqi9pq0NzCH/CdD/5Mf0DJ5eKSJXlyAPZXrvuADtZUT1n9IjRb++63J6oFNpjn7sPDyXL6A2o2Wj6ufDM2Epq3mGdh68L7TOhVb4v29R91pp127O4wGpqLTKQ/oMA0LaC0ytC8xR0HOrX8jVUFz/HGaq9OnK75f9BkZB79Ad31dy2g11pA83p5JLrTv7RVzEpoMjKP/oDeLz91u5S9sMIsaF4vA+roPFwQuwKajMwj8CwsMUSSpIoTsNcOvAzoxfHq8lzYIWgyMo/IcaDzu1ZsQS9x4GVAv4cnXJL33/c8OBuZRn9Aw7YbMa+3r8Zn9QmsEf5WrhFTkhkEXo2X6k8+AM/r/ekcuSfhycg8+gNy7IivL7WZDZ58U+IvMEVOYifV757UpGx3aF4GZBNiAZ1dNaBSWWheBmQTAgEdfrd7hfIPvJcNzVtcQLm7Vmbw1a8SQn9AraSKPd+H/7AtJqBdzVvFRrbBzxghX9Af0CNLzxgwr/uAjoZ9rCwX1z1rwCxkOt0B5Ry6ZMS87gOa6XyDcswiI6Yhs+kO6HLwaiPmdR/QkPe01YxxRkxDZtP/K2xMbyPeQeM+oDFvaKuxrxswC5lOf0Ar2tw5ftZbCmhe9wFtaayeC3Syrrc7RpbSH1CtQtC8xTwLG9o8efvCBi9DmydfMe/yLlf+LO65fnHHgdL6tu+/Ve90ZA2T3tZzeGL9AEmq2HCCu09D4pFomzDnbT27A+sOnZO8JPGFiGA3H4vNgGzCnLf13Bdd8Gau3Me7uv48/9esApVfz8r6TfnJ2SyuS/S6mSlv66mWUvjVl9Vdf741olCZGhERDQ7L8rMRXJfodYgpb+tpM6zwq1fcfKQff4XZhDlv60kJiF64ff+BnUv6lE0p+h6+CsjxYa+7+n7pm7lKJZPe1pPaWVIF3J/m5g4+Cij/kQ4ff7uoPo9qm8a0t/Wc/CE9fe8Jtzf7KKAV7dXz8Y/xEuSm8fO39aAGOacZ8k+fzFYaCR6JPrDJ/X8uXvFRQDHOP8HGTffJbKWR/oAOd4uXV5eVbnZzhNBLPgpokvMyVfeu98lspZH+gB4OXyu37H6wW09oXh8FdCg0VZbzZjQv+jNaCKc/oBrvyUekdHlZTWheXz2N39b0zv+5rRvfWm8a/QFVXy4vqnRJ3lAFmtdnBxJzv12Hv4+W3NIf0AMdPrujt3z+oTbQvDwSbRP6A9oTJlXbLd9eYQM0LwOyCYGn8RczlKfwK3/G5mVANiFyHOjEttXfYNep/0tA/4qJ7DaPn8ZcAukPKH9MJUmSqoOH5q4PaHrE4n2bunTjZYBKHv0BvRIwbt+pH18qg11p/LqADtVUP+krv8tCaJNkBf0BNXIe3B1n4HWiFz2prZY+Bm2SxOXM6dri8QyRkQIHEtdpq9Qgkemuui6gxOHaKi0a2iQJu9i+1+bvk2qLPLEReCnjOW31fGeB2a65LqB05zGll+OhTZKwhBh1mXmzwEVXdAa0Z8+ezSH91u5c+3iAu1PFvHNdQHltJuco/6eFunsPEJms6yZt9dA6/UN1BiS50D+Zi+ufhR3pHd6jeTOeeWqVKOeVd/t/qH+ozoAyXeifzMWNBxJ/3fQ9DwNZ5ul/qMu8ht/pH6r7b6C8959u37zXfPTR5pFof7Kn1lblL+nh3QSG6g1ob0spIrpXQ6kjeIYEA/IrnzZo2T3sCZGzTHUGdL5B08/V9dYmodiF7hiQf8nZ/clvQgN1BjQx6Ffnt78GTRSarxADsgmdAd3zQuH3I1pD8zIgm9AZUNCCwu/frwrNy4BsQmdAEVdfhH+zITQvAyrOyQtW74HXdAbUp0vh910fheZlQO4taxBcJeorq/fCSzoD+jRgjvPbxIDN0LwMyK23I3fJ+Sm1SkhBeo8DjZCiV//4w+oe0ghsXgbkTk7IL+rqo66e7ugfdB+J/rCW+jpYWDI4r+8COrmv5PxBofqpibY6gz1J8Rn9p3Pk7k9ZvR8++dRXAf2nW3BklSHgG/l96kAjbXUSO9/KZ8y7zG/xfBTQn7XfyZPPDetUgj49Kq/WXnW1AHvruF4/vzNlxRWRgTYPaJp29pujVbpPZjPGkvqz306aGirwyri4aWHPvxrdSOTDAWweUOxKbRU/0yezGeN0+0rlylcc5Mv/NNc3Va/3tKyRwDUobB5Q/w+01fNzfDKbMfoNzclxnLpntg+njHE+J4oS+HwAmwe0oId6Jb7zdff6ZDZDnArSnjV+08yHc7bbqa2eXKJ/qM0DyomK/e7klrbDPN/Tb3zvvPztpUo+nLP3Cm11j8CfijYPSL407Y7gdslGfMKZr/weqv31k1XXh3Muv+ucstxU77L+oXYPqARqr53xMGS0D6d0jKz3ytwna4u8esKA/M5P9QZ+tLBzlBEfbey9b6YMfUdoRgbkf87P6vv0kpJy6JMBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEcRfAjqxk5/tXiL5R0DHBwa3q9N2t7FTkC/4RUCOTiMvyI6lYWLXSicr+UVAW1to7z0eO8HYOcgH/CKgt53Xv9/4oLFzkA/4RUDzB2mr5THGzkE+4BcBZYapV8F0PJRk7BzkA34RkDyp+dr/7nikg9BVHslS/hGQvKFr7bYz2U8J5CcBUUnFgAjCgAhiZkD5mUfd3saAbMKcgAarFxyeWU2S6q1ycw8GZBPmBCQlyXKSNHDN+rgym4q+BwOyCfMCuiNO/eqFjkXfgwHZhHkBBa5Tv0qtVvQ9GJBNmBdQq0T1q9eaFn0PBmQTJgUU2iXugZAs2bE0dJTrz/8dUajMLD3bI79lTkBrEob1uL3CcjlDanfR9ee5WYVa7NCzPfJbJh4Hyr8iH0t3d73sKAZkD6YeiZ6Q7fYmBmQTpgYkHXB7EwOyCQZEEAZEEFMD+uyC25sYkE1YdToHA7IJBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQmwR0evYz47YaukXyjj0C+ir8qQXTGw1ydyUHMo8tAsq5LU1ZXu7Ei1b5ni0C+vJubZXa1cBtkndsEVCK83Oi9kUauE3yji0C+rq5tlrZ08BtkndsEVB+i7nKMrvZGgO3Sd6xRUDyf1rcOykubKqRmyTv2CMgOW/dtKQsQ7dI3rFJQGQVBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQBkQQywKaMO96vZ8YAHmgLzY+OgYb3+thbPwjPbHxsd2x8Y8/Ok9IbYsCmh93g0oNmkIC62Hjg8Kx8SGh2PhaNbDxdapi42+reuMj4p0RZ6wJ6C8a/YyN77wFGz8gGRs/4Q1sfOJwbPyaR7Dxu9pi4z1iQMVjQB4woOIxIA8YUPEYkAcMqHgMyAMGVDwG5AEDKh4D8oABFY8BeWB2QE0PYuO7fYmNH7QcG/9yAjb+n6Ox8RtisfHfdsDGe2R2QMfB8Scc2PgzOdj485ew8ZfPYePzTmHjHSew8R6ZHRDZHAMiCAMiCAMiCAMiCAMiCAMiCAMiCAMiCAMiCAMiCAMiiLkBpdwd1Pk78eFrJNVg0eHp67CdcI4X3Ym321VpPDNXfP7C8YLznxtVv3LrVbL4/F4yNaDUgKGrelQ+LDw+ITRJsVVwdH7beGgnCsYL7sRU6cXUCeUmC89/dbzg/AOqzk57SkrHHwQPTA2ocw9Zvlh3ovD4YV3E5/7t3b9J8cBOXB0vthNXqo1QlmNuyhOc/9p4sflPByyRZUfjQfiD4IGZAZ2UFirL5+oLb6BHnPjkaZ06VYoHdqJwvOBOZEmblWWKdFBw/qvjBef/+b5MZfm3x/AHwQMzA/pB2q4s5wRcEd3A7d3vqnznPOH5G8ZjO6GNF9yJy5mXleXomy4Jzn91PPCP4Ei7KRl/EDwwM6DPpP3KMln6U3B8foWQOWsHS8InlWoBADuhjUd24sNyY6F/BHW8+PxzKkmj8AfBEzMDSpcOKMsl0knB8VeWZynLgdXyBcdrAQA7oY0X34nsJ6VBucD8zvHi8x/8eGz5BPhB8MTMgPZKO5VlYkVsK2ukTMGRWgDATjh/hYnuxMbQ+muR+QvGC8+vGt3AoAfBPTMDOqE+E5CHNxAdfyxDPaN+vXRUcLwWALAT2njRndhYdph2Pr7o/IXjBedf9aA6bIF0AX0QPDH1afz9fWQ5N2K86PB0aamyfLae6Hjn/yDiO1HwK1BoJ3LDnyz4Smz+q+MF50+TvlaWz9SBHwRPTA0oreyr254IFn5nWF670KkbR5RZJTreGZD4TmjjBXfic2ncYtUlwfmvjhecP6dDxOJP4sskwQ+CJ+a+lLGqbVAX4Cj6xVFNqnbYJDy84G8Y4Z1wjhfbiSTJ6ajg/NfGC/4jnB3cuEob9T8v9EHwgC+mEoQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBEYQBefJowZuMpWlylQVW74z/YUCe7EhJSanbUVnsl/ulW70z/ocBeSOyn9V74LcYkDcKAqqu/AqrM7t9YMTcIz2r36pe+mJx68Dmi6zdN4sxIG+4BlR+yuc9A+rMTmtb6ZycWH5K2qiAuRbvnaUYkDdcA4qV5QPSSFlOlfadD5mq/HRIqLU7Zy0G5A3XgGbIcp60TK1ozy5p+/Hjx5Ml0z5HoARgQN5wDShBDShFC2hFwRP8vRbvnpUYkDfcBLRFyrZ4x6zHgLzhJqDsiurnUEwGPhKm5GNA3nATkDy+wrS0sQFzLN47SzEgb7gLyJEQGdgsyeKdsxYDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIggDIsj/A4NB6COn/MNOAAAAAElFTkSuQmCC" /><!-- --></p> <p>For comparison, the FOMC model is fitted as well, and the <span class="math inline"><em>χ</em><sup>2</sup></span> error level is checked.</p> -<pre class="r"><code>m.L1.FOMC <- mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet=TRUE) -plot(m.L1.FOMC, show_errmin = TRUE, main = "FOCUS L1 - FOMC")</code></pre> -<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAkAAAAHgCAMAAAB6sCJ3AAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+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////isF19AAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO2dCXgM5xvAvy0hkiDkkkSQCHFfiTvu+64IxT9UW9VSR6i7qJa2qNbVI1opgrjiJo44S4PW0YZWStwtSsQVOSTZ9z+zR2zWbjKz38zOzs77ex6zc3zzva/N79m5vvk+AghCAZE6AUTeoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQgUKhFCBAiFUoEAIFSgQQoX9C/Q+0XOUWdo7tL5zvYgNau029fJQN0f/4X+w8xMJiWM/4wiJZD7SZjX1dKkz8Nf8euYREmVY73fD/+YXriOzrh878yEzU8E4vFxRlkDPR+pmuz1jN93prF0qtgiMBbrmris5V1+PkUDnHTWG8AjXUe9NC91MgfByRRECjYrScAeGEuLSb3q4CyFvsptaEVJ+8MzwYoTsMBaoDSHhq1aPKUVeO62rx1CgK3vGuxCzApkJxwpErgNkO+oEKhBerihCoN262QRCqrPHnb+rEZIAsJGQWjeZxV2ENDASKMuBtGcX1hPymW5nQ4F65h+jeIRjBKpA1gGcJu7FWIEKhpcrihIohJDDmpkDhIQAhBKyX7P4RoNGLwoKdJ+QJrnMQtbatSd1OxsKtHT48GZFC2QUjhEojHwAsJj0Ls4KVDC8XFGEQGNWsPz5UEWq6lZWIapUcCPuBuUKHsICCAmctDlF/XK70TnQYvMCmQnHCPQVacj4QuZpBCoYXq4oQiAt35xnzmZ1K5nT13OpzK+MQbmCAiX6avZxDT+s385dINPhWIFOkmJPoRL5mRXIKLxcUZJA5wjprlvZnZBfmcNUC4NyRpfx6bHD6zmwuy3TbS9EoF2jGb4tIhwrUEoQOXiHOGSwAhmFlyuKEEh3UpL68phSlZD/oBzx0S69yMpSGwvEkv3LO8x1VJ52oRCBZrDCdC0qHCvQ2+TTLcwPj+YQVjC8XFGSQBCs/6MfJezJSHNCtNfoDYnDc5hDSDS7EE3IbLg0b5725Jk5+Pyj3Zm3QMbhWIGiSdeJjJ8agQqGlyuKEmg/ITWuMp9XahCyF2A1Ic0fMYu7VaQlwBZCerGFehGyGf4mpA17FaZuQpxztTtzPwcyE44VKJmUbU42aQUqGF6uKEogiCCkzKCPB5UmJAI0dpAq4xf+rxh5jfn7PvZj/qAzZzC/CxUfsb8KpPk3W5a3J2SAbl9GoAjN5dUKzc8FB4GMwrECgTtRkdtagQqGlyvKEuhpJ90Zbqen7OLNVtolh6/YpWNu2iX3Y8zCDW9dyTqPdfvOy39GoTmmcRGoYDiNQL0JqQRagYzCyxRlCQTq2MF1nesOjtWdtOYte71K6UbvXNEuPZrS3M2t+ZRHmoX0L0MrlgzosDJHv6sFAhUMpxFoASFv6AUyCi9P7F8gRFRQIIQKFAihAgVCqECBECpQIIQKFAihAgVCqECBGK72GRpVdCnEFCgQw6Rz0ErqHOQKCsSQob7dt4giR1/nssocadf5JSQjUCCWTX3uFlGCi0CTK3p/qJk519T7HYAfgrwiMm40C1oCMCpZmDxtEBSIYfcEoxaBL3T/Xi5xEGh/SMajoH3MTGaVpKxOe89Uufuk++z39uU2evL3KMFTthkULtCEaqdhc423eoRHaBaX+gdOyDsxeOAM9h/MDfCfpFkCjS2ajV22AYQc08waC/TLXoB+W5mZzQMAsjK3fgkQHTH56+u1Mob+a/X/mNVQuECwpUtG5fzmXEfr308fOPdE2WRg/+0LepIVGsPOsZte125cOQSuBhzRzGoE6lZZQ5K2su49s5iPBX0aeQ1KZ2YeNt725L2OB36bJsX/zEooXSB1tTd75S/MrNK2be0eJ9oAsP8mfwTw7TB2DliBtBsfe734/BPt7KtHtTub67Myfux99VnfWQBrA9dqVoc/2tllYpYV/jNSoHSBYH6Jl+255n/BnMA8OdGTEYj5N4k5dEUNYeeAFUi7Efrsb3BNO6sRqLOvBraDja2JAF+9z8x8/yYjz+t5Q3prz8z3L0htdHOmnFsdFobiBRro9vK34XyN+1nd4/UC7a35LLvVynyBtBthfdt2unLGv0DL2j150mGJOjnntu/VZz0XbO2gPTNX98682UG9abrV/kfWRekCnfDutfblUlRV3zGgFwjmBAZ8mJcvkHYjpDut0pUzFihvXOWKH+RkkRuwMbDiyOzxJVxcXN4E2LCCOajVbHXPav8l66JwgfKCVx4Jzim6HGIOhQsUHaKGiIrXpE5DxihcIIQWFAihAgVCqECBECpQIIQKFAihAgVCqECBECpQIIQKFAihAgVCqBBYoAnBiF0QyrX1gMACNVt5BrEHqnEdg0pogU4WXQaRAfVQIIQGFAihAgVCqECBECpQIIQKFAihAgVCqECBECpsRKBDvSu0/gZfu5IhtiHQ3GqxfYd2bSffIYmVi00IdN3zPtxwu9vpR2FjIFbAJgSKHspMxo2L7S9sDMQK2IRAS0czkwfu0d3MlUYKJ2mEFTA9WplNCHSoETv9uO5kYWMohxUtl4tO8CaToW1CoLwmU7MA1st5RFBpWfGO+DH627BAcC/cq31Qvche5kojhaN4gQD+OXQpN6vKCWGDKAYUKGff4k2pENNCbaY4UiiKFyipTui4MK+1ecFxwkZRCkoXKLPKOmb6l/eZwwHZwoZRCEoXaEcHzceXo6DLN8KGUQhKF2jJWM3H3q6Q5PVE2DjKQOkCrdU+w4iOABj2kbBxlIHSBXrgcZGZPq+/nbmcd7slbCBFoHSBYIPn9LjFgZqxjqa9LWwgRaB4geDG9H5jjmnmnvqcAXiUeDlX2Hj2DQpkQFS79A/KNQ+scVjYgHaNiAKd/003IxuBcuu2GPIYYL/X78JGtGdEFChCP9CHbASCmOLP2I/FQ4WNaM+8ItCBli4+Qwq/HFG/epKQXaBNsbaJugwF2lphTPdKjWb80kDYiPaMsUDLA3Y++3eut/mhPVLI3kqqGvOZuafjAh0D5qgBKmycWO4GrAp2qv0Ts7ritx1J2f4PoRkh5LFmF/kItNe/2I83z75dqamwEe0ZI4Ey3S6zH5+/aXaHFFLq3e0TVR8DhJVbuHMciWUECumzOXOpw6z4SNV3jEDl/nd0QfEP4GG/rve0T7jlI9BNVd9I5iOotbAR7RlWIPXq/MaD0yprPua6v2xPeLzgDimkHzOdVPop9ItiZoImMwLVU0O62xxm6V1PRqAGjDa9WsryEBZf0690wq9D/fEXiDOsQLmR+a2Xe3tpPgY7v2zQHFtwhxTCtnu4QE4z02fnVzhMZASaAnCaJKampq4ht6AiO/7vxGayFGhr319qlwqZnYjnQJwxOoQ9LvtAs/bVIez1pBC27d4jRqMT9VQ+3XxYgb4G2Ei0JEHFL0G2Al3xzsqpswsWvSlsRHvG+CR6eugV9mLkN9OlQf8L9Bf5Jc1h1H/MX4MVaDHAYXJfV6DiQpCtQPDGoLSEqtu9koSNaM8YC5S3yK1K+eBCGginkAHMdLrTowMkhTnp9tYJdL9kNLN6ZgeZC5Q+zrWJk9cxYQPaNa/eSFRfe1jYDinEcdSeaa9NhxsO/Y/vaOrcPVUjEEwpMTd+kmqJgUBv1T6jvWEkI4EAnpw+4n5H2IB2De870SlkWw/X6p8zV1obgpyb7FpdbpZWIPXCOk612OuyfIGO+JfWNtGyukDZD9IK2Vpk9y5ThvENqGAsEOgM3xjWFejWdH8VISUDp5q7GVqkQE99TvOKqGjsTqBzTn4jl6yJWTo6oJyZR6JFdzAVE4wtOrhidwK17ZahnckZ1NF0iaIFUrdZziekouEtUM4N3q+/WFWgMvnvdx13NV2CQxd3Fzwf8ImpZOyuQVnIKP3cbDMPJLj0kTh2BJ+YSsbuBIpTdYtOvJR8KiasmJl3TbkI9MT3FJ+gCsbuBILd7TTPVFTt480U4NRL66rgPF5RFYv9CQSQdjEhIcn4bujv/fU4r+NQhzqUPY/+o3/VhhMKva+qeOxRIJa8K0b9rj7cpMfxey4V/OF1H7Z6R6Wci6z0jwXxFYP9CbQ7vOdKWO5GHOea6bDFhVuPrB8Oy/XRPFOeZYWvSL7YnUCbSOPuDqNc5uyd4RBtugRHgZ5XWVFbM3O1Mp/4SsPuBGo4EuB78gUzN62h6RIcBYKtlRondvWqOuKCB5/4SmNFH/EHR+1oTYGc9gA80DR62+tsugRXgaCrg8+6eykzy7XkE19p/GyF4ZlDTD/8EEegavMBTpFVzNySGqZLcBboqqrxPYCTrp34xEeshjgCfeEY+YlvY58DqTvczHTXwlmg816B5ZrUqhzlxyc+YjXEEShnho/H+y+GE0J6PDddgrNApxrXXXn6z5x7nnziI1ZD1PtAl2LMtunhLNCTsvt90gBiu1gQHxEfqZq0chYIprYd+i4k+PwsbHxEIGxfoNx5rg4eQfuFDY8Ihe0LBPAiKiBT2OiIYMhBIICwWcJGRwRDHgL963FB2PCIUMhDIPixCbawt01kIpC680Jh4yMCIROB4LrnZWETQIRBLgLB4jY4FJQtIhuB8kI5tWFErIxsBIK/PHAQBBtEPgLB513Yg1hSiEvZ1ua7H0WsjIwEymmyHGBz8VpzZ1VxOF50ccQqyEgguOSeAqWHMDPqdtg6yFaQk0Awv/01labXhpsqfOPQRpCVQHmhE4tpZ8gzYdNBLEVWAkGym0rTtHtDcWGzQSxGXgLB16XdLwL84txX2GwQi5GZQHmtvIqXdy3eMkvYbBCLkZlAcLvCiokfYe+JtoPcBIJVtbF1oi0hO4Gg/xQh80AokZ9AD3yOCpkIQof8BIKdAU8FTAShQ4YCwXtvCpYGQoscBXpeI7boQoh1kKNAcNbzhmCJIHTIUiCY3wpf0rAR5ClQXofPhUoEoUOeAsE/Xng32jaQqUCwq+pjYRJB6JCrQPBBf0HyQCiRrUBZDVYIkghCh2wFgsseXFNHRES+AkF07QwBEkHokLFAMPg9+joQSuQs0LOaMfSVIHTIWSC44PGXALUgNMhaIPixjpluqBFrIW+BYCiOAiUxMhcIT4OkRuYCwUXPJGEqQixD7gLB+mr4UExKxBMo+0FaIVsFEwhG98K+7yREJIFuTfdXEVIycKq5rqCEE+hF6HyhqkL4I45A55z8Ri5ZE7N0dEC5302XEE4guOt7QLC6EL6II1DbbrrHVDmDOpouIaBAcMgHe0+UDHEEKhOnnzvuarqEkALBgob4WFUqxBEoZJR+bnZT0yUEFQgGDhGyNoQH4ggUp+oWnXgp+VRMWLE40yWEFSgjeJnmc9sHwxanC1kxUhQiXYXtbkdYVO3jzRQQViC44X2U0ahr42Urh/jhnUVrItp9oLSLCQlJD41WHi1fTodK4GvvA8yJ9Oz+bNebaxoIWzNSKKLeiX502agzVXWaHmdhf4EAFjZ6XvusZq4KDstiRUQSaE3zJLjXk5DSi80UEPgQxvBOP8//NDOh2Au5FRFHoCWk1T3o7Pv93knFzTwtF16g7Fa+R9jPHM/bQleNmEccgfynAdwj7EjdUxuaLiG8QHC3fGAaQN60boLXjJhHHIHctgBcIOwFdXxp0yVEEAh+c3IfPr5+2/+ErxkxizgCvd49G3JKH2TmplnlRqKODT4Lv8b+76yLOAJd8qyz4Mhcn5+OTHNYb7qEKALB7GBsI21lRLoKuzyyvOZOYu2NZgqII5A6IhxHYbEuot0HenHt5J5E8wPDiSMQZLeZIUq9iDlk36TViPsB2MreqtibQPCnF55GWxO7EwiO4Xsa1sQigW5doI4rnkCw3v+uaHUjxlgg0L4KhECHpXRxRRQIPg3GNkFWg79Aa4qPWENgpmo5VVwxBYJ3emEvwNaCv0C1xkMqszCpDlVcUQV60QV7DrIW/AVyitcIFO9MFVdUgeB589liVo+8hL9AjT7WCPRpPaq44goED4K+EbV+RA9/gVY6zEkk96NLfE0VV2SB4KrPFnEDIFosuApb6kYIKTmV7o10sQWCM54nRI6AsFhyHyj99MbDDyjjii4QHPDCboCtgP3dic5nm/ffosdA+AvUUw9VXCsIBKuq3hE/iNLhL9Awltc9XxtNFdcaAsFndYxfTEOExtJDWHrnxlRxxRKo4C3oCS3woYbIWHwOdJTcp4krikAp/VxLNd5lsEL9XvtMEeIgL7FYoFWlqK7jxRAo2WtRWva+at8brMob3CVL+EDIS/gLtFbDbO+2VHHFEGig5t7mFTdDZXLDw3KEj4Tkw18gRw1OzS5RxRVDID/tWM6NfjNcmdXpbeyEU0Ts6T6Qj/ad5ianCqxNb/U+GiQe9iRQH02d/7oaXXmltx6DBokGT4HcDaCKK4ZAZ7025EFSyGfG6580iRQ+GKKFp0BRUVGLXXwjv/4oyG83VVxRLuN/bVHa0/eHV39uHjeeIEI0hIX/Iez9NuxVTm5XukZ/It1ITDfds8LDBtNECYdYIJCvttvM7T5Uca3yKOMlj0I+tGo85WCBQNpOx5ZWooprZYHgUaPgKiXq/IBvzgsNf4HeLcs+LNhd1iYPYWZ5XtPz7exTLUdYN6oC4C/Qs7akfJ3ypD3dY0prC7RoQFpIpPp5xYvWDWv/WHIf6PC8yC9/poxrbYH6bYLHLd7NG/F90UURPtjTjcTC6LONuUTrMGgk5Qu1iDE8BVqRACv0UMW1tkCfjoR/frvbuTR2ASwwPAVyGQiueqjiWlugVK8q3iFl/D17Y/sgYVHKIexfj4bOVVwatRjU/ql1A9s7FgqUvJeytbG1BfpoAmRcz1XXPzIuhPaNJMQQ/gLd6jQRthQj7mbGsuSItQXqtUPzMXYJfFzzH+uGtm/4C/S6z3Zo0OVaJ9t/rceQ/rGaj7eXA3xRNcW6se0a/gKV/wHukASI9aCKa22BorqzD+nTKlxhpst9zlk3uD3DXyDXDfCTYybscqGKa22BspuHnbqzs472ofx2j/3WjW7H8Beoc4uDdftAeo8QqrjWFgiy5od4t9+uWzjqaa4DdIQn/AX63YuUOQfVS+wqrHSRWF2ggvzui/0HCYMFl/EZZ5hL+E2UwwJKLBBcrzEWm3YIgXjd/GY/SCtkq9QCQVqrcLwpLQAidfN7a7q/ipCSgVPNjZYhuUCQ9Ua7R1LnYAeI083vOSe/kUvWxCwdHVDOzP1G6QWCvLF1bkidg/wRp5vftt0ytDM5gzqaLmEDAgEsr/ir1CnIHnG6+S0Tp587buahvU0IBPs8zYyHh3BFnG5+Q0bp52Zbc8hL/pzzWyh1CjJHnG5+41TdohMvJZ+KCSsWZ7qEjQgEtxu8ky11DrJGpG5+d7fTjHipah9vpoCtCATpfVviMM8UiNbNb9rFhIQk40ZDOVf1ONuKQKD+uCq+qmE5fAV6tn/jTeZL//fPoz047LTdqBu841UDdLy2gE+W4rLGk+6xjKLhKdBffsyBaeQfldjDE5edDprbYjOHMJazlafgcw0L4SlQL7c1Ses9fJtvOHC6sBcLt0ZoIR0iIkyXsCmB4H6b/tidq2XwFMhjPjOZT4q6g7vHiQQ3YyA1mzUzXcK2BIKst+ubH6IcKQSeAhG2Hc2Wot/QSA5uwo4zIJNDGMsyr71SpyBL+ArE3tbZzuEVnxdTXL6XlUBwwhdPhCxALIEAjvp1uysngeB+ux74eJ434gkEj95wl5VA8GJ0Nbp3lZQIX4EcHB0dHYimq+ii91kfaXbAJVsUiDm781osdQpyg6dAMwygimubAkFy3SHPpc5BXijl3XiupP+vXrLUOcgKFMiYHzxipE5BTqBAr/BX3SF4W5ozKNCrpA+rlSR1DrIBBTLFGo+lOLwGN1Agk1xv2QkH7OUECmSaF1N9sQMGLqBA5jhcaWyG1DnIABTILI+H1DwrdQ62DwpUCOs85+UWXUrZoECFcbN9C8peSOweFKhQ1Mvd52ErocJAgYrg7+Ydrkudgy2DAhVF7jz3b/GuollQoKJJadPSbLsmxYMCcSD3K/dFeDlmGhSIE1faNsXnqyZBgbihXu01JUvqJGwRFIgrd8KCjkidgw2CAnFnV+X+94supTBQIB48i/RaiVf0BUGBeHG2SSsufWQrCBSIH+rV3mOfSJ2ELYEC8eXB2xVjzRzH7u3arLg+PlAg/iQGtzb5CvRcz87dvYcrbAAFFMgC8pZXGPXqmLHfV6/mUckrZKQECUkICmQRaaM9l7wwWlfZ6yDAhXqOynqpDAWykOQe1TcVXKPSdIl9tRjXb9Q+QIEsZndQ1wKX9Cp2PFZQF/tZmnQkAgWynBdLvd69+3LRIZKdxpa4KlU+koAC0fBoktsnz/QLHSsM3HlgnJu3su5Vo0B0XB/s/Z3ubDrJo3fndn29t0qbkLVBgWg516naJu2PzpUBnmU7JkqcjrVBgehJCG58QOocJAMFEoKEei2PSp2DRKBAgpC7snLPM1InIQkokEC8WO7bUYmv0qNAgpG52DtcWXehWVAgAXn+lXfYeamTsDIokKBkLPLpo6xzIRRIYDKWVeqqpKdhKJDgZEdXax2vmOcZ4gmU/SCtkK12LBBzUR9bv/66HKmzsA4iCXRrur+KkJKBU821EbZrgQDU8W38lymiZZk4Ap1z8hu5ZE3M0tEB5cwMoGTnAjGcDPOYcU/qJMRHHIHadtN1cJozqKPpEvYvEMCVkeXe/VPqJMRGHIHKxOnnjruaLqEEgQAefOrdZZ99n0+LI1DIKP3c7KamSyhDIICsVQ1qfm/PJ0PiCBSn6hadeCn5VExYsTjTJZQiEMORMLcP7beZq0hXYbvbERZV+3gzBRQkEMCNye699tppZ6+i3QdKu5iQkGT89t3d5XpKfsezPnmTER0cuPDVdxHtALEEyr2mfcU3s8CV7F8j9DgobnTb00Ndh5yQOgnhEUegnJmlSKnJbL+Uq8zsp6hDmI6HXwXVWWZvQ9OLI9CC4h/GjS/+FqBABVEfHuQ65JhdXdeLI1D1acxkLdmBAr1C6uLa1b/4V+oshEMcgZz3sNMI/0wUyASJw8v13JItdRYCIY5AjSew0/88P0CBTJK+uo3nWPtoeCaOQEvJ2ANZAHuKDZ2CApnm2ieBteffljoLekS6jJ9ThqQwH7t8CApkDvXxEW7tf5J7h4ti3QfKvqrp1z3nUJTp7SgQS9bWMNfwLbLuFA+btEpM2ooO5Yftk2/zRRRIeu4sbu7x3iGZDgeEAtkEN74M8Xz/oBwdQoFshavzG3u8u092t4dQIBvi+sKW5SO2PJc6DV6gQLbFne86lekdLaNBgVAgm+NR7ADXlvP+kjoNjqBAtkjWvlGVqo4/KIcTIhTIVjk/t5lrv+g7UqdRFCiQDXN/9RvlG0772abvMqJAtk3uiY+Cy4UtvyF1HmZBgWyf/9YO8ao+attjqfMwCQokC9Tnv+xcuvmMI7Y38jgKJBsyD01vWrrzFydt65QIBZIVj3dE1i/bY8Fp25EIBZIdqVvG1ivT9fMTtnE4Q4FkRcqIJh3nZwI83D6+sUvrGfHSn1ijQHJij2f4G4Pa19b2HfjswMx2peuNWivtQNEokIzI8fINW/JZULXxL9ec/jrcxzts4YkMqXJCgWTEmVI/MdPs5l4FV99YN6axU5Mxay5LkRMKJCNWO2o+tpZ4dVPGia8GVC7fbdbOu69uExUUSEZsdtAcqX4oaWb7vZ0zu7r59Z2774H1ckKBZMRlpyHPAC541yis0NVNkzu4Vnr9013WeQEfBZITTZt6dm9ZoZqZd+1eok7ZNKWzu1fXqZv+FrtjNBRITlxv0GzwQN8xHKW4vfOTvgEuTd/77riIt4tQIFmRu+3jr/iNSfbk+DcjmrlU6TUtNumFGBmhQEogL2Xrp/1rlqo94NMtyQI/RkOBlEP277HT+gSWqjtg9sY/BHuQhgIpjczzsTPCa5cK7DHxh2P/0VeHAimTnMs75r/Tws21ScScDWefUlSEAima1MSfpoXXd6rQ6q3PNp616FoNBUKYC/4jP07p16C0e9NBM346eovXrSMUCMnnv1/WfDI01Kdk9S7vz994mtsJEgqEGJN5ac+3E8ND3J1qdR/5xboT/xT6i4QCIeZ4dnHXN5MHtvApUTk0Yuo3O8+b7PIBBUKKIvvasZjPRvas7+ZYLSR0vdFYDSgQwpkbbXxadPBaVWAdCoRwps20HIBLlQ8ZrkOBEK5cCMiDLIBVYYYrUSCEK1v7rq5RynXQgTqGK1EghCuH/UJO5T763K3AOMooEMKVm68dYKbq6qGGK1EghCu7gz1nH9zQJqSR4UoUCOHKzh43IjsMWP1bQ8OVKBDClX/dNM/rPxlpuBIFQjgzvvUlyP6mwnXDdSgQwpncpd7uzt2TC6wTT6DsB2mFbEWB5MkD41c7RBLo1nR/FSElA6ea63sEBbITxBHonJPfyCVrYpaODij3u+kSKJCdII5Abbvp+qvJGdTRdAkUyE4QR6Aycfq5466G6//or8dhEZ/6EJtFHIFCRunnZhd4cPJwk55Kh/nUh9gs4ggUp+oWnXgp+VRMWLE40yWaneRTH2KziHQVtrsdYVG1jzdTAAWyE0S7D5R2MSEh6aHZzSiQnSDVnWgUyE5AgRAqUCCEChQIoQIFQqhAgRAqJBNo6vKC9BkcIQDh3YSoJaLjICFqGdBFiFoiOg4UopY3OgtRS0TnBUZ/N1+JBPpxhBGOVWsKgG8ZIWqpWTJAiFr8XISopWapKkLUUslZiFpqOnc2+ruNfSKNQK9QTZChQjYOEKIW7r/LhbKzlxC1CHS0P9hBiFqgncXPMFEg3qBAhqBAvEGBDEGBeIMCGYIC8QYFMgQF4g0KZAgKxBsUyBAUiDcokCFiC1RTkNGrtw4WohZo9KcQteztK/tKla0AAAR+SURBVEQtEPqbELUc6yJELdDpuKV7ii1QqiC15AgztJowyeQ+KroMBx6qhahFbb51KB8sT0ZsgRA7BwVCqECBECpQIIQKFAihAgVCqECBECpQIIQKFAihAgVCqECBECpQIIQKcQWKa1y23XnqWrZqOiUaTldJwg4hMtLWQpnQsqYuQV/m0Cajr4UumWeR/s7Bm8HyZEQVaLdq5Oauzrdoq1noGcVwlKqOvCYTBchIVwtdQnPIhN1Ti8+kTCa/FrpkIkovin+LJFiejKgCtesKkOE3nbaaUdSNpm5/25pMpM4ovxaqhLLLjGWmH5bKpUrmZS1UyTxWxQCog4ZZ/s2IKVAaiWam7/vT1tN1BG0N8aGhjhOpM9LXQpfQVcIOzxVHrlElk18LXTKX26Yw09YDLP9mxBToIklkpktU2ZT1VO/SyLn+cspKAicKkZGmFrqEslKymOn4UplUyeTXQv/tqONLrbH8mxFToIPkEjNdQx7QVZNXwm3J9uFkIV0tmj89dUaaWgRIaG3xSQJ8PWwt1MkscSSRFN+MmAIlEHYEmBhS2OgsHMjecJWZDi2TR1WL5k9PnZGmFuqE7g8hw3Kok9HWQp3MtW2THBZanoyYAiWRU8x0aUlBKttKUqj21/zpqTPSHsIoE9rj6b+dPhldLbTJsIyvankyYgr0kD3FhzFVKav57wzb4nsnuUdVi+ZPT52RphbKhPYUG5UJ1Mnoa6FLZnN3du8V5LnFyYh6Gd8+DCAnYAplLQlkHTN9rxJdLdrfDtqMdAdCmoRyfIbo5miSya+FLpl48iszfaei5cmIKlB8sU9ODC5H+2ZYblPPOXvGvraZrhatQLQZaWqhS+gQmbyKJZMqmfxa6JJ50SJg1b6Jr0VZ/s2I+yhjc5OyHegfZWRE1ijdYi9lJbqzF8qMtLVQJRRFtNyjSuZlLXTfztPhQS4h7G+Ypcngw1SEChQIoQIFQqhAgRAqUCCEChQIoQIFQqhAgRAqUCCEChQIoQIFQqhAgRAqUCCEChQIoQIFQqhAgRAqUCCEChQIoQIFQqhAgRAqUCCEChQIoQIFQqhAgRAqUCCEChQIoQIFKopw3UvEZC64rJA6GdsDBSqKk3FxcX4tmcklGJggdTK2BwrEhToDpc7AZkGBuKATyJU5hFVc1Nwp4Ls7PV0rs31arAp2qv2TtLlJDArEBUOBHGYd6qmquCi+ieMzWOowKz5S9Z3E2UkKCsQFQ4H6AySTcQC7yYV0tznM2nc9pU1OWlAgLhgKNB8gl8SyFv1+miSmpqauIdSDOcgYFIgLhgItZAWK0wi0UXeBnyRxelKCAnHBjECHyX2JE5MeFIgLZgS6X5IdYGIm9VgwcgYF4oIZgWBKibnxk1RLJM5OUlAgLpgTSL2wjlOtKImTkxYUCKECBUKoQIEQKlAghAoUCKECBUKoQIEQKlAghAoUCKECBUKoQIEQKlAghAoUCKECBUKoQIEQKlAghAoUCKECBUKoQIEQKlAghAoUCKECBUKoQIEQKlAghAoUCKHi/zXgn1DiioA0AAAAAElFTkSuQmCC" /><!-- --></p> +<pre class="r"><code>m.L1.FOMC <- mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet=TRUE)</code></pre> +<pre><code>## Warning in mkinfit("FOMC", FOCUS_2006_L1_mkin, quiet = TRUE): Optimisation did not converge: +## false convergence (8)</code></pre> +<pre class="r"><code>plot(m.L1.FOMC, show_errmin = TRUE, main = "FOCUS L1 - FOMC")</code></pre> +<p><img src="data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAkAAAAHgCAMAAAB6sCJ3AAADAFBMVEUAAAABAQECAgIDAwMEBAQFBQUGBgYHBwcICAgJCQkKCgoLCwsMDAwNDQ0ODg4PDw8QEBARERESEhITExMUFBQVFRUWFhYXFxcYGBgZGRkaGhobGxscHBwdHR0eHh4fHx8gICAhISEiIiIjIyMkJCQlJSUmJiYnJycoKCgpKSkqKiorKyssLCwtLS0uLi4vLy8wMDAxMTEyMjIzMzM0NDQ1NTU2NjY3Nzc4ODg5OTk6Ojo7Ozs8PDw9PT0+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////isF19AAAACXBIWXMAAA7DAAAOwwHHb6hkAAAgAElEQVR4nO2dB1gU1xaA7wZRBFQIvSkiFiJWsKAYxYJiiQa7wZJYosTeazRqjC0qJCYYxShWFDU2LBjFWKKJUYM+RcUeKwoWlL73zWzBZd2Fmb0zOzs75/8+Zmdm755zhN+dfi/CAEAAEroAQNyAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBABAgEEAECAUSAQAARIBBAhPkLNBypSaaWDgyoZ1M3Yqtc+Z58VbCDVdUh/9LzExFKoF8TEBpLvWR83cTZ1r/PX0VxFiIUoxn3pyHX2KVrS63rTs9MoGZctdOLFWkJ9GaEajbsNf3Ww1DlksVyrC3QLUdVy/nqOFoCXbBSGMIiXVu1N81UM8XSixVJCBQZo+AhHoCQbffpPWwRGki/1QKhD/vN6mGB0G5tgVoi1GPd+lHl0QdnVXE0Bbqxf5wt0iuQnnS0QOg2xrlWKoGKpRcrkhBon2o2CaEa9HbnWnWEkjCOR+iju9TiXoTqawmUY4la0wtbEPpW9WFNgToXbaNYpKMEckWbMD6LHC1ogYqnFyuSEigQoaOKmcMIBWIcjNAhxWLv+g3zigv0FKHGBdRCzsaNf6o+rClQ9JAhTUsXSCsdJVA4+grjFeiTMrRAxdOLFUkINGoNzf+ey1A11UpvJHuGHZCjRrvimzAfhHwnbU+Tv3tfax9ohX6B9KSjBPoeNaB8QQsVAhVPL1YkIZCSHy9Qe7OqldTu6/ln1LeMRrviAp32UHzGrsdR9fvMBdKdjhboT2TxCldGf9ACaaUXK1IS6DxCHVUrOyL0F7WZaqbRTuswPmvzkLqW9Md+UL1fgkB7R1KsLCUdLVBaTXTkIbJ8SwuklV6sSEIg1U7Js3fblGoIPcH2yF25lJeTI9cWiCb31GDqOKpQuVCCQDNpYTqUlo4W6As0dwf1xaPYhBVPL1akJBAOUP/RkxG9MxKEkPIYvQGyfIPnIRRLL8QiNAdfXbhQufNMbXz+U36YtUDa6WiBYlGHiZSfCoGKpxcrkhLoEEK1blKvN2ohdADj9QgFZVKL+2SoOcY7EOpCN+qC0HZ8DaGW9FGYvDGyKVB+mPk+kJ50tECpqFIQ2qYUqHh6sSIpgXAEQhX7zu5bAaEIrLADeY9b+pkF+oD6+77wov6gs2ZS3wuemfS3Agr6cceq1gj1Un2WEihCcXi1RvF1wUAgrXS0QNgRydB9pUDF04sVaQn0qp1qD7fdK3rxbgvlkuX39NJxB+WS43Fq4Y6bqqX/C9VnFxZdo1Bs05gIVDydQqBPEKqMlQJppRcp0hIIyzf3q2NTp99m1U5r4Q/dvCs0HHxDuZQ5JcjBIWhKpmIha0mwZzmfNr/mqz9qgEDF0ykEWoxQb7VAWunFifkLBPAKCAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEMXNrgNiSm8F6AIEoph0HrcQugaxAgJRvJXf/7SUJsndmKzSR8ZtdgWJCBCIZlvXR6W0YCLQZE+3CYqZ803cBmP8S02XiLd3mtaMwjgylZs6TRAQiGLfeK07AvNUP++WGAh0KPBtZs2D1Ey2d0pOuwPnvB+97Djny4MFDV9ei+S8ZJNB4gKNr34Wb6/1eaceEYrF6Kq+4wtP9uszk/7B832qTlIsYYUtijfb78I48LhiVlugUwcw7r6TmtneC+Oc7J1LMI6NmLzs9kdvBzww+j/MaEhcILyj/dsqRbdzJdd7mtVn/slKqZj+OVjzZU5wHD1Hv9VN+eav/fFNn2OKWYVAYVUUpCiDdeycQ70s7trQpW8WNfO80a6XX7Y9/Pc0If5lRkLqAsmrD+xStDDLu1Wr2p1OtsSY/pk8A+OVg+g5TAukfPOFS96Cb5Sz72/VHm6vR8s42+3m60+/xnij70bF6h6Ze9pPzDHCP0YIpC4QXlT23f1ci76jdmBenuxMCUT9TKI2XTH96TlMC6R8E3c9VP+WclYhUKiHArqDjZ2nMf5+ODXz80BKnm6F/T9R7pkfWvys4d1ZYr7rsCQkL1Afh3ffDRdqPc3pmKgW6IDf69wWvxYJpHwTb2kVomqn/Q30Q8jLl22i5Kn59z1uvu68eGcb5Z65/JPsu23k26Yb7V9kXKQu0Em3LhvfLcVU8xiF1QLheb4+EwqLBFK+ibOs16naaQtUOKaK51f5OegOjvf1HJE7rqytre1AjLeuoTZqfi0eG+2fZFwkLlBhwK/HAvJLbwfoQ+ICxQbKcYTnLaHLEDESFwggBQQCiACBACJAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACI4Fmh8AGAWBDO9e4BjgZr+eg4wB6ozHYOKa4H+LL0NIALqgkAACSAQQAQIBBABAgFEgEAAESAQQAQIBBABAgFEmIhASZ1dg6PhsSsRYhoCza0RHz4wrKV4hySWLiYh0C3ndHzH4VHoL9zmAIyASQgUO4CajBmzpSe3OQAjYBICRY+kJumOsWH6WgMlkzLMCOgercwkBPq9Ad3Byew6k7nNIR3WNF/FOwHbdKY2CYEKm07OwXiTmEcEFZY1g/nP0dOEBcJPejm38q03rou+1kDJSF4gjB8cu1aY432S2ySSAQTKP7AiPh3HNZPraQ6UiOQF+rd2izE9XDYUBiRwm0UqSF2gbO/N1PSK27mjPrncppEIUhdodxvFy5JI3P5HbtNIBKkLFDVa8XKwA05xecltHmkgdYE29VC8xPbHeNAMbvNIA6kLlO50mZq+qbsb4/8c7nGbSBJIXSAc7zxt+/JqX9Gz077gNpEkkLxA+O6MHqP/UMy9cj+H8fOTqXB3GQtAIA1iQl6PsG9eo8YRbhOaNTwKdOFv1YxoBCqoEzSQOhRLcr3AbUZzhkeBItQDfYhGILy+zGv6Jao/txnNmfcEOtzc1r1/yYcj8oL3VuUWu6dYuRMhQoF2uI4M82ow41R9bjOaM9oCrfLZ8/rBfDf9Q3ukoQOVZbUWUXOvxvha+cyTY+waP9H+Dl4XYF17LbXac2VbVKnnc9wUIfRC8RHxCHSgqsXqe/8MrtyY24zmjJZA2Q7X6ZcFA/V+IA2VH/rbRNlsjMPtl+4ZgzZTAgV23Z4dbfl14ljZT5RA9p8lLy7zFX7evcNj5RVu8Qh0V/bpWOqlxsfcZjRnaIHk64tuHpxWRfEy3/Hd/YQnin8gDXWnppMqvMLdY6iZmpMpgerKcZbDPGppqDMlUH1Kmy7NRbkJS/TzrHDoTIQPfAMxhhaoYGzR3cufuChe+tm8u6F5c/EPpCH6vodL6Cw1fX1hjeVESqApGJ9Fp589e7YB3cOe9Pi/E5uKUqCdn572L9943mnYB2KM1ibsRaV0xdr3h7BXk4boe/cyKY1O1pW5h7nTAi3DOB4pScGeS7BoBbrhlpPvvxcvH8htRnNGeyd6evAN6j+i69+6W2P1N9AVdCrDMvIJ9degBVqB8VH0VNXAcykWrUC4d9+MI9V2uaRwm9Gc0RaocLmD94cBJdwgnIZ6UdPp1pmHURq10+2mEuhpuVhq9aw2IhfozdhKjaxdjnOb0Kx5/0Si/Nbzkj6Qhqwi90/7YDq+Y9nzxO4mNh2fKQTCU8rOT5wki9IQ6PPa55QnjEQkEMav/kp2fMBtQrOG9ZnoNLSrk12NBdSR1taaNo33rrf/WimQfKm/9Uf0cVmRQMeqVlDeomV0gXLTM0p4t9TuXaYMYptQwhgg0Dm2OYwr0L3pVWUIlfOdqu9kaKkCvfI4wyqjpDE7gc5be42I2hAXPdLH/qLuFqV3MBUX8P7FGkA3ZidQq7C3ypn8vm11tyhdIHmrGDYpJQ1rgfLvsH78xagCVSx6vuuEne4WDLq4u+yczianlDG7G8oCI9Vzc5robsGkj8TRw9jklDJmJ1CCLCz29NXUM3HhFnqeNWUi0EvYj2aI2QmE94UorqnIWifqacCol9b1AYWsskoW8xMI44zLSUkp2mdDL/RUY7ORQQx58CpqerGHT/1xJZ5XlTzmKBBN4Q2tflcztqmx+plJgH9dnuIdbqtuXhjn9Z8B+SWD+Qm0r0fnX/EqB2Q1X0+HLbarGcWZMCjfXXFNefYQVvklhtkJtA016mgZaTvvwEzLWN0tGAr0xnt1bcXMzSps8ksNsxOowQiMf0bfUXPTGuhuwVAgvLNKo1PtnX2GpjixyS811nTlf3DUtsYUyHo/xumKm94O2OhuwVQg3MHSY/OTm1/bN2eTX2r8YYThmQN1X/zgR6DqizA+g9ZRc1G1dLdgLNBNWcAjjE9VascmP2A0+BHoO6ux33g0cj/8bLeDnu5aGAt0wcXXvlEt79VebPIDRoMfgfJnujsNzxuCEOr0RncLxgKdaVRn/V+pBY+d2eQHjAav54Guxp3V9xZjgV5VOuiegfGmDgbkB/hHqFtaGQuEp7UcMBQfdj9RektAAExfoILFdmUc/ZK4TQ9whekLRO1Rxfhkc5sd4AwxCIRx91ncZgc4QxwCPXC+xG16gCvEIRBe3QjusDdNRCKQPHQJt/kBjhCJQPi203VuCwC4QSwC4RUtYSgoU0Q0AhUG/8RtBQAniEYgfMXpLrclAFwgHoHwglB6I/ZvgE3F4DRuiwEMR0QCFTRZhXF8Gf8Fc6pa/sFtNYDBiEggfNUxDVcYSM3IW3tyWw1gMGISCC9ufUum6LXhrgyeODQRRCVQYYuJFsoZ9JrbcgBDEZVA+JqDTPGU2JYy3FYDGIy4BMLLKjhewvikTTi31QAGIzKBClu4lrGvVCYYRgc3FUQmEL7vunbSTP0daAPGRmwC4XW14e5EU0J0AuGeU7isAyBEfAKluydzWQhAhvgEwnt8XnFYCECGCAXCXw7krAyAFDEK9KbW5tIbAcZBjALhf5zvcFYIQIYoBcKLWsBDGiaCOAUqbLOAq0IAMsQpEP7PRW+/H4BREalAeG+1F9wUApAhVoHwVz05qQMgRLQC5dRfw0khABmiFQhfd2JaOsAj4hUIx9Z+y0EhABkiFgj3+5I8BkCImAV67RdHHgQgQ8wC4UtOVziIApAgaoHwan893VADxkLcAuEBRhimBigJkQuU5beekziAoYhcIHzZOYWbQIBhiF0gvMUXLooJCX8C5aZnlPAuZwLhUV2g7zsB4Umge9OryhAq5zv1lp4G3AmU12IhV6EA9vAj0HlrrxFRG+KiR/rYX9TdgjuB8CPPQ5zFAtjCj0CtwlSXqfL7ttXdgkOB8FH3e9wFA9jBj0AVE9RzJ+x0t+BSILykAZxPFAp+BAqMVM/NaaK7BacC4YERXEYDWMCPQAmysNjTV1PPxIVbJOhuwa1AbwOiFa87Iwcuh77LjApPR2H7QhCNrHWingbcCoTvuB2jNGrfeOW6AZ5wn5kx4e08UMblpKSU51ork+3VyBaxjFcKSdSO9Oye9BmhjfW5jQyUCK9nojOva3WmKs9QY8PtNxDG3zd8U/u8Ys4bhmUxIjwJtCEoBT/ujFCFFXoacLwJoxgS7vxEMRMM4/MaEX4EikItHuNQj58PTCqj56ZB7gXK/djzGP2a7/wf16EB/fAjUNVpGD9G9HgEUxvobsG9QPjxh74ZGBdODeM8MqAffgRy2IHxJZRFzSVW0N2CB4HwOWvHIWPrhjzhPjKgF34E6tYxF+dXOELNTTPKiUQV8R5LVxznIS6gH34Euursv/jYfPe1x6ZZbtHdgheB8DcBcE3DyPB0FHZ9xIeKM4m14/U04Ecgef/uMAqLceHtPFDerT/3n9Z3NxBfAuHcVjN4iQvoQ/S3tGrx1Afusjcq5iYQvuJyjKfIgC7MTiB8HJ7TMCYGCXTvEnFe/gTCW6o+4i02oI0BAh10RQi3iSbLy6NAeF4A3BNkNNgLtKHMsA0Iz5KtIsrLp0B4cBfoBdhYsBfoo3H4GbUwyZ8oL68C5XUYxmN0QBP2AlknKgRKtCHKy6tA+E3QHD7DA+9gL1DD2QqB5tYlysuvQDi95g+8xgfUsBfoV8t5p9HT2LLLiPLyLBC+5bGD3wSAEgOOwqIdEELlppI9kc63QPgf55M8ZwBoDDkPlHU2/mg6YV7eBcKHXeDxDCNgfmeii9jllsp7DoC9QJ3VEOU1gkB4XbWH/CeROuwFGkTTzfmDkUR5jSEQXuCv/WAawDWGbsKyQhsR5eVLoPxiSxOC4KIGzxi8D5SMnpLk5UWgG+F25QP3aKyQf9k6m4c8wDsMFmhdeaLjeD4EuuqyIjP3UI2fNFYV9gvN4T4R8A72Am1UMMetFVFePgTqs5yepn2oqUxBz/B8Pc0BLmAvkJUC66ZXifLyIZCXciznhn9rrsxp9zncZ88j5nQeyP2+4qXxmWJr33z8JXTjyh/mJFC3X+jpA/us4quzPh4KBvEGS4EcNSDKy4dA/7hsKcT/Brw3IvjLJmO4TwYoYSlQTEzMCluPsctm1PTaR5SXl8P4v5rbOnmufv/r5kWjcTxkA2jYb8KGt6SPcgo6kA0XyNOJxCzdJ6ee15/KSzrAAIE8lN1m/uZOlNcolzLekdloglHzSQcDBFJ2OhZdmSivkQXCmQ0DqljWXgVH9FzDXqChlfZS032VTHITppc3tZy/yDsbDDfbcw17gV63Qh/6f4haZ5XUulSMLdDy3hmNxsjfeJI/EgkUw5DzQEcXjl3yB2FeYwvUfRt+0Wxo4bCfjZvW/DGnE4kl0XUXdYjWtu8IeFiDY1gKtCYJr1FDlNfYAs0dju//9bhLBbjTnmNYCmTbB9upIcprbIGeuVZxa1TJ27kL3B/ELVLZhP3n1NCmim3DoH4hr4yb2NwxUKDUA4R3GxtboBnjcfadAnm9Y2MCSJ9IAjRhL9C9dhPxDgvkqGcsS4YYW6AuuxUvo6PwQr/7xk1t3rAXqJv7b7h++1vtTP+xHk16bla8fLEK44XV0oyb26xhL9CHv+CHKAlvdiLKa2yBYjrSF+kzXG9Q01/c/zFucnOGvUB2W/Faq2y815Yor7EFyg0KP/Nwj/80xcJuJxjnmSvYCxTa7EidrjirUyBRXmMLhHMWBbq1/k21kOy81cjpzRb2Al10QRXP4xpl9xLlNbpAxbnoAaekucGAw/i356hD+G2EwwIKLBC+7Tcabu3gAv66+c1NzyjhXaEFwhkfd4eT0hzAUze/96ZXlSFUzneqvtEyBBcI5/QJyRS6BjOAn25+z1t7jYjaEBc90sdez/lG4QXChWP87whdg/jhp5vfVmFvlTP5fdvqbmECAmG8yhUuzpPCTze/FRPUcyf0XLQ3CYHwQWc94+EBTOGnm9/ASPXcHGMOecme815LhC5B5PDTzW+CLCz29NXUM3HhFgm6W5iIQPi/Bl/kCl2DqOGpm999IYoRL2WtE/U0MBWBcHa/ZjDMMwG8dfObcTkpKUX7pqH8m2psTEUgLF9Y7bLQNYgYtgK9PhR/l/qlP/hfcicGH/pN60njE9V8VHywmE2V/LLJeU/pjQDdsBToihe1YRrxb2V688TkQ0f0vWMymzCa895T4LqGgbAUqIvDhpQtTh5BWw+fLenBwp0RSlCbiAjdLUxKIPy0ZZeXQtcgUlgK5LSImixCpZ3B3W+NAppSIL+mTXW3MC2BcM7gujeFrkGcsBQIxVOTHaU/oZEa0PgaFs0mjOZHF30HjEBJsBWIPq3zG4NHfPKm2P4sKoHwSQ/YETIAvgTCONkr7JGYBMJPW3cq6QYUQCf8CYQzezuKSiCcN8r3gtA1iA62AllaWVlZIkVX0aV/ZsvYa/reMkWBqKNHlxVClyA2WAo0UwOivKYpEL5Wpz9Zv0eSQyrPxjMlK6IOWRf8UgME0ma103qhSxATINB7XKnTH0YZYwwI9D5Zg/xgvF6mgEC62OgUBcNrMAME0smd5u1gwF5GgEC6yZvmflDoGkQBCKSPY1VGvRW6BhEAAunl5TA/6EeoVECgEtjs/F2B0DWYOiBQSdxrE0TYC4nZAwKViHyV00K4S6gkQKBSuN6s9W2hazBlQKDSKFjk+COcVdQLCFQ6aS2bpwpdg8kCAjGgcJnjMjgc0w0IxIgbIY3h+qpOQCBmyNe7TIE+FXUAAjHlYfcaR4WuwQQBgZizt0pP6AlGGxCIBa/HuayFI/rigECsON8kOEXoGkwLEIgd8vVuo6EjDw1AILakD/bcpGc79njvdn0dq5stIBB7/gxsofMR6PnOoR3dhkjsYB8EMoDCX1xHPHtv7c81qjtVdgkcIUBBAgICGUTGKOcVeVrrqrgcwfhSXStpPRsNAhlIauca24qvkSm6xL5pIa1rHiCQweyv1b7YsFcyejxWLLf4Q5hyBAIEMpy8H1yGaDw9ZjmWnm4uK63OFkEgEl5MdphT9Bx9W9c+ew6PcXCT1rlqEIiMO5+5rVTtTac4fRIa8qnbTmELMjYgECnnQ6vHK790bvRyrtT2tMDlGBsQiJwjgYHSHYceBOKCpHrNkoWuQSBAIE4oWOfd6ZzQRQgCCMQReas820pRIRCIM7Kj3LvrGaPajAGBOOTtMrdPzwtdhJEBgTjl7QqPT/4WugijAgJxTPaPVdofF7oIIwICcU7e2hot9kvmegZ/AuWmlzT2jRkLRB3Ub6lfd2O+0FUYB54Euje9qgyhcr5T9d0jbNYCURxo5R0tiTvL+BHovLXXiKgNcdEjfez1HNiau0AYn+nuNOOR0EXwDz8CtQpTdXCa37et7hbmLxDGNyLth5j9kPT8CFQxQT13wk53CykIhHH6PPfQRPPen+ZHoMBI9dycJrpbSEMg6lBifQO/n8x57BZ+BEqQhcWevpp6Ji7cIkF3C6kIRJHc3WF8mtBF8AZPR2H7QhCNrLW+sbQlJBDGd6c4dU40085eeTsPlHE5KSnludbKR6vUlPuJZTxx83ZtYLUl7z+LaAbwJVDBLeUjvtmPNddeGabGUnKj254dZBdxQugiuIcfgfJnlUflJ9P9Uq7T8zlJbcJUPF9Wq3a0uQ1Nz49Ai8tMSBhX5nMMAhVHfqyfXUSyWR3X8yNQjWnUZCPaDQK9x/MV/jW+eyB0FdzBj0A2++lpRNVsEEgHp4fad0rIFboKjuBHoEbj6ekT569AIJ28Wd/KaZR53HjGj0DRaPThHIz3WwyYAgLp5vbc6h8tvC90FeTwdBg/ryKiT77udUcgkD7kJ4c7hqx9IXQZhPB1Hij3Zg79kv97jO73QSCa3F3d7brvEHWneHBLq8BkxraxH3RQvLcvgkDC83BFkNOw30U6HBAIZBLcXdLI+cskMX4PgUCmwq3FjR2HHBDd6SEQyIS4832w/WcJ4roXHwQyLR7FhFb8JPap0GUwBwQyOV5s7mXXfOEVoctgCAhkiuQcjKxcbVySGHaIQCBT5cL8pnbhsQ9LbygsIJAJ8zSut0P9aX+Y9NE9CGTaFJyaEWAfHnNb6Dr0AgKZPk83DXCtHrnLNC+7gkCiQH5xSfsKQTOP5QhdyHuAQKIh+/cZTSu0W3DatHaJQCBR8XLPuPoVwxadMR2JQCDR8Wzn6HoV2397wjQ2ZyCQqEgb1rjtomyMM3aPb2z78cxE4XesQSAxsd+5R+++rWsrH058fXhWSIW6kRuFHSgaBBIR+S4e4VHf1qw+7t2as8t6uLuFLz35VqiaQCARca78WmqaG+RSfPWdTaMbWzcaueGaEI+8gkAiYr2V4mVn2fffyj61rLe3fYdZe4x98QwEEhHbLRVbql/K6Xn/8d6vwxw9u80/kG68mkAgEXHduv9rjC+51Sqp0a3tU9rYVe42d69xHsAHgcREkybOHZu7VtfzrN075GnbpoQ6urSfGp/K98MeIJCYuF2/ab8+HqMY9pZ3f+/c8Go2jYet/IPH00UgkKgo2DX7e6Z/MSUvT64c1tTWu8u0zSl5fFQEAkkBedrOuT39ytfu9U3CVY4vo4FA0iH34pbp3aqX9+85e+tFzp7HB4GkRs7FLbN6+pev1nHCL8kcjOUBAkmT/Ot7Fg9u7lip0Wdzt5x7SRAIBJI0z//8dXrP+tauLT7/Nv4fg47VQCCAOuA/tnpK9/oVHJv0nbk2+R6rPvVBIKCIJ6c2fDMg2L1cjfbDF8WffcLoMyAQoE321f0rJ/YIdCz/UdjwBRtP3C/xZDYIBOgj6/K+lVP6NvcoWyU4YuqPey7o7PIBBAJKI+/28bhvR3Su52DlGxi8JbP4myAQwJg7Ld2btXFZV2wdCAQwpuW0fIyvVvldcx0IBDDlkk8hzsF4XbjmShAIYMrOT9fXKm/X97C/5koQCGDKUa/AMwWZCxyKjaMMAgFMufvBYWoqrxGsuRIEApiyL8B5zpGtLQMbaq4EgQCm7Ol0Z2ybXuv/bqC5EgQCmPLAQXG9/psRmitBIIAx4z6+inN/dL2tuQ4EAhhTEO3maNMxtdg6/gTKTS9phGsQSJykaz/awZNA96ZXlSFUzneqvr5HQCAzgR+Bzlt7jYjaEBc90sf+ou4WIJCZwI9ArcJU/dXk922ruwUIZCbwI1DFBPXcCTvN9f/2VGO5nE08wGThR6DASPXcnGIXTp5vU1P5KJt4gMnCj0AJsrDY01dTz8SFWyTobtH0TzbxAJOFp6OwfSGIRtY6UU8DEMhM4O08UMblpKSU53rfBoHMBKHORINAZgIIBBABAgFEgEAAESAQQAQIBBAhmEBTVxWna78IDugRxkWUiLZ9uYjSqz0XUSLa9uEiSu9QLqJEhC7W+rt5CCTQ6mFaWFXz4wCPilxE8Svnw0UUL1suoviV9+YiSmUbLqL42YRq/d1GM+3ejGOB3qP6dS6ixPfiIgrz7+US2dOFiygcbe2PtOEiCg4x+BomCMQaEEgTEIg1IJAmIBBrQCBNQCDWgECagECsAYE0AYFYAwJpAgKxBgTShDdX4iwAAASRSURBVG+B/DgZvXpnPy6i4Ib/4yLKgU+5iIKD/+YiyvH2XETB7U4Y+km+BXrGSZR8boZW46aYgszS2zDgOSeDMcv13x3KBsOL4VsgwMwBgQAiQCCACBAIIAIEAogAgQAiQCCACBAIIAIEAogAgQAiQCCACBAIIIJfgRIaVQq5QBxlp6JToiFkQZJ2c1GRMgphQT80sa25JJ+0GHUUsmJej61qE7AdG14MrwLtk43Y3sHmHmmYpc4xFMlEMQobT+SgIlUUsoLmofH7ppaZRVhMURSyYiIqLE/8HCUZXgyvAoV0wPit13TSMJHEN03dX/kxmkhcUVEUooJyK46mphPKFxAV8y4KUTEvZHEYy2sOMvw3w6dAGSiWmg6vShqnwzDSCInBwVYTiStSRyEr6Caih+dKQLeIiimKQlbM9VZp1PTjXob/ZvgU6DI6TU2jZLmEcWq0b2hTbxVhEN+JXFSkiEJWUE5aDjUdVz6bqJiiKOS/HXli+Q2G/2b4FOgIukpNN6B0sjCFZR2ifhuClpJFUfzpiStSROGgoI1lJnHw66GjEBcTZYXGEvxm+BQoCdEjwMShkkZnYUDu1pvUdEDFQqIoij89cUWKKMQFPe2PBuUTF6OMQlzMrV2TLJcaXgyfAqWgM9Q0uhwnwXaiNKLPK/70xBUpN2GEBe13rvobeTGqKKTF0IyrZngxfAr0nN7Fx6OqEYZ5co6+43sPekwURfGnJ65IEYWwoP0WkdmYuBh1FLJitnekP70GvTG4GF4P41uHY5zvM4UwShLaRE2/rEwWRfndQVqRakNIUlC+e3/VHEkxRVHIiklEf1HTwZ6GF8OrQIkW35zsZ0/6ZFhBE+d5+0d/sJ0silIg0ooUUcgK+h1NXkeTTVRMURSyYvKa+aw7OPGDGMN/M/xeytjeuFIb8ksZb8fWqtDsAGEQ1d4LYUXKKEQFxSAlj4mKeReF7LfzakhN20D6O8zQYuBiKkAECAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEEAECAQQAQIBRIBAABEgEEAECAQQAQKVRg/VQ8RoPrZdI3QxpgcIVBp/JiQkeDWnJldxnyShizE9QCAm+PcRugKTBQRigkogO2oT5rk8yNrnp4ed7arQfVqsC7CuvVbY2gQGBGKCpkCWX//eWea5PLGx1Wscbfl14ljZTwJXJyggEBM0BeqJcSoag/E+dCnLYR61dqizsMUJCwjEBE2BFmFcgDbTFl08i04/e/ZsAyIezEHEgEBM0BRoKS1QgkKgeNUBforA5QkJCMQEPQIdRU8FLkx4QCAm6BHoaTl6gIlZxGPBiBkQiAl6BMJTys5PnCSLErg6QQGBmKBPIPlSf+uPYgQuTlhAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACJAIIAIEAggAgQCiACBACJAIICI/wPUXZ91xo7FDwAAAABJRU5ErkJggg==" /><!-- --></p> <pre class="r"><code>summary(m.L1.FOMC, data = FALSE)</code></pre> <pre><code>## Warning in sqrt(diag(covar)): NaNs produced</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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:39 2022 -## Date of summary: Fri Jul 8 15:44:40 2022 +## Date of fit: Wed Sep 14 22:28:35 2022 +## Date of summary: Wed Sep 14 22:28:35 2022 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted using 357 model solutions performed in 0.071 s +## Fitted using 369 model solutions performed in 0.081 s ## ## Error model: Constant variance ## @@ -1666,34 +1668,39 @@ plot(m.L1.FOMC, show_errmin = TRUE, main = "FOCUS L1 - FOMC")</code></ ## Fixed parameter values: ## None ## +## +## Warning(s): +## Optimisation did not converge: +## false convergence (8) +## ## Results: ## -## AIC BIC logLik -## 95.88804 99.44953 -43.94402 +## AIC BIC logLik +## 95.88781 99.44929 -43.9439 ## ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 92.47 1.2820 89.720 95.220 -## log_alpha 11.37 NaN NaN NaN -## log_beta 13.72 NaN NaN NaN -## sigma 2.78 0.4621 1.789 3.771 +## log_alpha 13.78 NaN NaN NaN +## log_beta 16.13 NaN NaN NaN +## sigma 2.78 0.4598 1.794 3.766 ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma -## parent_0 1.0000000 NaN NaN 0.0005548 +## parent_0 1.0000000 NaN NaN 0.0001671 ## log_alpha NaN 1 NaN NaN ## log_beta NaN NaN 1 NaN -## sigma 0.0005548 NaN NaN 1.0000000 +## sigma 0.0001671 NaN NaN 1.0000000 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. ## t-test (unrealistically) based on the assumption of normal distribution ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper -## parent_0 92.47 NA NA 89.720 95.220 -## alpha 87110.00 NA NA NA NA -## beta 911100.00 NA NA NA NA -## sigma 2.78 NA NA 1.789 3.771 +## parent_0 9.247e+01 NA NA 89.720 95.220 +## alpha 9.658e+05 NA NA NA NA +## beta 1.010e+07 NA NA NA NA +## sigma 2.780e+00 NA NA 1.794 3.766 ## ## FOCUS Chi2 error levels in percent: ## err.min n.optim df @@ -1701,8 +1708,8 @@ plot(m.L1.FOMC, show_errmin = TRUE, main = "FOCUS L1 - FOMC")</code></ ## parent 3.619 3 6 ## ## Estimated disappearance times: -## DT50 DT90 DT50back -## parent 7.249 24.08 7.249</code></pre> +## DT50 DT90 DT50back +## parent 7.25 24.08 7.25</code></pre> <p>We get a warning that the default optimisation algorithm <code>Port</code> did not converge, which is an indication that the model is overparameterised, <em>i.e.</em> contains too many parameters that are ill-defined as a consequence.</p> <p>And in fact, due to the higher number of parameters, and the lower number of degrees of freedom of the fit, the <span class="math inline"><em>χ</em><sup>2</sup></span> error level is actually higher for the FOMC model (3.6%) than for the SFO model (3.4%). Additionally, the parameters <code>log_alpha</code> and <code>log_beta</code> internally fitted in the model have excessive confidence intervals, that span more than 25 orders of magnitude (!) when backtransformed to the scale of <code>alpha</code> and <code>beta</code>. Also, the t-test for significant difference from zero does not indicate such a significant difference, with p-values greater than 0.1, and finally, the parameter correlation of <code>log_alpha</code> and <code>log_beta</code> is 1.000, clearly indicating that the model is overparameterised.</p> <p>The <span class="math inline"><em>χ</em><sup>2</sup></span> error levels reported in Appendix 3 and Appendix 7 to the FOCUS kinetics report are rounded to integer percentages and partly deviate by one percentage point from the results calculated by mkin. The reason for this is not known. However, mkin gives the same <span class="math inline"><em>χ</em><sup>2</sup></span> error levels as the kinfit package and the calculation routines of the kinfit package have been extensively compared to the results obtained by the KinGUI software, as documented in the kinfit package vignette. KinGUI was the first widely used standard package in this field. Also, the calculation of <span class="math inline"><em>χ</em><sup>2</sup></span> error levels was compared with KinGUII, CAKE and DegKin manager in a project sponsored by the German Umweltbundesamt <span class="citation">(Ranke 2014)</span>.</p> @@ -1735,17 +1742,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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:40 2022 -## Date of summary: Fri Jul 8 15:44:40 2022 +## Date of fit: Wed Sep 14 22:28:35 2022 +## Date of summary: Wed Sep 14 22:28:35 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.045 s +## Fitted using 239 model solutions performed in 0.049 s ## ## Error model: Constant variance ## @@ -1780,10 +1787,10 @@ plot(m.L2.FOMC, show_residuals = TRUE, ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma -## parent_0 1.000e+00 -1.151e-01 -2.085e-01 -7.637e-09 -## log_alpha -1.151e-01 1.000e+00 9.741e-01 -1.617e-07 -## log_beta -2.085e-01 9.741e-01 1.000e+00 -1.387e-07 -## sigma -7.637e-09 -1.617e-07 -1.387e-07 1.000e+00 +## parent_0 1.000e+00 -1.151e-01 -2.085e-01 -7.828e-09 +## log_alpha -1.151e-01 1.000e+00 9.741e-01 -1.602e-07 +## log_beta -2.085e-01 9.741e-01 1.000e+00 -1.372e-07 +## sigma -7.828e-09 -1.602e-07 -1.372e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. @@ -1813,10 +1820,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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:40 2022 -## Date of summary: Fri Jul 8 15:44:40 2022 +## Date of fit: Wed Sep 14 22:28:36 2022 +## Date of summary: Wed Sep 14 22:28:36 2022 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1825,7 +1832,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.119 s +## Fitted using 581 model solutions performed in 0.135 s ## ## Error model: Constant variance ## @@ -1856,18 +1863,18 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, ## Optimised, transformed parameters with symmetric confidence intervals: ## Estimate Std. Error Lower Upper ## parent_0 93.950 9.998e-01 91.5900 96.3100 -## log_k1 3.113 1.845e+03 -4360.0000 4367.0000 +## log_k1 3.112 1.842e+03 -4353.0000 4359.0000 ## log_k2 -1.088 6.285e-02 -1.2370 -0.9394 ## g_qlogis -0.399 9.946e-02 -0.6342 -0.1638 ## sigma 1.414 2.886e-01 0.7314 2.0960 ## ## Parameter correlation: ## parent_0 log_k1 log_k2 g_qlogis sigma -## parent_0 1.000e+00 6.784e-07 -5.188e-10 2.665e-01 -5.800e-10 -## log_k1 6.784e-07 1.000e+00 1.114e-04 -2.191e-04 -1.029e-05 -## log_k2 -5.188e-10 1.114e-04 1.000e+00 -7.903e-01 5.080e-09 -## g_qlogis 2.665e-01 -2.191e-04 -7.903e-01 1.000e+00 -7.991e-09 -## sigma -5.800e-10 -1.029e-05 5.080e-09 -7.991e-09 1.000e+00 +## parent_0 1.000e+00 6.783e-07 -3.390e-10 2.665e-01 -2.967e-10 +## log_k1 6.783e-07 1.000e+00 1.116e-04 -2.196e-04 -1.031e-05 +## log_k2 -3.390e-10 1.116e-04 1.000e+00 -7.903e-01 2.917e-09 +## g_qlogis 2.665e-01 -2.196e-04 -7.903e-01 1.000e+00 -4.408e-09 +## sigma -2.967e-10 -1.031e-05 2.917e-09 -4.408e-09 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. @@ -1875,7 +1882,7 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, ## for estimators of untransformed parameters. ## Estimate t value Pr(>t) Lower Upper ## parent_0 93.9500 9.397e+01 2.036e-12 91.5900 96.3100 -## k1 22.4800 5.544e-04 4.998e-01 0.0000 Inf +## k1 22.4800 5.553e-04 4.998e-01 0.0000 Inf ## k2 0.3369 1.591e+01 4.697e-07 0.2904 0.3909 ## g 0.4016 1.680e+01 3.238e-07 0.3466 0.4591 ## sigma 1.4140 4.899e+00 8.776e-04 0.7314 2.0960 @@ -1887,7 +1894,7 @@ plot(m.L2.DFOP, show_residuals = TRUE, show_errmin = TRUE, ## ## Estimated disappearance times: ## DT50 DT90 DT50back DT50_k1 DT50_k2 -## parent 0.5335 5.311 1.599 0.03083 2.058</code></pre> +## parent 0.5335 5.311 1.599 0.03084 2.058</code></pre> <p>Here, the DFOP model is clearly the best-fit model for dataset L2 based on the chi^2 error level criterion.</p> </div> </div> @@ -1913,10 +1920,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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:40 2022 -## Date of summary: Fri Jul 8 15:44:41 2022 +## Date of fit: Wed Sep 14 22:28:36 2022 +## Date of summary: Wed Sep 14 22:28:36 2022 ## ## Equations: ## d_parent/dt = - ((k1 * g * exp(-k1 * time) + k2 * (1 - g) * exp(-k2 * @@ -1925,7 +1932,7 @@ plot(mm.L3)</code></pre> ## ## Model predictions using solution type analytical ## -## Fitted using 376 model solutions performed in 0.072 s +## Fitted using 376 model solutions performed in 0.081 s ## ## Error model: Constant variance ## @@ -1963,11 +1970,11 @@ plot(mm.L3)</code></pre> ## ## Parameter correlation: ## parent_0 log_k1 log_k2 g_qlogis sigma -## parent_0 1.000e+00 1.732e-01 2.282e-02 4.009e-01 -9.632e-08 -## log_k1 1.732e-01 1.000e+00 4.945e-01 -5.809e-01 7.145e-07 -## log_k2 2.282e-02 4.945e-01 1.000e+00 -6.812e-01 1.021e-06 -## g_qlogis 4.009e-01 -5.809e-01 -6.812e-01 1.000e+00 -7.925e-07 -## sigma -9.632e-08 7.145e-07 1.021e-06 -7.925e-07 1.000e+00 +## parent_0 1.000e+00 1.732e-01 2.282e-02 4.009e-01 -9.664e-08 +## log_k1 1.732e-01 1.000e+00 4.945e-01 -5.809e-01 7.147e-07 +## log_k2 2.282e-02 4.945e-01 1.000e+00 -6.812e-01 1.022e-06 +## g_qlogis 4.009e-01 -5.809e-01 -6.812e-01 1.000e+00 -7.926e-07 +## sigma -9.664e-08 7.147e-07 1.022e-06 -7.926e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. @@ -2021,17 +2028,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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:41 2022 -## Date of summary: Fri Jul 8 15:44:41 2022 +## Date of fit: Wed Sep 14 22:28:36 2022 +## Date of summary: Wed Sep 14 22:28:37 2022 ## ## Equations: ## d_parent/dt = - k_parent * parent ## ## Model predictions using solution type analytical ## -## Fitted using 142 model solutions performed in 0.027 s +## Fitted using 142 model solutions performed in 0.034 s ## ## Error model: Constant variance ## @@ -2063,9 +2070,9 @@ plot(mm.L4)</code></pre> ## ## Parameter correlation: ## parent_0 log_k_parent sigma -## parent_0 1.000e+00 5.938e-01 3.440e-07 -## log_k_parent 5.938e-01 1.000e+00 5.885e-07 -## sigma 3.440e-07 5.885e-07 1.000e+00 +## parent_0 1.000e+00 5.938e-01 3.387e-07 +## log_k_parent 5.938e-01 1.000e+00 5.830e-07 +## sigma 3.387e-07 5.830e-07 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. @@ -2085,17 +2092,17 @@ 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.0 +<pre><code>## mkin version used for fitting: 1.1.2 ## R version used for fitting: 4.2.1 -## Date of fit: Fri Jul 8 15:44:41 2022 -## Date of summary: Fri Jul 8 15:44:41 2022 +## Date of fit: Wed Sep 14 22:28:37 2022 +## Date of summary: Wed Sep 14 22:28:37 2022 ## ## Equations: ## d_parent/dt = - (alpha/beta) * 1/((time/beta) + 1) * parent ## ## Model predictions using solution type analytical ## -## Fitted using 224 model solutions performed in 0.041 s +## Fitted using 224 model solutions performed in 0.045 s ## ## Error model: Constant variance ## @@ -2130,10 +2137,10 @@ plot(mm.L4)</code></pre> ## ## Parameter correlation: ## parent_0 log_alpha log_beta sigma -## parent_0 1.000e+00 -4.696e-01 -5.543e-01 -2.563e-07 -## log_alpha -4.696e-01 1.000e+00 9.889e-01 4.066e-08 -## log_beta -5.543e-01 9.889e-01 1.000e+00 6.818e-08 -## sigma -2.563e-07 4.066e-08 6.818e-08 1.000e+00 +## parent_0 1.000e+00 -4.696e-01 -5.543e-01 -2.468e-07 +## log_alpha -4.696e-01 1.000e+00 9.889e-01 2.478e-08 +## log_beta -5.543e-01 9.889e-01 1.000e+00 5.211e-08 +## sigma -2.468e-07 2.478e-08 5.211e-08 1.000e+00 ## ## Backtransformed parameters: ## Confidence intervals for internally transformed parameters are asymmetric. |