Documentation Example 2 for PROC GLIMMIX
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: gmxex02 */
/* TITLE: Documentation Example 2 for PROC GLIMMIX */
/* Mating Experiment with Crossed Random Effects */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: Generalized linear mixed models */
/* Binary data */
/* PROCS: GLIMMIX */
/* DATA: Salamander mating data */
/* */
/* SUPPORT: Oliver Schabenberger */
/* REF: McCullagh, P. and Nelder. J.A. (1989) */
/* Generalized Linear Models, Second Edition */
/* London: Chapman and Hall */
/* MISC: */
/****************************************************************/
data salamander;
input day fpop$ fnum mpop$ mnum mating @@;
datalines;
4 rb 1 rb 1 1 4 rb 2 rb 5 1
4 rb 3 rb 2 1 4 rb 4 rb 4 1
4 rb 5 rb 3 1 4 rb 6 ws 9 1
4 rb 7 ws 8 0 4 rb 8 ws 6 0
4 rb 9 ws 10 0 4 rb 10 ws 7 0
4 ws 1 rb 9 0 4 ws 2 rb 7 0
4 ws 3 rb 8 0 4 ws 4 rb 10 0
4 ws 5 rb 6 0 4 ws 6 ws 5 0
4 ws 7 ws 4 1 4 ws 8 ws 1 1
4 ws 9 ws 3 1 4 ws 10 ws 2 1
8 rb 1 ws 4 1 8 rb 2 ws 5 1
8 rb 3 ws 1 0 8 rb 4 ws 2 1
8 rb 5 ws 3 1 8 rb 6 rb 9 1
8 rb 7 rb 8 0 8 rb 8 rb 6 1
8 rb 9 rb 7 0 8 rb 10 rb 10 0
8 ws 1 ws 9 1 8 ws 2 ws 6 0
8 ws 3 ws 7 0 8 ws 4 ws 10 1
8 ws 5 ws 8 1 8 ws 6 rb 2 0
8 ws 7 rb 1 1 8 ws 8 rb 4 0
8 ws 9 rb 3 1 8 ws 10 rb 5 0
12 rb 1 rb 5 1 12 rb 2 rb 3 1
12 rb 3 rb 1 1 12 rb 4 rb 2 1
12 rb 5 rb 4 1 12 rb 6 ws 10 1
12 rb 7 ws 9 0 12 rb 8 ws 7 0
12 rb 9 ws 8 1 12 rb 10 ws 6 1
12 ws 1 rb 7 1 12 ws 2 rb 9 0
12 ws 3 rb 6 0 12 ws 4 rb 8 1
12 ws 5 rb 10 0 12 ws 6 ws 3 1
12 ws 7 ws 5 1 12 ws 8 ws 2 1
12 ws 9 ws 1 1 12 ws 10 ws 4 0
16 rb 1 ws 1 0 16 rb 2 ws 3 1
16 rb 3 ws 4 1 16 rb 4 ws 5 0
16 rb 5 ws 2 1 16 rb 6 rb 7 0
16 rb 7 rb 9 1 16 rb 8 rb 10 0
16 rb 9 rb 6 1 16 rb 10 rb 8 0
16 ws 1 ws 10 1 16 ws 2 ws 7 1
16 ws 3 ws 9 0 16 ws 4 ws 8 1
16 ws 5 ws 6 0 16 ws 6 rb 4 0
16 ws 7 rb 2 0 16 ws 8 rb 5 0
16 ws 9 rb 1 1 16 ws 10 rb 3 1
20 rb 1 rb 4 1 20 rb 2 rb 1 1
20 rb 3 rb 3 1 20 rb 4 rb 5 1
20 rb 5 rb 2 1 20 rb 6 ws 6 1
20 rb 7 ws 7 0 20 rb 8 ws 10 1
20 rb 9 ws 9 1 20 rb 10 ws 8 1
20 ws 1 rb 10 0 20 ws 2 rb 6 0
20 ws 3 rb 7 0 20 ws 4 rb 9 0
20 ws 5 rb 8 0 20 ws 6 ws 2 0
20 ws 7 ws 1 1 20 ws 8 ws 5 1
20 ws 9 ws 4 1 20 ws 10 ws 3 1
24 rb 1 ws 5 1 24 rb 2 ws 2 1
24 rb 3 ws 3 1 24 rb 4 ws 4 1
24 rb 5 ws 1 1 24 rb 6 rb 8 1
24 rb 7 rb 6 0 24 rb 8 rb 9 1
24 rb 9 rb 10 1 24 rb 10 rb 7 0
24 ws 1 ws 8 1 24 ws 2 ws 10 0
24 ws 3 ws 6 1 24 ws 4 ws 9 1
24 ws 5 ws 7 0 24 ws 6 rb 1 0
24 ws 7 rb 5 1 24 ws 8 rb 3 0
24 ws 9 rb 4 0 24 ws 10 rb 2 0
;
proc glimmix data=salamander;
class fpop fnum mpop mnum;
model mating(event='1') = fpop|mpop / dist=binary;
random fpop*fnum mpop*mnum;
lsmeans fpop*mpop / ilink;
run;
ods graphics on;
ods select DiffPlot SliceDiffs;
proc glimmix data=salamander;
class fpop fnum mpop mnum;
model mating(event='1') = fpop|mpop / dist=binary;
random fpop*fnum mpop*mnum;
lsmeans fpop*mpop / plots=diffplot;
lsmeans fpop*mpop / slicediff=(mpop fpop);
run;
ods graphics off;