生存時間の差に関する優越性試験のサンプルサイズ設計の公式

公開日:

【2022年11月3週】 【A000】生物統計学 【A051】コホート研究 【A074】サンプルサイズの設計 【A100】生存時間分析 【A102】ログランク検定

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

本稿では、生存時間の差の検定としてログランク検定を行う優越性試験のサンプルサイズ設計の公式行っています。本稿の方法は、生存時間の分布に指数分布モデルを仮定した方法であり、現実的にはこの仮定は必ずしも満たされるわけではありませんが、最も簡単に導出できる方法としてよく使用されています。

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

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

【公式】生存時間の差に関する優越性試験のサンプルサイズ設計の公式

【公式】
生存時間の差に関する優越性試験のサンプルサイズ設計の公式
Sample Size Determination for Log-Rank Test

生存時間に関するコホート研究について、
生存時間の分布に指数分布モデル λ(t)=λ0<t を仮定し、 各被験者の登録時期が一様分布 RiU[0,R] に従い、追跡調査のランダムな終了はないとする。

各群のサンプルサイズを n1=n0=N 帰無仮説と対立仮説を H0:λ1λ0=0H1:λ1λ0=Δ とする。

介入群と非介入群のハザード λ1λ0 を具体的定数として与え、 有意水準と検出力をそれぞれ、 α1β 帰無仮説の下での共通のハザードを λ¯=λ1+λ02 とするとき、 N=[Zαϕ0Z1βϕ1λ1λ0]2ϕ02=2λ¯2E(δ|λ¯)ϕ12=λ12E(δ|λ1)+λ02E(δ|λ1)E(δ|λ)=1eλ(SR)eλSλR

導出

導出

生存時間分布の指数分布モデルとイベント発生回数の二重同次ポアソンモデルの同等性から、生存時間におけるハザードを平均発生率と同じであると仮定する。

このとき、母集団の平均発生率の差と標本平均発生率の差をそれぞれ θ=λ1λ0θ^=λ^1λ^0 とおくと、 標本平均発生率の差の漸近分布は、対立仮説の下で θ^N[λ1λ0,λ12E(d1|λ1)+λ02E(d0|λ0)] 帰無仮説の下で θ^N[0,λ2E(d1|λ)+λ2E(d0|λ)] 帰無仮説の下での共通のハザードの推定値は、 λ¯=λ1+λ02 二重同次ポアソンモデルの下で、 E(dj|λj)=NE(δi|λj) また、定数ハザードが与えられた下で研究期間中にイベントが発生する確率について、 E(δi|λj)=1eλj(SR)eλjSλjR したがって、検定統計量の分散について、 ϕ02N=λ¯2NE(δi|λ)+λ¯2NE(δi|λ)ϕ02=2λ¯2E(δi|λ) ϕ12N=λ12NE(δi|λ1)+λ02NE(δi|λ1)ϕ12=λ12E(δi|λ1)+λ02E(δi|λ1) これをサンプルサイズの設計公式に代入すると、 N=[Zαϕ0Z1βϕ1λ1λ0]2

参考文献

  • 山口 拓洋 著. サンプルサイズの設計 後悔先に立たず. 健康医療評価研究機構, 2010, p.79-82
  • ジョン・ラチン 著, 宮岡 悦良 監訳, 遠藤 輝, 黒沢 健, 下川 朝有, 寒水 孝司 訳. 医薬データのための統計解析. 共立出版, 2020, p.520-524

関連記事

自己紹介

自分の写真

yama

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

このブログを検索

ブログ アーカイブ

QooQ