최근 수정 시각 : 2022-11-27 12:43:03

역학(의학)/관련도와 영향도의 측정

상대위험도에서 넘어옴

파일:상위 문서 아이콘.svg   상위 문서: 역학(의학)
예방의학
Preventive Medicine
{{{#!wiki style="margin:0 -10px -5px; min-height:calc(1.5em + 5px)"
{{{#!folding [ 펼치기 · 접기 ]
{{{#!wiki style="margin:-5px -1px -11px; word-break: keep-all"
<colbgcolor=#3c6><colcolor=#fff> 역학
(疫學)
감염 ( 돌파감염 · 수직감염) · 검역 · 면역 · 방역 · 병원체 · 봉쇄 · 사회적 거리두기( 사회적 거품) · 생물 안전도 · 생물재해 · 손 씻기 · 슈퍼전파자 · 역학조사 · 예방접종( 백신) · 음압병실 · 자가격리 · 전염병( 국제적 공중보건 비상사태 · 범유행전염병 · 법정 감염병) · 온실 면역(방역의 역설) / 혼합 면역(하이브리드 면역) · 집단 면역 · 코호트 격리 · 항원 결정기 / 항원 결합부 / 항원의 원죄
관련도와 영향도의 측정 · 질병과 사망의 측정 · 힐의 기준 · 연구에서 발생할 수 있는 편향
역사적 사건
14세기 흑사병 · 1918-19년 인플루엔자 범유행 · 2003년 SARS 유행 · 2009년 인플루엔자 범유행 · 2014년 에볼라 유행 · 2015년 MERS 유행 · 2019-23년 코로나-19 범유행
환경보건학
공해병 ( 미나마타병 · 온산병 · 이타이이타이병 · 카네미 유증사건)
의료관리학
간호관리학 }}}}}}}}}


1. 오즈비(OR)
1.1. 예시1.2. 로지스틱 회귀분석과 오즈비1.3. 짝지어진 자료에서 오즈비
2. 상대위험도(RR)
2.1. 예시2.2. 오즈비와 상대위험도의 관계
3. 코크란-멘텔-헨젤 검정
3.1. 방법3.2. 예시
3.2.1. 원 자료3.2.2. 층화 분석된 자료3.2.3. 코크란-멘텔-헨젤 검정 적용
4. 기여위험도(AR)
4.1. 예시4.2. 예방가능분율
5. 위험비(HR)

1. 오즈비(OR)

통계학에서 오즈(odds)는 사건이 일어날 확률을 사건이 일어나지 않을 확률로 나눈 값이다. 따라서, 오즈비(비차비, odds ratio) 또는 승산비는 두 집단간의 오즈값의 비이다. 일반적으로 역학적 연구에서 오즈비는 위험 요인에 노출된 집단과 그렇지 않은 집단의 질병 발생에 대한 오즈값의 비 또는 질병이 발생한 집단과 그렇지 않은 집단의 위험 요인 노출에 대한 오즈값의 비를 뜻하며, 두 방식으로 구한 오즈비값은 서로 같다. 약자로는 OR로 나타낸다.
질병 발생 질병 발생 안함
위험 요인 노출 a b
위험 요인 노출 안됨 c d
위 표에서 오즈비는 [math((a/c)/(b/d) = (a/b)/(c/d) = \dfrac{ad}{bc})]로 구해진다.
log(오즈비)의 표준오차(standard error)는 [math(\sqrt{\dfrac1a+\dfrac1b+\dfrac1c+\dfrac1d})]로 구할 수 있으며, 이를 이용하여 오즈비의 신뢰구간을 구할 수 있다.

1.1. 예시

폐암 발생 폐암 발생 안함
흡연 250 350
흡연 안함 100 300
위 예시에서 오즈비는 [math(\dfrac{ad}{bc}=\dfrac{250×300}{350×100})]으로, [math(\dfrac{15}{7})], 약 2.14이다.
log(오즈비)의 표준오차는 [math(\sqrt{\dfrac{1}{250}+\dfrac{1}{350}+\dfrac{1}{100}+\dfrac{1}{300}})]로, 약 0.14이다.
95% 신뢰수준에서 오즈비의 신뢰구간은 [math([exp(log2.14-1.96×0.14)), exp(log2.14+1.96×0.14))]]로, [1.63, 2.82]이다.

1.2. 로지스틱 회귀분석과 오즈비

로지스틱 회귀에서 모형은 [math(log(\dfrac{P}{1-P})=log(odds)=b_0+b_1x)]와 같이 나타낼 수 있다. 여기서 독립변수가 요인에 대한 노출 여부라면, [math(log(odds_{exposed})=b_0+b_1)], [math(log(odds_{unexposed})=b_0)]이다. 이때, [math(log(odds_{exposed})-log(odds_{unexposed})=log(\dfrac{odds_{exposed}}{odds_{unexposed}})=log(OR))]이고, 이는 [math((b_0+b_1)-b_0=b_1)]과 같다. 따라서, [math(log(OR)=b_1)]이므로, OR은 로지스틱 회귀분석에서 [math(e^{계수})]와 같다.

1.3. 짝지어진 자료에서 오즈비

Control
노출 노출 안됨
Case 노출 a b
노출 안됨 c d
위 표와 같이 짝지어진 자료에서 오즈비는 [math(\dfrac{b}{c})]로 구한다. 짝지어진 환자-대조군 연구와 같은 상황에서는 이 방법으로 오즈비를 구한다.

2. 상대위험도(RR)

상대위험도(relative risk)는 특정 요인에 노출된 사람들과 그렇지 않은 사람들의 질병 발생률의 비이다. 약자로는 RR로 나타낸다. 환자-대조군 연구의 경우, 발생률을 구할 수 없으므로, 상대위험도를 구할 수 없다.
질병 발생 질병 발생 안함 발생률
위험 요인 노출 a b [math(\dfrac{a}{a+b}=I_e)]
위험 요인 노출 안됨 c d [math(\dfrac{c}{c+d}=I_u)]
위 표에서 상대위험도는 [math(\dfrac{a}{a+b}/\dfrac{c}{c+d})] 또는 [math(\dfrac{I_e}{I_u})]로 구할 수 있다.
log(RR)의 표준오차는 [math(\sqrt{\dfrac{b}{a(a+b)}+\dfrac{d}{c(c+d)}})]로 구할 수 있으며, 이를 이용하여 상대위험도의 신뢰구간을 구할 수 있다.

2.1. 예시

폐암 발생 폐암 발생 안함 발생률
흡연 250 350 [math(\dfrac{250}{600}≈0.417)]
흡연 안함 100 300 [math(\dfrac{100}{400}=0.25)]
위 예시에서 상대위험도는 [math(\dfrac{250}{600}/\dfrac{100}{400}=\dfrac53)]으로, 약 1.67이다.

log(RR)의 표준오차는 [math(\sqrt{\dfrac{350}{250×600}+\dfrac{300}{100×400}})]으로, 약 0.099이다.

이를 이용해 구한 95% 신뢰수준에서 상대위험도의 신뢰구간은 [1.24,1.82]이다.

2.2. 오즈비와 상대위험도의 관계

오즈비와 상대위험도가 1보다 큰 경우에는 오즈비>상대위험도이고, 오즈비와 상대위험도가 1보다 작은 경우에는 오즈비<상대위험도이다. 즉, OR>RR>1 또는 OR<RR<1이다. 따라서, 오즈비는 상대위험도에 비하여 인과성을 더 크게 추정한다. 특히, 발생률이 커질수록 오즈비와 상대위험도의 차이가 커진다.
질병 발생 질병 발생 안함 발생률
위험 요인 노출 99 1 99%
위험 요인 노출 안됨 80 20 80%
위 표와 같이 흔한 질병의 경우 OR=24.75, RR≈1.24로, OR과 RR의 차이가 매우 크다.
질병 발생 질병 발생 안함 발생률
위험 요인 노출 20 80 20%
위험 요인 노출 안됨 1 99 1%
위 표와 같이 드문 질병의 경우 OR=24.75, RR=20로, OR과 RR의 차이가 크지 않다.

파일:OR과 RR의 관계.jpg

위 그래프는 발생률(I0)에 따른 오즈비(OR)과 상대위험도(RR)의 관계를 나타내고 있다. 발생률이 낮은 경우(I0=0.01) OR과 RR이 거의 1:1 대응을 보이고 있지만, 발생률이 높은 경우(I0=0.5)에는 OR이 RR을 과대추정함을 알 수 있다.

3. 코크란-멘텔-헨젤 검정

혼란변수가 있는 경우 연구하고자 하는 요인과 질병과의 관계가 과대 또는 과소추정될 수 있으므로, 이를 보정하는 것이 필요하다. 혼란변수가 있는 경우 오즈비(OR), 상대위험도(RR)를 보정하는 방법으로 코크란-멘텔-헨젤(Cochran-Mantel-Haenszel) 검정법이 있다.

3.1. 방법

층화분석 후, 오즈비와 상대위험도는 각각 다음의 방법으로 구할 수 있다.
질병 발생 질병 발생 안함 합계
위험 요인 노출 a b a+b
위험 요인 노출 안됨 c d c+d
합계 a+c b+d n
  • [math(OR=\dfrac{\sum \dfrac{a_id_i}{n_i}}{\sum \dfrac{b_ic_i}{n_i}})]
  • [math(RR=\dfrac{\sum \dfrac{a_i(c_i+d_i)}{n_i}}{\sum \dfrac{c_i(a_i+b_i)}{n_i}})]

3.2. 예시

3.2.1. 원 자료

폐암 발생 폐암 발생 안함 발생률
음주 250 400 약 0.385
음주 안함 80 270 약 0.229
(총 1000명)
위 표에서 OR은 약 2.11, RR은 약 1.68로 계산된다. 그런데, 폐암 발생에 영향을 미치는 중요한 요인은 흡연 여부이므로, 흡연 여부가 혼란 변수로 작용했을 수 있다.

3.2.2. 층화 분석된 자료

  • 흡연자
폐암 발생 폐암 발생 안함 발생률
음주 200 250 약 0.444
음주 안함 30 70 0.3
(총 550명)
흡연자에서 OR은 약 1.87, RR은 약 1.48로 계산된다.
  • 비흡연자
폐암 발생 폐암 발생 안함 발생률
음주 50 150 0.25
음주 안함 50 200 0.2
(총 450명)
비흡연자에서 OR은 약 1.33, RR은 약 1.25로 계산된다.

3.2.3. 코크란-멘텔-헨젤 검정 적용

코크란-멘텔-헨젤 방법으로 구한 OR과 RR의 값은 다음과 같다.
  • [math(OR=\dfrac{\dfrac{200×70}{550}+\dfrac{50×200}{450}}{\dfrac{250×30}{550}+\dfrac{150×50}{450}})], 약 1.57
  • [math(RR=\dfrac{\dfrac{200×(30+70)}{550}+\dfrac{50×(50+200)}{450}}{\dfrac{30×(200+250)}{550}+\dfrac{50×(50+150)}{450}})], 약 1.38

층화하기 전의 OR/RR은 각각 2.11/1.68이었으나, 흡연 여부로 층화하고 코크란-멘텔-헨젤 방법으로 보정한 OR/RR은 각각 1.57/1.38이므로, 흡연 여부가 혼란변수로 작용해 음주 여부와 폐암 발생의 관련성이 과대추정되도록 했음을 알 수 있다.

4. 기여위험도(AR)

질병 발생 질병 발생 안함 발생률
위험 요인 노출 a b [math(\dfrac{a}{a+b}=I_e)]
위험 요인 노출 안됨 c d [math(\dfrac{c}{c+d}=I_u)]
인구 집단 [math(I_t)]
위 표에서 기여위험도(attributable risk, AR)는 [math(I_e-I_u)]로, 두 집단의 발생률의 차이로 구한다. Risk difference(RD) 또는 absolute risk increase(ARI)라고도 부른다.
기여위험분율(attributable fraction, AF)는 [math(\dfrac{I_e-I_u}{I_e})]로 구한다. 이는 [math(1-\dfrac{I_u}{I_e})]와 같으므로, [math(1-\dfrac{1}{RR})]로도 구할 수 있다.
인구집단 기여위험도(population attributable risk, PAR)은 [math(I_t-I_u)]이다.
인구집단 기여위험분율(population attributable fraction, PAF)는 [math(\dfrac{I_t-I_u}{I_t})]이며, 인구 집단에서 요인에 노출된 사람의 비율이 [math(P_e)]일 때, [math(PAF=\dfrac{P_e(RR-1)}{P_e(RR-1)+1})]이다.

4.1. 예시

폐암 발생 폐암 발생 안함 발생률
흡연 250 350 [math(\dfrac{250}{600}≈0.417)]
흡연 안함 100 300 [math(\dfrac{100}{400}=0.25)]
인구 집단 0.35
위 표에서 기여위험도(AR)는 0.417-0.25=0.167이고, 기여위험분율(RR)은 0.4이다. 이를 흡연자의 폐암 중 40%가 흡연에 의한 것이라고 해석할 수 있다.
위 표에서 인구집단 기여위험도(PAR)는 0.35-0.25=0.1이고, 인구집단 기여위험분율(PAR)은 약 0.29이다. 이를 인구집단의 폐암 중 약 29%가 흡연에 의한 것이라고 해석할 수 있다.

4.2. 예방가능분율

예방가능분율(preventable fraction)은 기여위험분율과 거의 동일한 개념으로, 특정 요인에 노출됨으로써(예를 들어, 백신 접종) 질병이 예방되는 정도를 나타낸 지표이다.
질병 발생 질병 발생 안함 발생률
보호 요인 노출 a b [math(\dfrac{a}{a+b}=I_e)]
보호 요인 노출 안됨 c d [math(\dfrac{c}{c+d}=I_u)]
위 표에서 예방가능분율은 [math(\dfrac{I_u-I_e}{I_u})]로 구할 수 있다.
COVID-19 감염 COVID-19 감염 안됨 발생률
백신 접종 8 21712 약 0.0368%
위약 접종 162 21566 약 0.7456%
[1]
백신의 효과(vaccine efficacy)도 같은 방법으로 측정할 수 있다. 위 표의 예시의 경우 백신의 효과는 (0.7456-0.0368)/0.7456으로 약 95%이다.

5. 위험비(HR)


[1] Fernando P. Polack외, "Safety and Efficacy of the BNT162b2 mRNA Covid-19 Vaccine", N Engl J Med 2020; 383:2603-2615 원본 논문