#ROC曲线
TPR=rep(0,1000)
FPR=rep(0,1000)
for(i in 1:1000){
p0=i/1000;
ypred<-1*(fitt.pi>p0)
TPR[i]=sum(ypred*audit2$is_run)/sum(audit2$is_run)
FPR[i]=sum(ypred*(1-audit2$is_run))/sum(1-audit2$is_run)}
plot(FPR,TPR,type="l",col=2,main = "ROC Curve")
points(c(0,1),c(0,1),type="l",lty=2)
auc(glm.safe1)
TPR=rep(0,1000)
FPR=rep(0,1000)
for(i in 1:1000){
p0=i/1000;
ypred<-1*(fitt.pi>p0)
TPR[i]=sum(ypred*audit2$is_run)/sum(audit2$is_run)
FPR[i]=sum(ypred*(1-audit2$is_run))/sum(1-audit2$is_run)}
plot(FPR,TPR,type="l",col=2,main = "ROC Curve")
points(c(0,1),c(0,1),type="l",lty=2)
auc(glm.safe1)