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