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

predict(CR.rsm2, newdata=data.frame(x1=summary(CR.rsm2)$canonical$xs[1],Block = "B1",x2=summary(CR.rsm2)$canonical$xs[2]))

predict(CR.rsm2, newdata=data.frame(x1=summary(CR.rsm2)$canonical$xs[1],Block = "B2",x2=summary(CR.rsm2)$canonical$xs[2]))

   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.