tmpar <- par(mar=c(3.1,2.1,0.1,0),mfrow=c(1,2))

#tomato1 <- tomato[tomato$yr==1&!is.na(tomato$mfwlog)&!is.na(tomato$tg430),]
#tomato.aov <- aov(mfwlog~tg430,tomato1)

# See Data/Cloning/clonci.s

alpha <- .05
tmp <- ci.plot(tomato.aov,level=alpha,
	crit=qt(1-alpha/2,df.resid(tomato.aov))/sqrt(2),
	xlab="",width=.1, ylim=c(0.2,1.2))
mtext(paste("(a) ",100*(1-alpha),"% CI / sqrt(2)",sep=""),1,2)
#mtext("log flower time",2,1)

attach(tomato)
boxplot(split(mfwlog[yr==1],tg430[yr==1]),notch=T,xlab="",ylab="",
	ylim=c(0.2,1.2))
detach()
mtext("(b) notched box-plots",1,2)

par(tmpar)