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

发布于03月27日

我正在绘制一个适用于R中Loblolly数据的线性回归，并试图 for each 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，我不应该使用它.有什么帮助吗？

推荐答案

正如@AndroMine在 comments 中所说的，你可以用points中的aggregate来获得平均值. Select 一个好的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)