SOLVED:

 xs <- canonical(smell.rsm)$xs
> xs
       x1        x2        x3
0.1219125 0.1995746 1.7705249
> persp(smell.rsm,~x1+x2,at=xs,contour=TRUE)
>




On Tue, Dec 3, 2013 at 9:33 PM, Erin Hodgess <erinm.hodg...@gmail.com>wrote:

> Dear R People:
>
> I am using Response Surface Models and and also using the "persp" function
> to plot them.
>
> Today my particular model is y~ x1 + x2 + x3
>
> I can fit my RSM easily.  My question is:  how do I use the persp such
> that I can fix the x3 at a particular value, please?
>
> > smell.code <- coded.data(smell.df,x1~(T-80)/40,x2~(R-0.5)/0.2,x3~(H-4)/2)
> > smell.rsm <- rsm(Odor~SO(x1,x2,x3),data=smell.code)
> >
>
> > persp(smell.rsm,x1~x2,contour=TRUE)
> >
>
> I would like to set the x3 to be at the optimum value (available from the
> model)
>
> Thanks in advance for any help.
>
>
> Sincerely,
> Erin
>
> --
> Erin Hodgess
> Associate Professor
> Department of Computer and Mathematical Sciences
> University of Houston - Downtown
> mailto: erinm.hodg...@gmail.com
>



-- 
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
mailto: erinm.hodg...@gmail.com

        [[alternative HTML version deleted]]

______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.

Reply via email to