Part 1. Suppose a nominal variable with three categories (A,B,C). Proc Logistic, as does GLMMIX and others procs, solve that. However, the PI insists that in CATMOD an analysis can be constructed to compare A vs (B+C) and B vs C in a single proc statement. However, he doesn't recall how. Is this possible? Please, details.
Part 2. The covariate list will include continuous variables that, as I read the documentation, catmod doesn't work well with. Can the analysis be done in logistic or glmmix? Details please.
Part 3. The dataset is multilevel, occasions within persons, so if the analysis can done in sas, can random effects be included? Again, details please.
Thank you, Gene Maguin
1. This variable that has A B C as levels, is it a predictor variable or a response variable?
2. GLM or GLIMMIX can handle a mix of continuous and class predictor variables.
3. Both GLM and GLIMMIX allow random effects.
I think PROC LOGISTIC and PROC GLIMMIX allow three-level response variables and allows both continuous and categorical predictor variables.
I believe there are ways built into these procedures to compare A to B+C and then also B to C, but I have never done this.
The two logit response functions that you have defined are known as continuation ratio logits. The continuation ratio logit model is one type of model used for modeling an ordinal response and is mentioned in this note about the types of logistic models available in SAS. The model can be fit in PROC CATMOD, as you mentioned, using weighted least squares. A version of it can also be fit by maximum likelihood estimation in PROC LOGISTIC by fitting the two logits separately. Both are discussed and illustrated in this note. Since you have continuous predictors, the approach using separate binary logit models is probably the better approach. Since ordinary binary logistic models are used in that approach, you could instead use a PROC GLIMMIX to include random effects.
It should also be possible to fit the overall model in a single step by defining the model in PROC NLMIXED which also can incorporate random effects. This is probably the best and most general approach. However, an example of doing this is not available. You can see an example of how NLMIXED can similarly be used to fit a cumulative logit model in this note which shows the NLMIXED code in the Addendum at the end.
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
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.