生存時間分析 by R

忘れないためにとりあえずの方法を。

R studioを立ち上げて
>library(Rcmdr)
次に生存時間解析のパッケージsurvialを立ち上げる
>library(suvival)

Rcommanderにエクセルで作成した表をクリップボードからコピーする。
とりあえずDataset（ファイル名）として保存する。

> survdiff(Surv(time, status)~treat, data=Dataset)
Call:
survdiff(formula = Surv(time, status) ~ treat, data = Dataset)

N Observed Expected (O-E)^2/E (O-E)^2/V
treat=BS 12 10 9.66 0.0118 0.0375
treat=EK 12 10 10.34 0.0110 0.0375

Chisq= 0 on 1 degrees of freedom, p= 0.847

が出てくる。
つまりこの生存時間曲線は違いがない。