Présentation

Choix	Montant	Montant moyen
0	0	0
1	1-249	168,10
2	250-499	335,13
3	500 et plus	759,43

1. # ordinal regression
2. library(MASS)
3. dt<-brabant[samp,]
4. attach(dt)
5. ord<-ifelse(ca97<1,0,ifelse(ca97<250,1,ifelse(ca97<450,2,3)))
6. ordformula<-as.factor(ord) ~ r6 + rr6 + f6m + m6m + mm6m + caa96 + cda96
7. rfm.ordreg<-polr(ordformula)
8. dt<-brabant[-samp,]
9. attach(dt)
10. ord<-ifelse(ca97<1,0,ifelse(ca97<250,1,ifelse(ca97<450,2,3)))
11. predictv<-predict(rfm.ordreg, type="prob", newdata=dt)
12. rownames(predictv)<-1:nrow(predictv)
13. gc<-gainchart_ord(predictv,ord)
14. par(mfrow=c(2,2))
15. plot_gaincharts_ord(gc)