05-10-2018 01:12 PM
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;
model y/n=blu filk mak/dist=bin;
05-11-2018 01:44 PM
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;
Need further help from the community? Please ask a new question.