Dear Reader, In my model, where treatment is fixed and block is random, there seems to be a different treatment effect for a different set of blocks. Attached you see code where I estimate treatment effects per set of random blocks, and code where I estimate the difference for a set of blocks. The estimates do NOT coincide, meaning that the difference per treatment does not flow from the estimates of a treatment. Any help on explaining this difference would be much appreciated. proc glimmix data=test_Per plots=all;
class TT block;
model Per=TT / s cl ddfm=kenwardroger;
random block / s cl g gcorr;
estimate 'trt=1 @ blocks 1-6' intercept 6 TT 6 0 | block 1 1 1 1 1 1 0 0 0 0 0 0 / divisor=6;
estimate 'trt=2 @ blocks 1-6' intercept 6 TT 0 6 | block 1 1 1 1 1 1 0 0 0 0 0 0 / divisor=6;
estimate 'trt=1 vs trt=2 @ blocks 1-6' TT -6 6 0 0 0 0 0 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=3 @ blocks 1-6' TT -6 0 6 0 0 0 0 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=4 @ blocks 1-6' TT -6 0 0 6 0 0 0 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=5 @ blocks 1-6' TT -6 0 0 0 6 0 0 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=6 @ blocks 1-6' TT -6 0 0 0 0 6 0 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=7 @ blocks 1-6' TT -6 0 0 0 0 0 6 0 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=8 @ blocks 1-6' TT -6 0 0 0 0 0 0 6 0 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=9 @ blocks 1-6' TT -6 0 0 0 0 0 0 0 6 0 | block 1 1 1 1 1 1 0 0 0 0 0 0,
'trt=1 vs trt=10 @ blocks 1-6' TT -6 0 0 0 0 0 0 0 0 6 | block 1 1 1 1 1 1 0 0 0 0 0 0 / divisor=6 e adjdfe=row adjust=simulate(nsamp=100000);
estimate 'trt=1 vs trt=2 @ blocks 7-12' TT -6 6 0 0 0 0 0 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=3 @ blocks 7-12' TT -6 0 6 0 0 0 0 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=4 @ blocks 7-12' TT -6 0 0 6 0 0 0 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=5 @ blocks 7-12' TT -6 0 0 0 6 0 0 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=6 @ blocks 7-12' TT -6 0 0 0 0 6 0 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=7 @ blocks 7-12' TT -6 0 0 0 0 0 6 0 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=8 @ blocks 7-12' TT -6 0 0 0 0 0 0 6 0 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=9 @ blocks 7-12' TT -6 0 0 0 0 0 0 0 6 0 | block 0 0 0 0 0 0 1 1 1 1 1 1,
'trt=1 vs trt=10 @ blocks 7-12' TT -6 0 0 0 0 0 0 0 0 6 | block 0 0 0 0 0 0 1 1 1 1 1 1 / divisor=6 e adjdfe=row adjust=simulate(nsamp=100000);
store glimmix_Per;
run; Kind regards, Gavin
... View more