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
_______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
