[R-sig-ME] lmer formula

Ben Bolker bbo|ker @end|ng |rom gm@||@com
Sun Apr 10 01:07:26 CEST 2022

   This is not possible in lme4, which requires a *linear* model (some 
incorporation of constants is possible via offsets, but only additive 

   Using nlme() from the nlme package, something like

nlme(Reaction ~ A1*mu1/(B2*mu2) + mu3,
      fixed = list(mu1 ~ Days,
                   mu2 ~ C1),
      random = ~ 1 | mu3

might work.  But keep in mind that when you fit nonlinear mixed effects 
models you often have to be much clearer by what you mean by "estimating 
fixed effects for" something -- this will typically involve a sub-model 
that is a linear model for a parameter that enters the top-level 
nonlinear equation.

On 4/9/22 3:11 PM, Andras Farkas via R-sig-mixed-models wrote:
> Hello,
> could you please provide input on the following formula (please ignore the widely known example of the sleepstudy data and its analysis for a moment and assume we have this formula here) :
>   m1 <- lmer(Reaction ~ ((A1*Days)/(B2*C1)) + (1 | Subject), sleepstudy)
> this simple example - similar to my actual real life model of greeter complexity - what I would try to do is to estimate fixed effects for "Days" and "C1", but not for "A1"
> and "B2". In this hypothetical example A1, B2, AND C1 all have a column with values in the data sleepstudy, but instead of estimating a coefficient for A1 and B2 we would just use the actual values of them as they appear in the spreadsheet in the calculations.
> hope I am not way off with my question, but could you please advise as to how the model formula would need to be written in this case?
> much appreciate the help,
> thanks,
> Andras
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

More information about the R-sig-mixed-models mailing list