Residuals in logistic regression

[Walternate]

Hello,

I have created a multiple logistic regression model and am trying to look at the residuals. The only process I have found (iplots) prints residuals for about 100 participants at a time, which is not ideal since I have over 5000 study subjects. Any help is much appreciated!

Thanks.

[Reeza]

Show your work/code please. 

[AncaTilea]

How about

ODS GRAPHICS ON;

proc logistic data = your_data plots = (all);

     model y = x;

run;

In my case, the plots are saved on the C:\Documents and Settings\User

Good luck,

Anca.

[Walternate · Accepted solution]

proc logistic data=mydata;

class var1(reference=1) var2(reference=1) var3(reference=1);

model case(event='1')=var1 var2 var3/rsq ctable pprob=(0 to 1 by .1) lackfit risklimits iplots;

run;

[Reeza]

I think Anca solution is what you're looking for. In addition, I think you need to add the /param=ref option to your class statement.

Take a look at your output, specifically the Class level information and see if its what you want (basically your design matrix) or if you want the option I suggested.

[Walternate]

Could you please explain what the /param=ref option is for?

Thanks!

[AncaTilea]

In your class statement you use 'reference = 1', which tells SAS that for the variable var1 the reference category is '1'.

You need the param = ref to tell SAS to use that reference coding in parameter estimates.

----

[Reeza]

Not quite, it can still use that as the reference for Effect coding, but it won't be what you want necessarily.

[AncaTilea]

Thank you.

[Reeza]

It specifies how SAS creates the design matrix for class variables, the default is 'EFFECT CODING' when most people want 'REFERENCE CODING' instead.

SAS FAQ: In PROC LOGISTIC why aren't the coefficients consistent with the odds ratios?

[AncaTilea]

When I tried to run plots with Logistic with a data set that has more than 5000 observations I get this warning:

-------------

WARNING: Some plots have more than 5000 observations and are suppressed.  Specify the

         PLOTS(MAXPOINTS=NONE) option in the PROC LOGISTIC statement to display the plots.

Did you try this?

proc logistic data = my_data plots(maxpoints = none);

     model...;