[R-es] Predicción de valor máximo en superficie de respuesta, con paquete rsm

[Argel Gastélum Arellánez]

     Hola compañeros de la lista. Estoy aprendiendo a usar el paquete 
"rsm" para superficies de respuesta. Siguiendo este ejemplo todo va 
bien, hasta que trato de obtener el valor máximo predicho para la 
variable de respuesta en los valores de x1 y x2 estimados.

-------------------------------------------------------------
library("rsm")
ChemReact
CR <- coded.data(ChemReact, x1 ~ (Time - 85)/5, x2 ~ (Temp - 175)/5)
CR
# Realizando un ajuste de segundo orden:
CR.rsm2 <- rsm(Yield ~ Block + SO(x1, x2), data = CR)
CR.rsm2
summary(CR.rsm2)
# Graficando:
contour(CR.rsm2, ~ x1 + x2, at = summary(CR.rsm2)$canonical$xs, image = 
TRUE, img.col = terrain.colors(50))
xs <- canonical(CR.rsm2)$xs
points (xs[1], xs[2], pch=2, col="red")
points(CR[1:7,], col="black", pch=19)
points(CR[8:14,], col="blue", pch=19)
#Esta última línea da un error:
predict(CR.rsm2, 
newdata=data.frame(x1=summary(CR.rsm2)$canonical$xs[1],x2=summary(CR.rsm2)$canonical$xs[2]))
--------------------------------------------------------------

La última línea me arroja el siguiente error:

Error en eval(expr, envir, enclos) : objeto 'Block' no encontrado

     De antemano agradezco la ayuda.

     Saludos.

--
     Argel.