<html><head></head><body><div class="yahoo-style-wrap" style="font-family:Helvetica Neue, Helvetica, Arial, sans-serif;font-size:13px;"><div><div><br></div><div>Prezados colegas, preciso de ajuda.</div><div>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'.<br>Alguém sabe como é calculado. Eu li a documentação do R, referências e procurei na internet, mas não achei nada.</div><div>Desde já agradeço.</div><div>Luiz</div><div># Data<br> 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)<br> 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)</div><div># Design matrix<br> X = model.matrix(~x)</div><div># Package<br> library(nlme)</div><div># varFixed<br> modelo = gls(y~x, weights = varFixed(~x))<br> h = sqrt(1/x)<br> w = 1/h^2<br> W = diag(w)<br> b = solve(t(X)%*%solve(W)%*%X)%*%t(X)%*%solve(W)%*%y<br> b ; coef(modelo) # 'b' computed by hand and 'coef(modelo)' match</div><div># varExp<br> modelo = gls(y~x, weights = varExp())<br> h = attr(modelo$modelStruct$varStruct[], which="weights") # how are these values calculated?<br> w = 1/h^2<br> W = diag(w)<br> b = solve(t(X)%*%solve(W)%*%X)%*%t(X)%*%solve(W)%*%y<br></div></div><b></b><i></i><u></u><sub></sub><sup></sup><strike></strike><div><br></div></div></body></html>