Skip to content

Fixing residual variance in lme4

2 messages · Sijia Huang, Ben Bolker

#
Hi everyone,

I am trying to fit a crossed random effect model with lme4. I want the
residual variance fixed to be 1. Below is the code I use. I am wondering
why the sigma=1 argument does not work as it does in the nlme package. Or
if there is another to fix the residual variance. Thanks so much!
+                data = meta, weights=Precision,
+                control=lmerControl(optimizer="bobyqa",sigma=1))

Error in lmerControl(optimizer = "bobyqa", sigma = 1) :
  unused argument (sigma = 1)



Best,
Sijia
#
??? You need to use the blme package.? See 
https://bbolker.github.io/mixedmodels-misc/glmmFAQ.html#setting-residual-variances-to-a-fixed-value-zero-or-other
On 5/11/20 6:40 PM, Sijia Huang wrote: