안녕하세요 교수님
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 |
---|---|---|---|
754 | 반복 측정한 자료에서 증가한 비율의 차이가 두 집단간에 유의하게 다른지 알아보려고 합니다. | 렁닥 | 2018.09.20 |
753 | data 전처리 관련 질문입니다.. [1] | 망고링고 | 2018.09.22 |
752 | 통계 초보가 효과 크기에 대해 질문드립니다. | 상냥한펭귄씨 | 2018.10.01 |
751 | R studio에서 data import 관련하여 문의드립니다. | Jun | 2018.10.04 |
750 | Univariate and Multiple Logistic Regression: predictor의 선택 | 신희순 | 2018.10.04 |
749 | Enter calculcation equation을 넣는 형식 | 신희순 | 2018.10.04 |
748 | 로지스틱 회귀분석 | 김밥친구 | 2018.10.04 |
747 | 교수님, 질문 드립니다. | kjg | 2018.10.06 |
746 | 메타분석 출간오류에 대하여 문의드립니다. | 노을 | 2018.10.14 |
745 | 교수님! cox 회귀분석 큐빅 스플라인 그래프 작성 질문드립니다! [1] | susan | 2018.10.15 |
744 | 출간올유 funnel plot 에 대해서 [2] | 노을 | 2018.10.15 |
743 | boxplot 질문드립니다. [2] | jjunglammy | 2018.10.17 |
742 | 결측치 메우는 방법 문의 [1] | 김태오 | 2018.10.19 |
741 | 생존분석 error | gracehn | 2018.10.25 |
740 | 다변량 분석 결과 작성시 기술방법에 대하여 질문 드립니다. [2] | drted | 2018.10.26 |
739 | ggplot에서 y axis limit정하는 방법 여쭙니다. [2] | fascinate23 | 2018.10.26 |
738 | 메타분석 문의 [1] | Maruya | 2018.10.29 |
737 | 교수님..ggBar를 시행하면 계속 에러메시지가 뜹니다..? [4] | 스튜던트 | 2018.10.30 |
736 | 변수 넣고 빼기 [1] | Yoonida | 2018.11.01 |
» | ROC 커브 관련 질문 [2] | PedEndo | 2018.11.02 |