Updated iyprplot.R to use se=sum(residuals(lm(d[,pbindaff]~d[,predicted]))^2)