tmpar <- par(mar=c(4.1,4.1,0,0)) attach(tomato) boxplot(split(mfwlog[yr==1],tg430[yr==1]), xlab="Tomato Allele Type",ylab="Log Flower Time") detach() par(tmpar)