

El 01/04/11 13:07, Jorge Ivan Velez escribió:
> Hola Argel,
>
> El nuevo conjunto de datos con que debes alimentar la funcion 
> predict() debe ser un data.frame en el que el nombre de las variables 
> deben corresponder a los utilizados para ajustar el modelo.  Esa 
> informacion puedes obtenerla de ?predict. En tu caso, no 
> proporcionaste el "Block" para el cual quieres hacer la prediccion.
>
> Trata lo siguiente:
>
> > predict(CR.rsm2,
>      
> newdata=data.frame(x1=summary(CR.rsm2)$canonical$xs[1],x2=summary(CR.rsm2)$canonical$xs[2], 
> Block = 'B1'))
>       x1
> 84.36561
> > predict(CR.rsm2,
> newdata=data.frame(x1=summary(CR.rsm2)$canonical$xs[1],x2=summary(CR.rsm2)$canonical$xs[2], 
> Block = 'B2'))
>       x1
> 79.90808
>
> Un comentario adicional:    La respuesta que se obtiene en este caso 
> no es el valor de "x1" utilizando el modelo CR.rsm2 sino el valor de 
> "Yield".
>
>
> Saludos,
> Jorge Ivan Velez

     Hola de nuevo, muchas gracias a Jorge Ivan, Luciano y Carlos por 
sus respuestas. Ahora entiendo cuál fue el error al usar la función 
"predict".

     Como todavía estoy aprendiendo a trabajar con las superficies de 
respuesta, ahora me queda la duda de cuál de las dos predicciones para 
la respuesta "Yield" es la más adecuada, si la hecha con el Block "B1" o 
con el "B2". Por la figura obtenida con la función "contour" me parece 
que la hecha con el "B1". ¿Ustedes qué opinan?

     Saludos.

--
     Argel.

	[[alternative HTML version deleted]]

