Hello!

I have been thinking about an analysis for so much time, and I still have so many doubts, that I definitely need some help, confronting or validating.

I have four blocks. Within each block, we established five plots, each with a fertilization treatment. Each plot is divided in two subplots, one grazed and the other protected. My aim is to estimate the more accurately as possible the interactive effects of fertilization and grazing (and their uncertainty) using marginal means as implemented in the `emmeans`

package.

Now, we took measurement repeatidly (3-4 times) in each subplot during three years. Two blocks were generally measured the same day. After reflexion and reading, I produced the following formula (in `brms`

).

```
form <- bf(y ~ Fertilization + Grazing + Fertilization:Grazing +
Block + (1|Block:Fertilization:Grazing) +
Year + (1|Parcelle:Date)
)
```

Does that make sense?

Thank you very much!

Lucas