具有线性拟合的箱线图不使用系数

Question

我尝试了两个类别的简单箱线图（传输类型为 0 或 1），并添加了线性拟合，并希望它通过箱线图的方式，但这里出了点问题：

linearFit <- lm(mpg~am, data=mtcars)
boxplot(mpg~am, data=mtcars) 
abline(linearFit, lwd=2, col="blue")

linearFit$coefficients
# (Intercept)          am 
#   17.147368    7.244939

但是，这条线远远超出了两个箱线图，而不是通过“类别”0 的截距。如何告诉 R 在基本绘图系统中使用类别均值？

编辑澄清： 所以我只需要一条蓝线，它在这个箱线图中穿过 (0, 17.14) 和 (1, 17.14+7.244)，其中变量“am”（传输类型）是一个因子而不是数字变量。 在基础绘图系统中是否有一种简单的方法可以做到这一点？

Answer 1

由于am 被视为一个因子，因此am 的系数是向截距的偏移。因此，您必须将其添加到截距中才能获得 am=1 截距。

linearFit <- lm(mpg~am, data=mtcars)
boxplot(mpg~am, data=mtcars, border=1:2 )
abline(h=c(linearFit$coefficients[1], linearFit$coefficients[1]+linearFit$coefficients[2]), col=1:2, lty=2, lwd=2)

此外，am 术语的重要性可以通过与仅包含截距（即“1”）的模型进行比较来测试

linearFit0 <- lm(mpg~1, data=mtcars) # model with intercept only
anova(linearFit0, linearFit)

编辑 - 在箱线图外绘制

plot(x=c(0,1), res$stats[3,], ylim=range(res$stats))
segments(x0=c(0,1), x1=c(0,1), y0=res$stats[2,], y1=res$stats[4,])
abline(linearFit)

Edit2 - 箱线图的移位值

linearFit <- lm(mpg~I(am+1), data=mtcars)
boxplot(mpg~I(am+1), data=mtcars, border=1:2 )
abline(linearFit$coefficients, col=1:2, lty=2, lwd=2)