Hi all,

This is my first post here so I apologise in advance if I've broken any rules, and I also apologise if this is a very simple question, as I am new to Stata.

To give some background, I am attempting to see whether or not alcohol consumption effects the level of income one earns. Observations have been split into 3 groups within one variable: Heavy drinkers, moderate drinkers and non-drinkers.

My plan is to first run a multinomial logit with the drinking variable as the dependent variable and moderate drinker as the base. I am doing this in order to obtain the mills ratios for each of the types of drinker, which will then be used in 3 OLS regressions with log income as the dependent variable in order to correct for selection bias.

My question is how do I obtain these Mills ratios?