Documentation Example 7 for PROC GLIMMIX
/****************************************************************/
/* S A S S A M P L E L I B R A R Y */
/* */
/* NAME: gmxex07 */
/* TITLE: Documentation Example 7 for PROC GLIMMIX */
/* Isotonic Contrasts for Ordered Alternatives */
/* PRODUCT: STAT */
/* SYSTEM: ALL */
/* KEYS: Generalized linear mixed models */
/* Ordered hypothesis */
/* LSMESTIMATE */
/* PROCS: GLIMMIX */
/* DATA: */
/* */
/* SUPPORT: Oliver Schabenberger */
/* REF: */
/* MISC: */
/****************************************************************/
data FerriteCores;
do Temp = 1 to 4;
do rep = 1 to 5; drop rep;
input MagneticForce @@;
output;
end;
end;
datalines;
10.8 9.9 10.7 10.4 9.7
10.7 10.6 11.0 10.8 10.9
11.9 11.2 11.0 11.1 11.3
11.4 10.7 10.9 11.3 11.7
;
proc glimmix data=FerriteCores;
class Temp;
model MagneticForce = Temp;
lsmestimate Temp
'avg(1:1)<avg(2:4)' -3 1 1 1 divisor=3,
'avg(1:2)<avg(3:4)' -1 -1 1 1 divisor=2,
'avg(1:3)<avg(4:4)' -1 -1 -1 3 divisor=3
/ adjust=simulate(seed=1) cl upper;
ods select LSMestimates;
run;