proc plm effectplot contour: how to generate plot for one value of the outcome only?

[ChristosK]

proc logistic data=HEERDT.VPICU_label;
 model VP_ICU1_VasoNE_012(ref='0')= Ind_DBP1|Age/ nooddsratio link=glogit;
  store logimodel;
  run; 

proc plm source=logiModel ;
   effectplot contour(x= Ind_DBP1 y=Age plotby=VP_ICU1_VasoNE_012 );
run;

I'm generating the following image, but just want to get a publication quality image for the value of VP_ICU1_VasoNE_012 =2.

Since this is the outcome, I cant find an option using plotby = to just generate this image.

Sounds like a question for Rick:-)

[Rick_SAS · Accepted solution]

I can't open the pdf for some reason. However, I think the values in the EFFECTPLOT statement are restricted to the explanatory variables. It sounds like you want to plot the probability for a particular level of a multinomial response variable. If so, use the UNPACK option to put each response level in its own graph.  You might also want to use the ILINK option to get the contours on the probability scale. Here's an example:

proc logistic data=sashelp.cars;
  model origin(ref="Asia")= mpg_city|weight / nooddsratio link=glogit;
  store logimodel;
run; 

proc plm source=logiModel ;
   effectplot contour(x=mpg_city y=Weight) / ilink unpack;
run;

[ChristosK]

Thanks Rick.

Not quite sure I understand what the link function does. I tried it with and without, but the output looked the same.

This is the image I came up with: