pdf("rpd1_1.pdf")

# RPD Problem 1.1 


ozone <- c(0.02, 0.07, 0.11, 0.15)
yield <- c(242, 237, 231, 201)

ozone
yield

plot(ozone,yield)

ozyld.reg <- lm(yield ~ ozone)

summary(ozyld.reg)

aov(ozyld.reg)

predict(ozyld.reg)

resid(ozyld.reg)

predict(ozyld.reg,int="c")

predict(ozyld.reg,int="p")

dev.off()