예전에 R에서 ordinal regression analysis 결과를 SPSS 처럼 만들어내지 못해서 SPSS이용하여 통계 분석을 하였다. 이번 리비젼 과정에서 SPSS를 사용할 수 없었기 때문에 R에서 ordinal regression analysis를 시도할 수 밖에 없었고, 인터넷을 통하여 거의 같은 결과를 만들어내는데 성공하였다.
우선 MASS패키지를 사용한다. 이 패키지에 포함되어 있는 polr() 함수가 기본적으로 ordinal regression analysis 에 사용하는 함수이다.
AER 패키지를 사용한다. coef()를 통하여 hazard ratio 를 구할 수 있다. coefci()를 통하여 95% 신뢰 구간을 구할 수 있다. coeftest()를 통하여 p.value를 구할 수 있다.