今天在使用R绘制风险三联图时,我跟着教程的命令输进去,R确不断报错,教程的命令如下:
fit <- cph(Surv(time,status)~ANLN+CENPA+GPR182+BCO2,LIRI)
#使用数据中的ANLN、CENPA、GPR182和BCO2四个基因的表达值构建COX比例风险模型
我跟着输进去,结果出错了, Error in eval(predvars, data, callenv):没有发现XXX
我寻思着难道是我变量有问题?我索性删除没有找到的变量,继续执行这个命令,结果系统还是原提示,XX改为删除变量后排第一的YY了。那就是命令有问题
我改了,突然想到,人家的数据集名称就是 LIRI
我的数据集名称是? Mydata
于是更改后的命令如下:
fit <- cph(Surv(time,status)~XX+YY+CC+ZZ,data=Mydata)
#要限定下自己的数据集,不然识别不了,数据集也不要填错了
然后就没然后了
回归模型构建成功
Ps:在R语言中,cph函数是一个常用的用于实现Cox比例风险模型的函数