topic LSD test in PROC GLIMMIX in Statistical Procedures

LSD test in PROC GLIMMIX

Raman3 — Thu, 09 Apr 2026 14:35:12 GMT

Hi SAS community,

I am trying to analyze grain yield data which follows a normal distribution. I would like get lsmeans comparisons using LSD. But, in the lsmeans statement, adjust option does not provide LSD. Below is the code that I am using in Proc Glimmix (using tukey method)-

proc glimmix data=Yieldkgha plots=residualpanel;

class Treatment SiteYear Block Plot;

model Yieldkgha = Treatment|SiteYear / dist= gaussian;

random _Residual_/ group=Treatment ;

Random Block(SiteYear);

/*calculate lsmeans for all possible treatment and crop and adjust the values using Tukey*/

lsmeans Treatment|SiteYear / adjust=tukey lines cl;

/*covtest "common variance" homogeneity;*/

/*Checking assumptions of ANOVA */

output out = resids predicted=pred residual=res;

ods output tests3=anova

lsmeans=lsmeansYieldkgha

lsmlines=lslinesYieldkgha;

run;

Can someone inform me how to conduct LSD in Proc Glimmix?

Re: LSD test in PROC GLIMMIX

data_null__ — Thu, 09 Apr 2026 16:09:36 GMT

This suggests that you simply omit ADJUST

proc glimmix data=sashelp.class;
   class age;
   model weight = age;
   lsmeans age / lines;
   run;

Re: LSD test in PROC GLIMMIX

Ksharp — Fri, 10 Apr 2026 06:26:32 GMT

Plz post it at Statistical Forum:

https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures

Re: LSD test in PROC GLIMMIX

sbxkoenk — Fri, 10 Apr 2026 08:10:28 GMT

I moved this topic from "New Users" to "Statistical Procedures"-board.

Re: LSD test in PROC GLIMMIX

sbxkoenk — Fri, 10 Apr 2026 08:16:59 GMT

Can someone inform me how to conduct LSD in Proc Glimmix?

The short answer is: you cannot get the LSD value from PROC GLIMMIX or PROC MIXED. There is such a thing as a least significant difference (LSD) only if all the cells have the same size and the model is simple enough so that the significance of all comparisons can be judged on the same scale. This is often not the case for mixed models, as a result, no LSD value is available in PROC GLIMMIX or PROC MIXED.

You may specify the LSMEANS statement with ADJUST=TUKEY and the LINES option ...

... and that will be the closest method I could suggest.

BR, Koen