I am using Glimmix to run a mixed model with a Poisson distribution. Here's the basic model layout: proc glimmix noclprint=12 method=laplace; class site year sht tb poscat efnbin; model esm=efnbin tb efnbin*tb c_shtefn c_shtefn*tb c_shtefn*efnbin/dist=poisson link=log solution; random year site year*site sht*year*site; I want to estimate means for the "efnbin" effect with within-subject standard errors. The problem is specifying this in the ESTIMATE statement because there are over 800 levels of the "sht*year*site" random effect. So, my question is if there is a short-hand way to specify the same coefficient for all the >800 levels of "sht*year*site". Thanks.
... View more