There is no variable called x1 in the chem dataframe.
Did you mean
gg1=rsm(yield~ FO(x1,x2),data=ggg,subset = (Block == "B1"))
hth.
daniel
________________________________________
Felad?: r-help-bounces at r-project.org [r-help-bounces at r-project.org] ;
meghatalmazó: IZHAK shabsogh [ishaqbaba at yahoo.com]
K?ldve: 2013. december 15. 15:53
To: r-help at r-project.org
T?rgy: [R] RSM
kindly help find out why the following code is given me error
time<-c(80,80,90,90,85,85,85,85,85,85,92,77,85,85)> tem<-c(170,180,170,180,175,175,175,175,175,175,175,175,182.07,167.93)
>
yield<-c(80.5,81.5,82,83.5,83.9,84.3,84,79.7,79.8,79.5,78.4,75.6,78.6,77)
>
Block<-c("B1","B1","B1","B1","B1","B1","B1","B2","B2","B2","B2","B2","B2","B2")
> chem<-data.frame(time,tem,Block,yield)
> ggg<-coded.data(chem, x1 ~ (time - 85)/5, x2 ~ (tem - 175)/5)
> gg1=rsm(yield~ FO(x1,x2),data=chem,subset = (Block == "B1"))
Error in lapply(X = X, FUN = FUN, ...) : object 'x1' not
found> gg1=rsm(yield~ SO(x1,x2),data=chem,subset = (Block == "B1"))
Error in lapply(X = X, FUN = FUN, ...) : object 'x1' not found
[[alternative HTML version deleted]]
______________________________________________
R-help at 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.