st: adding mean plot to anovaplot

[John Novak <jnovak@math.usc.edu>]

I would like to add a plot of the treatment means collapsed across the by 
variable to an -anovaplot-.  I have done this:

#delimit ;
quietly anova y a b a*b;
anovaplot ,
        scatter(msymbol(i) xsize(3) ysize(3) name(by_b, replace))
        plot(mband y a) ;
delimit cr

It is almost what I want, but adds the median band instead of a mean 
band.  Does anyone know how I can accomplish the same effect, but with 
means instead of medians?

*
*   For searches and help try:
*   http://www.stata.com/support/faqs/res/findit.html
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/