YimingZhang — Sep 12, 2013, 2:31 PM
n <- 55 a <- 2 b <- -3 sigSq <- 0.5 x <- runif(n) y <- a+b*x+rnorm(n,sd=sqrt(sigSq)) (avgX <- mean(x))
[1] 0.5082
plot(x,y) abline(a,b,col="red")