使用MST包创建生存树:未定义的列错误?

时间:2019-06-30 13:23:33

标签: r survival-analysis

我正在尝试使用R中的MST包创建生存树。我一直在研究this论文。

我用随机生成的数据复制了他们的示例,并且效果很好。我调整了数据以适应相同的模型。我的数据具有相同的列和相同的数据类型。

我不断收到此错误:

Error in `[.data.frame`(mf_data[col.split.var], , 3) : undefined columns selected

具有以下代码行:

fit <- MST(formula = Surv(time,status)  ~ x1 + | id), data = data)

我已经浏览了所有文档,但没有找到任何东西,我也不明白为什么会出现此错误。 本文中的代码如下所示:

set.seed(186117)
data <- rmultime(N = 200, K = 4, beta = c(-1, 0.8, 0.8, 0, 0),cutoff = c(0.5, 0.3, 0, 0), model = "marginal.multivariate.exponential", rho = 0.65)$dat
test <- rmultime(N = 100, K = 4, beta = c(-1, 0.8, 0.8, 0, 0), cutoff = c(0.5, 0.3, 0, 0), model = "marginal.multivariate.exponential",rho = 0.65)$dat
fit <- MST(formula = Surv(time, status) ~ x1 + x2 + x3 + x4 | id,data, test, method = "marginal", minsplit = 100, minevents = 20,selection.method = "test.sample")

1 个答案:

答案 0 :(得分:0)

我尝试运行您的代码,但确实出现了一个错误,尽管不是您所得到的,并且在查看它后我很确定您需要使用SO的[edit]功能来修改您的问题。

> fit <- MST(formula = Surv(time,status)  ~ x1 + | id), data = data)
Error: unexpected '|' in "fit <- MST(formula = Surv(time,status)  ~ x1 + |"

给出的公式显然是错误的,并且有不必要的右括号。我可以通过以下方式报告您的错误:

> fit <- MST(formula = Surv(time,status)  ~ x1 | id, data = data)
[1] "No test sample supplied, changed selection.method = 'bootstrap'"
Error in `[.data.frame`(mf_data[col.split.var], , 3) : 
  undefined columns selected

....,但没有原始代码:

  

fit <-MST(公式= Surv(时间,状态)〜x1 + x2 + x3 + x4 | id,data,test,method =“ marginal”,minsplit = 100,minevents = 20,selection.method =“ test.sample”)

我还看到公式的RHS上带有x1+x2|id的错误,但没有三个变量:

> fit <- MST(formula = Surv(time, status) ~ x1 +x2 | id,data, test, method = "marginal", minsplit = 100, minevents = 20,selection.method = "test.sample")
Error in `[.data.frame`(mf_data[col.split.var], , 3) : 
  undefined columns selected
> fit <- MST(formula = Surv(time, status) ~ x1 +x2+x3| id,data, test, method = "marginal", minsplit = 100, minevents = 20,selection.method = "test.sample")

因此,我认为这是开发人员未曾预料到的错误。以下是获取报告所需电子邮件地址的方法:

> maintainer("MST")
[1] "Peter Calhoun <calhoun.peter@gmail.com>"