stat545a-2013-hw01_ni-jac.r

Jack — Sep 12, 2013, 5:59 PM

n <- 50
a <- 2
b <- -3
sigSq <- 0.5
x <- runif(n)
y <- a + b * x + rnorm(n, sd = sqrt(sigSq))
(avgX <- mean(x))

[1] 0.4113

write(avgX, "avgX.txt")
plot(x, y)
abline(a, b, col = "purple")

plot of chunk unnamed-chunk-1 

dev.print(pdf, "toylinePlot.pdf")

pdf 
  2