Prezados colegas, preciso de ajuda.Eu sei como os pesos são calculados pra o caso 'varFixed' conforme exemplo abaixo, mas não sei como os pesos são calculados em attr(modelo$modelStruct$varStruct[], which="weights") para 'varExp', 'varPower' e 'varConstPower'. Alguém sabe como é calculado. Eu li a documentação do R, referências e procurei na internet, mas não achei nada.Desde já agradeço.Luiz# Data   x <- c(20,  30,  30,  30,  40,  40,  50,  50,  50,  60,  70,  70,  70,  80,  80,  80,  90,  90,  90,  90, 100, 100, 110, 110, 120)   y <- c(113, 121, 212, 273, 160, 244, 157, 221, 268, 224, 252, 323, 361, 342, 352, 399, 376, 377, 389, 468, 353, 420, 421, 435, 546)# Design matrix   X = model.matrix(~x)# Package   library(nlme)# varFixed   modelo = gls(y~x, weights = varFixed(~x))   h = sqrt(1/x)   w = 1/h^2   W = diag(w)   b = solve(t(X)%*%solve(W)%*%X)%*%t(X)%*%solve(W)%*%y   b ; coef(modelo) # 'b' computed by hand and 'coef(modelo)' match# varExp   modelo = gls(y~x, weights = varExp())   h = attr(modelo$modelStruct$varStruct[], which="weights") # how are these values calculated?   w = 1/h^2   W = diag(w)   b = solve(t(X)%*%solve(W)%*%X)%*%t(X)%*%solve(W)%*%y