안녕하세요 교수님
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 |
---|---|---|---|
1193 | 첫번째 행을 라벨로 사용 버튼 클릭시 에러 문의 [3] | 수담이 | 2023.05.04 |
1192 | 안녕하세요? web-R table 그리고 관련하여 문의드립니다. [4] | mkw7654 | 2023.05.02 |
1191 | 메타분석 U3 그래프 관련 문의드립니다 | 화이팅 | 2023.05.02 |
1190 | 교수님 안녕하세요~~ | 메시12 | 2023.04.29 |
1189 | 연속형 변수인데 범주형으로 인지됩니다. [2] | shararaku | 2023.04.28 |
1188 | 통계마당 [2] | ejo | 2023.04.27 |
1187 | moonBook 오류 관련 재문의드립니다. [1] | holyhyun | 2023.04.24 |
1186 | moonBook 사용시 오류메시지로 문의드립니다. [1] | holyhyun | 2023.04.24 |
1185 | 에러 문의 [2] | 화이팅 | 2023.04.19 |
1184 | Cox regression 에 대한 문의입니다. [2] | saga | 2023.04.18 |
1183 | PSM 문의 [3] | starwars | 2023.04.13 |
1182 | ggPS, ggPSMSummary function [6] | 로렌쵸메디치 | 2023.04.12 |
1181 | Survival 분석에 대해서 문의 드립니다. | 정훈신 | 2023.04.05 |
1180 | 안녕하세요. 분석 툴에 대해 문의드립니다. [1] | 스매직맨 | 2023.04.04 |
1179 | 드디어 publish 되었네요. [2] | 정지윤 | 2023.04.04 |
1178 | 연구동향분석 시 논문 데이터 정제 관련 [5] | vtlbtl | 2023.03.31 |
1177 | 안녕하세요 교수님 에러 문의드립니다 [2] | 소소 | 2023.03.26 |
1176 | 코딩값이 여러개인 변수에 대한 전체 p값을 알 수 있는 방법이 궁금합니다 | 쥬쥬즈주 | 2023.03.22 |
1175 | PSM 매칭 버튼 보이지 않음. | swpapa | 2023.03.19 |
1174 | group-based trajectory modeling | swpapa | 2023.03.19 |