WFU

2012年5月27日 星期日

[R] How to plot a KM survival curve without Rcmdr


1. 先設定好變數兩組,當然,也可以直接帶入 database 而略過這步。

>survival=rnorm(20,100,60)
>status=c(1,0,1,1,1,0,0,0,0,0,1,1,1,0,0,1,1,0,0,1)

把資料合在一起,看看長度是否正確。
>data=cbind(survival, status)

2.記得安裝package:"survival"

3.
先來試試這個 formula(公式正不正確),此步可省略。
> Surv(survival,status)

重來了,要開始用 survfit 這個 function