# NOT RUN {
library(estimate)
library(report)
library(lme4)
data <- iris
data$Group <- ifelse(data$Sepal.Width > 3, "A","B")
model <- lm(Petal.Width ~ Species * Group, data = data)
report(estimate_contrasts(model))
# }
# NOT RUN {
library(rstanarm)
model <- stan_glm(Sepal.Width ~ Species * Petal.Width, data = iris)
report(estimate_contrasts(model))
report(estimate_contrasts(model, fixed = "Petal.Width"))
report(estimate_contrasts(model, modulate = "Petal.Width", length = 2))
# }
# NOT RUN {
# }
# NOT RUN {
library(estimate)
library(rstanarm)
model <- stan_glm(Sepal.Width ~ Species * Petal.Width, data = iris)
report(estimate_means(model))
report(estimate_means(model, modulate = "Petal.Width"))
# }
# NOT RUN {
library(rstanarm)
model <- stan_glm(Sepal.Width ~ Species * Petal.Width, data = iris)
report(estimate_slopes(model))
# }
# NOT RUN {
library(rstanarm)
model <- stan_glm(Sepal.Width ~ poly(Petal.Length, 2), data = iris)
report(estimate_smooth(model))
# }
Run the code above in your browser using DataLab