evidence (version 0.8.10)

looicplot: A dot plot is produced for several related models showing for each model its LOOIC-value with its credible interval.

Description

The LOOIC-value (like the non-Bayesian AIC-value) is a useful measure of model performance for model prediction.

Usage

looicplot(looiclist, modnames, perc = 90)

Arguments

looiclist

a list of character-valued names of rstanarm model objects

modnames

a character-valued vector of model names for each of the models

perc

the percentage credibility for the credible intervals (defaults to 90%)

Value

None provided, but a printed list of looic-values, their standard errors, and credible intervals, and a dot plot with the same information are produced.

References

van Hulst, R. 2018. Evaluating Scientific Evidence. ms.

See Also

dotchart

Examples

# NOT RUN {
data(budworm)
Mbudworm1 <- stan_glm(formula = cbind(dead, 20 - dead) ~ ldose,
                      family = binomial, data = budworm,
                      prior = student_t(df = 7),
                      prior_intercept = student_t(df = 7))
Mbudworm2 <- stan_glm(formula = cbind(dead, 20 - dead) ~ ldose * sex,
                      family = binomial, data = budworm,
                      prior = student_t(df = 7),
                      prior_intercept = student_t(df = 7))
Mbudworm3 <- stan_glm(formula = cbind(dead, 20 - dead) ~ ldose + sex,
                      family = binomial, data = budworm,
                      prior = student_t(df = 7),
                      prior_intercept = student_t(df = 7))
looicplot(looiclist = list("Mbudworm1", "Mbudworm2", "Mbudworm3"),
          modnames = c("~ ldose", "~ ldose + sex", "~ ldose * sex") )
# }