マンテル・ヘンツェル推定量

公開日: 更新日:

【2022年10月3週】 【A000】生物統計学 【A051】コホート研究 【A052】ケース・コントロール研究 【A062】層別解析

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

本稿では、コホート研究やケース・コントロール研究のデータを層別し、全層に共通した曝露効果の指標を求める際のマンテル・ヘンツェル推定量の公式を紹介しています。

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

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
  • 曝露(発症)状況を表す右下の添え字は、「0」である場合($n_0,\pi_0$ など)や「2」である場合($n_2,\pi_2$ など)がありますが、どちらも「非曝露群(コントロール群)」を表しています。

【公式】マンテル・ヘンツェル推定量(発生リスク差)

【公式】
マンテル・ヘンツェル推定量(発生リスク差)
Mantel-Haenszel Estimator for Common Risk Difference

マッチングなし・層化ありのコホート研究において、全層に共通した発生リスク差のマンテル・ヘンツェル推定量は、 \begin{gather} \mathrm{{\widehat{RD}}_{MH}}=\frac{\sum_{k=1}^{K}\frac{a_kn_{0k}-b_kn_{1k}}{N_k}}{\sum_{k=1}^{K}\frac{n_{0k}n_{1k}}{N_k}} \end{gather} 推定量の分散は、 \begin{gather} V \left(\mathrm{{\widehat{RD}}_{MH}}\right)=\frac{\sum_{k=1}^{K} \left[ \left(\frac{n_{0k}n_{1k}}{N_k}\right)^2 \left\{\frac{a_kc_k}{n_{1k}^2 \left(n_{1k}-1\right)}+\frac{b_kd_k}{n_{0k}^2 \left(n_{0k}-1\right)}\right\}\right]}{ \left(\sum_{k=1}^{K}\frac{n_{0k}n_{1k}}{N_k}\right)^2} \end{gather} で与えられる。

【公式】マンテル・ヘンツェル推定量(発生リスク比)

【公式】
マンテル・ヘンツェル推定量(発生リスク比)
Mantel-Haenszel Estimator for Common Risk Ratio

マッチングなし・層化ありのコホート研究において、全層に共通した発生リスク比のマンテル・ヘンツェル推定量は、 \begin{gather} \mathrm{{\widehat{RR}}_{MH}}=\frac{\sum_{k=1}^{K}\frac{a_kn_{0k}}{N_k}}{\sum_{k=1}^{K}\frac{b_kn_{1k}}{N_k}} \end{gather} 推定量の対数の分散は、 \begin{gather} V \left(\log{\mathrm{{\widehat{RR}}_{MH}}}\right)=\frac{\sum_{k=1}^{K} \left(\frac{n_{0k}n_{1k}m_{1k}}{N_k^2}-\frac{a_kb_k}{N_k}\right)}{ \left(\sum_{k=1}^{K}\frac{a_kn_{0k}}{N_k}\right) \left(\sum_{k=1}^{K}\frac{b_kn_{1k}}{N_k}\right)} \end{gather} で与えられる。

【公式】マンテル・ヘンツェル推定量(発生率差)

【公式】
マンテル・ヘンツェル推定量(発生率差)
Mantel-Haenszel Estimator for Common Incidence Rate Difference

平均発生率を指標とするマッチングなし・層化ありのコホート研究において、全層に共通した発生率差のマンテル・ヘンツェル推定量は、 \begin{gather} \mathrm{{\widehat{IRD}}_{MH}}=\frac{\sum_{k=1}^{K}\frac{a_kT_{0k}-b_kT_{1k}}{T_k}}{\sum_{k=1}^{K}\frac{T_{0k}T_{1k}}{T_k}} \end{gather} 推定量の分散は、 \begin{gather} V \left(\mathrm{{\widehat{IRD}}_{MH}}\right)=\frac{\sum_{k=1}^{K} \left[ \left(\frac{T_{0k}T_{1k}}{T_k}\right)^2 \left(\frac{a_k}{T_{1k}^2}+\frac{b_k}{T_{0k}^2}\right)\right]}{ \left(\sum_{k=1}^{K}\frac{T_{0k}T_{1k}}{T_k}\right)^2} \end{gather} で与えられる。

【公式】マンテル・ヘンツェル推定量(発生率比)

【公式】
マンテル・ヘンツェル推定量(発生率比)
Mantel-Haenszel Estimator for Common Incidence Rate Ratio

平均発生率を指標とするマッチングなし・層化ありのコホート研究において、全層に共通した発生率比のマンテル・ヘンツェル推定量は、 \begin{gather} \mathrm{{\widehat{IRR}}_{MH}}=\frac{\sum_{k=1}^{K}\frac{a_kT_{0k}}{T_k}}{\sum_{k=1}^{K}\frac{b_kT_{1k}}{T_k}} \end{gather} 推定量の対数の分散は、 \begin{gather} V \left(\log{\mathrm{{\widehat{IRR}}_{MH}}}\right)=\frac{\sum_{k=1}^{K} \left(\frac{m_{1k}T_{0k}T_{1k}}{T_k^2}\right)}{ \left(\sum_{k=1}^{K}\frac{a_kT_{0k}}{T_k}\right) \left(\sum_{k=1}^{K}\frac{b_kT_{1k}}{T_k}\right)} \end{gather} で与えられる。

【公式】マンテル・ヘンツェル推定量(曝露オッズ比)

【公式】
マンテル・ヘンツェル推定量(曝露オッズ比)
Mantel-Haenszel Estimator for Common Odds Ratio

マッチングなし・層化ありのケース・コントロール研究において、全層に共通した曝露オッズ比のマンテル・ヘンツェル推定量は、 \begin{gather} \mathrm{{\widehat{OR}}_{MH}}=\frac{\sum_{k=1}^{K}\frac{a_kd_k}{N_k}}{\sum_{k=1}^{K}\frac{b_kc_k}{N_k}} \end{gather} 推定量の対数の分散は、 \begin{gather} V \left(\log{\mathrm{{\widehat{OR}}_{MH}}}\right)=\frac{S_3}{2S_1^2}+\frac{S_5}{2S_1S_2}+\frac{S_4}{2S_2^2} \end{gather} で与えられる。 ただし、 \begin{gather} S_1=\sum_{k=1}^{K}\frac{a_kd_k}{N_k} \quad S_2=\sum_{k=1}^{K}\frac{b_kc_k}{N_k}\\ S_3=\sum_{k=1}^{K}\frac{a_kd_k \left(a_k+d_k\right)}{N_k^2} \quad S_4=\sum_{k=1}^{K}\frac{b_kc_k \left(b_k+c_k\right)}{N_k^2}\\ S_5=\sum_{k=1}^{K}\frac{ \left(a_k+d_k\right)b_kc_k+ \left(b_k+c_k\right)a_kd_k}{N_k^2} \end{gather}

参考文献

  • ケネス・ロスマン 著, 矢野 栄二, 橋本 英樹, 大脇 和浩 監訳. ロスマンの疫学:科学的思考への誘い 第2版. 篠原出版新社, 2013, p.248-262, p.268-269
  • 丹後 俊郎, 松井 茂之 編集. 医学統計学ハンドブック 新版. 朝倉書店, 2018, p.514-517
  • 丹後 俊郎, 小西 貞則 編集. 医学統計学の事典 新装版. 朝倉書店, 2018, p.99-100
  • ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.135-139
  • 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
  • Robins, J., Breslow, N. & Greenland, S.. Estimators of the Mantel-Haenszel variance consistent in both sparse data and large-strata limiting models. Biometrics. 1986, 42(2), p.311-323, doi: 10.2307/2531052
  • Robins, J., Greenland, S. & Breslow, N.. A general estimator for the variance of the Mantel-Haenszel odds ratio. Am J Epidemiol. 1986, 124(5), p.719-723, doi: 10.1093/oxfordjournals.aje.a114447

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ