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;
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
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.