尝试在R中使用混合模型(LMER)约束参数固定效应估计

时间:2020-09-24 20:55:40

标签: r lme4 effects mixed

我正在尝试使用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似乎是最好的选择,但是希望得到任何反馈!

谢谢!

0 个答案:

没有答案
相关问题