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

Tobias Verbeke tobias.verbeke at telenet.be
Sat Dec 9 22:42:38 CET 2006


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

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

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

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

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

Would it also be an improvement if for
a model like

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

abline(reg = mod2) would be equivalent to

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

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

abline(a = 0, b = 1)

Kind regards,
Tobias



More information about the R-devel mailing list