Hello,
I initially built a logistic model for my dataset, and used the score option within proc logistic to get a probability of the event for each observation. I would like to use proc genmod (with binomial distribution, logit link) so I can use the bayes option with coeffprior= in order to get coefficients based on a categorical variable.
The issue - I don't know how to get the probability of the event for each observation within proc genmod. How can I do something like the score option within proc logistic?
Thanks!
Hi,
You need to use ods with obstats option something liket this.
ods output obstats=obstats;
PROC GENMOD DATA=have;
MODEL Y=X / D=B OBSTATS;
RUN;
proc print data=obstats;
run;
Thank you for your response. Gave that a try, but the proc print tells me there's no dataset called obstats. In the line: ods output obstats=obstats, the obstats does not become blue as a recognized option.
Are you putting obstats in model options?
MODEL Y=X / D=B OBSTATS;
Yes, I put that in. Actually, I just realized I got a note too
NOTE: The OBSTATS, RESIDUALS, XVARS, PREDICTED, and DIAGNOSTICS options are ignored for a Bayesian analysis.
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.