aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2019-10-26 22:52:53 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2019-10-26 22:52:53 +0200
commit6e5630a0df7e857697ff2ce4730a5f7f45b67377 (patch)
tree1d65d0b5c0540b1ae1ab160ed3c6cd1c50a4446f
parentcb03bf1fb1e6326f06029969663f3399af96b05d (diff)
Fix checks and tests
Static documentation rebuilt by pkgdown
-rw-r--r--DESCRIPTION3
-rw-r--r--NAMESPACE1
-rw-r--r--R/confint.mkinfit.R4
-rw-r--r--R/parms.mkinfit.R1
-rw-r--r--_pkgdown.yml12
-rw-r--r--docs/news/index.html1
-rw-r--r--docs/reference/confint.mkinfit.html252
-rw-r--r--docs/reference/index.html24
-rw-r--r--docs/reference/parms.html200
-rw-r--r--docs/sitemap.xml6
-rw-r--r--man/confint.mkinfit.Rd6
-rw-r--r--man/parms.Rd4
-rw-r--r--test.log18
-rw-r--r--tests/figs/plotting/plot-errmod-with-sfo-lin-a-obs.svg12
-rw-r--r--tests/testthat/FOCUS_2006_D.csf2
15 files changed, 513 insertions, 33 deletions
diff --git a/DESCRIPTION b/DESCRIPTION
index 9efa2b89..88412fa8 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -17,8 +17,7 @@ Description: Calculation routines based on the FOCUS Kinetics Report (2006,
equation models are solved using compiled C functions. Please note that no
warranty is implied for correctness of results or fitness for a particular
purpose.
-Imports: stats, graphics, methods, deSolve, R6, inline, parallel, plyr,
- numDeriv
+Imports: stats, graphics, methods, deSolve, R6, inline, parallel, numDeriv
Suggests: knitr, rbenchmark, tikzDevice, testthat, rmarkdown, covr, vdiffr,
benchmarkme, tibble
License: GPL
diff --git a/NAMESPACE b/NAMESPACE
index b4b781ac..3718f45f 100644
--- a/NAMESPACE
+++ b/NAMESPACE
@@ -71,6 +71,7 @@ importFrom(stats,optimize)
importFrom(stats,pt)
importFrom(stats,qchisq)
importFrom(stats,qf)
+importFrom(stats,qnorm)
importFrom(stats,qt)
importFrom(stats,rnorm)
importFrom(utils,write.table)
diff --git a/R/confint.mkinfit.R b/R/confint.mkinfit.R
index 58918321..58a3c8db 100644
--- a/R/confint.mkinfit.R
+++ b/R/confint.mkinfit.R
@@ -28,7 +28,9 @@
#' @examples
#' f <- mkinfit("SFO", FOCUS_2006_C, quiet = TRUE)
#' confint(f, method = "quadratic")
-#' confint(f, method = "profile")
+#' \dontrun{
+#' confint(f, method = "profile")
+#' }
#' @export
confint.mkinfit <- function(object, parm,
level = 0.95, alpha = 1 - level,
diff --git a/R/parms.mkinfit.R b/R/parms.mkinfit.R
index 8f9b4136..0628cb92 100644
--- a/R/parms.mkinfit.R
+++ b/R/parms.mkinfit.R
@@ -5,6 +5,7 @@
#' considering the error structure that was assumed for the fit.
#'
#' @param object A fitted model object
+#' @param \dots Not used
#' @return A numeric vector of fitted model parameters
#' @export
parms <- function(object, ...)
diff --git a/_pkgdown.yml b/_pkgdown.yml
index d51a7f27..bbf63301 100644
--- a/_pkgdown.yml
+++ b/_pkgdown.yml
@@ -17,13 +17,15 @@ reference:
desc: Functions working with mkinfit objects
contents:
- plot.mkinfit
+ - parms.mkinfit
+ - confint.mkinfit
- summary.mkinfit
- - mkinresplot
- - mkinparplot
- - mkinerrplot
- - endpoints
- - mkinerrmin
- logLik.mkinfit
+ - mkinresplot
+ - mkinparplot
+ - mkinerrplot
+ - endpoints
+ - mkinerrmin
- title: Work with mmkin objects
desc: Functions working with aggregated results
contents:
diff --git a/docs/news/index.html b/docs/news/index.html
index 67f64150..1322f75d 100644
--- a/docs/news/index.html
+++ b/docs/news/index.html
@@ -134,6 +134,7 @@
<a href="#mkin-0-9-49-6-unreleased" class="anchor"></a>mkin 0.9.49.6 (unreleased)<small> Unreleased </small>
</h1>
<ul>
+<li><p>Add ‘parms’ and ‘confint’ methods for mkinfit objects. Confidence intervals based on the quadratic approximation as in the summary, and based on the profile likelihood</p></li>
<li><p>Move long-running tests to tests/testthat/slow with a separate test log. They currently take around 7 minutes on my system</p></li>
<li><p>‘mkinfit’: Clean the code and return functions to calculate the log-likelihood and the sum of squared residuals</p></li>
<li><p>Vignette ‘twa.html’: Add the maximum time weighted average formulas for the hockey stick model</p></li>
diff --git a/docs/reference/confint.mkinfit.html b/docs/reference/confint.mkinfit.html
new file mode 100644
index 00000000..64ea0d4c
--- /dev/null
+++ b/docs/reference/confint.mkinfit.html
@@ -0,0 +1,252 @@
+<!-- Generated by pkgdown: do not edit by hand -->
+<!DOCTYPE html>
+<html lang="en">
+ <head>
+ <meta charset="utf-8">
+<meta http-equiv="X-UA-Compatible" content="IE=edge">
+<meta name="viewport" content="width=device-width, initial-scale=1.0">
+
+<title>Confidence intervals for parameters of mkinfit objects — confint.mkinfit • mkin</title>
+
+
+<!-- jquery -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.3.1/jquery.min.js" integrity="sha256-FgpCb/KJQlLNfOu91ta32o/NMZxltwRo8QtmkMRdAu8=" crossorigin="anonymous"></script>
+<!-- Bootstrap -->
+
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/css/bootstrap.min.css" integrity="sha256-916EbMg70RQy9LHiGkXzG8hSg9EdNy97GazNG/aiY1w=" crossorigin="anonymous" />
+
+<script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/js/bootstrap.min.js" integrity="sha256-U5ZEeKfGNOja007MMD3YBI0A3OSZOQbeG6z2f2Y0hu8=" crossorigin="anonymous"></script>
+
+<!-- Font Awesome icons -->
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.7.1/css/all.min.css" integrity="sha256-nAmazAk6vS34Xqo0BSrTb+abbtFlgsFK7NKSi6o7Y78=" crossorigin="anonymous" />
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.7.1/css/v4-shims.min.css" integrity="sha256-6qHlizsOWFskGlwVOKuns+D1nB6ssZrHQrNj1wGplHc=" crossorigin="anonymous" />
+
+<!-- clipboard.js -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.4/clipboard.min.js" integrity="sha256-FiZwavyI2V6+EXO1U+xzLG3IKldpiTFf3153ea9zikQ=" crossorigin="anonymous"></script>
+
+<!-- headroom.js -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.9.4/headroom.min.js" integrity="sha256-DJFC1kqIhelURkuza0AvYal5RxMtpzLjFhsnVIeuk+U=" crossorigin="anonymous"></script>
+<script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.9.4/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script>
+
+<!-- pkgdown -->
+<link href="../pkgdown.css" rel="stylesheet">
+<script src="../pkgdown.js"></script>
+
+
+
+
+<meta property="og:title" content="Confidence intervals for parameters of mkinfit objects — confint.mkinfit" />
+<meta property="og:description" content="Confidence intervals for parameters of mkinfit objects" />
+<meta name="twitter:card" content="summary" />
+
+
+
+
+<!-- mathjax -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script>
+<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script>
+
+<!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]-->
+
+
+
+ </head>
+
+ <body>
+ <div class="container template-reference-topic">
+ <header>
+ <div class="navbar navbar-default navbar-fixed-top" role="navigation">
+ <div class="container">
+ <div class="navbar-header">
+ <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
+ <span class="sr-only">Toggle navigation</span>
+ <span class="icon-bar"></span>
+ <span class="icon-bar"></span>
+ <span class="icon-bar"></span>
+ </button>
+ <span class="navbar-brand">
+ <a class="navbar-link" href="../index.html">mkin</a>
+ <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">0.9.49.6</span>
+ </span>
+ </div>
+
+ <div id="navbar" class="navbar-collapse collapse">
+ <ul class="nav navbar-nav">
+ <li>
+ <a href="../reference/index.html">Functions and data</a>
+</li>
+<li class="dropdown">
+ <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" aria-expanded="false">
+ Articles
+
+ <span class="caret"></span>
+ </a>
+ <ul class="dropdown-menu" role="menu">
+ <li>
+ <a href="../articles/mkin.html">Introduction to mkin</a>
+ </li>
+ <li>
+ <a href="../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a>
+ </li>
+ <li>
+ <a href="../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a>
+ </li>
+ <li>
+ <a href="../articles/twa.html">Calculation of time weighted average concentrations with mkin</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a>
+ </li>
+ </ul>
+</li>
+<li>
+ <a href="../news/index.html">News</a>
+</li>
+ </ul>
+ <ul class="nav navbar-nav navbar-right">
+
+ </ul>
+
+ </div><!--/.nav-collapse -->
+ </div><!--/.container -->
+</div><!--/.navbar -->
+
+
+
+ </header>
+
+<div class="row">
+ <div class="col-md-9 contents">
+ <div class="page-header">
+ <h1>Confidence intervals for parameters of mkinfit objects</h1>
+
+ <div class="hidden name"><code>confint.mkinfit.Rd</code></div>
+ </div>
+
+ <div class="ref-description">
+ <p>Confidence intervals for parameters of mkinfit objects</p>
+ </div>
+
+ <pre class="usage"><span class='co'># S3 method for mkinfit</span>
+<span class='fu'><a href='https://rdrr.io/r/stats/confint.html'>confint</a></span>(<span class='no'>object</span>, <span class='no'>parm</span>, <span class='kw'>level</span> <span class='kw'>=</span> <span class='fl'>0.95</span>, <span class='kw'>alpha</span> <span class='kw'>=</span> <span class='fl'>1</span> -
+ <span class='no'>level</span>, <span class='kw'>method</span> <span class='kw'>=</span> <span class='fu'><a href='https://rdrr.io/r/base/c.html'>c</a></span>(<span class='st'>"profile"</span>, <span class='st'>"quadratic"</span>), <span class='kw'>transformed</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>,
+ <span class='kw'>backtransform</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>, <span class='kw'>distribution</span> <span class='kw'>=</span> <span class='fu'><a href='https://rdrr.io/r/base/c.html'>c</a></span>(<span class='st'>"student_t"</span>, <span class='st'>"normal"</span>),
+ <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='no'>...</span>)</pre>
+
+ <h2 class="hasAnchor" id="arguments"><a class="anchor" href="#arguments"></a>Arguments</h2>
+ <table class="ref-arguments">
+ <colgroup><col class="name" /><col class="desc" /></colgroup>
+ <tr>
+ <th>object</th>
+ <td><p>An <code><a href='mkinfit.html'>mkinfit</a></code> object</p></td>
+ </tr>
+ <tr>
+ <th>parm</th>
+ <td><p>A vector of names of the parameters which are to be given
+confidence intervals. If missing, all parameters are considered.</p></td>
+ </tr>
+ <tr>
+ <th>level</th>
+ <td><p>The confidence level required</p></td>
+ </tr>
+ <tr>
+ <th>alpha</th>
+ <td><p>The allowed error probability, overrides 'level' if specified.</p></td>
+ </tr>
+ <tr>
+ <th>method</th>
+ <td><p>The 'profile' method searches the parameter space for the
+cutoff of the confidence intervals by means of a likelihood ratio test.
+The 'quadratic' method approximates the likelihood function at the
+optimised parameters using the second term of the Taylor expansion, using
+a second derivative (hessian) contained in the object.</p></td>
+ </tr>
+ <tr>
+ <th>transformed</th>
+ <td><p>If the quadratic approximation is used, should it be
+applied to the likelihood based on the transformed parameters?</p></td>
+ </tr>
+ <tr>
+ <th>backtransform</th>
+ <td><p>If we approximate the likelihood in terms of the
+transformed parameters, should we backtransform the parameters with
+their confidence intervals?</p></td>
+ </tr>
+ <tr>
+ <th>distribution</th>
+ <td><p>For the quadratic approximation, should we use
+the student t distribution or assume normal distribution for
+the parameter estimate</p></td>
+ </tr>
+ <tr>
+ <th>quiet</th>
+ <td><p>Should we suppress messages?</p></td>
+ </tr>
+ <tr>
+ <th>...</th>
+ <td><p>Not used</p></td>
+ </tr>
+ </table>
+
+ <h2 class="hasAnchor" id="value"><a class="anchor" href="#value"></a>Value</h2>
+
+ <p>A matrix with columns giving lower and upper confidence limits for
+ each parameter.</p>
+ <h2 class="hasAnchor" id="references"><a class="anchor" href="#references"></a>References</h2>
+
+ <p>Pawitan Y (2013) In all likelihood - Statistical modelling and
+ inference using likelihood. Clarendon Press, Oxford.</p>
+
+ <h2 class="hasAnchor" id="examples"><a class="anchor" href="#examples"></a>Examples</h2>
+ <pre class="examples"><div class='input'><span class='no'>f</span> <span class='kw'>&lt;-</span> <span class='fu'><a href='mkinfit.html'>mkinfit</a></span>(<span class='st'>"SFO"</span>, <span class='no'>FOCUS_2006_C</span>, <span class='kw'>quiet</span> <span class='kw'>=</span> <span class='fl'>TRUE</span>)
+<span class='fu'><a href='https://rdrr.io/r/stats/confint.html'>confint</a></span>(<span class='no'>f</span>, <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"quadratic"</span>)</div><div class='output co'>#&gt; 2.5% 97.5%
+#&gt; parent_0 71.8242430 93.1600766
+#&gt; k_parent_sink 0.2109541 0.4440528
+#&gt; sigma 1.9778868 7.3681380</div><div class='input'><span class='co'># \dontrun{</span>
+ <span class='fu'><a href='https://rdrr.io/r/stats/confint.html'>confint</a></span>(<span class='no'>f</span>, <span class='kw'>method</span> <span class='kw'>=</span> <span class='st'>"profile"</span>)</div><div class='output co'>#&gt; <span class='message'>Profiling the likelihood</span></div><div class='output co'>#&gt; 2.5% 97.5%
+#&gt; parent_0 71.3471007 93.9447024
+#&gt; k_parent_sink 0.2030765 0.4491067
+#&gt; sigma 2.9810656 8.8633278</div><div class='input'># }
+</div></pre>
+ </div>
+ <div class="col-md-3 hidden-xs hidden-sm" id="sidebar">
+ <h2>Contents</h2>
+ <ul class="nav nav-pills nav-stacked">
+ <li><a href="#arguments">Arguments</a></li>
+ <li><a href="#value">Value</a></li>
+ <li><a href="#references">References</a></li>
+ <li><a href="#examples">Examples</a></li>
+ </ul>
+
+ </div>
+</div>
+
+
+ <footer>
+ <div class="copyright">
+ <p>Developed by Johannes Ranke.</p>
+</div>
+
+<div class="pkgdown">
+ <p>Site built with <a href="https://pkgdown.r-lib.org/">pkgdown</a> 1.4.1.</p>
+</div>
+
+ </footer>
+ </div>
+
+
+
+
+ </body>
+</html>
+
+
diff --git a/docs/reference/index.html b/docs/reference/index.html
index 328f5aca..fa424060 100644
--- a/docs/reference/index.html
+++ b/docs/reference/index.html
@@ -185,12 +185,30 @@ more datasets</p></td>
</tr><tr>
<td>
+ <p><code><a href="parms.html">parms()</a></code> </p>
+ </td>
+ <td><p>Extract model parameters from mkinfit models</p></td>
+ </tr><tr>
+
+ <td>
+ <p><code><a href="confint.mkinfit.html">confint(<i>&lt;mkinfit&gt;</i>)</a></code> </p>
+ </td>
+ <td><p>Confidence intervals for parameters of mkinfit objects</p></td>
+ </tr><tr>
+
+ <td>
<p><code><a href="summary.mkinfit.html">summary(<i>&lt;mkinfit&gt;</i>)</a></code> <code><a href="summary.mkinfit.html">print(<i>&lt;summary.mkinfit&gt;</i>)</a></code> </p>
</td>
<td><p>Summary method for class "mkinfit"</p></td>
</tr><tr>
<td>
+ <p><code><a href="logLik.mkinfit.html">logLik(<i>&lt;mkinfit&gt;</i>)</a></code> </p>
+ </td>
+ <td><p>Calculated the log-likelihood of a fitted mkinfit object</p></td>
+ </tr><tr>
+
+ <td>
<p><code><a href="mkinresplot.html">mkinresplot()</a></code> </p>
</td>
<td><p>Function to plot residuals stored in an mkin object</p></td>
@@ -219,12 +237,6 @@ with mkinfit</p></td>
<p><code><a href="mkinerrmin.html">mkinerrmin()</a></code> </p>
</td>
<td><p>Calculate the minimum error to assume in order to pass the variance test</p></td>
- </tr><tr>
-
- <td>
- <p><code><a href="logLik.mkinfit.html">logLik(<i>&lt;mkinfit&gt;</i>)</a></code> </p>
- </td>
- <td><p>Calculated the log-likelihood of a fitted mkinfit object</p></td>
</tr>
</tbody><tbody>
<tr>
diff --git a/docs/reference/parms.html b/docs/reference/parms.html
new file mode 100644
index 00000000..8ab26240
--- /dev/null
+++ b/docs/reference/parms.html
@@ -0,0 +1,200 @@
+<!-- Generated by pkgdown: do not edit by hand -->
+<!DOCTYPE html>
+<html lang="en">
+ <head>
+ <meta charset="utf-8">
+<meta http-equiv="X-UA-Compatible" content="IE=edge">
+<meta name="viewport" content="width=device-width, initial-scale=1.0">
+
+<title>Extract model parameters from mkinfit models — parms • mkin</title>
+
+
+<!-- jquery -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/jquery/3.3.1/jquery.min.js" integrity="sha256-FgpCb/KJQlLNfOu91ta32o/NMZxltwRo8QtmkMRdAu8=" crossorigin="anonymous"></script>
+<!-- Bootstrap -->
+
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/css/bootstrap.min.css" integrity="sha256-916EbMg70RQy9LHiGkXzG8hSg9EdNy97GazNG/aiY1w=" crossorigin="anonymous" />
+
+<script src="https://cdnjs.cloudflare.com/ajax/libs/twitter-bootstrap/3.3.7/js/bootstrap.min.js" integrity="sha256-U5ZEeKfGNOja007MMD3YBI0A3OSZOQbeG6z2f2Y0hu8=" crossorigin="anonymous"></script>
+
+<!-- Font Awesome icons -->
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.7.1/css/all.min.css" integrity="sha256-nAmazAk6vS34Xqo0BSrTb+abbtFlgsFK7NKSi6o7Y78=" crossorigin="anonymous" />
+<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/5.7.1/css/v4-shims.min.css" integrity="sha256-6qHlizsOWFskGlwVOKuns+D1nB6ssZrHQrNj1wGplHc=" crossorigin="anonymous" />
+
+<!-- clipboard.js -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/clipboard.js/2.0.4/clipboard.min.js" integrity="sha256-FiZwavyI2V6+EXO1U+xzLG3IKldpiTFf3153ea9zikQ=" crossorigin="anonymous"></script>
+
+<!-- headroom.js -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.9.4/headroom.min.js" integrity="sha256-DJFC1kqIhelURkuza0AvYal5RxMtpzLjFhsnVIeuk+U=" crossorigin="anonymous"></script>
+<script src="https://cdnjs.cloudflare.com/ajax/libs/headroom/0.9.4/jQuery.headroom.min.js" integrity="sha256-ZX/yNShbjqsohH1k95liqY9Gd8uOiE1S4vZc+9KQ1K4=" crossorigin="anonymous"></script>
+
+<!-- pkgdown -->
+<link href="../pkgdown.css" rel="stylesheet">
+<script src="../pkgdown.js"></script>
+
+
+
+
+<meta property="og:title" content="Extract model parameters from mkinfit models — parms" />
+<meta property="og:description" content="This function always returns degradation model parameters as well as error
+model parameters, in order to avoid working with a fitted model without
+considering the error structure that was assumed for the fit." />
+<meta name="twitter:card" content="summary" />
+
+
+
+
+<!-- mathjax -->
+<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/MathJax.js" integrity="sha256-nvJJv9wWKEm88qvoQl9ekL2J+k/RWIsaSScxxlsrv8k=" crossorigin="anonymous"></script>
+<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.5/config/TeX-AMS-MML_HTMLorMML.js" integrity="sha256-84DKXVJXs0/F8OTMzX4UR909+jtl4G7SPypPavF+GfA=" crossorigin="anonymous"></script>
+
+<!--[if lt IE 9]>
+<script src="https://oss.maxcdn.com/html5shiv/3.7.3/html5shiv.min.js"></script>
+<script src="https://oss.maxcdn.com/respond/1.4.2/respond.min.js"></script>
+<![endif]-->
+
+
+
+ </head>
+
+ <body>
+ <div class="container template-reference-topic">
+ <header>
+ <div class="navbar navbar-default navbar-fixed-top" role="navigation">
+ <div class="container">
+ <div class="navbar-header">
+ <button type="button" class="navbar-toggle collapsed" data-toggle="collapse" data-target="#navbar" aria-expanded="false">
+ <span class="sr-only">Toggle navigation</span>
+ <span class="icon-bar"></span>
+ <span class="icon-bar"></span>
+ <span class="icon-bar"></span>
+ </button>
+ <span class="navbar-brand">
+ <a class="navbar-link" href="../index.html">mkin</a>
+ <span class="version label label-default" data-toggle="tooltip" data-placement="bottom" title="Released version">0.9.49.6</span>
+ </span>
+ </div>
+
+ <div id="navbar" class="navbar-collapse collapse">
+ <ul class="nav navbar-nav">
+ <li>
+ <a href="../reference/index.html">Functions and data</a>
+</li>
+<li class="dropdown">
+ <a href="#" class="dropdown-toggle" data-toggle="dropdown" role="button" aria-expanded="false">
+ Articles
+
+ <span class="caret"></span>
+ </a>
+ <ul class="dropdown-menu" role="menu">
+ <li>
+ <a href="../articles/mkin.html">Introduction to mkin</a>
+ </li>
+ <li>
+ <a href="../articles/FOCUS_D.html">Example evaluation of FOCUS Example Dataset D</a>
+ </li>
+ <li>
+ <a href="../articles/FOCUS_L.html">Example evaluation of FOCUS Laboratory Data L1 to L3</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/FOCUS_Z.html">Example evaluation of FOCUS Example Dataset Z</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/compiled_models.html">Performance benefit by using compiled model definitions in mkin</a>
+ </li>
+ <li>
+ <a href="../articles/twa.html">Calculation of time weighted average concentrations with mkin</a>
+ </li>
+ <li>
+ <a href="../articles/web_only/NAFTA_examples.html">Example evaluation of NAFTA SOP Attachment examples</a>
+ </li>
+ </ul>
+</li>
+<li>
+ <a href="../news/index.html">News</a>
+</li>
+ </ul>
+ <ul class="nav navbar-nav navbar-right">
+
+ </ul>
+
+ </div><!--/.nav-collapse -->
+ </div><!--/.container -->
+</div><!--/.navbar -->
+
+
+
+ </header>
+
+<div class="row">
+ <div class="col-md-9 contents">
+ <div class="page-header">
+ <h1>Extract model parameters from mkinfit models</h1>
+
+ <div class="hidden name"><code>parms.Rd</code></div>
+ </div>
+
+ <div class="ref-description">
+ <p>This function always returns degradation model parameters as well as error
+model parameters, in order to avoid working with a fitted model without
+considering the error structure that was assumed for the fit.</p>
+ </div>
+
+ <pre class="usage"><span class='fu'>parms</span>(<span class='no'>object</span>, <span class='no'>...</span>)
+
+<span class='co'># S3 method for mkinfit</span>
+<span class='fu'>parms</span>(<span class='no'>object</span>, <span class='kw'>transformed</span> <span class='kw'>=</span> <span class='fl'>FALSE</span>, <span class='no'>...</span>)</pre>
+
+ <h2 class="hasAnchor" id="arguments"><a class="anchor" href="#arguments"></a>Arguments</h2>
+ <table class="ref-arguments">
+ <colgroup><col class="name" /><col class="desc" /></colgroup>
+ <tr>
+ <th>object</th>
+ <td><p>A fitted model object</p></td>
+ </tr>
+ <tr>
+ <th>...</th>
+ <td><p>Not used</p></td>
+ </tr>
+ <tr>
+ <th>transformed</th>
+ <td><p>Should the parameters be returned
+as used internally during the optimisation?</p></td>
+ </tr>
+ </table>
+
+ <h2 class="hasAnchor" id="value"><a class="anchor" href="#value"></a>Value</h2>
+
+ <p>A numeric vector of fitted model parameters</p>
+
+ </div>
+ <div class="col-md-3 hidden-xs hidden-sm" id="sidebar">
+ <h2>Contents</h2>
+ <ul class="nav nav-pills nav-stacked">
+ <li><a href="#arguments">Arguments</a></li>
+ <li><a href="#value">Value</a></li>
+ </ul>
+
+ </div>
+</div>
+
+
+ <footer>
+ <div class="copyright">
+ <p>Developed by Johannes Ranke.</p>
+</div>
+
+<div class="pkgdown">
+ <p>Site built with <a href="https://pkgdown.r-lib.org/">pkgdown</a> 1.4.1.</p>
+</div>
+
+ </footer>
+ </div>
+
+
+
+
+ </body>
+</html>
+
+
diff --git a/docs/sitemap.xml b/docs/sitemap.xml
index b62b9a44..9e3363ba 100644
--- a/docs/sitemap.xml
+++ b/docs/sitemap.xml
@@ -55,6 +55,9 @@
<loc>https://pkgdown.jrwb.de/mkin/reference/add_err.html</loc>
</url>
<url>
+ <loc>https://pkgdown.jrwb.de/mkin/reference/confint.mkinfit.html</loc>
+ </url>
+ <url>
<loc>https://pkgdown.jrwb.de/mkin/reference/endpoints.html</loc>
</url>
<url>
@@ -118,6 +121,9 @@
<loc>https://pkgdown.jrwb.de/mkin/reference/nafta.html</loc>
</url>
<url>
+ <loc>https://pkgdown.jrwb.de/mkin/reference/parms.html</loc>
+ </url>
+ <url>
<loc>https://pkgdown.jrwb.de/mkin/reference/plot.mkinfit.html</loc>
</url>
<url>
diff --git a/man/confint.mkinfit.Rd b/man/confint.mkinfit.Rd
index 94f55d14..943904b9 100644
--- a/man/confint.mkinfit.Rd
+++ b/man/confint.mkinfit.Rd
@@ -37,6 +37,8 @@ the student t distribution or assume normal distribution for
the parameter estimate}
\item{quiet}{Should we suppress messages?}
+
+\item{\dots}{Not used}
}
\value{
A matrix with columns giving lower and upper confidence limits for
@@ -48,7 +50,9 @@ Confidence intervals for parameters of mkinfit objects
\examples{
f <- mkinfit("SFO", FOCUS_2006_C, quiet = TRUE)
confint(f, method = "quadratic")
-confint(f, method = "profile")
+\dontrun{
+ confint(f, method = "profile")
+}
}
\references{
Pawitan Y (2013) In all likelihood - Statistical modelling and
diff --git a/man/parms.Rd b/man/parms.Rd
index 6de52557..73cb23cd 100644
--- a/man/parms.Rd
+++ b/man/parms.Rd
@@ -12,10 +12,10 @@ parms(object, ...)
\arguments{
\item{object}{A fitted model object}
+\item{\dots}{Not used}
+
\item{transformed}{Should the parameters be returned
as used internally during the optimisation?}
-
-\item{complete}{Unused argument for compatibility with the generic coef function from base R}
}
\value{
A numeric vector of fitted model parameters
diff --git a/test.log b/test.log
index 7e3f1c4e..bc3b3c8c 100644
--- a/test.log
+++ b/test.log
@@ -2,30 +2,30 @@ Loading mkin
Testing mkin
✔ | OK F W S | Context
✔ | 2 | Export dataset for reading into CAKE
-✔ | 3 | Confidence intervals and p-values
-✔ | 10 | Error model fitting [37.1 s]
+✔ | 9 | Confidence intervals and p-values [2.3 s]
+✔ | 10 | Error model fitting [37.2 s]
✔ | 5 | Calculation of FOCUS chi2 error levels [3.5 s]
-✔ | 13 | Results for FOCUS D established in expertise for UBA (Ranke 2014) [3.5 s]
+✔ | 13 | Results for FOCUS D established in expertise for UBA (Ranke 2014) [3.3 s]
✔ | 6 | Test fitting the decline of metabolites from their maximum [0.9 s]
✔ | 1 | Fitting the logistic model [0.9 s]
✔ | 1 | Test dataset class mkinds used in gmkin
✔ | 12 | Special cases of mkinfit calls [2.7 s]
✔ | 9 | mkinmod model generation and printing [0.2 s]
✔ | 3 | Model predictions with mkinpredict [0.3 s]
-✔ | 16 | Evaluations according to 2015 NAFTA guidance [4.1 s]
-✔ | 4 | Calculation of maximum time weighted average concentrations (TWAs) [2.2 s]
+✔ | 16 | Evaluations according to 2015 NAFTA guidance [4.0 s]
+✔ | 4 | Calculation of maximum time weighted average concentrations (TWAs) [2.3 s]
✔ | 3 | Summary
✔ | 11 | Plotting [0.6 s]
✔ | 3 | AIC calculation
-✔ | 2 | Complex test case from Schaefer et al. (2007) Piacenza paper [5.3 s]
+✔ | 2 | Complex test case from Schaefer et al. (2007) Piacenza paper [5.2 s]
✔ | 4 | Fitting the SFORB model [1.7 s]
✔ | 1 | Summaries of old mkinfit objects
-✔ | 4 | Results for synthetic data established in expertise for UBA (Ranke 2014) [7.2 s]
+✔ | 4 | Results for synthetic data established in expertise for UBA (Ranke 2014) [7.0 s]
══ Results ═════════════════════════════════════════════════════════════════════
-Duration: 70.2 s
+Duration: 72.1 s
-OK: 113
+OK: 119
Failed: 0
Warnings: 0
Skipped: 0
diff --git a/tests/figs/plotting/plot-errmod-with-sfo-lin-a-obs.svg b/tests/figs/plotting/plot-errmod-with-sfo-lin-a-obs.svg
index 920f70cd..2b4930ba 100644
--- a/tests/figs/plotting/plot-errmod-with-sfo-lin-a-obs.svg
+++ b/tests/figs/plotting/plot-errmod-with-sfo-lin-a-obs.svg
@@ -108,7 +108,7 @@
<line x1='315.91' y1='472.85' x2='323.55' y2='472.85' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
<line x1='319.73' y1='476.67' x2='319.73' y2='469.04' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
<polyline points='69.07,75.47 71.60,318.44 74.13,417.65 76.66,458.17 79.19,474.71 81.73,481.47 84.26,484.23 86.79,485.36 89.32,485.82 91.85,486.00 94.39,486.08 96.92,486.11 99.45,486.12 101.98,486.13 104.51,486.13 107.05,486.13 109.58,486.13 112.11,486.13 114.64,486.13 117.17,486.13 119.71,486.13 122.24,486.13 124.77,486.13 127.30,486.13 129.83,486.13 132.37,486.13 134.90,486.13 137.43,486.13 139.96,486.13 142.49,486.13 145.03,486.13 147.56,486.13 150.09,486.13 152.62,486.13 155.15,486.13 157.69,486.13 160.22,486.13 162.75,486.13 165.28,486.13 167.81,486.13 170.35,486.13 172.88,486.13 175.41,486.13 177.94,486.13 180.47,486.13 183.01,486.13 185.54,486.13 188.07,486.13 190.60,486.13 193.13,486.13 195.67,486.13 198.20,486.13 200.73,486.13 203.26,486.13 205.79,486.13 208.33,486.13 210.86,486.13 213.39,486.13 215.92,486.13 218.45,486.13 220.99,486.13 223.52,486.13 226.05,486.13 228.58,486.13 231.11,486.13 233.65,486.13 236.18,486.13 238.71,486.13 241.24,486.13 243.77,486.13 246.31,486.13 248.84,486.13 251.37,486.13 253.90,486.13 256.43,486.13 258.97,486.13 261.50,486.13 264.03,486.13 266.56,486.13 269.09,486.13 271.63,486.13 274.16,486.13 276.69,486.13 279.22,486.13 281.75,486.13 284.29,486.13 286.82,486.13 289.35,486.13 291.88,486.13 294.41,486.13 296.95,486.13 299.48,486.13 302.01,486.13 304.54,486.13 307.07,486.13 309.61,486.13 312.14,486.13 314.67,486.13 317.20,486.13 319.73,486.13 ' style='stroke-width: 0.75;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
-<polyline points='69.07,486.13 71.60,333.73 74.13,317.72 76.66,343.39 79.19,376.31 81.73,405.38 84.26,428.15 86.79,445.03 89.32,457.21 91.85,465.86 94.39,471.96 96.92,476.24 99.45,479.23 101.98,481.32 104.51,482.78 107.05,483.80 109.58,484.51 112.11,485.00 114.64,485.34 117.17,485.58 119.71,485.75 122.24,485.87 124.77,485.95 127.30,486.00 129.83,486.04 132.37,486.07 134.90,486.09 137.43,486.10 139.96,486.11 142.49,486.12 145.03,486.12 147.56,486.13 150.09,486.13 152.62,486.13 155.15,486.13 157.69,486.13 160.22,486.13 162.75,486.13 165.28,486.13 167.81,486.13 170.35,486.13 172.88,486.13 175.41,486.13 177.94,486.13 180.47,486.13 183.01,486.13 185.54,486.13 188.07,486.13 190.60,486.13 193.13,486.13 195.67,486.13 198.20,486.13 200.73,486.13 203.26,486.13 205.79,486.13 208.33,486.13 210.86,486.13 213.39,486.13 215.92,486.13 218.45,486.13 220.99,486.13 223.52,486.13 226.05,486.13 228.58,486.13 231.11,486.13 233.65,486.13 236.18,486.13 238.71,486.13 241.24,486.13 243.77,486.13 246.31,486.13 248.84,486.13 251.37,486.13 253.90,486.13 256.43,486.13 258.97,486.13 261.50,486.13 264.03,486.13 266.56,486.13 269.09,486.13 271.63,486.13 274.16,486.13 276.69,486.13 279.22,486.13 281.75,486.13 284.29,486.13 286.82,486.13 289.35,486.13 291.88,486.13 294.41,486.13 296.95,486.13 299.48,486.13 302.01,486.13 304.54,486.13 307.07,486.13 309.61,486.13 312.14,486.13 314.67,486.13 317.20,486.13 319.73,486.13 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
+<polyline points='69.07,486.13 71.60,333.73 74.13,317.73 76.66,343.39 79.19,376.31 81.73,405.38 84.26,428.15 86.79,445.03 89.32,457.21 91.85,465.86 94.39,471.96 96.92,476.24 99.45,479.23 101.98,481.32 104.51,482.78 107.05,483.80 109.58,484.51 112.11,485.00 114.64,485.34 117.17,485.58 119.71,485.75 122.24,485.87 124.77,485.95 127.30,486.00 129.83,486.04 132.37,486.07 134.90,486.09 137.43,486.10 139.96,486.11 142.49,486.12 145.03,486.12 147.56,486.13 150.09,486.13 152.62,486.13 155.15,486.13 157.69,486.13 160.22,486.13 162.75,486.13 165.28,486.13 167.81,486.13 170.35,486.13 172.88,486.13 175.41,486.13 177.94,486.13 180.47,486.13 183.01,486.13 185.54,486.13 188.07,486.13 190.60,486.13 193.13,486.13 195.67,486.13 198.20,486.13 200.73,486.13 203.26,486.13 205.79,486.13 208.33,486.13 210.86,486.13 213.39,486.13 215.92,486.13 218.45,486.13 220.99,486.13 223.52,486.13 226.05,486.13 228.58,486.13 231.11,486.13 233.65,486.13 236.18,486.13 238.71,486.13 241.24,486.13 243.77,486.13 246.31,486.13 248.84,486.13 251.37,486.13 253.90,486.13 256.43,486.13 258.97,486.13 261.50,486.13 264.03,486.13 266.56,486.13 269.09,486.13 271.63,486.13 274.16,486.13 276.69,486.13 279.22,486.13 281.75,486.13 284.29,486.13 286.82,486.13 289.35,486.13 291.88,486.13 294.41,486.13 296.95,486.13 299.48,486.13 302.01,486.13 304.54,486.13 307.07,486.13 309.61,486.13 312.14,486.13 314.67,486.13 317.20,486.13 319.73,486.13 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
<polyline points='69.07,486.13 71.60,461.63 74.13,419.15 76.66,380.16 79.19,350.09 81.73,328.88 84.26,314.90 86.79,306.36 89.32,301.75 91.85,299.93 94.39,300.04 96.92,301.48 99.45,303.81 101.98,306.73 104.51,310.04 107.05,313.59 109.58,317.26 112.11,321.01 114.64,324.76 117.17,328.51 119.71,332.21 122.24,335.86 124.77,339.45 127.30,342.96 129.83,346.41 132.37,349.78 134.90,353.07 137.43,356.29 139.96,359.43 142.49,362.50 145.03,365.49 147.56,368.42 150.09,371.27 152.62,374.06 155.15,376.77 157.69,379.43 160.22,382.01 162.75,384.54 165.28,387.00 167.81,389.41 170.35,391.75 172.88,394.04 175.41,396.27 177.94,398.45 180.47,400.58 183.01,402.66 185.54,404.68 188.07,406.66 190.60,408.58 193.13,410.46 195.67,412.30 198.20,414.09 200.73,415.84 203.26,417.54 205.79,419.21 208.33,420.83 210.86,422.41 213.39,423.96 215.92,425.47 218.45,426.94 220.99,428.37 223.52,429.77 226.05,431.14 228.58,432.47 231.11,433.78 233.65,435.05 236.18,436.28 238.71,437.49 241.24,438.67 243.77,439.82 246.31,440.95 248.84,442.04 251.37,443.11 253.90,444.16 256.43,445.17 258.97,446.17 261.50,447.14 264.03,448.08 266.56,449.01 269.09,449.91 271.63,450.78 274.16,451.64 276.69,452.48 279.22,453.29 281.75,454.09 284.29,454.87 286.82,455.63 289.35,456.37 291.88,457.09 294.41,457.79 296.95,458.48 299.48,459.15 302.01,459.81 304.54,460.44 307.07,461.07 309.61,461.67 312.14,462.27 314.67,462.85 317.20,463.41 319.73,463.96 ' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
<rect x='241.17' y='81.22' width='75.06' height='57.60' style='stroke-width: 0.75; fill: #FFFFFF;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
<line x1='244.41' y1='95.62' x2='266.01' y2='95.62' style='stroke-width: 0.75;' clip-path='url(#cpNTkuMDR8MzI5Ljc2fDUwMi41Nnw1OS4wNA==)' />
@@ -150,7 +150,7 @@
<line x1='419.04' y1='357.34' x2='411.84' y2='357.34' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
<line x1='419.04' y1='292.94' x2='411.84' y2='292.94' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
<line x1='419.04' y1='228.54' x2='411.84' y2='228.54' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
-<line x1='419.04' y1='164.14' x2='411.84' y2='164.14' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
+<line x1='419.04' y1='164.15' x2='411.84' y2='164.15' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
<line x1='419.04' y1='99.75' x2='411.84' y2='99.75' style='stroke-width: 0.75;' clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)' />
<g clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)'><text transform='translate(401.76,489.47) rotate(-90)' style='font-size: 12.00px; font-family: Liberation Sans;' textLength='6.67px' lengthAdjust='spacingAndGlyphs'>0</text></g>
<g clip-path='url(#cpMC4wMHw3MjAuMDB8NTc2LjAwfDAuMDA=)'><text transform='translate(401.76,425.07) rotate(-90)' style='font-size: 12.00px; font-family: Liberation Sans;' textLength='6.67px' lengthAdjust='spacingAndGlyphs'>5</text></g>
@@ -180,13 +180,13 @@
<circle cx='453.90' cy='326.46' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<circle cx='430.36' cy='482.64' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<circle cx='430.36' cy='485.13' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
-<circle cx='429.07' cy='328.65' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
+<circle cx='429.07' cy='328.66' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<circle cx='429.07' cy='354.25' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<circle cx='429.07' cy='481.50' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<circle cx='429.07' cy='328.36' r='2.70pt' style='stroke-width: 0.75;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<polyline points='506.65,446.70 510.29,453.00 503.02,453.00 506.65,446.70 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
-<polyline points='506.65,391.22 510.29,397.52 503.02,397.52 506.65,391.22 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
-<polyline points='517.01,149.16 520.64,155.46 513.37,155.46 517.01,149.16 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
+<polyline points='506.65,391.22 510.29,397.51 503.02,397.51 506.65,391.22 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
+<polyline points='517.01,149.17 520.64,155.46 513.37,155.46 517.01,149.17 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<polyline points='517.01,479.70 520.64,485.99 513.37,485.99 517.01,479.70 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<polyline points='463.78,479.37 467.42,485.67 460.14,485.67 463.78,479.37 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<polyline points='463.78,416.53 467.42,422.82 460.14,422.82 463.78,416.53 ' style='stroke-width: 0.75; stroke: #FF0000;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
@@ -223,7 +223,7 @@
<line x1='447.84' y1='480.11' x2='455.47' y2='480.11' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<line x1='451.66' y1='483.93' x2='451.66' y2='476.29' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<line x1='437.55' y1='249.04' x2='445.19' y2='249.04' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
-<line x1='441.37' y1='252.86' x2='441.37' y2='245.22' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
+<line x1='441.37' y1='252.86' x2='441.37' y2='245.23' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<line x1='437.55' y1='423.26' x2='445.19' y2='423.26' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<line x1='441.37' y1='427.07' x2='441.37' y2='419.44' style='stroke-width: 0.75; stroke: #00CD00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
<line x1='419.04' y1='402.95' x2='689.76' y2='402.95' style='stroke-width: 0.75; stroke-dasharray: 4.00,4.00;' clip-path='url(#cpNDE5LjA0fDY4OS43Nnw1MDIuNTZ8NTkuMDQ=)' />
diff --git a/tests/testthat/FOCUS_2006_D.csf b/tests/testthat/FOCUS_2006_D.csf
index e3d2a303..22a4b125 100644
--- a/tests/testthat/FOCUS_2006_D.csf
+++ b/tests/testthat/FOCUS_2006_D.csf
@@ -5,7 +5,7 @@ Description:
MeasurementUnits: % AR
TimeUnits: days
Comments: Created using mkin::CAKE_export
-Date: 2019-10-25
+Date: 2019-10-26
Optimiser: IRLS
[Data]

Contact - Imprint