Hi,
The below code gives me exactly what I need. I wonder if there is any way to store the 2.5% and 97.5% centiles as well.
Thanks,
NM

statsby "ci means ei_x`i', level(95) " N=r(N) mean=r(mean) ub=r(ub) lb=r(lb) se=r(se) , by(profile_new age )