[R-br] Dúvida sobre a função gls do pacote nlme

[Luiz Leal]

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

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20200625/d0463ac4/attachment.html>