<div dir="ltr"><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Luiz,</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Você terá que mudar a função que está usando.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default"><span style="font-family:monospace">library(nlme)<br><br>table(x)<br><br>lattice::xyplot(y ~ x)<br><br>modelo0 <- aov(y ~ x)<br>a <- anova(modelo0)<br>a<br><br># Componentes de variância pelo método ANOVA.<br>c("sigma^2_b" = (a["x", "Mean Sq"] - a["Residuals", "Mean Sq"])/6,<br>  "sigma^2_e" = a["Residuals", "Mean Sq"])<br><br>modelo1 <- lme(y ~ 1,<br>               random = ~1 | x,<br>               method = "REML")<br>modelo1<br><br># Componentes de variância pelo método REML.<br>VarCorr(modelo1)<br><br>modelo2 <- lme(y ~ 1,<br>               random = ~1 | x,<br>               weights = varExp(),<br>               method = "ML")<br>modelo2<br><br># Componentes de variância pelo método REML.<br># ATTENTION: não são imediatamente comparáveis com os anteriores.<br>VarCorr(modelo2)<br></span></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif"><br></div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">À disposição.</div><div class="gmail_default" style="font-family:trebuchet ms,sans-serif">Walmes.<br></div></div>