[R-sig-ME] Mixed effect model with lme4

Thierry Onkelinx th|erry@onke||nx @end|ng |rom |nbo@be
Wed Jan 19 08:50:31 CET 2022

Dear Sijia,

I think you want (1 + d | delta). Keep in mind that this will fit a
different alpha_i and beta_i for every level of the random effect.

Best regards,

ir. Thierry Onkelinx
Statisticus / Statistician

Vlaamse Overheid / Government of Flanders
Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance
thierry.onkelinx using inbo.be
Havenlaan 88 bus 73, 1000 Brussel

To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to say
what the experiment died of. ~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data. ~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of data.
~ John Tukey


Op di 18 jan. 2022 om 17:10 schreef Sijia Huang <huangsjcc using gmail.com>:

> Hi everyone,
> I am trying to fit a mixed effect model using lme4. The random effect part
> of the model is kind of tricky:
> delta_{1s}*(alpha_i + beta_i * d_{t=0,i=1})
> In which *delta_1s* is the random effect and are multiplied by the
> quantities in the parentheses. The *alpha_i* and *beta_i* are parameters
> that need to be estimated, while the d_t=0,i=1 is a dummy variable.
> I wonder if this model can be specified with lme4 syntax. Thank you!
> Best,
> Sijia
>         [[alternative HTML version deleted]]
> _______________________________________________
> R-sig-mixed-models using r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models

	[[alternative HTML version deleted]]

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