stat545a-2013-hw01_wang-ton.R

tongwang — Sep 12, 2013, 12:03 PM

a <- 2
b <- 7
sigSq <- 0.5
n <- 400

set.seed(1234)
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 = "blue", lwd = 2)