目前有不少粉丝私信我说使用ggrcs包出现如下错误:
Error in reformulate(attr(termobj, “term.labels”)[-dropx], response = if (keep.response) termobj[[2L]], :
'termlabels’必需是长度至少为一的字节矢量
我查看了一下,目前报错来源于新版本后的RMS包,主要是预测函数的报错,这个只能等R包作者来修复这个错误。目前需要急用的话,我提供了一个方案,请看下面视频操作
关于目前ggrcs包的报错解决方案
代码
library(ggrcs)
library(rms)
library(ggplot2)
library(scales)
library(cowplot)dt<-smoke
dd<-datadist(dt)
options(datadist='dd')
fit<- cph(Surv(time,status==1) ~ rcs(age,4)+gender, x=TRUE, y=TRUE,data=dt)ggrcs(data=dt,fit=fit,x="age"