pdf("F:\\sta6208\\sc_wine.pdf") sc_wine <- read.table("http://www.stat.ufl.edu/~winner/data/sugarcane_wine.dat",header=F, col.names=c("x1","x2","x3","sugar","pH","temp","solids","eth","col","aro","flav","overall")) attach(sc_wine) library(rsm) sc_wine.rsm1 <- rsm(overall ~ SO(sugar, pH, temp)) summary(sc_wine.rsm1) par(mfrow=c(1,3)) contour(sc_wine.rsm1, ~ sugar + pH + temp, at=summary(sc_wine.rsm1)$canonical$xs) dev.off()