Skip to content
Snippets Groups Projects
Commit 1c2cc7b8 authored by Walmes Marques Zeviani's avatar Walmes Marques Zeviani
Browse files

Adds comparison between aov and lmer classes.

parent dd0df83f
Branches
No related tags found
No related merge requests found
......@@ -579,3 +579,97 @@ do.call(rbind, lht)
# Pooled Error 30.74135 687742.4 22371.90
#--------------------------------------------
```
# AOV vs LMER
```{r}
#-----------------------------------------------------------------------
# {emmeans} aceita `aov` e `lmer`.
data(warpbreaks)
str(warpbreaks) # Analisar como parcela subdividida para demonstração.
library(lme4)
library(lmerTest)
library(emmeans)
data(Soybean, package = "nlme")
str(Soybean)
# Um recorte dos dados com todos fatores categóricos.
my_data <- subset(Soybean, Year == 1990 & Time > 50)
my_data <- transform(my_data, Time = factor(Time))
lattice::xyplot(weight ~ Time | Variety,
groups = Plot,
data = my_data,
type = c("p", "a"))
# Modelo de efeitos aleatórios por REML.
m0 <- lmer(weight ~ Variety * Time + (1 | Plot),
data = my_data)
anova(m0)
# Mesmo modelo mas por mínimos quadrados.
m1 <- aov(weight ~ Variety * Time + Error(Plot),
data = my_data)
summary(m1)
#--------------------------------------------
# Efeito principal de variedade.
emm0 <- emmeans(m0, spec = ~Variety)
emm1 <- emmeans(m1, spec = ~Variety)
emm0
emm1
contrast(emm0, method = "pairwise")
contrast(emm1, method = "pairwise")
multcomp::cld(emm0, Letters = letters)
multcomp::cld(emm1, Letters = letters)
#--------------------------------------------
# Efeito principal de tempo.
emm0 <- emmeans(m0, spec = ~Time)
emm1 <- emmeans(m1, spec = ~Time)
emm0
emm1
contrast(emm0, method = "pairwise")
contrast(emm1, method = "pairwise")
multcomp::cld(emm0, Letters = letters)
multcomp::cld(emm1, Letters = letters)
#--------------------------------------------
# Desdobramento de tempo dentro de variedade.
emm0 <- emmeans(m0, spec = ~Time | Variety)
emm1 <- emmeans(m1, spec = ~Time | Variety)
emm0
emm1
contrast(emm0, method = "pairwise")
contrast(emm1, method = "pairwise")
multcomp::cld(emm0, Letters = letters)
multcomp::cld(emm1, Letters = letters)
#--------------------------------------------
# Desdobramento de variedade dentro de tempo.
emm0 <- emmeans(m0, spec = ~Variety | Time)
emm1 <- emmeans(m1, spec = ~Variety | Time)
emm0
emm1
contrast(emm0, method = "pairwise")
contrast(emm1, method = "pairwise")
multcomp::cld(emm0, Letters = letters)
multcomp::cld(emm1, Letters = letters)
#-----------------------------------------------------------------------
```
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment