topic Re: How to output the LSMEANS differences? in SAS Programming

How to output the LSMEANS differences?

Hello_there — Fri, 25 Mar 2022 16:32:52 GMT

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?

Re: How to output the LSMEANS differences?

svh — Fri, 25 Mar 2022 17:52:31 GMT

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;

Re: How to output the LSMEANS differences?

sbxkoenk — Fri, 25 Mar 2022 17:54:36 GMT

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

Re: How to output the LSMEANS differences?

Hello_there — Sat, 26 Mar 2022 14:10:13 GMT

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

Re: How to output the LSMEANS differences?

Hello_there — Sat, 26 Mar 2022 14:10:36 GMT

Thanks, this worked!