我从 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
我的真实数据也会发生这种情况。
谁能看到我做错了什么？