问题描述
因此,我一直在尝试在林区图中可视化我的Cox模型结果,但始终收到错误Error in
[。data.frame (data,var) : undefined columns selected
。这在我看来有些奇怪。定义模型并提供摘要,结果如常打印,似乎没有问题。此外,模型不包含strata()
和tt()
对象。
经过一番谷歌搜索,这似乎是不时发生的错误,但是我找不到任何令人满意的有效解决方案来管理它。
鉴于所有其他功能都能正常工作,我很确定这实际上与数据本身无关。
或者,对于任何建议如何对结果进行信息性概述的建议,我也感到非常高兴,尤其是当我的其他一些模型包含层次结构时,ggforest似乎不支持它。
编辑:我发现了问题:在coxph函数中指定ID
参数时,我得到未定义的列选择错误。
duration_analysis_10 <- coxph(formula = Surv(start1,stop1,event) ~ Var1 +
Var2 +
log(Var3)+
Var4 + Var5 +
Var6 +
Var7 + Var8 +
Var9 + Var10 +
var11 + Var12 +
Var13,ties = "efron",na.action = na.exclude,id = CountryName,data = df)
来自eha软件包的可复制样本:
oldmort <- arrange(oldmort,id)
oldmort$event <- ifelse(oldmort$event == "TRUE",1,0)
oldmort$sex <- ifelse(oldmort$sex == "female",0)
oldmort$id <- as.character(oldmort$id)
oldmorttest <- coxph(formula = Surv(enter,exit,event) ~ imr.birth + sex,id = id,data = oldmort)
summary(oldmorttest)
ggforest(oldmorttest,data = oldmort)
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)