jenny — Sep 6, 2013, 3:15 PM
a <- 2
b <- 7
sigSq <- 0.5
n <- 400
set.seed(1234)
x <- runif(n)
y <- a + b * x + rnorm(n, sd = sqrt(sigSq))
(avgX <- mean(x))
[1] 0.4969
plot(x, y)
abline(a, b, col = "blue", lwd = 2)
sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-apple-darwin10.8.0 (64-bit)
locale:
[1] en_CA.UTF-8/en_CA.UTF-8/en_CA.UTF-8/C/en_CA.UTF-8/en_CA.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] knitr_1.4.1 roxygen2_2.2.2 digest_0.6.3 lattice_0.20-15
loaded via a namespace (and not attached):
[1] brew_1.0-6 evaluate_0.4.7 formatR_0.9 grid_3.0.1
[5] stringr_0.6.2 tools_3.0.1