No tengo mucha idea del paquete rsm. Pero lo que probablemente te este
diciendo el error es que te falta la columna Block en los nuevos datos para
predecir.

Espero que mi respuesta te sea de ayuda.

Luciano


El 1 de abril de 2011 15:50, Argel Gastélum Arellánez <
argel.gastelum@gmail.com> escribió:

>    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.
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>

	[[alternative HTML version deleted]]

