Hello, I am running a multilevel Poisson model with students within classes within schools, and I would like to estimate the sum (and its SE) of the fixed and random effects for exposure at the school level. PROC GLIMMIX data = dat7 noclprint ; class school class / ref=first ; model O = exposure / solution ; random exposure / subject=school solution cl ; random intercept / subject=class (school) ; estimate 'school 1' intercept 0 exposure 1 | intercept 0 exposure 1 / subject 1; estimate 'school 2' intercept 0 exposure 1 | intercept 0 exposure 1 / subject 0 1; estimate 'school 3' intercept 0 exposure 1 | intercept 0 exposure 1 / subject 0 0 1; estimate 'school 4' intercept 0 exposure 1 | intercept 0 exposure 1 / subject 0 0 0 1; estimate 'school 5' intercept 0 exposure 1 | intercept 0 exposure 1 / subject 0 0 0 0 1; run; I get the following error message : ERROR: Contrasts between subjects require that all random effects have the same subject; the CONTRAST statement is ignored. Glimmix does not allow the estimate statement when there are more than two subject levems. I tried by curiosity with Proc HPMIXED and it does work (but not with Proc Mixed). Is there a way to do this ? Thank you ! Alexis
... View more