topic PROC GLIMMIX Contrasts in Statistical Procedures

PROC GLIMMIX Contrasts

mthorne — Tue, 17 Nov 2020 20:27:47 GMT

Hi,

I am using proc glimmix to analyze count data. I notice that the pvalues for a contrast between two simple effects are the same as the pvalues for lsmeans diff for the same simple effect. If I use the contrast comparison, do I report the actual means or the lsmeans?

Re: PROC GLIMMIX Contrasts

PaigeMiller — Tue, 17 Nov 2020 20:57:51 GMT

Depending on how you did the contrast, it can indeed be the same as the LSMEANs difference. Not just the p-values should match, but the actual difference should match.

Again, depending on what you did, I would report the LSMEANs and their differences as a general rule, since I have seen only very rare situations where the contrast is actually needed. I would report Means instead of LSMEANS only if you have a completely balanced and orthogonal experiment, or you have the equivalent of a one-way ANOVA.

Re: PROC GLIMMIX Contrasts

mthorne — Tue, 17 Nov 2020 21:19:29 GMT

Thanks for the reply!

This is my model.

ods graphics on; 
proc glimmix data=counts method=laplace nobound plot=studentpanel; by loc time;
class loc time herb rep;
tc2=log(c2+1);
model tc2 = herb / dist=normal link=identity;
random intercept / subject=rep;
lsmeans herb / ilink diff; 
contrast '0 vs 1' herb 1 -1 0;
contrast '0 vs 2' herb 1 0 -1;
contrast '1 vs 2' herb 0 1 -1;
*output out=allresid pred=p resid=student;
run; quit;
ods graphics off;

The design is a balanced factorial with no missing data. Three locations, four application times, three herbicide treatments, and four replicate blocks at each location. Glimmix couldn't handle all zeros from some treatments, so I had to transform and compare to a normal distribution.

Re: PROC GLIMMIX Contrasts

PaigeMiller — Tue, 17 Nov 2020 22:31:55 GMT

contrast '0 vs 1' herb 1 -1 0;
contrast '0 vs 2' herb 1 0 -1;
contrast '1 vs 2' herb 0 1 -1;

There's no need to type all this out, the LSMEANS statement with the DIFF produces the same results.

Re: PROC GLIMMIX Contrasts

mthorne — Wed, 18 Nov 2020 00:18:02 GMT

Great, thanks Paige!

Mark