Dear Dr. Rick, Thank you for your code and clear explanation. I am trying to use your code in PROC LOGISTIC, but I don't know how to add other variables to adjusted (like gender, education...) . In the last example, we can used ADDINPUTVARS in GLMSELECT and output the SPL_ variables to PROC REG, but I can't find the similar option in PROC LOGISTIC statement (I need to add other variables). Is there any other way to solve it? I am appreciated for your reply!