# WLoss3 Example -- source code file wloss3.r wloss = read.table('wloss.txt', header=T) model = lm(loss ~ A + B + C + D, data=wloss) summary(model) p = predict(model) r = residuals(model) print(cbind(p, r)) pdf("wloss3.pdf") plot(p, r, main="Residual Plot", xlab="Predicted Values", ylab="Residuals") abline(h=0, col="red") qqnorm(r, xlab="Predicted Normal Scores", ylab="Residuals") dev.off( )