Luiz,
Você terá que mudar a função que está usando.
library(nlme)
table(x)
lattice::xyplot(y ~ x)
modelo0 <- aov(y ~ x)
a <- anova(modelo0)
a
# Componentes de variância pelo método ANOVA.
c("sigma^2_b" = (a["x", "Mean Sq"] - a["Residuals", "Mean Sq"])/6,
"sigma^2_e" = a["Residuals", "Mean
Prezados, boa tarde.
No modelo de efeitos aleatórios estimamos a componente da variância conforme
abaixo.
y <-c(2370, 1687, 2592, 2283, 2910, 3020, 1282, 1527, 871, 1025, 825, 920,
562, 321, 636, 317, 485, 842, 173, 127, 132, 150, 129, 227, 193,
71, 82, 62, 96, 44)x