And after some google later, I found the link = log is a model called "log-binomial regression" (https://stats.stackexchange.com/questions/581678/what-is-the-difference-between-a-logistic-regression-and-log-binomial-regression )   The thread says, Logistic regression always produces estimates of 𝑝𝑖 that are between 0 and 1. This is not true for log-binomial regression; I guess this is the reason why few people use the log-binomial regression to model their data because the risk of deriving a risk higher than 100%. Thanks. ... View more

Exactly! This is the issue. I specified serumALB instead of SplALB in the MODEL statement. After correcting this error, everything followed as expected. Thank you. ... View more

Hello, Dave.   Wow! That's a great reference. Thanks. I copied and edited the code to suit my project. But, so far, I have never get the Odd Ratio plot as shown in the first example. My code only generated one column for the variable of serumAlb. I am just not able to make it. Your investigation would be of great important. Thanks.   My code: %macro Logistic2; proc logistic data= &raw_data_name plots = oddsratio(type=vertical); class &class_var; effect SplALB=spline(serumAlb/naturalcubic); model &censor(event = '1') = &indep_list; effectplot fit(x=serumAlb); effectplot fit(x=serumAlb) / link; oddsratio serumAlb / at(serumAlb=34.8125 36.0425 37.35 39.0 40.92); store kmod; run; %mend; %Logistic2 ... View more

Thanks for these advice, ballardw.   Yep, I will follow your guide. ... View more

I also want to learn if there any statement to generate spline for odd ratios in SAS? ... View more

Hmmm, thanks. The EFFECTPLOT statement is powerful. However, I don't understand one thing.   proc logistic data=d202ps_3; class obesity_status(REF='0') genderC(REF='0') AF(REF='0') Hypertension(REF='0') CHF(REF='0') akiC(ref='0') ddiC(REF='0') indicationC(REF='0') CYP2C9code(REF='0') VKORC1Code(REF='0'); model firstInrAttain(event='1') = obesity_status Weight__kg_ ageY genderC AF Hypertension CHF serumAlb akiC ddiC indicationC CYP2C9code VKORC1Code / clodds=both clparm=both; effectplot slicefit(x=ageY sliceby=obesity_status=0 1 2) / clm noobs; run; and proc logistic data=d202ps_3; class /*obesity_status(REF='0')*/ genderC(REF='0') AF(REF='0') Hypertension(REF='0') CHF(REF='0') akiC(ref='0') ddiC(REF='0') indicationC(REF='0') CYP2C9code(REF='0') VKORC1Code(REF='0'); model firstInrAttain(event='1') = obesity_status Weight__kg_ ageY genderC AF Hypertension CHF serumAlb akiC ddiC indicationC CYP2C9code VKORC1Code / clodds=both clparm=both; effectplot slicefit(x=ageY sliceby=obesity_status=0 1 2) / clm noobs; run; The two produced obviously different graphs. The only difference in the code is whether obesity_status was listed as class/category variable. Also, if I treat obesity_status as a class variable, I cannot define its levels in the EFFECTPLOT analysis (i.e., it has 4 levels but I prefer to use only three levels.) Thanks. ... View more

Hello, everyone   Today I learned the PROC TRANSREG could be used to draw spline, e.g.,   ods graphics on; proc transreg data=d202ps_4; model identity(addMg) = spline(serumAlb); run; However, my linear regression model should include many other independent variables. I tried to add them to the PROC TRANSREG but it did not work as I expected. Does anyone know how to solve this issue (spline with other variables adjusted)?   Thanks ... View more

Oh, great. Let met study it. Thank you, sir.   Hmmm, it seems that the example in the thread used simple poisson regression with only one independent variable. But in reality, we often have multiple independent variables.   Tom ... View more