#ps("/tmp/p",he=4.5,wi=4.5,po=10)

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

tmp <- c(40,20,40)
names(tmp) <- c("A","H","B")
attach(tomatos)
for (i in levels(tg430[yr==1]))
{
	print(i)
	hist(mfwlog[yr==1][tg430[yr==1]==i],xlab="",yaxt="n",xlim=c(-.3,1.74),
		nclass=tmp[i])
	mtext(paste("Allele",i),1,2,at=mean(par("usr")[1:2]),cex=1)
}

par(tmpar)