st: Mimicking graph bar's multiple over options in ciplot

[Venable <venablito@gmail.com>]

Dear Statalist,

I would like to create a graph of estimates of means with standard
error bars. ciplot seems like an excellent tool for this. I would like
to make one small change, though: I would like to arrange the
estimates in a format similar to that produced by graph bar with
multiple over() options.

As an example:

sysuse citytemp, clear

// ciplot produces the content I want
//  mean plotted with a point
//  confidence interval around point estimate plotted with capped bar
ciplot tempjan, by(division)

// graph bar produces the formatting I want
//  separating first by region and then by division within region
gr bar (mean) tempjan, over(division) over(region) nofill


How could I produce a graph that combines these?

Many thanks.
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/