session01_toy-line.r

jenny — Jan 27, 2014, 11:08 AM

a <- 2
b <- 4
sigSq <- 0.5
n <- 40

x <- runif(n)
(avgX <- mean(x))

[1] 0.4957

write(avgX, "avgX.txt")

y <- a + b * x + rnorm(n, sd = sqrt(sigSq))

plot(x, y)
abline(a, b, col = "green", lwd = 2)

plot of chunk unnamed-chunk-1


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

pdf 
  2