Presentation

1. # Tobit regression
2. library(survival)
3. dt<-brabant[samp,]
4. attach(dt)
5. rfm.tobit<-survreg(Surv(ca97, ca97>0, type="left") ~r6 + rr6 + f6m + m6m + mm6m + caa96 + cda96, dist='gaussian')
6. dt<-brabant[-samp,]
7. attach(dt)
8. predictv<-predict(rfm.tobit, type="response", newdata=dt)
9. #predictv<-predictv+rfm.tobit$scale # a verifier et justifier
10. predictv<-ifelse(predictv>0,predictv,0)
11. predictv<-predictv + 110 # a verifier et justifier
12. names(predictv)<-1:length(predictv)
13. gc<-gainchart_cont(predictv,ca97)
14. par(mfrow=c(2,2))
15. plot_gaincharts_cont(gc)