Hoofdstuk 6 Les 6: Two-way ANOVA

Lees Chapter 18 (Multiple explanatory variables)

Stel, je wilt het effect van drie soorten bemesting op plantengroei uittesten op een groot proefveld. Het proefveld heeft een duidelijk gradient van droog naar vochtig. Tijdens de lessen over experimental design heb je gehoord over het randomized block design (gewarde blokproef). Deze opzet past precies bij dit onderzoek. Je verdeelt dan het proefveld in blokken, dwars op de richting van de vochtgradiënt. Binnen de blokken verdeel je dan random de plots voor iedere soort bemesting (NB: volgend blok gaan we het in meer detail hebben over verschil tussen fixed en random effects). In je statistisch model moet je dan het blokeffect meenemen. Dat kan met een two-way ANOVA.

Je gebruikt dezelfde functie lm() als voor de regressie en de one-way ANOVA:

fit <- lm(groei ~ bemesting + blok)

Daarna kan je met een posthoctoets checken welke bemesting verschilt van elkaar. De posthoctoets houdt rekening met de verschillen tussen de blokken.

emmeans(fit, specs = pairwise ~ bemesting)