回归线到数据点/如何创建垂直线？

Question

如何在R中获得以下可视化（见下文）： 让我们考虑一个简单的三点案例。

# Define two vectors
x <- c(12,21,54)
y <- c(2, 7, 11)

# OLS regression
ols <- lm(y ~ x)

# Visualisation
plot(x,y, xlim = c(0,60), ylim =c(0,15))    
abline(ols, col="red")

我想要的是，从OLS线（红线）到点绘制垂直距离线。

Answer 1

您可以使用ggplot2

非常好地完成此操作

library(ggplot2)
set.seed(1)     
x<-1:10
y<-3*x + 2 + rnorm(10)
m<-lm(y ~ x)
yhat<-m$fitted.values
diff<-y-yhat     
qplot(x=x, y=y)+geom_line(y=yhat)+
       geom_segment(aes(x=x, xend=x, y=y, yend=yhat, color="error"))+
       labs(title="regression errors", color="series")

Answer 2

有一个更简单的解决方案：

segments(x, y, x, predict(ols))

Answer 3

如果你构造一个点矩阵，你可以使用apply来绘制这样的行：

创建坐标矩阵：

cbind(x,x,y,predict(ols))
#   x  x  y          
#1 12 12  2  3.450920
#2 21 21  7  5.153374
#3 54 54 11 11.395706

这可以绘制为：

apply(cbind(x,x,y,predict(ols)),1,function(coords){lines(coords[1:2],coords[3:4])})

有效地在矩阵的行上运行for循环，并为每一行绘制一行。