2 weeks ago
I am trying to run multiple comparison testing to see whether the value of a categorical variable education (4 levels) differs significantly depending on which one of 5 groups a person is in. This is weighted survey data so I am trying to test these multiple comparisons by running a multinomial logistic regression using education level as the dependent variable and using the 5 group variable as the independent variable, then using lsmeans to check for significant pairwise differences using Tukey procedure.
My syntax is below:
PROC SURVEYLOGISTIC DATA=data1;
CLASS EDUCATION_LEVEL GROUP / PARAM=GLM;
MODEL EDUCATION_LEVEL = GROUP / LINK=GLOGIT;
LSMEANS GROUP / ADJUST=TUKEY PDIFF=ALL;
However, the differences table in my output doesn't show the differences between groups for the reference category of education. How do I change the code to show whether the reference group of education varies by the 5 groups?
Should I be using LSESTIMATE and if so, how do I specify that? I'm using SAS version 9.4
Thanks in advance!