stat545a-2013-hw01_gao-wen.html

Sibyl — Sep 14, 2013, 3:22 PM

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

[1] 0.4511

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

plot of chunk unnamed-chunk-1

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