R 我看不出。。。当模型函数包装在另一个函数中时

R 我看不出。。。当模型函数包装在另一个函数中时,r,R,如果我在另一个函数中包装一个模型拟合函数,例如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, "fa

如果我在另一个函数中包装一个模型拟合函数,例如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)
库(lme4)
x=rnorm(100)
beta0=1
β1=2
事实=系数(代表(1:4,各=25))
y=beta0+beta1*x+rep(形式(4),每个=25)+rnorm(100)

dat你是说这个代码会为你返回一个错误?这对我来说似乎很管用。这里到底有什么问题?使用
lme4_1.1-21进行测试
否,但控制参数未显示,出于某种原因,这很重要?这个问题只是为了获得更漂亮的输出?如果您试图构建完全相同的函数调用,那么如果您将f的输出传递给另一个使用getCall的函数,则需要执行类似于
f的操作。这可能会导致错误