[Rd] Fitted values with small weights in lm.wfit (PR#11979)

ablocker at gmail.com ablocker at gmail.com
Thu Aug 7 11:28:35 CEST 2008

Full_Name: Alexander Blocker
Version: 2.7.1
OS: Ubuntu 8.04 / Windows XP
Submission from: (NULL) (

When running lm(modeleq, weights=wt, data=dataset) with small weights (<1e-10),
I have encountered an odd phenomenon with fitted values. Due to numerical
precision issues, the fitted values and residuals returned by lm.wfit (from its
.Fortran call to dqrls) can differ greatly from those returned by running
predict on the resulting lm object. This is completely attributable to the
numerical precision passed to the given function, but I wonder if a warning
message for weights below as certain threshold may be in order.

More information about the R-devel mailing list