R：来自Zelig'的错误；s coxph

我从泽利格提供的考克斯生存模型中得到了一个错误

下面是一个使用虚拟数据的示例：

require(survival)
require(Zelig)
dat <- data.frame(survival = rnorm(1000, 100, 20), event = c(rep(1, 200), rep(0, 800)), ind = sample(rep(c("A","B","C"), 1000 * c(0.3, 0.3, 0.4)), 1000))
zelig(Surv(survival, event) ~ ind, model = "coxph", data = dat)

我的真实数据也是如此

---

有人知道我做错了什么吗？

这是个棘手的问题。显然，这甚至让哈佛感到困惑：他们为此开了一张GitHub票，但那是4年前的事了，现在票不见了。虽然在这种情况下，他们使用了

model=“relogit”

，他们猜测这是relogit中的一个bug……谢谢。。。我想我可能不得不使用一种不同的方法！Zelig很有吸引力，因为它需要多个插补数据集，但有更多的手动选项……是的，我发现如果我更改模型类型，我不会出现错误，但我可以使用相同的公式运行

survival:：coxph（）

，运行良好。所以这一定是软件包中的一些深层次的技术缺陷，但是如果你打开GitHub票据，他们可能会很快做出响应。好的，太好了。我开了一张票。我也会在这里发布任何更新。得到了一个答案--“很遗憾，Zelig 5不支持coxph。我们将致力于重新引入支持。这是一个棘手的问题。显然，这甚至让哈佛感到困惑：他们为此开了一张GitHub票，但那是4年前的事了，现在票不见了。虽然在这种情况下，他们使用了

model=“relogit”

，他们猜测这是relogit中的一个bug……谢谢。。。我想我可能不得不使用一种不同的方法！Zelig很有吸引力，因为它需要多个插补数据集，但有更多的手动选项……是的，我发现如果我更改模型类型，我不会出现错误，但我可以使用相同的公式运行

survival:：coxph（）

，运行良好。所以这一定是软件包中的一些深层次的技术缺陷，但是如果你打开GitHub票据，他们可能会很快做出响应。好的，太好了。我开了一张票。我也会在这里发布任何更新。得到了一个答案--“很遗憾，Zelig 5不支持coxph。我们将致力于重新引入“支持”。

Error in eval(expr, envir, enclos) : attempt to apply non-function