We are utilizing a complex survey (NHATS) with stratification, clustering, and weighting. Our final models are Poisson with GEE. We understand that SURVEY procedures (SURVEYMEANS, SURVEYLOGISTIC, etc.) are necessary to use in order to appropriately account for this sample design and obtain accurate standard errors. We also know that PROC GENMOD does not have a direct corrollary with a SAS SURVEY statement.
We have had the suggestion to use SUDAAN, STATA, or some other program, but have used SAS for the rest of the analysis and would like to continue to do so. We found what appeared to be a very helpful article outlining with sample code how to use a non-survey procedure (PROC GENMOD) to obtain valid point estimates and extended programming to obtain valid standard errors for our model's parameter estimates (Poisson Regressions for Complex Surveys--attached here) and have followed the steps outlined therein carefully. Unfortunately Step 2 will not run due to the GEE portion of our model (if we take out the GEE coding lines starting with "repeated subjects" then the code runs, but with them included we get errors).
Here is one of our original models as an example (we did change to proc genmod data=jkweights appropriately when attempting to run the suggested code):
proc genmod data = analysis descending;
class SPID frailty time/param=ref ref=first;
weight weight;
model drive = frailty time frailty*time/type3 dist = poisson link = log;
repeated subject = SPID/ type=cs;
estimate "frailty t1" frailty 1 / exp;
estimate "frailty t2" frailty 1 frailty*time 1 0 0/ exp;
estimate "frailty t3" frailty 1 frailty*time 0 1 0/ exp;
estimate "frailty t4" frailty 1 frailty*time 0 0 1/ exp;
run;
In comparison, here is the Poisson model (without GEE) that is presented in the paper we are attaching:
proc genmod data=jkweights;
class sex race marital private education;
weight SamplingWeight;
model visits = sex race marital private education / dist=poisson;
run;
Is anybody aware of how to take the information presented in this article and apply it to a Poisson model including GEE? Any help is greatly appreciated.
