About Ksharp

Ksharp · ‎02-13-2026

@Rick_SAS wrote a blog about this: https://blogs.sas.com/content/iml/2014/02/19/scoring-a-regression-model-in-sas.html

Ksharp · ‎02-13-2026

https://support.sas.com/kb/22/912.html

Ksharp · ‎02-13-2026

Firstly using ODS TRACE to check whether the output object name contains "COEF" or not. You need to use "lsmeans" to get "COEF" object. Here is an example: ods trace on; proc logistic data=sashelp.heart(obs=1000); class sex bp_status/ param=glm; model status=sex bp_status weight height; lsmeans sex / e ilink; ods output coef=coeffs; store out=ques; run; ods trace off; And better to post your question at Stat Forum: https://communities.sas.com/t5/Statistical-Procedures/bd-p/statistical_procedures

Ksharp · ‎02-12-2026

You mean @Rick_SAS or @StatDave ?

Ksharp · ‎02-12-2026

https://blogs.sas.com/content/graphicallyspeaking/2015/11/30/7180/ https://blogs.sas.com/content/graphicallyspeaking/2015/09/21/response-colors-and-thickness/

Ksharp · ‎02-12-2026

If I am using @StatDave 's code, I could get the same OR. Check this: 1) Main Effect --> Same OR data have; set sashelp.heart(obs=1000 keep=status sex ageatstart rename=(ageatstart=age sex=gender)); run; data test; gender='Male';age=65;output; gender='Male';age=75;output; run; data mydata; set have test indsname=indsname; dsn=indsname; run; proc logistic data=mydata; class gender / param=ref; effect spla=spline(age/naturalcubic basis=tpf(noint)); model status= gender spla ; oddsratio age / at(age=65); units age=10; /* estimate 'For M: 45 vs 50(ref)' gender 2 splage[1,45][-1,50] gender*splage [1,2 45] [-1,2 50] / E EXP CL;*/ output out=out xbeta=logodds; run; proc transpose data=out(where=(dsn='WORK.TEST')) out=logodds; var logodds; run; proc report data=logodds nowd split='*'; column col1 col2 or; define col1/display 'log(odds) for age=65'; define col2/display 'log(odds) for age=75'; define or/computed 'when gender="Male"*OR=odds(age=75)/odds(age=65)'; compute or; or=exp(col2-col1); endcomp; run; 2) Interact Effect --> Also same OR proc logistic data=mydata; class gender / param=ref; effect spla=spline(age/naturalcubic basis=tpf(noint)); model status= gender spla gender*spla; oddsratio age / at(age=65); units age=10; /* estimate 'For M: 45 vs 50(ref)' gender 2 splage[1,45][-1,50] gender*splage [1,2 45] [-1,2 50] / E EXP CL;*/ output out=out xbeta=logodds; run;

Ksharp · ‎02-11-2026

I think that is due to treatnp=3 appeared before treatnp=2 in the last visit. Here is an example: data maze3 ; input avisitn avisit $ treatnp chg; cards; 1 visit1 2 23 1 visit1 2 23 1 visit1 2 23 1 visit1 3 123 1 visit1 3 123 1 visit1 3 123 2 visit2 3 23 2 visit2 3 23 2 visit2 3 23 2 visit2 2 123 2 visit2 2 123 2 visit2 2 123 ; proc npar1way data = maze3 wilcoxon HL ALPHA = 0.05; by avisitn avisit; class treatnp; var chg; EXACT WILCOXON; ods output wilcoxontest =WT23 HodgesLehmann = HL23; run; If you changed the order of 2 and 3 in the last visit, you got this: data maze3 ; input avisitn avisit $ treatnp chg; cards; 1 visit1 2 23 1 visit1 2 23 1 visit1 2 23 1 visit1 3 123 1 visit1 3 123 1 visit1 3 123 2 visit2 2 23 2 visit2 2 23 2 visit2 2 23 2 visit2 3 123 2 visit2 3 123 2 visit2 3 123 ;

Ksharp · ‎02-11-2026

Thanks

Ksharp · ‎02-10-2026

If I understood @StatDave 's code ,you can report OR at age=65 to age=75 . a.k.a OR= odds(age=65)/odds(age=75) data have; set sashelp.heart(obs=1000 keep=status sex ageatstart rename=(ageatstart=age sex=gender)); run; data test; gender='Male';age=65;output; gender='Male';age=75;output; run; data mydata; set have test indsname=indsname; dsn=indsname; run; proc gampl data=mydata; class gender / param=ref; id dsn; model status= param(gender) spline(age) / dist=binary link=logit; output out=out xbeta=logodds; run; proc transpose data=out(where=(dsn='WORK.TEST')) out=logodds; run; proc report data=logodds nowd split='*'; column col1 col2 or; define col1/display 'log(odds) for age=65'; define col2/display 'log(odds) for age=75'; define or/computed 'OR=odds(age=65)/odds(age=75)'; compute or; or=exp(col1-col2); endcomp; run; NOTE: as @StatDave indicated the Male and Female have the same OR since model did not include interact effect “gender*age” . You can test it with the following data. Thanks Dave. data test; /*gender='Male ';age=65;output;*/ /*gender='Male ';age=75;output;*/ gender='Female';age=65;output; gender='Female';age=75;output; run;

Ksharp · ‎02-07-2026

OK. This one ? /*Example 2*/ ods listing gpath="%sysfunc(pathname(work))" image_dpi=300 style=htmlblue; ods graphics / reset=all imagename='sepallength' noborder width=77px height=77px; ods html exclude sgplot; /*generate a histogram graph*/ title; proc sgplot data=sashelp.iris pad=0 noborder noautolegend; histogram sepallength ; xaxis display=none; yaxis display=none; run; ods listing gpath="%sysfunc(pathname(work))" image_dpi=300 style=htmlblue; ods graphics / reset=all imagename='petallength' noborder width=77px height=77px; ods html exclude sgplot; /*generate a histogram graph*/ title; proc sgplot data=sashelp.iris pad=0 noborder noautolegend; histogram petallength ; xaxis display=none; yaxis display=none; run; %sganno data sganno; %SGIMAGE(IMAGE="%sysfunc(pathname(work))\sepallength.png",ANCHOR="topleft",BORDER="FALSE",DRAWSPACE="LAYOUTPERCENT" ,x1=5,y1=94) %SGIMAGE(IMAGE="%sysfunc(pathname(work))\petallength.png",ANCHOR="topleft",BORDER="FALSE",DRAWSPACE="LAYOUTPERCENT" ,x1=51,y1=48) /*get rid of the tick markers around histogram*/ %SGRECTANGLE(X1=10, Y1=91.9, HEIGHT=2, WIDTH=40, ANCHOR="TOPLEFT" , DISPLAY="FILL" , /*DISPLAY="ALL"*/ DRAWSPACE="GRAPHPERCENT", FILLCOLOR="white" , HEIGHTUNIT="PERCENT" , WIDTHUNIT="PERCENT" ) %SGRECTANGLE(X1=8.1, Y1=90, HEIGHT=40, WIDTH=2, ANCHOR="TOPLEFT" , DISPLAY="FILL" , /*DISPLAY="ALL"*/ DRAWSPACE="GRAPHPERCENT", FILLCOLOR="white" , HEIGHTUNIT="PERCENT" , WIDTHUNIT="PERCENT" ) %SGRECTANGLE(X1=50, Y1=10.2, HEIGHT=2, WIDTH=40, ANCHOR="TOPLEFT" , DISPLAY="FILL" , /*DISPLAY="ALL"*/ DRAWSPACE="GRAPHPERCENT", FILLCOLOR="white" , HEIGHTUNIT="PERCENT" , WIDTHUNIT="PERCENT" ) %SGRECTANGLE(X1=89.8, Y1=50, HEIGHT=40, WIDTH=2, ANCHOR="TOPLEFT" , DISPLAY="FILL" , /*DISPLAY="ALL"*/ DRAWSPACE="GRAPHPERCENT", FILLCOLOR="white" , HEIGHTUNIT="PERCENT" , WIDTHUNIT="PERCENT" ) /*Add X an Y variable label*/ %SGTEXT( LABEL="Sepal Length (mm)", WIDTH=40, WIDTHUNIT= "PERCENT" , ANCHOR="CENTER" , BORDER="FALSE", DRAWSPACE="GRAPHPERCENT", TEXTSIZE=12, X1=30, Y1=95 ) %SGTEXT( LABEL="Petal Length (mm)", WIDTH=40, WIDTHUNIT= "PERCENT" , ANCHOR="CENTER" , BORDER="FALSE", DRAWSPACE="GRAPHPERCENT", TEXTSIZE=12, X1=70, Y1=95 ) %SGTEXT( ROTATE=90, LABEL="Sepal Length (mm)", WIDTH=40, WIDTHUNIT= "PERCENT" , ANCHOR="CENTER" , BORDER="FALSE", DRAWSPACE="GRAPHPERCENT", TEXTSIZE=12, X1=4, Y1=75 ) %SGTEXT( ROTATE=90, LABEL="Petal Length (mm)", WIDTH=40, WIDTHUNIT= "PERCENT" , ANCHOR="CENTER" , BORDER="FALSE", DRAWSPACE="GRAPHPERCENT", TEXTSIZE=12, X1=4, Y1=25 ) run; ods graphics/reset; title; proc sgscatter data=sashelp.iris sganno=sganno pad=40px; matrix sepallength petallength; run;

Ksharp · ‎02-07-2026

It is a nonparametric regression unlike logistic model. Maybe @Rick_SAS could give you a hand .

Ksharp · ‎02-06-2026

You want this ? /*Example 1*/ ods listing gpath="%sysfunc(pathname(work))" image_dpi=300 style=htmlblue; ods graphics / reset=all imagename='sepallength' noborder width=80px height=80px; ods html exclude sgplot; /*generate a histogram graph*/ title; proc sgplot data=sashelp.iris pad=0 noborder noautolegend; histogram sepallength ; xaxis display=none; yaxis display=none; run; ods listing gpath="%sysfunc(pathname(work))" image_dpi=300 style=htmlblue; ods graphics / reset=all imagename='petallength' noborder width=80px height=80px; ods html exclude sgplot; /*generate a histogram graph*/ title; proc sgplot data=sashelp.iris pad=0 noborder noautolegend; histogram petallength ; xaxis display=none; yaxis display=none; run; %sganno data sganno; %SGIMAGE(IMAGE="%sysfunc(pathname(work))\sepallength.png",ANCHOR="topleft",BORDER="FALSE",DRAWSPACE="LAYOUTPERCENT" ,x1=6,y1=93) %SGIMAGE(IMAGE="%sysfunc(pathname(work))\petallength.png",ANCHOR="topleft",BORDER="FALSE",DRAWSPACE="LAYOUTPERCENT" ,x1=54,y1=48) run; ods graphics/reset; proc sgscatter data=sashelp.iris sganno=sganno; title "Scatterplot Matrix for Iris Data"; matrix sepallength petallength; run;

Ksharp · ‎02-06-2026

That was supposed to be if you are using "proportional" option . You can not apply OFFSETMIN/MAX to the different the panel with using PROC SGPANEL. 1) You can plot these bar separated one by one and combine them together: https://blogs.sas.com/content/graphicallyspeaking/2022/09/10/complex-layouts-using-the-sg-procedures/ 2)You can change the bar width to make it look better: title1 "Plot with different x-axes across panels using PROPORTIONAL"; title2 "Oops! The bars are not squeezed together!"; proc sgpanel data=plotdata; where not (panel = 1 and group = 'C'); panelby panel / rows=1 onepanel proportional uniscale=row ; vbarparm response=y category=group /barwidth=0.4; run; 3) Or make some dummy group values as quickbluefish did: data plotdata; call streaminit(123); do panel = 1, 2; group='_';y=0;output; do group = 'A', 'B', 'C'; y = rand("Uniform", 0, 100); output; end; group='-';y=0;output; end; run; title1 "Plot with different x-axes across panels using PROPORTIONAL"; title2 "Oops! The bars are not squeezed together!"; proc sgpanel data=plotdata; panelby panel / rows=1 onepanel proportional uniscale=row ; vbarparm response=y category=group /grouporder=data nooutline; colaxis values=('_' 'A' 'B' 'C' '-') valuesdisplay=(' ' 'A' 'B' 'C' ' '); run;

Ksharp · ‎02-06-2026

Since you are asking OR for a SPLINE effect(a.k.a a non-linear effect), that is not easy to get it. I think you need help from macro %margin , %nlmeans: https://support.sas.com/kb/35/189.html

Ksharp · ‎02-06-2026

Or using STRIP() function to get rid of these heading/trailing blanks . data have; string='A '; run; proc sql; /*select * from sashelp.class where name contains 'A';*/ select * from sashelp.class where name contains strip((select string from have)); quit;

Online Status	Offline
Date Last Visited	yesterday