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
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
