How to output the LSMEANS differences?

[Hello_there]

Hi, 

I'm running a model in proc GLM. 

That looks something like this: 

PROC GLM data=a;
  by subject;
class activity;
model cont_var=activity;
lSMEANS activity / adjust = Tukey out=b;
run; quit;

The activities are a. (reference), b., c., d.

The current code only outputs cont_var LSMEANS and the p values for the difference from a.

How would I be able to output the LSMEANS difference from a?

[svh]

This sounds like a job for the ODS tables, which allow you to export certain statistics from PROC GLM to an output data set.

https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.4/statug/statug_glm_details70.htm

E.g., you might be looking to add a line like this to your code:
ods output LSMeans=LSMEANS_diffs;

[Hello_there]

Thank you for giving me this link, it was helpful!

[sbxkoenk · Accepted solution]

Hello,

Is this what you are after?

*ods trace on;
*ods output Diff=work.Difference_Matrix;
ods output  LSMeans=work.LSMeans;

PROC GLM data=sashelp.class;
*by subject;
class age;
model height=age;
LSMEANS age / adjust = Tukey out=b;
run; 
quit;

data work.LSMeans_DIFF;
 set work.LSMeans;
 retain ref .;
 if _N_=1 then do; ref=LSMean; end;
 Diff = LSMean - ref;
run;
/* end of program */

Koen

[Hello_there]

Thanks, this worked!