統計のコツのこつ(24)

 

前回に引き続き生存分析のお話です。

それでは、「すぐに役立つ統計のコツ」第８章(134ページ)を開いて下さい。

「K-M法による Logrank test」を「R」で行ってみましょう。

前回のデータであれば、次のコマンド（関数）だけでＯＫです。

 

***
survdiff(Surv(観察期間, イベント発生 == 1)~ グループ, data=dat)

出力結果：第８章(136ページ参照)

> survdiff(Surv(観察期間, イベント発生 == 1)~ グループ, data=dat)
Call:
survdiff(formula = Surv(観察期間, イベント発生 == 1) ~ グループ,
    data = dat)

            N Observed Expected (O-E)^2/E (O-E)^2/V
グループ=A 10        3     6.13      1.60      4.03
グループ=B 11        8     4.87      2.02      4.03

 Chisq= 4  on 1 degrees of freedom, p= 0.0448
***

 

この結果は、「すぐに役立つ統計のコツ」第８章(137ページ)を見て下さい。

 

次回は、「コックス比例ハザード・モデル」をご紹介します。

 

情報統計研究所はここから！