R：来自Zelig的coxph的错误

Question

我从Zelig提供的Cox生存模型中得到错误。

以下是使用虚拟数据的示例：

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)

这给出了错误：

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

这也适用于我的真实数据。

任何人都可以看到我做错了吗？