03-23-2016 05:44 PM
I am running a poisson regression model for the number of admissions for a group of patients (identified by "ID"). The explanatory variables are post (pre vs post period) and case (case vs control) with an interaction term post*case. As each patient can have more than one admission, I have added a repeated subject statement to the model. However, after adding the repeated subject statement to the code the output for the model no longer gives the deviance. How can I specify that the output for the model include the deviance?
proc genmod data=alldata;
model admit = post case post*case / dist=poisson link=log offset=logtime;
Any suggestions would be appreciated!
03-24-2016 04:41 PM
When you add a repated statement, you are perform GEE, which involves a quasi-likeliood (not a true likelihood or true distribution). Thus, the deviance (which is based on a true likelihood) is not printed.
03-29-2016 12:23 PM
Thank you for your reply. Do you know how I can run this type of model with the repeated subject statement and force SAS to output the deviance? Or do you have a suggestion of another model that would be better for this situation? Thank you.