Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- set.seed(1234)
- x = runif(100, 0, 1)
- y = rbinom(100, size=1, prob = x)
- model = glm(y ~ x, family = binomial("logit"))
- newdata = data.frame(x = c(.25, .75))
- predicted.probs = predict(model, newdata, type="response", se.fit = T)
- upper.interval = predicted.probs$fit + 1.96 * predicted.probs$se.fit
- lower.interval = predicted.probs$fit - 1.96 * predicted.probs$se.fit
- errbar(newdata$x, predicted.probs$fit, upper.interval, lower.interval)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement