本稿では、標本発症リスク差の漸近分布の導出を行っています。この漸近分布は、母集団発症リスク差の検定の実施や信頼区間の導出の基礎となります。
なお、閲覧にあたっては、以下の点にご注意ください。
- スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
- 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。
- 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
【定理】標本リスク差の漸近分布
【定理】
標本リスク差の漸近分布
Asymptotic Distribution of Sample Risk Difference
マッチングなしのコホート研究における発症リスク差と標本発症リスク差を
\begin{gather}
\delta=\mathrm{RD}=\pi_1-\pi_0\\
\hat{\delta}=\mathrm{\widehat{RD}}={\hat{\pi}}_1-{\hat{\pi}}_0
\end{gather}
帰無仮説と対立仮説をそれぞれ
\begin{gather}
H_0:\pi_1=\pi_0 \left(=\pi\right) \quad H_1:\pi_1 \neq \pi_0
\end{gather}
とするとき、
標本リスク差の漸近分布は、
〔1〕対立仮説
\begin{align}
H_1:\hat{\delta}\xrightarrow[]{d}\mathrm{N} \left[\pi_1-\pi_0,\frac{\pi_1 \left(1-\pi_1\right)}{n_1}+\frac{\pi_0 \left(1-\pi_0\right)}{n_0}\right]
\end{align}
標本リスク差の漸近分散の一致推定量は、
\begin{align}
{\hat{\sigma}}_1^2=\frac{{\hat{\pi}}_1 \left(1-{\hat{\pi}}_1\right)}{n_1}+\frac{{\hat{\pi}}_0 \left(1-{\hat{\pi}}_0\right)}{n_0}
\end{align}
〔2〕帰無仮説 \begin{align} H_0:\hat{\delta}\xrightarrow[]{d}\mathrm{N} \left[0,\pi \left(1-\pi\right) \left(\frac{1}{n_1}+\frac{1}{n_0}\right)\right] \end{align} 共通の母比率の推定量は、 \begin{align} \hat{\pi}=\frac{m_1}{N} \end{align} 標本リスク差の漸近分散の一致推定量は、 \begin{align} {\hat{\sigma}}_0^2=\frac{m_1m_0}{Nn_1n_0} \end{align}
導出
〔1〕対立仮説の下での漸近分布
二項分布の正規近似により、標本比率は漸近的に
\begin{gather}
{\hat{\pi}}_i \sim \mathrm{N} \left[\pi_i,\frac{\pi_i \left(1-\pi_i\right)}{n_i}\right]
\end{gather}
標本リスク差を $\hat{\delta}={\hat{\pi}}_1-{\hat{\pi}}_0$ とすると、正規分布の再生性より、
\begin{align}
\hat{\delta} \sim \mathrm{N} \left\{\pi_1-\pi_0,\frac{\pi_1 \left(1-\pi_1\right)}{n_1}+\frac{\pi_0 \left(1-\pi_0\right)}{n_0}\right\}
\end{align}
母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、
\begin{align}
{\hat{\sigma}}_1^2=\frac{{\hat{\pi}}_1 \left(1-{\hat{\pi}}_1\right)}{n_1}+\frac{{\hat{\pi}}_0 \left(1-{\hat{\pi}}_0\right)}{n_0}
\end{align}
〔2〕帰無仮説の下での漸近分布
$\pi_1=\pi_0=\pi$ を代入すると、帰無仮説のもとで、
\begin{align}
\hat{\delta}\xrightarrow[]{d}\mathrm{N} \left[0,\pi \left(1-\pi\right) \left(\frac{1}{n_1}+\frac{1}{n_0}\right)\right]
\end{align}
共通の母比率 $\pi$ の推定量は、
\begin{align}
\hat{\pi}=\frac{m_1}{N}
\end{align}
これを用いて、母比率を一致推定量である標本比率で置き換えると、漸近分散の一致推定量は、
\begin{align}
{\hat{\sigma}}_0^2&=\frac{m_1}{N} \cdot \frac{m_0}{N} \cdot \frac{n_1+n_0}{n_1n_0}\\
&=\frac{m_1}{N} \cdot \frac{m_0}{N} \cdot \frac{N}{n_1n_0}\\
&=\frac{m_1m_0}{Nn_1n_0}
\end{align}
$\blacksquare$
参考文献
- ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.23-24
0 件のコメント:
コメントを投稿