메뉴 건너뛰기

웹에서 하는 R 통계

Risk scoring???

2019.01.28 16:04

홍성엽 조회 수:16

전에 외부에서 하는 문교수님 워크샵을 듣고, web-r을 사용하면서 항상 감사하게 생각하면서 도움을 받고 있습니다. 감사합니다.

최근에 논문을 하나 쓰고 있는데 힘든 문제가 있어 도움을 받고자 글을 남깁니다.


제가 하고자 하는 분석은 약물중독환자에서 사망 및 중환자실 입원 등에 대한 예후인자를 분석하고 합니다.

이미 많은 연구가 진행된 분야이지만 logistic regression Neural network를 이용한 분석을 비교하여 논문을 쓰고자 합니다.

 

  천건 정도의 환자자료를 모았습니다. 첨부된 엑셀파일은 약물의 종류를 ICD-10 중독분류에 따라 분류한 표와 사망여부, ICU 재원일수를 표시한 자료입니다.

제가 분석을 시작하고자 하는데 가장 난관이 되는 문제가 있어 질문드리고자 합니다.

 

약물중독의 예후인자는 나이, 성별, 질환 등 여러가지가 있겠지만 약물의 종류와 용량도 크게 좌우합니다.

용량은 정확히 계량이 불가능한 경우가 많아 bias가 너무 심할 것 같습니다

한 환자가 한가지 성분의 약을 먹는 것도 아니고 여러가지 섞어서 먹는 경우도 있어서 분석이 더 힘듭니다.

 

그래서 여러 방법 중 약물종류를 ICD-10 질병분류에 기준해서 중독약물별로 분류하기로 하였습니다. 

첨부하여 드린 엑셀표의 샘플자료에서 toxcatg1 toxcatg2 toxcatg3 toxcatg4 toxcatg5가 먹은 약물의 ICD 카테고리 분류입니다.

 

문제는 여기서 약물카테고리 별로 분류는 했으나 이걸 독립변수로 넣을 방법이 막막하여 

약물 종별로 사망  또는 중환자실 입원일수에 대한 risk score를 환자별로 계산하고 이를 독립변수의 하나로 사용해서 분석할 예정입니다.


Propensity score matching처럼 될 수 있을 거라 생각했지만 propensity는 결과를 보여주지 않고 매칭만 하기 때문에 저런 risk score를 계산할 수 있을지 해결책을 알아보려고 열심히 찾아 보았으나 찾을 수가 없어서 염치불구하고 질문드리게 되었습니다. 1중독약물이 각각 다르고 먹은 사람마다 주요약물이 달라서 여러 종의 약물을 동시에 먹은 경우 복합적인 hazards를 모두 고려하여야 하는데 해결할 방법이 있을지 질문드립니다.


CRAN에 riskRegression이라는 패키지가 있던데 이걸 사용하면 가능할까요?



어떠한 의견이든 저에겐 큰 도움이 될 것 같습니다.

 

감사합니다.

번호 제목 글쓴이 날짜 조회 수
공지 묻고답하기 게시판을 운영합니다 [1] cardiomoon 2015.04.20 622
420 Survdiff로 log-rank test 시 P-value 소수점 begabung 2019.02.15 3
419 mytable 관련 질문드립니다 [3] file 낡은찻잔 2019.02.12 18
418 web-R 기반 분석 후 다운로드 관련 문의점 kato7 2019.02.11 5
417 stepwise logistic 회귀에서 결측치 문제 질문드립니다. 매직하마 2019.02.10 12
416 문의 드립니다. [1] 김날 2019.01.31 20
» Risk scoring??? file 홍성엽 2019.01.28 16
414 깃허브 적용 bnmy6581 2019.01.24 14
413 짝지은 자료 (PSM 등 매칭을 통한 자료)에서의 COX(생존분석)방법 의학통계 2019.01.23 76
412 통계 초보 질문있습니다(저에게는 매우 어려운 질문입니다.) file 이스라 2019.01.22 94
411 변수를 변경하는데 질문이 있습니다. [2] 최00 2019.01.21 267
410 NRI, IDI... 쑤니 2019.01.14 29
409 안녕하세요. K-M 곡선에서 censored mark에 대해서 문의드립니다. 간간해요 2019.01.14 40
408 로지스틱 회귀분석의 표본 추출에 관하여 문의드립니다 TAOKAKA 2019.01.12 24
407 안녕하세요~ 통계 질문드립니다. [3] 김두영 2019.01.04 59
406 ggplot2에 관하여 쑤니 2018.12.30 110
405 ios bnmy6581 2018.12.28 94
404 Incidence rate 그리는 방법 ? [1] 심심심심 2018.12.25 49
403 Box plot 그래프를 그리고 싶은데요. [2] file 알은어려워 2018.12.25 54
402 knit PDF 및 ztable에서 표 출력에 대하여 질문드립니다. [5] file zepiros15 2018.12.22 35
401 생존분석에서 연속형변수의 cutpoint에 대해 질문드립니다. [2] zepiros15 2018.12.17 48