標本対数平均発生率比の漸近分布

公開日:

【2022年10月4週】 【A000】生物統計学 【A051】コホート研究 【A071】標本分布

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

本稿では、平均発生率を指標とするコホート研究における標本対数平均発生率比の漸近分布の導出を行っています。この漸近分布は、母集団における平均発生率比の信頼区間を導出するうえでの基礎となります。

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

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

【定理】標本対数平均発生率比の漸近分布

【定理】
標本対数平均発生率比の漸近分布
Asymptotic Distribution of Sample Log Incidence Rate Ratios

平均発生率を指標とするマッチングなしのコホート研究における対数平均発生率比と標本対数平均発生率比を δ=logIRR=logλ1λ0δ^=logIRR^=logλ^1λ^0 とするとき、 標本対数平均発生率比の漸近分布は、 δ^dN[logλ1λ0,1λ1T1+1λ0T0] 漸近分散の一致推定量は、 σ^12=1d1+1d0

証明

証明

それぞれの標本平均発生率を以下のようにおくと、 λ^1=d1T1λ^0=d0T0 ポアソン分布の正規近似により、漸近的に λ^iN(λi,λiTi) ここで、 g(λi)=logλig(λ^i)=logλ^i と変数変換する。 デルタ法を用いて、g(λ^i) を期待値 E(λ^i)=λi まわりでテイラー展開すると、g(λ^i) の1階微分は、 g(λ^i)=1λ^i よって、デルタ法における期待値と分散の公式より、 E{g(λ^i)}E[g(λi)]=logλi V[g(λ^i)]{g(πi)}2V(λ^i)=1λi2λiTi=1λiTi スラツキーの定理より、 logλ^idN[logλi,1λiTi]

母集団の対数平均発生率比と標本対数平均発生率比をそれぞれ δ=logIRR=logλ1λ0δ^=logIRR^=logλ^1λ^0 とおくと、 正規分布の再生性より、 δ^N(logλ1λ0,1λ1T1+1λ0T0) 母集団の平均発生率を一致推定量である標本平均発生率で置き換えると、漸近分散の一致推定量は、 σ^12=1λ^1T1+1λ^0T0=1d1+1d0

参考文献

  • ケネス・ロスマン 著, 矢野 栄二, 橋本 英樹, 大脇 和浩 監訳. ロスマンの疫学. 篠原出版新社, 2013, p.236-238
  • 丹後 俊郎, 松井 茂之 編集. 医学統計学ハンドブック. 朝倉書店, 2018, p.509

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ