Amanda — Sep 12, 2013, 5:24 PM
a <- 2 b <- -3 n <- 40 sigSq <- 0.5 x <- runif(n) y <- a + b * x + rnorm(n, sd = sqrt(sigSq)) (avgX <- mean(x))
[1] 0.5453
write(avgX, "avgX.txt") plot(x,y) abline(a, b, col = "purple")
