stat545a-2013-hw01_ma-hui.R

user — Sep 13, 2013, 9:29 PM

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

[1] 0.4353

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