x<-matrix(scan(file='alloy.dat'),ncol=3,byrow=T)
strength<-x[,3]
fnames<-list(lab=LETTERS[1:6],bar=paste('Bar ',1:4))
alloy.design<-fac.design(c(6,4),fnames)
alloy.df<-data.frame(alloy.design,strength)
aov2<-aov(strength~bar+lab,alloy.df)
