我正在尝试使用lme4软件包中的lmer函数运行混合模型,并希望在我的参数上包括变化的边界(例如,参数a和b的固定效果估计值为0.4,而我希望限制为0.3)。 / p>
model <- lme4::lmer(y ~ a + b + c + d + e +
(a + b + c + d | t),
data = data,
control = lmerControl(optimizer = "nloptwrap",
optCtrl=list(method="nlminb",starttests = FALSE)),
upper = c(0.3, 0.3, 500, 100, 100))
当我包含“ upper”语句时,出现以下错误:
Error in lme4::lmer(y ~ a + b + c + :
unused argument (upper = c(0.3, 0.4, 500, 100, 100))
是否可以将边界应用于变量? 我一直在试验lmer vs. glmer vs. nlme vs. glmmPQL,并认为lmer似乎是最好的选择,但是希望得到任何反馈!
谢谢!