'oddsrisk'에 해당되는 글 1건

  1. 2011.10.20 [Stata] oddsrisk 명령어
관련성을 평가할 때 많이 쓰는 교차비(odds ratio)는 엄밀하게는 위험비(risk ratio)와 다르다. 교차비와 위험비가 비슷하려면 (1) 관심 질환의 유병률이 낮거나(희귀 질환 가정), (2) 잘 설계된 환자-대조군 연구에서 산출됐을 경우이다. 실제 데이터에서는 이 두 가지 가정 어느 쪽도 만족하지 않는데도 교차비를 제시한 경우가 흔하다. 의학보건학 연구에서 교차비가 많이 쓰이는 이유는 범주형 결과 변수를 분석할 때 로지스틱 회귀분석을 주로 이용하고, 그 결과로 교차비를 쉽게 산출할 수 있기 때문이다. 위험비는 푸아송(Poisson)[각주:1] 회귀분석을 적용해서 계산할 수 있는데 많은 연구자들이 익숙치 않다. Joseph HilbeLogistic Regression Models 책 에서 제시한 oddsrisk 명령어로 교차비와 위험비를 둘다 산출해보고 비교해보도록 하자. (명령어는 findit oddsrisk로 찾아서 설치할 수 있고, 데이터셋은 링크한 CRC Press 책 소개 사이트에서 다운로드할 수 있다.)

. use anterior

. list

     +--------------------------+
     | count   death   anterior |
     |--------------------------|
  1. |   120       1          1 |
  2. |    67       1          0 |
  3. |  2005       0          1 |
  4. |  2504       0          0 |
     +--------------------------+

. oddsrisk death anterior [fw=count]

---------------------------------------------------------------------
Incidence for unexposed risk group =     0.0261
---------------------------------------------------------------------
Predictor    Odds Ratio   Risk Ratio     [95% Conf. Interval]
---------------------------------------------------------------------
anterior        2.2368       2.1670       1.6220       2.8807
---------------------------------------------------------------------

비노출군에서 발생 위험이 2.6%=67/(67+2504)로 낮으므로 교차비와 위험비가 비슷할 것이라는 예상을 해볼 수 있고, 실제 결과도 교차비가 2.24, 위험비가 2.17로 큰 차이가 없다.

poisson 명령어를 이용하여 푸아송 회귀분석을 수행하여 산출한 위험비와 비교해보자.

. poisson death anterior [fw=count], nolog irr robust

Poisson regression                                Number of obs   =       4696
                                                  Wald chi2(1)    =      26.69
                                                  Prob > chi2     =     0.0000
Log pseudolikelihood =  -776.2572                 Pseudo R2       =     0.0171

------------------------------------------------------------------------------
             |               Robust
       death |        IRR   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
    anterior |   2.166953   .3243484     5.17   0.000     1.616003    2.905741
       _cons |   .0260599   .0031423   -30.25   0.000     .0205748    .0330073
------------------------------------------------------------------------------

계산된 위험비가 2.17로 같음을 알 수 있다.

  1. 지금까지 통계학 관련 서적에서는 포아송으로 번역해 왔으나 국립국어원 외국 인명 표기에 따르면 '푸아송'이 맞다. [본문으로]
Posted by cyberdoc
: