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
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;
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;
Don't miss out on SAS Innovate - Register now for the FREE Livestream!
Can't make it to Vegas? No problem! Watch our general sessions LIVE or on-demand starting April 17th. Hear from SAS execs, best-selling author Adam Grant, Hot Ones host Sean Evans, top tech journalist Kara Swisher, AI expert Cassie Kozyrkov, and the mind-blowing dance crew iLuminate! Plus, get access to over 20 breakout sessions.
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.