CoxPH 와 model selection

없을리가 없다고 생각했지만 그 동안 R에서 못썼던 기능 중 하나는 Cox PH 분석에서 forward 혹은 backward 방식의 model selection 이다. 너무 쉽게 찾아낸 방법은 step() 함수이다.

step() 괄호 안에 coxph 명령어를 넣고 실행하면 SPSS 랑 똑같이 나온다. 다만, forward 순서는 조금 다르게 나오는 듯함. 내가 잘못 사용하고 있는지도 모르겠다.

SPSS 에서는 결측값이 있는 항목이 들어갈 경우 자동으로 그 부분을 인식하지만, R에서는 에러가 나기 때문에 결측값이 없는 자료만을 입력해야 한다.

trace=0 을 사용하면 최종 결과만을 확인할 수 있다.