x <- 1:30 y <- seq(1, 300, 10) + rnorm(30)
y[5] <- y[5] + 10 y[10] <- y[10] + 5 y[25] <- y[25] - 10
fit <- lm(y ~ x) plot(fit)