R包survival拟合cox模型模型出现“Error in `[.default`(y, , 3) : 下标出界”的错误

今天利用R拟合Cox回归模型时候，意外发现模型出错？

下面是我的代码：

fit<-cph(Surv(生存时间,心血管) ~ rcs(bmi,3)+地区+sex+age,data=mydata)

因为我做的是限制性立方样条，所以有rcs。下面是我出错的情况。

经过分析后发现，原因是我的因变量 “心血管”为字符变量导致出错。

这是因为我在前面用如下代码，为 心血管（它的值是1和0）设置了参照组

##设定因子

mydata$心血管<-as.factor(mydata$心血管)

##设置分类变量的参照组

mydata$心血管<-relevel(mydata$心血管, ref="0")

当我删除上面的代码，并重新导入数据，心血管为数值变量时，模型运行成功了。

总结：经过对比，我发现cph拟合模型时，要求结局变量在数据库中必须是数值变量，否则会出错

希望我的经验可以帮助到大家。