Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- set.seed(1)
- probLower = vector(length=1000)
- for(i in 1:1000){
- x = rnorm(1000)
- y = rbinom(n=1000, size=1, prob=pnorm(x))
- logitModel = glm(y~x, family=binomial(link="logit"))
- probitModel = glm(y~x, family=binomial(link="probit"))
- probLower[i] = deviance(probitModel)<deviance(logitModel)
- }
- sum(probLower)/1000
- [1] 0.695
- deviance(probitModel)
- [1] 1025.759
- deviance(logitModel)
- [1] 1026.366
- deviance(logitModel)-deviance(probitModel)
- [1] 0.6076806
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement