Storing estimates of proc Glimmix

[ammarhm]

Hi all,

I am using a proc glimmix to perform a poison regression, and using the "store" statement to store the output parameters

proc glimmix data=Poisson_sim;

class prog ;

model num_awards = prog math / dist=poisson  s cl;

store p1;

run;

Then using the following to estimate the IRR, and lower and upper CI

ods output ParameterEstimates = est;

proc plm source = p1;

  show parameters;

run;

data est_exp (drop=stderr llci luci estimate );

  set est;

  IRR = exp(estimate);

  llci=estimate - 1.96*stderr;

  luci=estimate + 1.96*stderr;

  UpperCI=exp (luci);

  LowerCI= exp (llci);

run;

proc print data = est_exp;

run;

The problem is: I am not managing to get the p-values from the proc Glimmix to be displayed in the output

Could anyone please advice me on how to do it?

Best wishes

[ammarhm · Accepted solution]

Solved!

proc glimmix data=Poisson_sim;

class prog ;

model num_awards = prog math / dist=poisson  s cl;

ods output parameterestimates = p1;

run;

data p2;

      set p1;

      IRR = exp(estimate);

      LL = exp(lower);

      UL = exp(upper);

proc print data = p2;

      var Effect OR LL UL probt;

run;