# fit model, collecting information on projections
# caution: qr=T option for proj() may take a long time
tomato.fit <- aov(mfwlog~tg430+yr+yr:tg430+Error(entry/tg430),tomatos,qr=T)

tmpar <- par(mfcol=c(1,2),mar=c(4,4,.5,.5))

effect.plot(tomato.fit,stratum="entry",
   xlab="(a) entry effects",ylab="")
effect.plot(tomato.fit,stratum="Within",
   xlab="(b) year effects")

par(tmpar)