マンテル・ヘンツェル検定

公開日:

【2022年11月1週】 【A000】生物統計学 【A051】コホート研究 【A052】ケース・コントロール研究 【A073】統計的仮説検定

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

本稿では、単層の2×2分割表に対するマンテル・ヘンツェル検定統計量の導出を行っています。この検定統計量は、プール化による層別解析を行う際に用いられるマンテル・ヘンツェル検定統計量に拡張されます。

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

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

【定理】マンテル・ヘンツェル検定

【定理】
マンテル・ヘンツェル検定
Mantel-Haenszel Test

マッチングなしのコホート研究(超幾何分布モデル)における曝露群と非曝露群の発症確率 π1π0 に関する 帰無仮説と対立仮説を H0:π1=π0(=π)H1:π1π0 とするとき、 マンテル・ヘンツェル検定統計量は、漸近的に χMH2=(an1m1N)2m1m0n1n0N2(N1) で与えられる。

導出

導出

超幾何分布モデルを仮定するとき、曝露群の発症者数について、 aHG(N,n1,m1) 超幾何分布の期待値と分散は、 E(a)=n1m1NV(a)=m1m0n1n0N2(N1) 超幾何分布の正規近似(超幾何分布の二項近似→二項分布の正規近似)により、漸近的に aN[E(a),V(a)] 標準化変換の性質より、 aE(a)V(a)=Z0N(0,1) したがって、χ2分布の定義より、 χMH2={aE(a)}2V(a)=(an1m1N)2m1m0n1n0N2(N1) これは、 χMH2χ2(1)

参考文献

  • ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.41
  • Mantel, N. & Haenszel, W.. Statistical Aspects of the Analysis of Data From Retrospective Studies of Disease. Journal of the National Cancer Institute. 1959, 22(4), p.719-748, doi: 10.1093/jnci/22.4.719
  • Mantel, N.. Chi-square tests with one degree of freedom: Extensions of the Mantel-Haenszel procedure. Journal of the American Statistical Association. 1963, 58, p.690-700, doi: 10.2307/2282717

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ