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

tmpr <- is.na( wheatfix$species )
tmpr <- range( wheatfix$estimate[tmpr] + wheatfix$ci[tmpr],
   wheatfix$estimate[tmpr] - wheatfix$ci[tmpr], wheatfix$estimate[!tmpr] )
nested.tree(means=wheatfix.means,bar="ci",small=c(0,0,.04),
   ylim=tmpr, xlim=c(-.16,4.5), xlab = "", ylab = "" )
mtext( "yield", 2, 2 )
mtext( "(a) fixed effects estimates", 1, 1.5,
   at = mean( par( "usr" )[1:2] ) )

######################################################################

tmps <- names( wheatran.means )[1]

tmpr <- is.na( wheatran$access )
tmpr <- range( wheatran$estimate[tmpr] + wheatran$ci[tmpr],
   wheatran$estimate[tmpr] - wheatran$ci[tmpr], wheatran$estimate[!tmpr] )
nested.tree( means = wheatran.means, bar = "ci", small = c(0,0.1,0),
   ylim = tmpr, xlab = "", ylab = "" )
#axis( 2, seq( 900, 1100, by = 100 ) )
#axis( 2, seq( 850, 1150, by = 100 ), lab = F )
mtext( paste("(b)", tmps, "random predictions" ), 1, 1.5,
   at = mean( par( "usr" )[1:2] ) )

se.bar( .95, wheatran.means[[tmps]]$estimate,
   2 * wheatfix.means$species[wheatfix.means$species$species==tmps,"ci"],
   cap = "" )
se.bar( 3.5, wheatran.means[[tmps]]$estimate,
   2 * wheatran.means$access$ci[1], cap = "" )

par( tmpar )
rm( tmpr, tmps, tmpar )