データ

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)