# NOT RUN {
library(psycho)
library(ggplot2)
fit <- glm(Sex ~ Adjusting, data = affective, family = "binomial")
refgrid <- psycho::refdata(affective, "Adjusting")
predicted <- get_predicted(fit, newdata = refgrid)
ggplot(predicted, aes(x = Adjusting, y = Sex_Predicted)) +
geom_line() +
geom_ribbon(aes(
ymin = Sex_CI_2.5,
ymax = Sex_CI_97.5
),
alpha = 0.1
)
# }
Run the code above in your browser using DataLab