Mina — Sep 12, 2013, 11:05 AM
a <- 2
b <- -3
n <- 400
sigSq <- 0.5
set.seed(1234) #set the seed for random number generator
x <- runif(n)
y <- a + b * x + rnorm(n, sd = sqrt(sigSq))
(avgX <- mean(x))
[1] 0.4969
plot(x, y)
abline(a, b, col = "purple")