母比率の滑らかな関数にもとづく検定

公開日: 更新日:

【2022年11月1週】 【A000】生物統計学 【A051】コホート研究 【A073】統計的仮説検定

この記事をシェアする
  • B!
サムネイル画像

本稿では、対数リスク比や対数オッズ比といった母比率の滑らかな関数にもとづく検定に用いる検定統計量の導出を行っています。

なお、閲覧にあたっては、以下の点にご注意ください。

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
  • 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。
  • 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
  • デルタ法を用いる際、剰余項(2次の項)が漸近的に無視できる($0$に確率収束する)と仮定しています。

【定理】母比率の滑らかな関数にもとづく検定

【定理】
母比率の滑らかな関数にもとづく検定
Z-Test Based on Smooth Function of Population Proportions

対数リスク比や対数オッズ比といった母比率 $\pi_i$ の滑らかな関数 $g \left(\pi_i\right)$ に関する帰無仮説 \begin{align} H_0:\pi_1=\pi_0 \left(=\pi\right)\Leftrightarrow H_0:g \left(\pi_1\right)=g \left(\pi_0\right) \left\{=g \left(\pi\right)\right\} \end{align} に対する $\mathrm{Z}$検定統計量は、 \begin{align} Z_g=\frac{g \left(p_1\right)-g \left(p_2\right)}{\sqrt{\hat{V} \left\{g \left(p_1\right)-g \left(p_2\right)\middle| H_0\right\}}} \end{align} で与えられる。

証明

証明

二項分布の正規近似により、標本比率は漸近的に \begin{gather} p_i \sim \mathrm{N} \left[\pi_i,\frac{\pi_i \left(1-\pi_i\right)}{n_i}\right] \end{gather} ここで、 \begin{gather} \pi_i\rightarrow g \left(\pi_i\right) \quad p_i\rightarrow g \left(p_i\right) \end{gather} と変数変換する。 デルタ法における期待値と分散の公式より、 \begin{align} E \left\{g \left(p_i\right)\right\}\cong E \left[g \left(\pi_i\right)\right] \end{align} \begin{align} V \left[g \left(p_i\right)\right]\cong \left\{g^\prime \left(\pi_i\right)\right\}^2V \left(p_i\right) \end{align} スラツキーの定理より、 \begin{align} g \left(p_i\right)\xrightarrow[]{d}\mathrm{N} \left[g \left(\pi_i\right), \left\{g^\prime \left(\pi_i\right)\right\}^2V \left(p_i\right)\right] \end{align} したがって、正規分布の再生性より、 \begin{align} g \left(p_1\right)-g \left(p_0\right)\xrightarrow[]{d}\mathrm{N} \left[g \left(\pi_1\right)-g \left(\pi_0\right), \left\{g^\prime \left(\pi_1\right)\right\}^2V \left(p_1\right)+ \left\{g^\prime \left(\pi_0\right)\right\}^2V \left(p_0\right)\right] \end{align} 帰無仮説 $H_0:\pi_1=\pi_0 \left(=\pi\right)$ のもとでは、 \begin{align} g \left(p_1\right)-g \left(p_0\right)\xrightarrow[]{d}\mathrm{N} \left[0, \left\{g^\prime \left(\pi\right)\right\}^2 \left\{V \left(p_1\right)+V \left(p_0\right)\right\}\right] \end{align} ここで、帰無仮説のもとでの漸近分散を以下のようにおくと、 \begin{align} \sigma_0^2= \left\{g^\prime \left(\pi\right)\right\}^2 \left\{V \left(p_1\right)+V \left(p_0\right)\right\} \end{align} これを用いて標準化した値は、 \begin{align} \frac{g \left(p_1\right)-g \left(p_2\right)}{{\hat{\sigma}}_0}=Z_g \sim \mathrm{N} \left(0,1\right) \end{align} したがって、母比率の滑らかな関数にもとづく検定統計量は、漸近的に通常の $\mathrm{Z}$検定を行うための検定統計量として用いることができる。 $\blacksquare$

参考文献

  • ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.38-40

関連記事

自己紹介

自分の写真

yama

大学時代に読書の面白さに気づいて以来、読書や勉強を通じて、興味をもったことや新しいことを学ぶことが生きる原動力。そんな人間が、その時々に学んだことを備忘録兼人生の軌跡として記録しているブログです。

このブログを検索

ブログ アーカイブ

QooQ