###################################################
### chunk number 1: sim1 eval=FALSE
###################################################
## #line 52 "/home/bates/st849/Fall10/hwk4.Rnw"
## underfit <- lm(Ysim ~ 1 + x1 + x2)


###################################################
### chunk number 2: sim2 eval=FALSE
###################################################
## #line 58 "/home/bates/st849/Fall10/hwk4.Rnw"
## overfit <- lm(Ysim ~ 1 + x1 + x2 + x3 + x4)


###################################################
### chunk number 3: job
###################################################
#line 109 "/home/bates/st849/Fall10/hwk4.Rnw"
DataURL <- "http://www.stat.wisc.edu/~st849-1/data/"
job <- read.table(paste(DataURL, "jobdata.txt", sep=''),
                  col.names=c("Y","X1","X2","X3","X4"))
str(job)