x <- runif(100, -5, 5)
y <- 0.5 + 0.75 * x + rnorm(100)
data <- data.frame(x,y)



Y0 <- 0.5 + 0.75 * 1.5 + rnorm(20)