stat545a-2013-hw01_okkels-chr.R

Christian — Sep 14, 2013, 5:01 PM

# STAT 545A - Homework 01 - Christian Okkels
n <- 500
a <- 1
b <- 3
sigSq <- 0.5
x <- runif(n)
y <- a + b * x + rnorm(n, sd = sqrt(sigSq))
(avgX <- mean(x))

[1] 0.5192

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

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

pdf 
  2