pdf("rpd11_1r.pdf") rad1 <- matrix(c( 0, 118.4, 0.7, 0, 215.2, 2.9, 0, 283.9, 5.6, 0, 387.9, 8.7, 0, 451.5, 12.4, 0, 515.6, 17.4, 1, 109.1, 1.3, 1, 199.6, 4.8, 1, 264.2, 6.5, 1, 358.2, 9.4, 1, 413.2, 12.9, 1, 452.5, 12.3), byrow=T,ncol=3) level <- rad1[,1] rad <- rad1[,2] seedwt <- rad1[,3] sqrtsw <- sqrt(seedwt) reg11_1.reg <- lm(sqrtsw ~ rad + level) reg11_1.rstudent <- rstudent(reg11_1.reg) reg11_1.yhat <- predict(reg11_1.reg) plot(reg11_1.yhat,reg11_1.rstudent,pch=level) abline(h=0) segments(reg11_1.yhat,0,reg11_1.yhat,reg11_1.rstudent) reg11_3.reg <- lm(seedwt ~ rad + level) reg11_3.rstudent <- rstudent(reg11_3.reg) reg11_3.yhat <- predict(reg11_3.reg) plot(reg11_3.yhat,reg11_3.rstudent,pch=level) abline(h=0) segments(reg11_3.yhat,0,reg11_3.yhat,reg11_3.rstudent) reg11_10.inf <- influence.measures(reg11_3.reg) reg11_10.inf dev.off()