Presentation

1. # Linear discriminant analysis
2. library(MASS)
3. dt<-brabant[samp,]
4. attach(dt)
5. formula<-if_cde97 ~ r6 + rr6 + f6m + m6m + mm6m + caa96 + cda96
6. rfm.lda<-lda(formula)
7. dt<-brabant[-samp,]
8. attach(dt)
9. predictv<-pnorm(predict(rfm.lda, newdata=dt)$x)[,1]
10. names(predictv)<-1:length(predictv)
11. gc<-gainchart(predictv,if_cde97)
12. par(mfrow=c(2,2))
13. plot_gaincharts(gc)