DiagnosisMed

DiagnosisMedパッケージを使ってROCカーブを描いてみる

こちらを参考 library(DiagnosisMed) Ncase<-100 Ncont<-100 mean.case<-40 sd.case<-3 mean.cont<-35 sd.cont<-2 Case<-rnorm(Ncase,mean.case,sd.case) Cont<-rnorm(Ncont,mean.cont,sd.cont) boxplot(cbind(Case,Cont)) ra<-range(c(Case,Cont)) breaks<-…