Advertisement
Guest User

Untitled

a guest
Jun 15th, 2019
78
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.44 KB | None | 0 0
  1. set.seed(1234)
  2.  
  3. x = runif(100, 0, 1)
  4. y = rbinom(100, size=1, prob = x)
  5.  
  6. model = glm(y ~ x, family = binomial("logit"))
  7.  
  8. newdata = data.frame(x = c(.25, .75))
  9.  
  10. predicted.probs = predict(model, newdata, type="response", se.fit = T)
  11.  
  12. upper.interval = predicted.probs$fit + 1.96 * predicted.probs$se.fit
  13. lower.interval = predicted.probs$fit - 1.96 * predicted.probs$se.fit
  14.  
  15. errbar(newdata$x, predicted.probs$fit, upper.interval, lower.interval)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement