BookmarkSubscribeRSS Feed
🔒 This topic is solved and locked. Need further help from the community? Please sign in and ask a new question.
amng
Calcite | Level 5

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_0-1606600475086.png

 

1 ACCEPTED SOLUTION

Accepted Solutions
FreelanceReinh
Jade | Level 19

@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;

 

amng_0-1606600475086.png

 


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.

View solution in original post

6 REPLIES 6
PaigeMiller
Diamond | Level 26

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.

--
Paige Miller
amng
Calcite | Level 5

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.

amng
Calcite | Level 5
Sorry - the dataset is too large to be posted. What is the best way to upload it in that case?
PaigeMiller
Diamond | Level 26

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.

--
Paige Miller
FreelanceReinh
Jade | Level 19

@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;

 

amng_0-1606600475086.png

 


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.

sas-innovate-2024.png

Available on demand!

Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.

 

Register now!

What is ANOVA?

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.

Discussion stats
  • 6 replies
  • 630 views
  • 1 like
  • 3 in conversation