[SAS] 생존분석 (Survival Analysis) - Time-dependent Cox Proportional Hazard Model

생존분석을 진행하다 보면 때때로 time dependent (또는 time varying) cox proportional hazard model을 정의해야 할 때가 있다. ​Time dependent cox ph model이라고 하면 말 그대로, 시간에 따라 변화하는 공변량(ex. 이식 여부 등 0과 1로 나뉘는 이변수)이 존재하는 경우를 말할 수도 있지만 조금 더 복잡하게는 반복적으로 발생하는 event가 존재하는 경우를 말할 수도 있다.​​​1. 시간에 따라 변화하는 공변량이 존재하는 경우​간단한 코드로 분석할 수 있다.예를 들어 심장 이식 여부를 생각해 보면, 우리는 심장 이식이 생존 시간에 미치는 영향을 알고 싶을 것이다.하지만 심장 이식은 생존한 시간 이후에 이루어질 수 없다. ㅋㅋ그래서 심장 ..

  • format_list_bulleted SAS
  • · 2024. 11. 22.
  • textsms

[SAS] 생존분석 (Survival Analysis) - 다중공선성 및 비례위험가정, Kaplan-Meier Plot

1. 다중공선성다중공선성은 SAS에서 3가지 옵션을 통해 계산할 수 있다.주석에 적힌 다중공선성의 기준은 사람마다 다르게 받아들여지기 때문에 참고만 하자. PROC REG DATA=MAIN_TABLE;MODEL SURVIVED_TIME=AGE SEX DM / VIF; /*VIF, 기준 5 이상*/RUN;PROC REG DATA=MAIN_TABLE;MODEL SURVIVED_TIME=AGE SEX DM / TOL; /*TOL, 기준 0.2 이하*/RUN;PROC REG DATA=MAIN_TABLE;MODEL SURVIVED_TIME=AGE SEX DM / COLLIN; /*CI, 15 이상일 경우 의심*/RUN; ​2. 비례위험가정Cox Proportional Hazard Model을 생성하기 위해 선행되..

  • format_list_bulleted SAS
  • · 2024. 11. 22.
  • textsms

[SAS] 생존분석 (Survival Analysis) - 사망연계자료 JOIN 및 중도절단

아래 코드는 건강보험공단 DB에 사용하였던 코드를 표본코호트 DB에 맞게 변수명을 변형한 것으로, 실제로 적용하기에는 호환이 잘 안 될 수도 있으며 추후 (내가 다시) 참고하려고 작성하였으니 너무 믿지 마시길 바랍니다. (저는 제가 아직도 왜 전공자인지 의심이 됨ㅋㅋ)​코드에 나와 있는 추적관찰기간은 전부 제가 한 실제 분석과 관계 없는 가짜 데이터입니다. ​​1. 사망연계자료 JOIN하여 전처리 테이블 만들기PROC SQL;CREATE TABLE DEATH_TABLE ASSELECT A.*, B.DTH_YYYYMMFROM MAIN_DATA_TABLE AS AINNER JOIN NSC2_BND AS BON A.RN_INDI=B.RN_INDI;QUIT;PROC SQL;CREATE TABLE DEATH_TA..

  • format_list_bulleted SAS
  • · 2024. 11. 22.
  • textsms