在Byar Dataset上使用幸存函数

Question

我必须为Byar Dataset的2个死亡原因绘制分布函数。在clustMD的Byar数据集中，事件的范围从1到8，有几种死因。我使用以下代码：

library(clustMD)
data(Byar)
A<-Byar

library(survival)
attach(A)
time<-Observation
event<-SurvStat

B<-survfit(Surv(time,event)~1, data=A)
plot(B,xlab="Time",ylab="Survival Probability")

这给了我一条曲线，但它并不对应于两种死亡原因。将 ＆＃39; survfit（监测网（时间，事件==＆＃34; 1＆＃34; |事件==＆＃34; 2＆＃34; |事件==＆＃34; 3＆＃34）〜1，数据= A）＆＃39; Plot obtained给我2个死亡原因的情节，其中事件= 1是活着的，其他两个是死亡原因。

Answer 1

请注意，SurvStat中有超过2或3个级别。如果多个级别描述相同的死因，则需要重新编码。

pacman::p_load(clustMD)
data(Byar)
A<-Byar

library(survival)
attach(Byar)
time<-Observation
event<-SurvStat

B<-survfit(Surv(Observation,event=SurvStat, type="mstate")~1, data=Byar)
plot(B,xlab="Time",ylab="Survival Probability")