Leah — Sep 13, 2013, 10:21 AM
a <- 5 b <- 10 n <- 100 sigSq <- 0.5 x <- runif(n) y <- a + b * x + rnorm(n,sd=sqrt(sigSq)) (avgX<-mean(x))
[1] 0.5114
plot(x,y) abline(a,b, col="purple")