[R프로그래밍] rms 패키지, 로지스틱회귀 오즈비 (OR, Odds Ratio), 신뢰구간 (CI, confint) 구하기

R 에서 로지스틱 회귀분석을 위해서 주로 사용하는 패키지는 크게 2가지가 있습니다.

하나는 rms 패키지의 lrm 명령어, 

나머지는 glm 명령어 입니다.

보통은 glm 을 많이 쓰게 되는데...

RCS (restricted cubic spline curve) 를 사용할 때는 어쩔 수 없이 rms 패키지를 사용하게 됩니다. 

저는 glm 이 더 익숙해서... rms 패키지의 lrm 보다 더 편하더군요. 

어쨌든... 

lrm 명령어에서 OR, 오즈비 (Odds ratio) 및 95% 신뢰구간 (CI, confidence interval) 을 구하는 방법입니다.

코드는 아래와 같습니다.

lrmmodel <- lrm (y ~ x, data=df)
exp(coef(lrmmodel)
exp(confint.default(lrmmodel)) 

lrmmodel 이라는 함수에서, 

OR 을 구하기 위해서는 coefficient 값을 승산비 (exp) 해주면 됩니다.

그 코드가 2번째 줄입니다.

그리고 신뢰구간을 구하기 위해서는...

confint.default 명령어로 구한 값을 승산비 (exp) 해주면 됩니다.

glm 에서는 그냥 confint 만 입력하면 되는데...

lrm 에서는 confint.default 를 입력해야 되어 처음에 이걸 몰라서 헤맸던 기억이 있습니다.

혹시 glm 패키지를 이용해서 RCS 를 그릴 수 있는 방법을 아시는 분이 계시다면 알려주시면 감사하겠습니다. ^^

로지스틱회귀분석을 restricted cubic spline curve 로 표현하시고 싶으시다면 아래의 포스팅을 참고하여 주세요.

https://bpapa.tistory.com/26

도움이 되셨다면, '공감' 버튼 부탁드립니다^^