横断研究・コホート研究【有病率・発生割合】(マッチングあり・層化あり)

公開日: 更新日:

【2022年10月1週】 【A000】生物統計学 【A050】研究デザイン 【A051】コホート研究

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

本稿では、横断研究・コホート研究の研究デザインのうち、①有病率(横断研究)や発生割合(コホート研究)を曝露効果の指標とする、②マッチングあり、③層化ありのデザイン・パターンについて、その分割表の形式、統計モデル、曝露効果の指標の定義をまとめています。

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

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

分割表の形式

指示変数 j を任意の被験者の曝露状況を表す変数 j={1Exposed(E)0Unexposed(E¯) とし、 曝露者1人に対し、背景因子の水準が同程度の非曝露者を1人マッチングし、計2人のペアを作る。そして、調整したい交絡因子の水準にもとづいて、互いに独立な K 個の層に層化する。
このペアの総数(サンプルサイズ)を Nk とし、発症状況について調べる。 指示変数 Yi 番目のペアの j 番目のメンバーの発症状況を表す変数 yijk={1Disease(D)0Not disease(D¯) とする。

このとき、各ペアの曝露・発症状況は、 曝露者・非曝露者=①発症あり・発症あり、②発症あり・発症なし、
③発症なし・発症あり、④発症なし・発症なし
(yi1k,yi0k)=(1,1),(1,0),(0,1),(0,0) のいずれかに分類される。 それぞれの曝露・発症状況に該当するペアの数を ek(=n11k)fk(=n12k)gk(=n21k)hk(=n22k) とする。

また、周辺度数として、 ①曝露者が発症したペア、②曝露者が発症しなかったペア、
③非曝露者が発症したペア、④非曝露者が発症しなかったペア
が得られる。 それぞれの合計ペア数を nEDk(=n1k)nED¯k(=n0k)nE¯Dk(=n1k)nE¯D¯k(=n0k) とする。

表1 ペア・マッチングを行った横断研究・コホート研究に関する
2×2 分割表(第 k 層の観測値)
非曝露者
(E¯)
合計
発症あり
(D)
発症なし
(D¯)
曝露者
(E)
発症あり
(D)
ek
(=n11k)
fk
(=n12k)
nEDk
(=n1k)
発症なし
(D¯)
gk
(=n21k)
hk
(=n22k)
nED¯k
(=n0k)
合計 nE¯Dk
(=n1k)
nE¯D¯k
(=n0k)
Nk

統計モデル

k 層における各セルの観測値 nk=(ekfkgkhk) が四項分布

nkMN(Nk,πk)πk=(π11kπ12kπ21kπ22k) に従うとする。

表2 マッチングを行った横断研究・コホート研究に関する
2×2 分割表(第 k 層の統計モデル)
非曝露者
(E¯)
合計
発症あり
(D)
発症なし
(D¯)
曝露者
(E)
発症あり
(D)
π11k π12k πEDk
(=π1k)
発症なし
(D¯)
π21k π22k πED¯k
(=π0k)
合計 πE¯Dk
(=π1k)
πE¯D¯k
(=π0k)
1

四項尤度

L(πk)=Nk!ek!fk!gk!hk!π11kekπ12kfkπ21kgkπ22khk

曝露効果の指標

標本比率

π^11k=ekNkπ^12k=fkNkπ^21k=gkNkπ^22k=hkNk

条件付き周辺発症オッズ比

ORzk=π12k|zπ21k|z

条件付き発症オッズ比

ORCk=π12kπ21kOR^Ck=π^12kπ^21k=fkgk

母集団平均発症リスク比

RRAk=π1kπ1k=π11k+π12kπ11k+π21kRR^Ak=π^1kπ^1k=π^11k+π^12kπ^11k+π^21k=ek+fkek+gk=n1kn1k

参考文献

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

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ