ロジスティック・モデル

公開日:

【2022年11月2週】 【A000】生物統計学 【A051】コホート研究 【A090】多変量回帰モデル 【A092】ロジスティック回帰分析

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

本稿では、関心のあるイベントの発生確率の予測などに用いられるロジスティック回帰モデルのうち、最も単純(共変量が1つの場合)なロジスティック・モデルについて、定義式の紹介、発症確率とロジスティック関数の関係、対数オッズ比と回帰係数の関係の証明などを行っています。

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

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

ロジスティック・モデル

マッチングなしのコホート研究において、対数発症オッズ(発症確率のロジット)について、 log(π11π1)=α+βπ11π1=eα+βlog(π01π0)=απ01π0=eα が成り立つとする 統計モデルをロジット連結モデル logit link model、または、ロジスティック・モデル logistice model という。

このとき、曝露群・非曝露群の発症確率は、ロジスティック関数として、 π1=eα+β1+eα+β1π1=11+eα+βπ0=eα1+eα1π0=11+eα で与えられる。 また、発症オッズ比 φ=OR は、 logφ=βφ=eβ で与えられる。

証明:ロジスティック関数

非曝露群のオッズの定義式より、 π01π0=eαπ0=eα(1π0)(1+eα)π0=eαπ0=eα1+eα 1π0=1eα1+eα=1+eαeα1+eα=11+eα 同様に、曝露群のオッズの定義式より、 π11π1=eα+βπ1=eα+β(1π1)(1+eα+β)π1=eα+βπ1=eα+β1+eα+β 1π1=1eα+β1+eα+β=1+eα+βeα+β1+eα+β=11+eα+β

証明:対数オッズ比と回帰係数の関係

対数オッズ比の定義式より、 logφ=log(π11π11π0π0)=logπ11π1logπ01π0=logeα+βlogeα=α+βα=β 両辺の指数を取ると、 β=logφ

対数リスクモデル

マッチングなしのコホート研究において、曝露群・非曝露群の発症確率について、 logπ1=α+βlogπ0=α が成り立つとする統計モデルを対数リスクモデルという。

このとき、曝露群・非曝露群の発症確率は、 π1=eα+β1π0=1eα+βπ2=eα1π0=1eα で与えられる。 また、発症リスク比 φ=RR は、 logRR=βRR=eβ で与えられる。

参考文献

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

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ