> You would get better advice, I bet, if you told us exactly what LSMEANS
does. 
LSMEANS produces the adjusted mean for each class of a variable of interest.
Exemple : 
proc glm;
 class age Cas sexec; 
 model vitA  = age Cas sexec;
 lsmeans Cas / stderr;
run;
In this case, lsmean give the mean ajusted by sexe and age of vitamine A for each class of the variable Cas : case and control.
Equivalent in STATA : ???
Thank's for your help !
Cecile
*
*   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/