我有以下内容:
splits<-split(df,disease)
我用column1分割了一个数据帧。但我想使用Surv功能:
library(survival)
with(df,Surv(TIME,STATUS)~GENDER)
如何将其扩展到我的分割数据集,以便在每个分割数据集上单独运行Surv()函数?我尝试过使用sapply()函数:
surv.obj<-sapply(splits, function(x) Surv(TIME,STATUS)~GENDER)
这似乎工作正常,因为我得到:
surv.obj
$heart_disease
Surv(TIME,STATUS)~GENDER
<environment: 0X0000000000000019cac10>
$cancer
Surv(TIME,STATUS)~GENDER
<environment: 0X0000000000000019cac11>
.
.
.
,但是在尝试使用幸存函数时:survfit(surv.obj)
,我尝试绘制surv.obj $ heart_disease我得到了
plot(suv.obj$heart_disease)
Error in Surv(TIME,STATUS) : object TIME not found
这显然意味着公开声明没有完成它的工作。有人可以告诉我正确的sapply语法应该是什么?