r - 从 R 中的实际点到建模点的线

Question

昨天我worked up an example of普通最小二乘法 (OLS) 与主成分分析 (PCA) 之间的差异。对于那个插图，我想展示 OLS 和 PCA 最小化的误差，所以我绘制了实际值、预测线，然后我手动(使用 GIMP)绘制了一条垂线来说明几个误差项。如何在 R 中编写错误行的代码？这是我用于示例的代码:

set.seed(2)
x <- 1:100

y <- 20 + 3 * x
e <- rnorm(100, 0, 60)
y <- 20 + 3 * x + e

plot(x,y)
yx.lm <- lm(y ~ x)
lines(x, predict(yx.lm), col="red")

然后我手动添加了黄线以产生以下内容:

Answer 1

?segments
我会提供一个例子，但我今天很忙，选择要点并没有那么复杂。 ;-)

好吧，我没那么忙……

n=58; segments(x[n],y[n],x[n],predict(yx.lm)[n])
n=65; segments(x[n],y[n],x[n],predict(yx.lm)[n])