Storing estimates of proc Glimmix

ammarhm — Fri, 20 Feb 2015 03:28:03 GMT

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

Re: Storing estimates of proc Glimmix

ammarhm — Fri, 20 Feb 2015 04:22:23 GMT

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;