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

No.21154　KMのY軸スケール　　【赤岳】　2014/07/07(Mon) 16:19

いつもご教示ありがとうございます。
今回，カプランマイヤー曲線をRを使って描きたいのですが，Rのデフォルトでは，Y軸の生存率が0～1の表示になっているかと思います。それを，％表示にするのはどのように指定すればよいかわかりません。
ご教示よろしくお願いします。

現在のプログラムは次のとおりです。
> library(survival)
> NP <-read.csv("NP.csv")
> km1<-survfit(Surv(TIME, EVENT) ~TREATMENT, data=NP)
> plot(km1, lty=1:3, main="KM", xlab="Time", ylab="Survival rate")

No.21155　Re: KMのY軸スケール　　【青木繁伸】　2014/07/07(Mon) 22:06

こんなことできないのかな？と思ったら，まずはオンラインヘルプを見るべき。
methods 関数により，plot は plot.survfit を呼ぶというのがわかるので，? plot.survfit で，引数の yscale であなたのやりたいことができることがわかる。
library(survival)
a <- survfit(Surv(time, status) ~ 1, data=kidney)
plot(a, xlab="Time", ylab="Survival rate", yscale=100)
今回はまあいいですが，質問の際に添付するプログラムは，コピー・ペーストだけで動くものにするのがよいですね。

No.21156　Re: KMのY軸スケール　　【赤岳】　2014/07/08(Tue) 05:43

青木先生，
数冊のR関連の本やWEBで探しましたが，見つからずこの質問をさせていただきました。
プログラムの件，大変失礼しました。以後気をつけます。
本当に，ありがとうございました。