# predict
rf.pred = predict(rf, newdata = testGroup, type="response")
bg.pred = predict(bagg, newdata = testGroup, type="prob")
ev.pred = predict(ev, newdata = testGroup, type="response")
ct.pred = predict(tree, newdata = testGroup, type="response")
rp.pred = predict(rp, newdata = testGroup, type="response")

plot(rf.pred)
plot(bg.pred)
plot(ev.pred)
plot(ct.pred)
plot(rp.pred)

#library(ROCR)
#library(pROC)
library(sensitivity)
library(AUC)

# predictions
bagging.pred = prediction(bg.pred, testGroup$`whether he/she donated blood in March 2007`)
randomForest.pred = prediction(rf.pred, testGroup$`Numeric whether he/she donated blood in March 2007`)
evtree.pred = prediction(ev.pred, testGroup$`Numeric whether he/she donated blood in March 2007`)
ctree.pred = prediction(ct.pred, testGroup$`whether he/she donated blood in March 2007`)
rpart.pred = prediction(rp.pred, testGroup$`whether he/she donated blood in March 2007`)

bagging.specificity = auc(specificity(factor(bg.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
bagging.sensitivity = auc(sensitivity(factor(bg.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
bagging.roc = roc(factor(bg.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`)))
bagging.auc = auc(bagging.roc)
plot(performance(bagging.pred, "sens", "spec"))

randomForest.specificity = auc(specificity(factor(rf.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
randomForest.sensitivity = auc(sensitivity(factor(rf.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
randomForest.roc = roc(factor(rf.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`)))
randomForest.auc = auc(randomForest.roc)
plot(performance(randomForest.pred, "sens", "spec"))

evtree.specificity = auc(specificity(factor(ev.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
evtree.sensitivity = auc(sensitivity(factor(ev.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
evtree.roc = roc(factor(ev.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`)))
evtree.auc = auc(evtree.roc)
plot(performance(evtree.pred, "sens", "spec"))

ctree.specificity = auc(specificity(factor(ct.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
ctree.sensitivity = auc(sensitivity(factor(ct.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
ctree.roc = roc(factor(ct.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`)))
ctree.auc = auc(ctree.roc)
plot(performance(ctree.pred, "sens", "spec"))

rpart.specificity = auc(specificity(factor(rp.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
rpart.sensitivity = auc(sensitivity(factor(rp.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`))))
rpart.roc = roc(factor(rp.pred), factor(as.numeric(testGroup$`whether he/she donated blood in March 2007`)))
rpart.auc = auc(rpart.roc)
plot(performance(rpart.pred, "sens", "spec"))