How to get fixed effects estimates glimmix poisson

[rfriesdo]

Hi there, I am a SAS novice running a mixed model using PROC glimmix, to analyze longitudinal count data (Poisson distribution, autoregressive covariance structure), and would like to produce fixed effect estimates for the continuous predictors in my model. I am working with the following code:

proc glimmix data = activedata;

class id time;

model DV = time var1 var2 var3/ dist=POISSON;

random _residual_;

random intercept var1 var2  / subject=id type=un;

random time / subject=id type=ar(1)  residual;

run;

This code runs without error and produces the expected "Type III tests of Fixed Effects" table (which gives F- and p-values), but I also need to estimate coefficients, standard errors, t- (or Wald), and p-values, as usually produced in the "Solutions for Fixed Effects" table. Any help for code that would help to produce these estimates (or why I will/should not be able to produce them) would be greatly appreciated!

[PGStats · Accepted solution]

Add option SOLUTION to the MODEL statement to get the estimates.

[rfriesdo]

PGSTATS, thank you for this answer! I had previously tried the “solution” option in the model statement but it did not work. I may have had an extra / in there. It seems to be working now though 🙂