Вчера я разработал пример разницы между обычным методом наименьших квадратов (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")
Затем я вручную добавил желтые линии, чтобы получить следующее: