Ich versuche, eine Formel für ein lineares Mischeffektmodell (mit lme4
) für mein experimentelles Design anzugeben , bin mir aber nicht sicher, ob ich es richtig mache.
Das Design: Im Grunde messe ich einen Antwortparameter an Pflanzen. Ich habe 4 Behandlungsstufen und 2 Bewässerungsstufen. Die Pflanzen sind in 16 Parzellen gruppiert, wobei in jeder Parzelle I 4 Unterparzellen beprobt werden. In jeder Teilfläche nehme ich zwischen 15 und 30 Beobachtungen (abhängig von der Anzahl der gefundenen Pflanzen). Das heißt, es gibt insgesamt 1500 Zeilen.
Ursprünglich war die Subplot-Ebene nur zu Stichprobenzwecken hier, aber ich dachte, ich würde sie gerne im Modell berücksichtigen (als Variable mit 64 Ebenen), da ich sah, dass es von einer Subplot-zur anderen eine große Variabilität gibt , sogar innerhalb derselben Parzelle (größer als die Variabilität zwischen ganzen Parzellen).
Meine erste Idee war zu schreiben:
library(lme4)
fit <- lmer(y ~ treatment*irrigation + (1|subplot/plot), data=mydata)
oder
fit <- lmer(y ~ treatment*irrigation + (1|subplot) + (1|plot), data=mydata)
Ist das korrekt? Ich bin mir nicht sicher, ob ich beide Plot- / Subplot-Ebenen in meiner Formel behalten muss. Kein fester Effekt ist signifikant, aber die zufälligen Effekte sind sehr signifikant.
quelle