[R] lm methods inconsistency?
Martin Maechler
maechler at stat.math.ethz.ch
Mon Jan 10 17:39:33 CET 2000
>>>>> "Tom" == Tomislav Goles <tom at ait-tech.com> writes:
Tom> I'm very new to R so maybe i'm doing something wrong so
Tom> please let me know it that is that case.
Tom> Here is an example where the summary() and residuals()
Tom> applied to lm object produce different results
Tom> (I think the residuals() results is correct since SAS produces
Tom> those numbers - the second residual corresponding to observation
Tom> with weight 4 is wrong in summary()).
Tom> x <- c( 10, 20, 30, 40 )
Tom> y <- c( 33.8, 62.2, 92.0, 122.4 )
Tom> w <- c( 1, 4, 1, 1 )
Tom> lm.obj <- lm( y ~ x, weights=w )
Tom> print( summary( lm.obj ) )
(showing the weighted residuals)
Tom> print( residuals( lm.obj ) )
(showing the raw residuals)
No real bug here.
One could argue that summmary() should say
"Weighted residuals:" instead of
"Residuals:"
(but note that R's behavior here is really like S').
Tom> Best regards,
Tom> Tomislav Goles
Tom> PS. my version or R is 0.64.2 on Linux.
which is *really* too old.
but it doesn't matter here...
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list