Multinomial Ordinal Models
Multinomial Ordinal Models
DUSTIN C. BROWN
POPULATION RESEARCH CENTER
Objectives
Examples:
Smoking status – never, current, former smoker
Marital status – married, divorced, widowed, never
married
Multinomial Logit (Probit) Model
The data are from the NHIS Adult Sample Files (2009)
Predictors:
Education: <High School, High School, Some College, College (Ref.)
Race-Ethnicity: NH White (Ref.), NH Black, Hispanic,
Age in years
Stata Code: “mlogit smk3 lths hs scol nhb hispanic age, base(0) rrr”
Examples:
Self-Rated Health – poor, fair, good, very good, excellent
Ordered Logit (Probit) Models
The data are from the NHIS Adult Sample Files (2009)
Predictors:
Education: <High School, High School, Some College, College (Ref.)
Race-Ethnicity: NH White (Ref.), NH Black, Hispanic,
Age in years
Stata Code: “ologit health lths hs scol nhb hispanic age, or”
The odds of reporting poor vs. fair, good, very good, and excellent
health are 3.97 times greater for persons who did not graduate high
school in comparison to persons with a college degree net of race-
ethnicity and age.
Once again, the ordered logit (probit) model assumes that the
distance between each category of the outcome is proportional.
A Brant test can be used to test whether the proportional odds (i.e.,
parallel lines) assumption holds.
This is available as a user-added post-estimation command in Stata.
The Brant test indicates that the influence of education and race-ethnicity
are not proportional across each category of self-rated health. Note, that
the association between age and self-rated health is proportional though.
When the Proportionality Assumption is Violated…
Option 2: Use a multinomial logit model. This frees you of the proportionality
assumption, but it is less parsimonious and often dubious on substantive grounds.
Option 3: Dichotomize the outcome and use binary logistic regression. This is
common, but you lose information and it could alter your substantive conclusions.
Option 4: Use a model that does not assume proportionality. Increasingly, this is
common. Two user-submitted Stata commands fit these kinds of models:
“gologit2” – generalized ordered logit models (see Williams 2007, Stata Jn.)
“oglm” – heterogeneous choice models (see Williams 2010, Stata Jn.)
Recommendation: Try all the above and decide what to do based on your results.