tmpar <- par( mar = c(3,3,0,0) )
on.exit( { par( tmpar ) ; rm( tmpar ) } )

tmpy <- cor( wasp[,2:14], wasp$DA )[ names( wasp.t[-1] ), 1 ]

mplot( c(1,-1) * wasp.t["DA"], c(1,-1), rep( "DA", 2), fulltext = T,
   xlab ="", ylab = "" )

lines( wasp.t[-1], tmpy )
points( wasp.t[-1], tmpy )
text( wasp.t[-1] + 5 * sign( wasp.t[-1] ), uncollide(tmpy,.05), names( wasp.t[-1] ) )

abline(h=0,v=0,lty=2)
mtext("queen larger",1,2,adj=1)
mtext("worker larger",1,2,adj=0)
mtext("t statistics",1,2)
mtext("correlation with DA",2,2)