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

Carlos J. Gil Bellosta cgb en datanalytics.com
Vie Abr 1 21:07:59 CEST 2011


Pues... ¡añádelo!

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

Un saludo,

Carlos J. Gil Bellosta
http://www.datanalytics.com

El día 1 de abril de 2011 20:50, Argel Gastélum Arellánez
<argel.gastelum en 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 en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>



Más información sobre la lista de distribución R-help-es