[R] Response surface methodology

Peter.Watkins@foodscience.afisc.csiro.au Peter.Watkins at foodscience.afisc.csiro.au
Fri Sep 27 05:25:40 CEST 2002


Dear all,

I am attempting to use RSM in order to optimise the analysis of volatile
compounds. After changing two variables and performing a series of
measurements, I get the following:

    X1   X2   resp
1  -0.8 -1.0 1101039  
2  -0.8  0.0  809947   
3  -0.8  1.0  814664 
4   0.0 -1.4  759484  
5   0.0  0.0  665916  
6   0.0  0.0  654054  
7   0.0  0.0  641134
8   0.0  0.0  673143 
9   0.0  0.0  705040
10  0.0  1.4  793333
11  1.0 -1.0  650786
12  1.0  1.0  525182 

and analyse the data with 

resp.lm <- lm(resp ~ X1 + X2 + X1*X2+I(X1^2) + I(X2^2)) where X1 and X2 are
the encoded variables.

I'd like to see a plot of the data - could someone suggest how to do so? I
think I could use 'contour' but I'm not certain of the syntax for the
command.

TIA, Peter

Peter Watkins
Research Scientist

Food Science Australia
Private Bag 16
(671 Sneydes Road)
Werribee
Victoria 3030
Australia

Tel: +61-3-9731 3467   Fax: +61-3-9731 3250
E-mail: Peter.Watkins at foodscience.afisc.csiro.au
Web: www.foodscience.afisc.csiro.au

Food Science Australia is a joint venture of CSIRO and Afisc.
	


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._



More information about the R-help mailing list