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

Question

我正在为 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)

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

Answer 1

正如@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)