有病率・発生割合の信頼区間(補対数対数変換)

公開日:

【2022年10月5週】 【A000】生物統計学 【A051】コホート研究 【A072】統計的推定

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

本稿では、横断研究やコホート研究における有病率・発生割合の信頼区間の導出を行っています。本稿における方法は、母比率の補対数対数変換にもとづく方法で、非対称な信頼区間とはなりますが、上下信頼限界が区間 [0,1] に収まる方法です。

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

  • スマートフォンやタブレット端末でご覧の際、数式が見切れている場合は、横にスクロールすることができます。
  • 漸近的な性質を用いる際は、①中心極限定理が成り立つ、②漸近分散を推定する際に、母数をその一致推定量で置き換えることができるということが成り立つと仮定しています。
  • デルタ法を用いる際、剰余項(2次の項)が漸近的に無視できる(0に確率収束する)と仮定しています。

【定理】有病率・発生割合の信頼区間(補対数対数変換)

【定理】
有病率・発生割合の信頼区間(補対数対数変換)
Confidence Interval for Prevalence or Incidence Proportion Based on Complementary log-log-Transformation

マッチングなしのコホート研究における有病率・発生割合を補対数対数変換したものを θ=log(logπ)θ^=log(logπ^) とするとき、 補対数対数変換にもとづく母集団の有病率・発生割合の非対称な 100(1α)% 信頼区間は、漸近的に 100(1α)% C.I.=[πL,πU]πL=exp{exp(θU)}πU=exp{exp(θL)} で与えられる。 ただし、 θL=θ^σ^Z0.5αθU=θ^+σ^Z0.5ασ^2=1πnπ(logπ)2σ^=S.E.^(θ^)

なお、θ の上側信頼限界 θU から π の下側信頼限界 πL が得られ、θ の下側信頼限界 θL から π の上側信頼限界 πU が得られる点に注意する。

導出:母比率の補対数対数変換にもとづく方法

導出

標本補対数対数比率の漸近分布は、 θ^dN[log(logπ),1πnπ(logπ)2] 標本補対数対数比率の漸近分散の一致推定量は、 σ^12=1π^nπ^(logπ^)2 標本補対数対数比率の漸近的な標準誤差 S.E.^(θ^)=σ^ は、 σ^=1π^nπ^|logπ^| 母集団の対数オッズの 100(1α)% 信頼区間は、 θ^σ^Z0.5αlog(logπ)θ^+σ^Z0.5α ここで、スペースの節約のために以下のようにおいて、 θL=θ^σ^Z0.5αθU=θ^+σ^Z0.5α 逆変換を行うと、 θLlog(logπ)θUexp(θL)logπexp(θU)exp(θU)logπexp(θL)exp{exp(θU)}πexp{exp(θL)}

θ の上側信頼限界 θU から π の下側信頼限界 πL が得られ、θ の下側信頼限界 θL から π の上側信頼限界 πU が得られる点に注意する。

参考文献

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

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ