교수님, 안녕하십니까? mytable과 multipleROC를 특히 유용하게 잘 쓰고 있는 연구자입니다.
다름이 아니라, multivariable logistic regression analysis를 시행하여 dichotomous한 outcome에 대해 6개의 factor가 유의미하다는 사실을 확인하였고, 이에 대한 ROC curve를 6개 factor를 모두 넣어 확인해 본 결과 AUC값은 0.981 로 확인되는 상태입니다. 하나하나 했을 때는 0.8~0.9 사이인 것에 비해 월등한 결과입니다.
여기서 요지는, 이렇게 나온 ROC curve에서 6개 factor 모두의 optimal cutoff를 찾고 싶은데, multipleROC package에서 automatic 하게 계산되어 나오는 값은 Youden's index (sensitivity+specificity가 최대값인 포인트) 결과대로 나오는 것 같은데 이 값은 제가 보여주고자 하는 값이 아닙니다. 그래서 sensitivity와 specificity를 곱한 값이 최대인 포인트를 확인하고 싶은데 multivariable 이 입력된 ROC curve에서는 이 값을 어떻게 찾을 수 있는지 잘 모르겠습니다. 구글링을 해보면
y= b0 + b1X1 + b2X2 + b3X3 ... 이렇게 수식이 만들어지기 때문에 optimal 한 X1, X2, X3 (multivariables) 의 값을 구할 수가 없다고 하는데 multipleROC에서는 그 값이 나오더라구요.
혹시 single variable에 대한 ROC curve에서 optimal cutoff를 찾을때처럼 lr.eta 값으로 raw variable의 값을 추정할 수 있는 수식이 있을까요? 아래와 같이 말입니다.
#######
x<- multipleROC(shock_sepsis3~PCT, data=train)
x$fit
write.csv(x$df,
"PCT.csv")
fit <- glm(shock_sepsis3~PCT, data=train, family= binomial)
(qlogis(x$cutpoint)-unname(coef(x$fit)[1]))/unname(coef(x$fit)[2])
(qlogis(0.607635)-unname(coef(x$fit)[1]))/unname(coef(x$fit)[2]) ## multiplication
########
github에 올려주신 코드를 봐도 어떻게 계산되어 나오는 것인지 이해가 되지 않아 여쭤봅니다.
고진선처 부탁드립니다. 감사합니다.
Comment 0
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 |