diff options
| author | Johannes Ranke <jranke@uni-bremen.de> | 2022-10-28 13:39:15 +0200 |
|---|---|---|
| committer | Johannes Ranke <jranke@uni-bremen.de> | 2022-10-28 13:39:15 +0200 |
| commit | f820bf5b91be0f589de16c3e3250f5f79672df75 (patch) | |
| tree | 2b1406e1c9286634ca017db586e09e2299dec048 /vignettes | |
| parent | b1740ade9a1746ccdb325b95915ef88872489f03 (diff) | |
Rename parhist to parplot and make it generic
That parhist name was not the brightest idea, as it does
not show histograms.
Diffstat (limited to 'vignettes')
| -rw-r--r-- | vignettes/web_only/multistart.rmd | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/vignettes/web_only/multistart.rmd b/vignettes/web_only/multistart.rmd index 94f4ef98..27a8a96a 100644 --- a/vignettes/web_only/multistart.rmd +++ b/vignettes/web_only/multistart.rmd @@ -40,7 +40,7 @@ with different starting values. ```{r} f_saem_full_multi <- multistart(f_saem_full, n = 16, cores = 16) -parhist(f_saem_full_multi) +parplot(f_saem_full_multi) ``` This confirms that the variance of k2 is the most problematic parameter, so we @@ -51,7 +51,7 @@ for k2. f_saem_reduced <- update(f_saem_full, no_random_effect = "log_k2") illparms(f_saem_reduced) f_saem_reduced_multi <- multistart(f_saem_reduced, n = 16, cores = 16) -parhist(f_saem_reduced_multi, lpos = "topright") +parplot(f_saem_reduced_multi, lpos = "topright") ``` The results confirm that all remaining parameters can be determined with sufficient @@ -67,7 +67,7 @@ The parameter histograms can be further improved by excluding the result with the low likelihood. ```{r} -parhist(f_saem_reduced_multi, lpos = "topright", llmin = -326, ylim = c(0.5, 2)) +parplot(f_saem_reduced_multi, lpos = "topright", llmin = -326, ylim = c(0.5, 2)) ``` We can use the `anova` method to compare the models, including a likelihood ratio |