我从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
这也适用于我的真实数据。
任何人都可以看到我做错了吗？