[prev in list] [next in list] [prev in thread] [next in thread] 

List:       r-sig-mixed-models
Subject:    Re: [R-sig-ME] Fixing residual variance in lme4
From:       Ben Bolker <bbolker () gmail ! com>
Date:       2020-05-11 22:56:39
Message-ID: d81dedc7-e08e-c185-9242-361b4e26bbfe () gmail ! com
[Download RAW message or body]

     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:
> 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!
> 
> > example<- lmer(formula = g ~ 0 + (1|Study) + (1|Subscale),
> +                data = meta, weights=Precision,
> +                control=lmerControl(optimizer="bobyqa",sigma=1))
> 
> Error in lmerControl(optimizer = "bobyqa", sigma = 1) :
> unused argument (sigma = 1)
> 
> 
> 
> Best,
> Sijia
> 
> 	[[alternative HTML version deleted]]
> 
> _______________________________________________
> R-sig-mixed-models@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

_______________________________________________
R-sig-mixed-models@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models


[prev in list] [next in list] [prev in thread] [next in thread]