안녕하세요? 저는 현재 박사과정에 재학중인 학생입니다.
다름이 아니라 현재 건강보험공단 데이터로 분석을 하고 있고,
ARX tool을 이용하여 데이터를 비식별화하여 original 데이터와 비식별화된 데이터 간의 비교를 하고 있습니다.
regression 결과를 비교를 하려고 하는데,
ROC curve는 선생님께서 만들어주신 패키지를 이용하였고, 결과가 아주 잘 나왔습니다. 감사합니다.
여기서 질문드리고자 하는 것은,
여러 모델에 대한 auc를 delong's test로 test하고싶은데,
이전에 답해주신 내용을 보면
multipleROC를 이용하여 답변주신 내용이 있습니다.
a=multipleROC(male~height,data=radial)
plot(a)
b=multipleROC(male~age,data=radial)
plot_ROC(list(a,b))
따라해보니 위에 내용 자체도 실행이 잘 안되고, multipleROC 함수 자체가 로딩이 잘 안되서 설치가 잘못되었는지,
또 roc.test 예시도 그렇고 보통 predictor 1~2개에, outcome 1개 이렇게 놓고 확인하는데,
제가 test하고자하는 모델은 dental_caries~age+gender+bmi+cholesterol+cr+gtp +smoking_status+alcohol+plaque
이렇게 길어서 auc 비교검정하기가 어려운 부분일지 질문드립니다.
감사합니다.
아래처럼 해보시겠어요?
roc1 <- pRoc::roc(...)
roc2 <- pRoc::roc(...)
pROC::roc.test(roc1, roc2)$p.value