0

我正在尝试使用 forestmodel 的 forest_model 函数绘制一个用生存的 coxph() 计算的 cox 比例风险比模型。

协变量是因子，它们的一些名称中有一个空格，即。“手术切除”。当我将这些输入 coxph 模型时，我用反引号引用它们，以便 R 知道忽略空格：

  res.coxm <- coxph(Surv(Survival_Overall, Death == 1) ~ Age+`Surgical Resection`
               +`Intracranial BoD`+`Systemic BoD`+`Treatment prior to CNS involvement`+
                 `Systemic treatment after CNS involvlement`+`Brain-directed radiotherapy`, data = df)

模型运行良好。当我尝试绘制此模型时，出现错误：

  forest_model(res.coxm, format_options = list(color="black", text_size=4.5))

  Error in grouped_df_impl(data, unname(vars), drop) : Column `Surgical Resection` is unknown

有没有人遇到过这个问题？

如果： - 我将变量名称更改为它们的“无空格”版本，即“surgical_resection”而不是“Surgical Resection”，或者 - 如果我将变量从因子更改为整数，但将名称保留为空格

我尝试使用各种不同的列命名函数（包括名称（sd））重新编码我的原始 df

有任何想法吗？谢谢！！

4

1 回答 1

0

这似乎是重命名列的阻力最小的路径，因此它们没有空格。如果成功，那么您还可以向维护人员发送功能请求，要求他/她使包符合 R 回归函数中的典型编程实践。如果它不成功，那么您仍然会发送内容有所不同的请求。

我也无法用一个简单的测试用例重现错误：

test1 <- list(time=c(4,3,1,1,2,2,3), 
               status=c(1,1,1,0,1,1,0), 
               x=c(0,2,1,1,1,0,0), 
               'Sex m f'=c(0,0,0,0,1,1,1)) 

 f <- coxph(Surv(time, status) ~ x + `Sex m f`, test1)
 forestmodel::forest_model(f)

  # no error, plot as expected.
于 2018-01-04T00:09:37.407 回答