From 062a1be18677a4e83a64d23622b5214258d9bc7d Mon Sep 17 00:00:00 2001 From: Johannes Ranke Date: Fri, 16 Sep 2022 21:31:32 +0200 Subject: Add llhist to the multistart vignette --- vignettes/web_only/multistart.html | 3 +++ vignettes/web_only/multistart.rmd | 11 +++++++++++ 2 files changed, 14 insertions(+) (limited to 'vignettes') diff --git a/vignettes/web_only/multistart.html b/vignettes/web_only/multistart.html index ab1b768d..004b62be 100644 --- a/vignettes/web_only/multistart.html +++ b/vignettes/web_only/multistart.html @@ -391,6 +391,9 @@ parhist(f_saem_full_multi, lpos = "bottomright") 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

+
llhist(f_saem_reduced_multi)
+

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