ログランク

• こちらはリンク切れ
  • 目次のページ
    • 前篇(カプラン・マイヤー法とログランク検定，人年法について)
    • 中篇(Cox 回帰分析について)
    • 後篇(競合リスクや再発事象に関する解析方法について)

library(survival)
# survdiff(Surv(時間変数, イベント変数) ~ 群変数, データセット名)
> survdiff(Surv(time, censor) ~ group, Event) # ログランク検定
N Observed Expected (O-E)^2/E (O-E)^2/V
group=A 9 3 3.5 0.0714 0.146
group=B 9 4 3.5 0.0714 0.146
Chisq= 0.1 on 1 degrees of freedom, p= 0.702 ← 差なし！
（総イベント数が少なすぎ…）
> (1-9/18)+(0-8/16)+(1-2*7/14)+(1-5/10)+(0-4/8)+(0-2/4) # スコアU1
[1] -0.5
> (9*9)/18^2+(8*8)/16^2+(7*7*2*(14-2))/(13*14^2)+(5*5)/10^2+
+ (4*4)/8^2+(2*2)/4^2 # 分散V
[1] 1.711538
> 0.5^2/1.711538 # χ2 = U1
2/V
[1] 0.1460675
> pchisq(0.146, df=1) # p 値

• こちらも
• こちらも
• 検定に当たってログランクならsurvdiff()関数、一般化Wilcoxsonなら(こちら)