Prezados, bom dia.
Alguém sabe como estimar "variance between groups" e "variance within groups" no modelo abaixo?

y <- c(1,2,3, 5,9,8, 2,3,4, 9,7,8 )
x <- as.factor(rep(letters[1:4], each=3))
library(nlme)
model <- gls(y~x, weights = varExp())

Desde já agradeço

Luiz