Re: [Rd] abline for intercept-only simple lm models (with and without offset)

From: Martin Maechler <maechler_at_stat.math.ethz.ch>
Date: Mon 11 Dec 2006 - 13:17:48 GMT

Yes, I think all your propositions would be improvements. OTOH, I don't think the improvements warrant a big increase in code (complexity), nor do I think the improvements are crucial for R's integrity.

So...
If you (our someone else) provides a patch {against R-devel, as always} which keeps the code simple,
I'd strongly consider adding that to R.

Regards,
Martin Maechler

>>>>> "Tobias" == Tobias Verbeke <tobias.verbeke@telenet.be>
>>>>> on Sat, 09 Dec 2006 22:42:38 +0100 writes:

    Tobias> The abline function can be used to draw the
    Tobias> regression line when one passes the lm object as an
    Tobias> argument.

    Tobias> However, if it's an intercept-only model, it appears     Tobias> to use the intercept as the slope of the abline:

    Tobias> mod <- lm(dist ~ 1, data = cars) plot(dist ~ speed,     Tobias> data = cars) abline(reg = mod) # nothing appears

    Tobias> This behaves as documented, but might catch
    Tobias> someone. Would it be an improvement if this
    Tobias> situation was detected so as to plot the appropriate
    Tobias> horizontal line, i.e.

    Tobias> abline(a = coef(mod), b = 0) ?

    Tobias> Would it also be an improvement if for a model like

    Tobias> mod2 <- lm(dist ~ 1 + offset(speed), data = cars)

    Tobias> abline(reg = mod2) would be equivalent to

    Tobias> abline(a = coef(mod2), b = 1) ?

    Tobias> For models through the origin, the current function
    Tobias> works fine, but one might even consider models
    Tobias> through the origin and having the independent
    Tobias> variable in an offset() to be fully fool-proof, i.e.

    Tobias> abline(a = 0, b = 1)

    Tobias> Kind regards, Tobias

    Tobias> ______________________________________________
    Tobias> R-devel@r-project.org mailing list     Tobias>
https://stat.ethz.ch/mailman/listinfo/r-devel

R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Tue Dec 12 08:30:47 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Mon 18 Dec 2006 - 14:00:56 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.