par(mfrow=c(3,2),bty='l',las=1,oma=c(0,0,0,0))
x<-matrix(scan(file='fisher.dat'),ncol=4,byrow=T)
yield<-x[,4]
fnames<-list(var=LETTERS[1:5],year=paste('Year ',1:2),
place=paste('Pl ',1:6))
fisher.design<-fac.design(c(5,2,6),fnames)
fisher.df<-data.frame(fisher.design,yield)
aov3<-aov(yield~var*place,fisher.df)
aov4<-aov(yield~year+var+place,fisher.df)
plot.design(fisher.df,main='(a)')
plot.factor(fisher.df,main='(bbb)')
attach(fisher.df)
interaction.plot(year,place,yield,main='(e)')
plot(fitted.values(aov4),residuals(aov4),xlab='Fitted Values',
ylab='Residuals',main='(f)')
