# NOT RUN {
library(lme4)
lmer1 <- lmer(Sepal.Length ~ Petal.Length + (1 | Species), data = iris)
lmer2 <- lmer(Sepal.Length ~ Petal.Length + (Petal.Length | Species), data = iris)
lmer3 <- lmer(
Sepal.Length ~ Petal.Length + (Petal.Length | Species) + (1 | Petal.Width),
data = iris
)
m <- bayesfactor_models(lmer1, lmer2, lmer3, denominator = 1)
m
update(m, reference = "bottom")
as.matrix(m)
# }
Run the code above in your browser using DataLab