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(x1,x2,x3)) summary(sc_wine.rsm1) par(mfrow=c(1,3)) contour(sc_wine.rsm1, ~ x1 + x2 + x3, at=summary(sc_wine.rsm1)$canonical$xs) dev.off() > 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) The following object(s) are masked from 'sc_wine (position 3)': aro, col, eth, flav, overall, pH, solids, sugar, temp, x1, x2, x3 > > library(rsm) > > sc_wine.rsm1 <- rsm(overall ~ SO(x1,x2,x3)) > summary(sc_wine.rsm1) Call: rsm(formula = overall ~ SO(x1, x2, x3)) Estimate Std. Error t value Pr(>|t|) (Intercept) 3.546667 0.098714 35.9287 3.144e-07 *** x1 0.232500 0.060450 3.8462 0.012048 * x2 -0.132500 0.060450 -2.1919 0.079905 . x3 -0.110000 0.060450 -1.8197 0.128460 x1:x2 -0.045000 0.085489 -0.5264 0.621127 x1:x3 -0.090000 0.085489 -1.0528 0.340641 x2:x3 0.030000 0.085489 0.3509 0.739948 x1^2 -0.165833 0.088980 -1.8637 0.121388 x2^2 -0.375833 0.088980 -4.2238 0.008297 ** x3^2 -0.150833 0.088980 -1.6951 0.150818 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 Analysis of Variance Table Response: overall Df Sum Sq Mean Sq F value Pr(>F) FO(x1, x2, x3) 3 0.66970 0.223233 7.6363 0.02584 TWI(x1, x2, x3) 3 0.04410 0.014700 0.5029 0.69678 PQ(x1, x2, x3) 3 0.63813 0.212709 7.2762 0.02843 Residuals 5 0.14617 0.029233 Lack of fit 3 0.14110 0.047033 18.5658 0.05154 Pure error 2 0.00507 0.002533 Stationary point of response surface: x1 x2 x3 0.9160124 -0.2575967 -0.6635437 Eigenanalysis: $values [1] -0.1101646 -0.2036432 -0.3786922 $vectors [,1] [,2] [,3] x1 0.64920255 0.75462858 -0.09524577 x2 -0.09757273 -0.04156234 -0.99436016 x3 -0.75433123 0.65483454 0.04664880