beta0 = 0.57
beta1 = 5
n = 50
x1 = rnorm(n,2.45,1.85)
epsilon = rnorm(n,mean = 0, sd=sqrt(1.64))
y = beta0 + beta1*x1+ epsilon
data = data.frame(x1,y)
head(data)