如何在回归图中添加每个 x 水平的均值?

How to add means of each x level in a regression plot?

我正在为 R 中的 Loblolly 数据绘制线性回归拟合,并尝试为每个 x 水平添加每个平均 y。我的代码如下:

 plot(Loblolly$height~Loblolly$age, xlab = "Age", ylab = "Height", main = "Loblolly Data")
    lob.fit <- lm(Loblolly$height~Loblolly$age)
    lob.sum <- summary(lob.fit)
    str(lob.sum)
    abline(lob.fit)

这可能是一个非常容易回答的问题,但我不知道如何添加方法。当我 google 时,所有显示的都是 ggplot,我不应该使用它。有帮助吗?

正如@AndroMine 在评论中所述,您可以在 points 中使用 aggregate 来获取方法。选择一个不错的 pch=,例如 "-",给它一个 col=or 并增加 cex= 的大小。请注意,您使用的函数有一个 data= 参数,因此您不必继续输入 Loblolly$.

plot(height ~ age, Loblolly, xlab="Age", ylab="Height", main="Loblolly Data")
abline(lm(height ~ age, Loblolly))
points(aggregate(height ~ age, Loblolly, mean), pch='-', cex=4, col=4)