set.seed(123)
x <- rnorm(100)
y <- 2 * x + rnorm(100)
model <- lm(y ~ x)
x_new <- seq(min(x), max(x), length.out = 100) 
y_pred <- predict(model, newdata = data.frame(x = x_new))
plot(x, y, main = "Regresi Linear Gaussian", xlab = "x", ylab = "y")
lines(x_new, y_pred, col = "red")