diff --git a/50-parcela-subdividida.Rmd b/50-parcela-subdividida.Rmd
index 563468a0c46a09ca4e9211f2810daec638b5c0ff..755753b3272bd7655f6fd503afb37f602154fe16 100644
--- a/50-parcela-subdividida.Rmd
+++ b/50-parcela-subdividida.Rmd
@@ -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)
+
+#-----------------------------------------------------------------------
+```