pdf("rpd1_12r.pdf")

#  Example 1.12  in RPD


dust <- c(75,100,150,350,600,900,1300,1650,2250)
rr <- c(1.10,1.05,0.97,1.90,1.83,2.45,3.70,3.52,4.16)

rr_1 <- rr-1

rr_1
dust

plot(dust,rr_1)
abline(lm(rr_1 ~ dust -1))

dustrr.reg <- lm(rr_1 ~ dust -1)

summary(dustrr.reg)

aov(dustrr.reg)

predict(dustrr.reg)

resid(dustrr.reg)

dev.off()