Présentation

1. #Probit Regression
2. dt<-brabant[samp,]
3. attach(dt)
4. formula<-if_cde97 ~ r6 + rr6 + f6m + m6m + mm6m + caa96 + cda96
5. rfm.probit<-glm(formula,family=binomial(link=probit))
6. par(mfrow=c(2,2))
7. plot(rfm.probit)
8. dt<-brabant[-samp,]
9. attach(dt)
10. predictv<-predict(rfm.probit, type="response", newdata=dt)
11. names(predictv)<-1:length(predictv)
12. gc<-gainchart(predictv,dt$if_cde97)
13. par(mfrow=c(2,2))
14. plot_gaincharts(gc)