Displaying 1 result from an estimated 1 matches for "newqr".
Did you mean:
newer
2003 Jun 04
1
standard errors from glm (PR#3180)
...ut don't we want
fitted values from parameter estimates at the final iteration?
I tried inserting the following lines upon completion of the iterative
weighted least squares loop:
mu.eta.val <- mu.eta(eta)
w <- sqrt((weights[good] * mu.eta.val[good]^2)/variance(mu)[good])
newQR<-qr(x[good, ] * w)
My standard errors based on newQR$qr agree with the standard
errors I had computed manually from the final fits.
Has anyone else encountered this issue? Thanks in advance
for any clarifications.
Brad
--please do not edit the information below--
Version:
platform = i686-p...