The GLIMMIX Procedure

Odds Ratios in Multinomial Models

The GLIMMIX procedure fits two kinds of models to multinomial data. Models with cumulative link functions apply to ordinal data, and generalized logit models are fit to nominal data. If you model a multinomial response with LINK=CUMLOGIT or LINK=GLOGIT, odds ratio results are available for these models.

In the generalized logit model, you model baseline category logits. By default, the GLIMMIX procedure chooses the last category as the reference category. If your nominal response has J categories, the baseline logit for category j is

\[ \log \left\{ \pi _ j / \pi _ J \right\} = \eta _ j = \mb{x}’\bbeta _ j + \mb{z}’\mb{u}_ j \]

and

\begin{align*} \pi _ j & = \frac{ \exp \{ \eta _ j \} }{\sum _{k=1}^ J \exp \{ \eta _ k\} } \\ \eta _ J & = 0 \end{align*}

As before, suppose that the two conditions to be compared are identified with subscripts 1 and 0. The log odds ratio of outcome j versus J for the two conditions is then

\begin{align*} \log \left\{ \psi \left(\eta _{j1},\eta _{j0}\right)\right\} & = \log \left\{ \frac{\pi _{j1} / \pi _{J1}}{ \pi _{j0}/ \pi _{J0} } \right\} = \log \left\{ \frac{ \exp \{ \eta _{j1}\} }{\exp \{ \eta _{j0}\} } \right\} \\ & = \eta _{j1} - \eta _{j0} \end{align*}

Note that the log odds ratios are again differences on the scale of the linear predictor, but they depend on the response category. The GLIMMIX procedure determines the estimable functions whose differences represent log odds ratios as discussed previously but produces separate estimates for each nonreference response category.

In models for ordinal data, PROC GLIMMIX models the logits of cumulative probabilities. Thus, the estimates on the linear scale represent log cumulative odds. The cumulative logits are formed as

\[ \log \left\{ \frac{\mr{Pr}(Y \leq j)}{\mr{Pr}(Y > j)} \right\} = \eta _ j = \alpha _ j + \mb{x}’\bbeta + \mb{z}’\bgamma = \alpha _ j + \tilde{\eta } \]

so that the linear predictor depends on the response category only through the intercepts (cutoffs) $\alpha _1, \ldots ,\alpha _{J-1}$. The odds ratio comparing two conditions represented by linear predictors $\eta _{j1}$ and $\eta _{j0}$ is then

\begin{align*} \psi \left(\eta _{j1},\eta _{j0}\right) & = \exp \left\{ \eta _{j1}-\eta _{j0} \right\} \\ & = \exp \left\{ \tilde{\eta }_1 - \tilde{\eta }_0 \right\} \end{align*}

and is independent of category.