使用R进行生存分析

Question

在使用生存包运行生存模型时，我遇到了一条奇怪的错误消息，我无法理解。

错误消息是

Error in Surv(time, event) : Time and status are different lengths

model1中使用的每列中的观察数量完全相同。

> mydata3<-read.csv(file.choose())
> attach(mydata3)
> event<-Event
> time<-time
> x<-cbind(WDCL,Age,SEX)
> ch<-coxph(Surv(time,event)~x,method="breslow")

谢谢和问候。

Answer 1

您需要通过添加所有参数来正确指定公式。

试试这个：

coxph(Surv(time, Event) ~ WDCL + Age + SEX, method="breslow", data = mydata3)