Dear All:
I want to plot observed vs predicted data, that estimated using PROC GEE , to check whether my model is
good fit of the data or not.
How can i access these predicted data from the estimated model using proc GEE? the code am using is given below.
Proc gee data=halu;
class blu;
model y/n=blu filk mak/dist=bin;
repeated subject=id/type=cs;
run;
Thanks alot.
use the OUTPUT statement to store the predicted probabilities in a data set. Since you are fitting a binomial model with event/Trials syntax, I assume you want to use the observed proportions as the "observed data"?
Proc gee data=halu;
class blu;
model y/n=blu filk mak/dist=bin;
repeated subject=id/type=cs;
output out=GEEOut predicted=Pred;
run;
data Score;
set GEEOut;
ObsProb = y / n; /* form observed proportions */
run;
proc sgplot data=Score;
scatter x=Pred y=ObsProb;
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.