Skip to menu

Risk scoring???

홍성엽 2019.01.28 16:04 Views : 181

전에 외부에서 하는 문교수님 워크샵을 듣고, 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이라는 패키지가 있던데 이걸 사용하면 가능할까요?



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

 

감사합니다.