Hierarchical structure split-plot design with repeated measures

General
2

Ok, writing the above post made me change my keywords, and it seems the correct structure should be

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

Following Gumperts and Brown 1993.

EDIT : Finally, after re-reading the above paper and thinking about some particularity of my data (blocks note measured the same day or with the same interval), I am now convinced this is good solution :

bf(Albedo ~ 
                    ### Between plot analysis
                    Parcelle + 
                    Fertilization + (1 | Parcelle:Fertilization) +
                    Grazing + Fertilization:Grazing + (1 | Parcelle:Fertilization:Grazing) +
                    ### Within plot analysis
                    #### Between year analysis
                    Year + Year:Parcelle + 
                    Year:Fertilization + (1 | Year:Parcelle:Fertilization) + 
                    Year:Grazing + Year:Fertilization:Grazing + (1|Year:Parcelle :Fertilization:Grazing) +
                    #### Within year analysis
                    Year:scale_DOY + Year:scale_DOY:Parcelle + 
                    Year:scale_DOY:Fertilization + 
                    Year:scale_DOY:Grazing + Year:scale_DOY:Fertilization:Grazing +
                    (0 + scale_DOY | Year:Parcelle:Fertilization) + (0 + scale_DOY | Year:Parcelle:Fertilization:Grazing)
)
1 Like