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

List:       r-sig-mixed-models
Subject:    Re: [R-sig-ME] lmer formula
From:       Ben Bolker <bbolker () gmail ! com>
Date:       2022-04-09 23:07:26
Message-ID: c2272f41-e82d-83ee-552d-6dddfd2ff4dd () gmail ! com
[Download RAW message or body]


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

   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@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] 

Configure | About | News | Add a list | Sponsored by KoreLogic