topic Re: Score a genmod model in SAS Procedures

Score a genmod model

BN_RN17 — Tue, 09 Sep 2014 15:04:47 GMT

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!

Re: Score a genmod model

stat_sas — Tue, 09 Sep 2014 15:26:49 GMT

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;

Re: Score a genmod model

BN_RN17 — Tue, 09 Sep 2014 17:25:24 GMT

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.

Re: Score a genmod model

stat_sas — Tue, 09 Sep 2014 17:34:28 GMT

Are you putting obstats in model options?

MODEL Y=X / D=B OBSTATS;

Re: Score a genmod model

BN_RN17 — Tue, 09 Sep 2014 17:58:20 GMT

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.