ログランク検定のソースを表示

'''ログランク検定'''（ログランクけんてい、{{Lang-en-short|logrank test, or log-rank test}}）は、2つの[[標本 (統計学)|標本]]の[[生存関数|生存分布]]を比較する[[仮説検定]]である。これは[[ノンパラメトリック手法|ノンパラメトリック]]検定で、データが右に歪んで[[打ち切り (統計学)|打ち切られている]]場合に使用するのに適している（技術的には、打ち切りは情報を与えないものでなければならない）。この検定は、[[臨床試験]]において、新しい治療法の有効性を対照群と比較して確立するため広く使用されており、測定する対象は事象発生までの時間（初回治療から心臓発作までの時間など）である。この検定は、{{仮リンク|ネイサン＝マンテル|en|Nathan Mantel}}と[[デイヴィッド・コックス]]にちなんで'''マンテル＝コックス検定'''と呼ばれることもある。ログランク検定は、時間層別化された[[コクラン・マンテル・ヘンツェルの統計量|コクラン＝マンテル＝ヘンツェル検定]]<!-- time-stratified Cochran–Mantel–Haenszel test -->と見なすこともできる。

この検定は、ネイサン＝マンテルによって最初に提案され、[[リチャード・ピート]]と{{仮リンク|ジュリアン＝ピート|en|Julian Peto}}によってログランク検定と命名された<ref name="Mantel1966">{{cite journal|author=Mantel, Nathan|author-link=Nathan Mantel|year=1966|title=Evaluation of survival data and two new rank order statistics arising in its consideration.|journal=Cancer Chemotherapy Reports|volume=50|issue=3|pages=163–70|pmid=5910392}}</ref><ref name="Peto1972">{{cite journal|author=Peto, Richard|author-link=Richard Peto|year=1972|title=Asymptotically Efficient Rank Invariant Test Procedures|journal=Journal of the Royal Statistical Society, Series A|volume=135|issue=2|pages=185–207|publisher=Blackwell Publishing|doi=10.2307/2344317|jstor=2344317|author2=Peto, Julian|hdl=10338.dmlcz/103602|hdl-access=free}}</ref><ref>
{{Cite book|title=Encyclopedia of Biostatistics|chapter=Linear Rank Tests in Survival Analysis|first=David|last=Harrington|doi=10.1002/0470011815.b2a11047|year=2005|publisher=Wiley Interscience|isbn=047084907X}}</ref>。

== 定義 ==
ログランク検定統計量は、観察された各事象の時間における2つのグループの[[ハザード関数]]の推定値を比較する。これは、観察された事象の時間ごとに、いずれかのグループでの事象の観測数と期待数を計算し、これらを加算して、事象があったすべての時点にわたる全体的な要約を得ることによって構築される。

患者の2つのグループ、たとえば治療群と対照群を考える。 どちらかのグループで観察された事象の明確な時間を <math>1, \ldots, J</math>とする。<math>N_{1,j}</math> および <math>N_{2,j}</math> を、それぞれのグループ内における、時間 <math>j</math> の開始時点での「リスクがある」（まだ事象が発生していない、または打ち切られていない）被験者の数とする。また、<math>O_{1,j}</math> および <math>O_{2,j}</math> を、時間 <math>j</math> における各群で観測された事象の数とする。最後に、 <math>N_j = N_{1,j} + N_{2,j}</math> と <math>O_j = O_{1,j} + O_{2,j}</math> を定義する。

[[帰無仮説]]は、2つのグループのハザード関数が同一であるというもので、<math> H_0 : h_1(t) = h_2(t)</math> となる。したがって、<math>H_0</math> の下では、各グループ <math>i = 1, 2</math> に対して、<math>O_{i,j}</math> はパラメータ <math>N_j</math>, <math>N_{i,j}</math>, <math>O_j</math> を持つ[[超幾何分布]]に従う。この分布は、期待値が <math>E_{i,j} = N_{i,j} \frac{O_j}{N_j}</math>、分散が <math>V_{i,j} = E_{i,j} \left( \frac{N_j - O_j}{N_j} \right) \left( \frac{N_j - N_{i,j}}{N_j - 1} \right)</math> である。

ログランク統計量は、すべての <math>j = 1, \ldots, J</math> について、<math>O_{i,j}</math> と <math>H_0</math> のもとでの期待値 <math>E_{i,j}</math> と比較するものである。これは <math>Z = \frac {\sum_{j=1}^J (O_{i,j} - E_{i,j})} {\sqrt {\sum_{j=1}^J V_{i,j}}}\ \xrightarrow{d}\ \mathcal N(0,1)</math> &nbsp;&nbsp;&nbsp;&nbsp; (<math>i=1</math> または <math>2</math> の場合) として定義されている。

[[中心極限定理]]により、<math>Z</math> の分布は、<math>J</math> が無限に近づくにつれて標準正規分布の分布に収束するため、十分に大きな <math>J</math> に対しては標準正規分布で近似することができる。Peto and Petoの論文の付録Bで記述されているように、この量を、ピアソンの第1種ベータ分布または第2種ベータ分布（最初の4つのモーメントを一致させる）と等しくすることにより、より良い近似が得られる<ref name="Peto1972" />。

== 漸近分布 ==
2つのグループが同じ[[生存関数]]を持つ場合、ログランク統計量は近似的に標準正規分布となる。片側レベル <math>\alpha</math> 検定は、<math>Z>z_\alpha</math> ならば帰無仮説を棄却する。ここで <math>z_\alpha</math> は、標準正規分布の上位 <math>\alpha</math> 分位点である。[[ハザード比]]を <math>\lambda</math>、被験者総数を <math>n</math> 人、どちらかの群の被験者が最終的に事象を起こす確率を <math>d</math> （したがって、<math>nd</math> は分析時の事象の期待数）、各群に無作為に割り振られた被験者の割合を50%とすると、ログランク統計量は平均 <math> (\log{\lambda}) \, \sqrt {\frac {n \, d} {4}} </math>、分散 1 の近似正規分布となる<ref>{{cite journal|last=Schoenfeld|first=D|year=1981|title=The asymptotic properties of nonparametric tests for comparing survival distributions|journal=Biometrika|volume=68|issue=1|pages=316–319|doi=10.1093/biomet/68.1.316|jstor=2335833}}</ref>。検出力 <math>1-\beta</math> の片側レベル <math>\alpha</math> 検定の場合、必要な標本サイズは <math> n = \frac {4 \, (z_\alpha + z_\beta)^2 } {d\log^2{\lambda}}</math> となり、ここに <math>z_\alpha</math> と <math>z_\beta</math> は標準正規分布の分位数である。

== 同時分布 ==

<math> Z_1 </math> および <math> Z_2 </math> を、同じ検定の2つの異なる時点でのログランク統計量であるとする（<math> Z_1 </math> が先）。ここでも、2つのグループのハザード関数がハザード比 <math>\lambda</math> に比例し、<math> d_1 </math> と <math> d_2 </math> は <math> d_1  \leq d_2 </math> の2つの時点で被験者が事象を起こす確率であると仮定する。<math> Z_1 </math> および <math> Z_2 </math> は、平均 <math> \log{\lambda} \, \sqrt {\frac {n \, d_1} {4}} </math> と <math> \log{\lambda} \, \sqrt {\frac {n \, d_2} {4}} </math> 、相関 <math>\sqrt {\frac {d_1} {d_2}} </math> を持つ近似二変量正規分布である。{{仮リンク|データモニタリング委員会|en|Data monitoring committee}}による1つの検査でデータが複数回が調査された場合、エラー率を正しく維持するためには、同時分布を含む計算が必要となる。

== 他の統計との関係 ==

* ログランク統計量は、2つのグループを比較する[[Cox比例ハザードモデル]]の{{仮リンク|スコア検定|en|Score test}}として導出できる。したがって、その統計量は、そのモデルに基づく[[尤度比検定]]統計量と漸近的に等価である。
* ログランク統計量は、比例ハザード代替性<!-- proportional hazard alternative -->{{訳語疑問点|date=2021年10月}}を持つ任意の分布族の尤度比検定統計量と漸近的に等価である。たとえば、2つの標本からのデータが[[指数分布]]を持つ場合がある。
* <math> Z </math> をログランク統計量、<math> D </math> を観察された事象の数、<math>\hat {\lambda} </math> をハザード比の推定値とすると、<math> \log{\hat {\lambda}} \approx Z \, \sqrt{4/D} </math> である。この関係は、2つの量が分かっている場合（たとえば、発表された論文から）、3つ目の量が必要な場合に有用である。
* ログランク統計量は、観測が打ち切られている場合に使用できる。データに打ち切られた観測が存在しない場合、[[ウィルコクソンの順位和検定]]が適切である。
* ログランク統計量は、事象が発生した時間にかかわらず、すべての計算に同じ重みを与える。ピートログランク検定統計量<!-- Peto logrank test -->は、観測値の数が多い場合、初期の事象により多くの重みを与える。

== 検定の仮定 ==
ログランク検定は、[[カプラン＝マイヤー推定量|カプラン＝マイヤー生存曲線]]と同じ仮定に基づいている。すなわち、[[打ち切り (統計学)|打ち切り]]は予後とは無関係であり、生存確率は研究の初期と後期に募集された被験者で同じであり、事象は指定された時間に起こったという仮定である。これらの仮定からの逸脱が問題となるのは、比較されるグループ間で充足の度合いが異なる場合である。たとえば、あるグループでは打ち切りが他のグループよりも起こりやすいなどである<ref>{{Cite journal|last1=Bland|first1=J. M.|last2=Altman|first2=D. G.|year=2004|title=The logrank test|journal=BMJ|volume=328|issue=7447|pages=1073|doi=10.1136/bmj.328.7447.1073|pmid=15117797|pmc=403858|author-link1=Martin Bland|author-link2=Doug Altman}}</ref>。

== 参照項目 ==
{{Portal|数学}}

* [[生存分析]]
* [[カプラン＝マイヤー推定量]]
* [[ハザード比]]

== 脚注 ==
{{Reflist}}

{{Statistics|analysis}}

{{DEFAULTSORT:ろくらんくけんてい}}
[[Category:生存分析]]
[[Category:統計検定]]