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

[Thierry Onkelinx]

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
INSTITUUT VOOR NATUUR- EN BOSONDERZOEK / RESEARCH INSTITUTE FOR NATURE AND
FOREST
Team Biometrie & Kwaliteitszorg / Team Biometrics & Quality Assurance
thierry.onkelinx using inbo.be
Havenlaan 88 bus 73, 1000 Brussel
www.inbo.be

///////////////////////////////////////////////////////////////////////////////////////////
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
///////////////////////////////////////////////////////////////////////////////////////////

<https://www.inbo.be>


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