Skip to menu

안녕하세요. 교수님과 통계 잘하시는 선생들께 여쭤보고 싶습니다.

ROC 비교하기 위해 교수님께서 만드신 library도 사용해보고, 다른 코드도 이용해보았는데,

다른 코드에서는 에러가 발생해서 이렇게 글을 올리게 되었습니다.

 

비교 모델은 데이터 셋 자체는 같으나 필요에 의해서 데이터셋이 변형이 되어 

데이터셋 길이가 다릅니다.

교수님 library로 이용했을 땐 특별히 오류는 나오지 않고, 그림에서 두 모델 간 차이가 없는지 p-value도 뜨지 않습니다.

그래서 다른 코드로 이용하여 비교했을 때는 아래와 같은 에러가 뜨나, 결과는 또 나옵니다.

Warning message:

In x1 + x2 :

  longer object length is not a multiple of shorter object length

 

에러를 무시할 수 있는 부분인지 궁금합니다.

감사합니다.

 

코드 참고하시라고 올려봅니다.

 

pred.ori<-predict(log_ori2, type="response")

pred.D1<-predict(log_D1, type="response")

 

library(Hmisc)

library(plyr)

improveProb(pred.ori, pred.D1, y=as.double(as.character(ori2$los_d)))

 

roc.ori<-roc(ori2$los_d~pred.ori)

roc.D1<-roc(D1$los_d~pred.D1)

 

plot(roc.ori)

plot(roc.D1, add=T, lty=2)

 

roc.test(roc.ori, roc.D1, method="delong")