使用forestmodel软件包的森林图对于多变量cox回归非常有用。但是,我在将原始变量的名称和数据框的因子替换为最终展示标签时遇到问题(即variable:age2代表“ Age”; factors代表“ <60”; 1代表“≥60”)。有限的R编码知识,但是我尝试使用expss包为变量和因子添加标签。但是,coxph()不适用于标签,而适用于因子。
这是我的编码:
Cox比例模型:
mcox<-coxph(pblsurv~age2+sex1+origin,data = pbl)
使用forestmodel包的森林图:
print(forest_model(mcox))
最终情节中的变量名称为age2,sex1,origen;因此,我使用expss包添加标签:
pbl <- apply_labels(pbl,
age2 = "Age",age2 = c("<60"=0,"≥60"=1),
sex1 = "Gender",sex1 = c("Female"=0,"Male"=1),
origin = "Ethnicity",origin =c("Non=hispanic"=0, "Hispanic"=1))
但是,在应用标签后,coxph无效:
mcox<-coxph(pblsurv~age2+sex1+origin,data = pbl)
Error in coxph(pblsurv ~ age2 + sex1 + origin, data = pbl) :
data contains an infinite predictor
您知道要在print(forest_model(mcox))
中使用什么其他代码来进行最终期刊介绍吗?