今天利用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拟合模型时,要求结局变量在数据库中必须是数值变量,否则会出错
希望我的经验可以帮助到大家。