How can access predicted data from Proc GEE

[Lijuu]

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.

[Rick_SAS]

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;

[Lijuu · Accepted solution]

Dear Rick!
Thank you much!!