「統計学関連なんでもあり」の過去ログ--- 045

No.16778　km.surv　　【赤岳】　2012/04/17(Tue) 19:25

RでKaplan-Meyer 法による生存曲線を描きたくて，本サイトの統計学自習ノートから
source("http://aoki2.si.gunma-u.ac.jp/R/src/km_surv.R", encoding="euc-jp")
を指定して描こうとしました。
関数定義をみると，km.surv（time，event）とあり，groupを指定できません。
Groupを指定する方法はどのようにすればよろしいでしょうか。
一つの図の中に，Group A, Group B, Group Cのそれぞれの曲線が入るようにしたいと思っています。
ご指導よろしくお願いします。

No.16779　Re: km.surv　　【青木繁伸】　2012/04/17(Tue) 21:21

種々の生存曲線を描くには，survival パッケージにある，Surv と survfit を組み合わせて使うという王道を行く方がよいでしょう。
library(survival) # splines パッケージもインストールしておくこと
? survfit.formula
example(survfit.formula)
をお試しください。

No.16780　Re: km.surv　　【赤岳】　2012/04/18(Wed) 10:12

青木先生，承知しました。勉強してみます。
ありがとうございました。