本稿では、希少疾患の仮定の下で、発症リスク比が曝露オッズ比で近似できることを証明しています。この関係は、ケース・コントロール研究の妥当性を保証するものとして非常に重要な基本事項です。
なお、閲覧にあたっては、以下の点にご注意ください。
- スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
- 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。
発症リスク比と曝露リスク比の関係
ケース・コントロール研究における、後ろ向き曝露リスク比は定義より、 \begin{gather} \mathrm{{\rm RR}_{Retro}}=\frac{\phi_1}{\phi_0} \end{gather} いっぽう、関心のある前向き発症リスク比の定義式は、 \begin{align} \mathrm{{\rm RR}_{Pro}}=\frac{P \left(D\middle| E\right)}{P \left(D\middle|\bar{E}\right)}=\frac{\pi_1}{\pi_0} \neq \frac{\phi_1}{\phi_0}=\frac{P \left(E\middle| D\right)}{P \left(E\middle|\bar{D}\right)}=\mathrm{{\rm RR}_{Retro}} \end{align} したがって、基本的に、後ろ向き曝露リスク比を用いて前向き発症リスク比を評価することはできない。
【定理】曝露オッズ比による発症リスク比の近似
【定理】
曝露オッズ比による発症リスク比の近似
Retrospective Odds Ratio Approximate the Risk Ratio under Rare Event Assumption
マッチングなしのケース・コントロール研究において、全体の有病割合 $\delta$ が分かっているとき、前向き発症リスク比の母数とその一致推定量は、 \begin{gather} \mathrm{RR}=\frac{\phi_1}{1-\phi_1} \cdot \frac{\delta \left(1-\phi_1\right)+ \left(1-\delta\right) \left(1-\phi_0\right)}{\delta\phi_1+ \left(1-\delta\right)\phi_0}\\ \mathrm{\widehat{RR}}=\frac{{\hat{\phi}}_1}{1-{\hat{\phi}}_1} \cdot \frac{\delta \left(1-{\hat{\phi}}_1\right)+ \left(1-\delta\right) \left(1-{\hat{\phi}}_0\right)}{\delta{\hat{\phi}}_1+ \left(1-\delta\right){\hat{\phi}}_0} \end{gather} で与えられる。 特に、希少疾患の仮定 \begin{gather} \delta\rightarrow0 \end{gather} の下では、 前向き発症リスク比は、後ろ向き曝露オッズ比で近似できる、すなわち、 \begin{gather} \lim_{\delta\rightarrow0}{\mathrm{RR}}=\mathrm{{\rm OR}_{Retro}} \end{gather} が成り立つ。
証明
全体の有病割合 $\delta$ を用いて得られた分割表を書き直すと、以下のようになる。
曝露あり $ \left(E\right)$ |
曝露なし $(\bar{E})$ | 合計 | |
---|---|---|---|
ケース群 $ \left(D\right)$ | $\delta\phi_1$ | $\delta \left(1-\phi_1\right)$ | $\delta$ |
コントロール群 $(\bar{D})$ | $ \left(1-\delta\right)\phi_0$ | $ \left(1-\delta\right) \left(1-\phi_0\right)$ | $1-\delta$ |
合計 |
$\delta\phi_1$ $+$ $ \left(1-\delta\right)\phi_0$ |
$\delta \left(1-\phi_1\right)$ $+$ $ \left(1-\delta\right) \left(1-\phi_0\right)$ | $1$ |
このデータがマッチングなしのコホート研究で得られたものであると仮定すると、曝露群と非曝露群の発症確率はそれぞれ \begin{gather} \pi_1=\frac{ \left(1-\delta\right)\phi_0}{\delta\phi_1+ \left(1-\delta\right)\phi_0}\\ \pi_0=\frac{\delta \left(1-\phi_1\right)}{\delta \left(1-\phi_1\right)+ \left(1-\delta\right) \left(1-\phi_0\right)} \end{gather} 前向き発症リスク比の定義式 $\mathrm{RR}=\frac{\pi_1}{\pi_0}$ より、 \begin{align} \mathrm{RR}&=\frac{\delta\phi_1}{\delta\phi_1+ \left(1-\delta\right)\phi_0} \cdot \frac{\delta \left(1-\phi_1\right)+ \left(1-\delta\right) \left(1-\phi_0\right)}{\delta \left(1-\phi_1\right)}\\ &=\frac{\phi_1}{1-\phi_1} \cdot \frac{\delta \left(1-\phi_1\right)+ \left(1-\delta\right) \left(1-\phi_0\right)}{\delta\phi_1+ \left(1-\delta\right)\phi_0} \end{align} この一致推定量は、 \begin{align} \mathrm{\widehat{RR}}=\frac{{\hat{\phi}}_1}{1-{\hat{\phi}}_1} \cdot \frac{\delta \left(1-{\hat{\phi}}_1\right)+ \left(1-\delta\right) \left(1-{\hat{\phi}}_0\right)}{\delta{\hat{\phi}}_1+ \left(1-\delta\right){\hat{\phi}}_0} \end{align}
多くの場合、全体の有病割合が事前には分からず、通常ケース・コントロール研究からは直接推定することができないため、このアプローチはほとんど実用的ではないが、希少疾患の仮定($\delta\rightarrow0$)の下では、 \begin{align} \lim_{\delta\rightarrow0}{\mathrm{RR}}&=\lim_{\delta\rightarrow0}{\frac{\phi_1}{1-\phi_1} \cdot \frac{\delta \left(1-\phi_1\right)+ \left(1-\delta\right) \left(1-\phi_0\right)}{\delta\phi_1+ \left(1-\delta\right)\phi_0}}\\ &=\frac{\phi_1}{1-\phi_1} \cdot \frac{1-\phi_0}{\phi_0}\\ &=\mathrm{{\rm OR}_{Retro}} \end{align} $\blacksquare$
参考文献
- ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.218-219
- Cornfield, J.A.. A method of estimating comparative rates from clinical data; applications to cancer of the lung, breast, and cervix. Journal of the National Cancer Institute. 1951, 11(6), p.1269-1275, PMID: 14861651
0 件のコメント:
コメントを投稿