Hi all,
For a class project I am modeling the impact of a binary health literacy variable (getadvice2) on the likelihood of completing up to 5 self-management behaviors for diabetic patients (manage5). I originally tried to run an ordinal logistic regression, but following the violation of the proportional odds assumption I instead modeled a multinomial regression. However, the odds ratios under "point estimate" are not showing up properly. I wanted to ask whether this is some kind of SAS Error? There is nothing out of the ordinary in the log.
PROC LOGISTIC DATA=diabetes;
class getadvice2(ref=first) manage5(ref="0") / param=ref;
MODEL manage5(Desc) = getadvice2 / link=glogit;
title "Model 1";
RUN;
@amng wrote:
PROC LOGISTIC DATA=diabetes;
class getadvice2(ref=first) manage5(ref="0") / param=ref;
MODEL manage5(Desc) = getadvice2 / link=glogit;
title "Model 1";
RUN;
Hi @amng,
Look at a cross tabulation of getadvice2 and manage5 (PROC FREQ). If the cell getadvice2=0 & manage5=0 is empty, the estimated odds of manage5=5 (and likewise 4, 3, 2, 1) vs. manage5=0, which you chose as the reference category, will have a zero in the denominator for getadvice2=0, leading to the "zero" odds ratio estimates in your table. In this case, pick a different reference category for manage5, i.e., one with non-empty cells for both getadvice2 categories, so that non-trivial odds ratio estimates (for manage5 ne 0) have a chance to occur.
When a SAS user says that there is an error in SAS, I generally believe SAS, and that the error is the user's fault, not SAS's fault.
But unless you can share your data with us, there's no way we can resolve this.
Does this help? https://support.sas.com/kb/22/954.html
I also think you should not use your response variable MANAGE5 in the CLASS statement, but I don't know if that makes any difference here.
Hi thanks. I wasn't trying to suggest there was an error in SAS itself - just trying to understand why it was giving me these results. The data is attached here.
We need you to show us a portion of the data (not the whole thing) by providing the data in this format and not any other format.
@amng wrote:
PROC LOGISTIC DATA=diabetes;
class getadvice2(ref=first) manage5(ref="0") / param=ref;
MODEL manage5(Desc) = getadvice2 / link=glogit;
title "Model 1";
RUN;
Hi @amng,
Look at a cross tabulation of getadvice2 and manage5 (PROC FREQ). If the cell getadvice2=0 & manage5=0 is empty, the estimated odds of manage5=5 (and likewise 4, 3, 2, 1) vs. manage5=0, which you chose as the reference category, will have a zero in the denominator for getadvice2=0, leading to the "zero" odds ratio estimates in your table. In this case, pick a different reference category for manage5, i.e., one with non-empty cells for both getadvice2 categories, so that non-trivial odds ratio estimates (for manage5 ne 0) have a chance to occur.
Available on demand!
Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.