如果我将一个模型拟合函数(例如 lme4 包中的 glmer )包装在另一个函数中并允许一个 ... 参数,则 getCall 不会识别它。这是一个错误吗?
library(lme4)
x = rnorm(100)
beta0 = 1
beta1 = 2
fact = factor(rep(1:4, each = 25))
y = beta0 + beta1 * x + rep(rnorm(4), each = 25) + rnorm(100)
dat <- data.frame("y" = y, "x" = x, "fact" = fact)
f <- function(data, formula, ...) {
lmer(formula, data, ...)
}
mod1 <- f(data = dat, formula = formula(y ~ x + (1|fact)),
control = lmerControl(optimizer="Nelder_Mead"))
mod2 <- lmer(data = dat, formula = y ~ x + (1|fact),
control = lmerControl(optimizer="Nelder_Mead"))
getCall(mod1)
getCall(mod2)