# NOT RUN {
library(dplyr)
qi.central.interval <- zelig(Petal.Width ~ Petal.Length + Species,
data = iris, model = "ls") %>%
setx(Petal.Length = 2:4, Species = "setosa") %>%
sim() %>%
zelig_qi_to_df() %>%
qi_slimmer()
# }
Run the code above in your browser using DataLab