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.
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;
Show your work/code please.
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.
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;
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.
Could you please explain what the /param=ref option is for?
Thanks!
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.
----
Not quite, it can still use that as the reference for Effect coding, but it won't be what you want necessarily.
Thank you.
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?
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...;
Join us for SAS Innovate April 16-19 at the Aria in Las Vegas. Bring the team and save big with our group pricing for a limited time only.
Pre-conference courses and tutorials are filling up fast and are always a sellout. Register today to reserve your seat.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.