[R-br] variance components no modelo gls

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çoLuiz