aboutsummaryrefslogtreecommitdiff
path: root/vignettes/web_only/dimethenamid_2018.rmd
diff options
context:
space:
mode:
authorJohannes Ranke <jranke@uni-bremen.de>2022-07-01 13:52:52 +0200
committerJohannes Ranke <jranke@uni-bremen.de>2022-07-01 13:52:52 +0200
commit16a7ed9548b37fe3c68c993651226fdc2dda6402 (patch)
treeefbc8d6c4b89737b337d91cf543b161af74202de /vignettes/web_only/dimethenamid_2018.rmd
parent7bd35821d04e104ad2fb2ed582c2c766964ffcf6 (diff)
Sensitivity of AIC values in DMTA vignette to BLAS version
Diffstat (limited to 'vignettes/web_only/dimethenamid_2018.rmd')
-rw-r--r--vignettes/web_only/dimethenamid_2018.rmd18
1 files changed, 16 insertions, 2 deletions
diff --git a/vignettes/web_only/dimethenamid_2018.rmd b/vignettes/web_only/dimethenamid_2018.rmd
index 41681d45..047cb6db 100644
--- a/vignettes/web_only/dimethenamid_2018.rmd
+++ b/vignettes/web_only/dimethenamid_2018.rmd
@@ -1,7 +1,7 @@
---
title: Example evaluations of the dimethenamid data from 2018
author: Johannes Ranke
-date: Last change 7 March 2022, built on `r Sys.setlocale("LC_TIME", "C"); format(Sys.Date(), format = "%d %b %Y")`
+date: Last change 1 July 2022, built on `r Sys.setlocale("LC_TIME", "C"); format(Sys.Date(), format = "%d %b %Y")`
output:
html_document:
toc: true
@@ -341,11 +341,19 @@ print(AIC_parent_saemix_methods)
```
The AIC values based on importance sampling and Gaussian quadrature are very
similar. Using linearisation is known to be less accurate, but still gives a
-similar value. In order to illustrate that the comparison of the three method
+similar value.
+
+In order to illustrate that the comparison of the three method
depends on the degree of convergence obtained in the fit, the same comparison
is shown below for the fit using the defaults for the number of iterations and
the number of MCMC chains.
+When using OpenBlas for linear algebra, there is a large difference in the
+values obtained with Gaussian quadrature, so the larger number of iterations
+makes a lot of difference. When using the LAPACK version coming with Debian
+Bullseye, the AIC based on Gaussian quadrature is almost the same as the one obtained
+with the other methods, also when using defaults for the fit.
+
```{r AIC_parent_saemix_methods_defaults, cache = FALSE}
f_parent_saemix_dfop_tc_defaults <- mkin::saem(f_parent_mkin_tc["DFOP", ])
f_parent_saemix_dfop_tc_defaults$so <-
@@ -394,6 +402,12 @@ attempts to quantify the variability of the slower rate constant of the
biphasic decline of dimethenamid indicate that the data are not sufficient
to characterise this variability to a satisfactory precision.
+# Session Info
+
+```{r sessionInfo, cache = FALSE}
+sessionInfo()
+```
+
# References
<!-- vim: set foldmethod=syntax: -->

Contact - Imprint