안녕하세요 교수님
roc 커브에 대해 공부하고 있습니다.
교수님 자료로 아주 좋은 공부가 되었습니다. 언제나 감사드립니다.
다름이 아니라 변수 1개씩 넣어서 2가지 그래프를 비교할 때 DeLong's test for two correlated ROC curves를 이용하셨고
이 값이 plot_ROC (a1, a2) 에 고정되어 그림에 표기 되는 것 같습니다.
단, 3개 이상 비교를 할떄는 이 테스트가 그림에 표기 되지 않는것 같습니다. plot_ROC (a1, a2, a3)
2개를 비교한 그림에서도 DeLong's test for two correlated ROC curves 결과가 표기 되지 않도록 하는 방법은 없을까요?
교수님 자료에서 show.AUC = FALSE or show.AUC=TRUE 와 같은 조절 방법이 가능할까요?
감사합니다.
Comment 2
-
cardiomoon
2018.11.07 11:57
다음과 같이 DeLong's test를 출력하는 부분을 주석처리하면 될 것 같습니다. 감사합니다.plot_ROC=function(x,y=NULL,z=NULL,type=0,col=c("black","red","blue"),lwd=2,show.points=TRUE,show.line=FALSE,show.lr.eta=TRUE,show.sens=TRUE,show.AUC=TRUE){plotROCcount<<-1plot(0:1,0:1,xlim=0:1,ylim=0:1,type="n",ylab="Sensitivity",xlab="1-Specificity")if(type==1) rect(par("usr")[1],par("usr")[3],par("usr")[2],par("usr")[4],border=FALSE,col=gray(0.90))abline(h=seq(0,1,0.1),col=ifelse(type,"white",gray(0.9)))abline(v=seq(0,1,0.1),col=ifelse(type,"white",gray(0.9)))if(!is.null(x)) add_plot_ROC(x,lwd=lwd,col=col[1],show.points=show.points,show.line=show.line,show.lr.eta=show.lr.eta,show.sens=show.sens,show.AUC=show.AUC)segments(0,0,1,1,lty=2)if(!is.null(y)) add_plot_ROC(y,lwd=lwd,col=col[2],show.points=show.points,show.line=show.line,show.lr.eta=show.lr.eta,show.sens=show.sens,show.AUC=show.AUC)if(!is.null(z)) add_plot_ROC(z,lwd=lwd,col=col[3],show.points=show.points,show.line=show.line,show.lr.eta=show.lr.eta,show.sens=show.sens,show.AUC=show.AUC)# if((!is.null(x)) & (!is.null(y)) & is.null(z)) {# if((ncol(x$lr$model)==2)&(ncol(y$lr$model)==2)) {# result=roc.test(ROC2roc(x),ROC2roc(y),plot=T)# if(result$p.value <0.001) {# temp="p < 0.001"# } else temp=paste("p = ",round(result$p.value,3),sep="")## text(0.45,plotROCcount*0.11-0.05,cex=0.8,pos=4,col=col,# paste("DeLong's test for two correlated ROC curves\n",# "Z = ",round(result$statistic,3),", ",temp,sep=""))# }# }} -
PedEndo
2018.11.12 18:40
네 감사합니다.
교수님
덕분에 좋은 것을 배워갑니다!!!
No. | Subject | Author | Date |
---|---|---|---|
1253 | logistic regression 명령어 문의드립니다. [2] | 로렌쵸메디치 | 2023.11.22 |
1252 | spss 파일 불러오기. [1] | 통계린 | 2023.11.22 |
1251 | PSM 매칭변수 숫자 관련 문의입니다. [1] | 오하나 | 2023.11.22 |
1250 | Multiple imputation 관련 질문 및, age adjusted incidence rate 관련 질문을 드립니다. | YSKBH | 2023.11.15 |
1249 | semMediation 설치 재문의드립니다 [3] | 마술샘 | 2023.11.11 |
1248 | 안녕하세요. 정회원 문의 [1] | lionking | 2023.11.06 |
1247 | PSM 관련 질문드립니다. | dr. | 2023.11.05 |
1246 | ggplot [1] | 통계린 | 2023.10.30 |
1245 | 정회원 [1] | 통계린 | 2023.10.29 |
1244 | PSM 관련 질문 드립니다. [2] | 로렌쵸메디치 | 2023.10.25 |
1243 | forest plot 그림이 전체가 다 나오지 않습니다. [1] | joyguni | 2023.10.08 |
1242 | PSM caliper default 값 관련 문의드립니다. [1] | Ooppps | 2023.10.05 |
1241 | reliability test tool,,, [1] | FAdavid | 2023.10.05 |
1240 | 정회원 가입 확인 요청드립니다. [1] | osman | 2023.10.03 |
1239 | 생존분석, cox 회귀 관련 기초적인 질문 입니다. | 만니톨 | 2023.09.27 |
1238 | Twang-mnps에서 매칭된 데이터를 추출하는 방법 좀 알려주세요 | Dorocy | 2023.09.26 |
1237 | 영수증 발급을 부탁드립니다. [1] | 규빈어멈 | 2023.09.26 |
1236 | 정회원 확인 [2] | lionking | 2023.09.24 |
1235 | autoReg 에러 문의드립니다 [1] | 소소 | 2023.09.23 |
1234 | 안녕하세요 에러 관련 문의드립니다 [1] | 화이팅 | 2023.09.19 |