looicplot function

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

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

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%)

Returns

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.

Author(s)

Robert van Hulst

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") )
## End(Not run)