You ought to be able to use the reshape command, with basic syntax
reshape long stubnames, i(varlist) [options]
to convert your data which is now in the wide format into the long format. Then you should be able to generate the plots that you want and get the means by using the bysort prefix.
N.B. you also can always type "help reshape" or the equivalent into stata and it will give you the syntax and often example(s). Also, of course you can always write back to the forum.