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;
class ID;
model admit = post case post*case / dist=poisson link=log offset=logtime;
repeated subject=ID;
run;
Any suggestions would be appreciated!
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.
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.
Save $250 on SAS Innovate and get a free advance copy of the new SAS For Dummies book! Use the code "SASforDummies" to register. Don't miss out, May 6-9, in Orlando, Florida.
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.