[Stata] oddsrisk 명령어
통계 연습 2011. 10. 20. 17:54 |
관련성을 평가할 때 많이 쓰는 교차비(odds ratio)는 엄밀하게는 위험비(risk ratio)와 다르다. 교차비와 위험비가 비슷하려면 (1) 관심 질환의 유병률이 낮거나(희귀 질환 가정), (2) 잘 설계된 환자-대조군 연구에서 산출됐을 경우이다. 실제 데이터에서는 이 두 가지 가정 어느 쪽도 만족하지 않는데도 교차비를 제시한 경우가 흔하다. 의학보건학 연구에서 교차비가 많이 쓰이는 이유는 범주형 결과 변수를 분석할 때 로지스틱 회귀분석을 주로 이용하고, 그 결과로 교차비를 쉽게 산출할 수 있기 때문이다. 위험비는 푸아송(Poisson) 회귀분석을 적용해서 계산할 수 있는데 많은 연구자들이 익숙치 않다. 1Joseph Hilbe가 Logistic 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 명령어를 이용하여 푸아송 회귀분석을 수행하여 산출한 위험비와 비교해보자.
. 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로 같음을 알 수 있다.
- 지금까지 통계학 관련 서적에서는 포아송으로 번역해 왔으나 국립국어원 외국 인명 표기에 따르면 '푸아송'이 맞다. [본문으로]
'통계 연습' 카테고리의 다른 글
[Stata] 다변수 분석에서 위험비(RR) 구하기 (0) | 2012.04.23 |
---|---|
[Stata] batplot 명령어 (0) | 2011.12.26 |
[Stata] 다변수 로지스틱 분석에서 범주형 설명 변수의 빈도표 (0) | 2011.10.18 |
[Stata] 표준화비 신뢰구간 계산 (0) | 2011.10.14 |
[Stata] 위험 차(risk difference) 계산 (2) | 2011.10.14 |