Hi R users:
How can I obtain the same "anova" table
for the effects for a 2^k experiment design that MINITAB shows
(and authors recommends Box, Hunter, & Hunter).
http://www.stat.psu.edu/online/development/stat503/06_2k/04_2k_unreplicate.html
Here is the code that I use for this case:
D<-C<-B<-A<-c("-","+")
design<-expand.grid(A=A,B=B,C=C,D=D)
design$yield<-c(45,71,48,65,68,60,80,65,43,100,45,104,75,86,70,96)
model1<-aov(yield~A+C+D+A:C+A:D,data=design)
summary(model1)
model.tables(model1,se=T)
Thank you for your help.
Kenneth
