10-13-2014 08:42 AM
I have fitted a model to some data using a Poisson distribution (see below). From this, I have a dataset containing the standarised resisduals. What I would also like is the 2.5% and the 97.5% confidence band for each standarised resisdual. I see in the dataset from using 'obstats' that it gives the lower and upper bounds for each predicted value but I dont see an option for the standardised resisduals.
Is there such an output or do I need to use (and find) a formula to work this out?
proc genmod data=mydata;
model count=trtcd / offset=logt dist=poisson link=log scale=Pearson obstats;
lsmeans trtcd / diff CL;
estimate 'rate ratio A vs P' trtcd 1 -1 / exp;
ods output obstats=stdres#
Thank you for reading my question.