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 ______________________________________________ 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.
