R help - Nov 2008 - 2^k designs "anova"

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