################################################### ### code chunk number 28: reg31b ################################################### fit1 <- plotPlane(m2centered, plotx1="x1c", plotx2="x2c", ticktype="detailed")