Follow up re #2:   I would like to do the use the "ODS output partialplotdata=PD" with PROC LOGISTIC, since in this case too I have multiple predictors, but in the ODS table for PROC LOGISTIC I cannot find an equivalent. Is there a way to do this?   My ultimate goal is to be able to produce a graph in which I can represent, using different colours, the partial effect of the same predictor on two different DVs.    Specifically, I am testing two PROC LOGISTIC:   proc logistic data=covid ; model close2 (event='1') = policy close1 identif ctrlvar2; proc logistic data=covid ; model distant2 (event='1') = policy distant1 identif ctrlvar2;   Using the "plots=effect (clband=yes)" I can obtain the graph attached.  What I would like to do, however, is to produce one single graph displaying the the partial effect of policy on close2 and distant2 using, say the color RED for close2 and the color BLUE for distant2.    Am I asking too much?       ... View more

Thank you! ... View more

Thank you! both are good. I may take advantage of your expertise for one more step. The design is in fact more complex, with a nested factor that is random. This is the syntax:    proc glm data=ebal_st plots=(meanplot(cl ascending)); class genre author; model lessico = genre author(genre) /effectsize; output out=GLMOut Predicted=Pred lclm=lclm uclm=uclm; lsmeans genre / cl; random author(genre); ods output LSMeanCL=LS; run; quit; run; quit;   SAS automatically outputs an additional graph, to plot the random factor, AUTHOR, see below. I am not sure why for some of the levels of the random factor AUTHOR it puts what seem to be error bars. I cannot find the documentation where this graph is explained. Any idea? The graph is below.       ... View more

Thank you Rick. For a simple plot of the two predictors on the same graph, the solution with SGPLOT reported above works. Of course, as noted, this is not the same as reporting the results of a multiple regression, for which I will be reading the documents you shared. Eman  ... View more

Thank you! ... View more

Dear all, I found out that the EFFECTSIZE option in the MODEL of PROC GLM outputs, for each predictor, Semi-partial eta-square, which are the same as squared semi-partial correlations (sr2). If that's the case, I can simply compute the root-square of those to obtain semi-partial correlations. Am I missing something? Eman  ... View more

That is perfect! And using sgplot indeed gives me a chance to make a better graph. Thank you very much, really appreciated. Eman ... View more