0

如果我将一个模型拟合函数(例如 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)
4

0 回答 0