🔒 This topic is solved and locked.
Need further help from the community? Please
sign in and ask a new question.
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Posted 05-10-2018 01:12 PM
(1777 views)
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.
1 ACCEPTED SOLUTION
Accepted Solutions
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Dear Rick!
Thank you much!!
Thank you much!!
2 REPLIES 2
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
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;
- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content
Dear Rick!
Thank you much!!
Thank you much!!