[R-br] componente da variância no modelo de efeitos aleatórios

[Luiz Leal]

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 <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))summary(aov(y~x))


Mas não sei como calcular quando utilizo mínimos quadrados generalizados.
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 <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))library(nlme)
modelo <- gls(y~x, weights = varExp()           , method = "ML")
Alguém sabe como calculo as componentes da variância neste caso?
Desde já agradeço
Luiz

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20191127/9bf1298e/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: 1574880488033blob.jpg
Tipo: image/png
Tamanho: 8829 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20191127/9bf1298e/attachment.png>