From 062a1be18677a4e83a64d23622b5214258d9bc7d Mon Sep 17 00:00:00 2001
From: Johannes Ranke <jranke@uni-bremen.de>
Date: Fri, 16 Sep 2022 21:31:32 +0200
Subject: Add llhist to the multistart vignette

---
 vignettes/web_only/multistart.rmd | 11 +++++++++++
 1 file changed, 11 insertions(+)

(limited to 'vignettes/web_only/multistart.rmd')

diff --git a/vignettes/web_only/multistart.rmd b/vignettes/web_only/multistart.rmd
index 5f109f3d..a25b73ce 100644
--- a/vignettes/web_only/multistart.rmd
+++ b/vignettes/web_only/multistart.rmd
@@ -46,3 +46,14 @@ f_saem_reduced <- update(f_saem_full, covariance.model = diag(c(1, 1, 0, 1)))
 f_saem_reduced_multi <- multistart(f_saem_reduced, n = 16, cores = 16)
 parhist(f_saem_reduced_multi, lpos = "bottomright")
 ```
+
+We can also analyse the log-likelihoods obtained in the multiple runs:
+
+```{r}
+llhist(f_saem_reduced_multi)
+```
+
+The one run with the lower likelihood is probably responsible for the outliers
+in the boxplots above, and caused by some weird starting value combination. In
+any case, the converged values from the original fit (red circles) appear
+perfectly acceptable, supporting the choice of starting values made by mkin.
-- 
cgit v1.2.1

