Mixed logit

Regression analysis
Part of a series on Statistics

Models
Linear regression Simple regression Ordinary least squares Polynomial regression General linear model
Generalized linear model Discrete choice Logistic regression Multinomial logit Mixed logit Probit Multinomial probit Ordered logit Ordered probit Poisson
Multilevel model Fixed effects Random effects Mixed model
Nonlinear regression Nonparametric Semiparametric Robust Quantile Isotonic Principal components Least angle Local Segmented
Errors-in-variables
Estimation
Least squares Ordinary least squares Linear (math) Partial Total Generalized Weighted Non-linear Non-negative Iteratively reweighted Ridge regression
Least absolute deviations Bayesian Bayesian multivariate
Background
Regression model validation Mean and predicted response Errors and residuals Goodness of fit Studentized residual Gauss–Markov theorem
Statistics portal

Mixed logit is a fully general statistical model for examining discrete choices. The motivation for the mixed logit model arises from the limitations of the standard logit model. The standard logit model has three primary limitations, which mixed logit solves: "It obviates the three limitations of standard logit by allowing for random taste variation, unrestricted substitution patterns, and correlation in unobserved factors over time."^[1] Mixed logit can also utilize any distribution for the random coefficients, unlike probit which is limited to the normal distribution. It has been shown that a mixed logit model can approximate to any degree of accuracy any true random utility model of discrete choice, given an appropriate specification of variables and distribution of coefficients."^[2] The following discussion draws from Ch. 6 of Discrete Choice Methods with Simulation, by Kenneth Train (Cambridge University Press), to which the reader is referred for more details and citations. See also the article on discrete choice for information on how the mixed logit relates to discrete choice analysis in general and to other specific types of choice models.

Random taste variation

The standard logit model's "taste" coefficients, or $\beta$ 's, are fixed, which means the $\beta$ 's are the same for everyone. Mixed logit has different $\beta$ 's for each person (i.e., each decision maker.)

In the standard logit model, the utility of person n for alternative i is:

U_{{ni}}=\beta x_{{ni}}+\varepsilon _{{ni}}

with

\varepsilon_{ni}

~ iid extreme value

For the mixed logit model, this specification is generalized by allowing $\beta _{n}$ to be random. The utility of person n for alternative i in the mixed logit model is:

U_{{ni}}=\beta _{n}x_{{ni}}+\varepsilon _{{ni}}

with

\varepsilon_{ni}

~ iid extreme value

\quad \beta _{n}\sim f(\beta |\theta )

where θ are the parameters of the distribution of $\beta _{n}$ 's over the population, such as the mean and variance of $\beta _{n}$ .

Conditional on $\beta _{n}$ , the probability that person n chooses alternative i is the standard logit formula:

L_{{ni}}(\beta _{{n}})={\frac {e^{{\beta _{{n}}X_{{ni}}}}}{\sum _{{j}}e^{{\beta _{{n}}X_{{nj}}}}}}

However, since $\beta _{n}$ is random and not known, the (unconditional) choice probability is the integral of this logit formula over the density of $\beta _{n}$ .

P_{{ni}}=\int L_{{ni}}(\beta )f(\beta |\theta )d\beta

This model is also called the random coefficient logit model since $\beta _{n}$ is a random variable. It allows the slopes of utility (i.e., the marginal utility) to be random, which is an extension of the random effects model where only the intercept was stochastic.

Any probability density function can be specified for the distribution of the coefficients in the population, i.e., for $f(\beta |\theta )$ . The most widely used distribution is normal, mainly for its simplicity. For coefficients that take the same sign for all people, such as a price coefficient that is necessarily negative or the coefficient of a desirable attribute, distributions with support on only one side of zero, like the lognormal, are used.^[3]^[4] When coefficients cannot logically be unboundedly large or small, then bounded distributions are often used, such as the $S_{b}$ or triangular distributions.

Unrestricted substitution patterns

The mixed logit model can represent general substitution pattern because it does not exhibit logit's restrictive independence of irrelevant alternatives (IIA) property. The percentage change in the probability for one alternative given a percentage change in the mth attribute of another alternative is

E_{{nix_{{nj}}^{m}}}=-{\frac {x_{{nj}}^{m}}{P_{{ni}}}}\int \beta ^{m}L_{{ni}}(\beta )L_{{nj}}(\beta )f(\beta )d\beta =-x_{{nj}}^{m}\int \beta ^{m}L_{{nj}}(\beta ){\frac {L_{{ni}}(\beta )}{P_{{ni}}}}f(\beta )d\beta

where β ^m is the mth element of $\beta$ .^[1]^[4] It can be seen from this formula that "A ten-percent reduction for one alternative need not imply (as with logit) a ten-percent reduction in each other alternative."^[1] The relative percentages depend on correlation between the likelihood that respondent n will choose alternative i, L _ni, and the likelihood that respondent n will choose alternative j, L _nj, over various draws of β.

Correlation in unobserved factors over time

Standard logit does not take into account any unobserved factors that persist over time for a given decision maker. This can be a problem if you are using panel data, which represent repeated choices over time. By applying a standard logit model to panel data you are making the assumption that the unobserved factors that affect a person's choice are new every time the person makes the choice. That is a very unlikely assumption. To take into account both random taste variation and correlation in unobserved factors over time, the utility for respondent n for alternative i at time t is specified as follows:

U_{{nit}}=\beta _{{n}}X_{{nit}}+\varepsilon _{{nit}}

where the subscript t is the time dimension. We still make the logit assumption which is that $\varepsilon$ is i.i.d extreme value. That means that $\varepsilon$ is independent over time, people, and alternatives. $\varepsilon$ is essentially just white noise. However, correlation over time and over alternatives arises from the common effect of the $\beta$ 's, which enter utility in each time period and each alternative.

To examine the correlation explicitly, assume that the β 's are normally distributed with mean ${\bar {\beta }}$ and variance $\sigma ^{2}$ . Then the utility equation becomes:

U_{{nit}}=({\bar {\beta }}+\sigma \eta _{{n}})X_{{nit}}+\varepsilon _{{nit}}

and η is a draw from the standard normal density. Rearranging, the equation becomes:

U_{{nit}}={\bar {\beta }}X_{{nit}}+(\sigma \eta _{{n}}X_{{nit}}+\varepsilon _{{nit}})

U_{{nit}}={\bar {\beta }}X_{{nit}}+e_{{nit}}

where the unobserved factors are collected in $e_{{nit}}=\sigma \eta _{{n}}X_{{nit}}+\varepsilon _{{nit}}$ . Of the unobserved factors, $\varepsilon _{{nit}}$ is independent over time, and $\sigma \eta _{{n}}X_{{nit}}$ is not independent over time or alternatives.

Then the covariance between alternatives $i$ and $j$ is,

Cov(e_{{nit}},e_{{njt}})=\sigma ^{2}(X_{{nit}}X_{{njt}})

and the covariance between time $t$ and $q$ is

Cov(e_{{nit}},e_{{niq}})=\sigma ^{2}(X_{{nit}}X_{{niq}})

By specifying the X's appropriately, one can obtain any pattern of covariance over time and alternatives.

Conditional on $\beta _{n}$ , the probability of the sequence of choices by a person is simply the product of the logit probability of each individual choice by that person:

L_{{n}}(\beta _{{n}})=\prod _{{t}}{\frac {e^{{\beta _{{n}}X_{{nit}}}}}{\sum _{{j}}e^{{\beta _{{n}}X_{{njt}}}}}}

since $\varepsilon _{{nit}}$ is independent over time. Then the (unconditional) probability of the sequence of choices is simply the integral of this product of logits over the density of $\beta$ .

P_{{ni}}=\int L_{{n}}(\beta )f(\beta |\theta )d\beta

Simulation

Unfortunately there is no closed form for the integral that enters the choice probability, and so the researcher must simulate P_n. Fortunately for the researcher, simulating P_n can be very simple. There are four basic steps to follow

1. Take a draw from the probability density function that you specified for the 'taste' coefficients. That is, take a draw from $f(\beta |\theta )$ and label the draw $\beta ^{r}$ , for $r=1$ representing the first draw.

2. Calculate $L_{n}(\beta ^{r})$ . (The conditional probability.)

3. Repeat many times, for $r=2,...,R$ .

4. Average the results

Then the formula for the simulation look like the following,

${\tilde {P}}_{{ni}}={\frac {\sum _{{r}}L_{{n}}(\beta ^{r})}{R}}$

where R is the total number of draws taken from the distribution, and r is one draw.

Once this is done you will have a value for the probability of each alternative i for each respondent n.

References

1 2 3 Train, K. (2003) Discrete Choice Methods with Simulation
↑ McFadden, D. and Train, K. (2000). “Mixed MNL Models for Discrete Response,” Journal of Applied Econometrics, Vol. 15, No. 5, pp. 447-470,
↑ David Revelt and Train, K (1998). "Mixed Logit with Repeated Choices: Households' Choices of Appliance Efficiency Level," Review of Economics and Statistics, Vol. 80, No. 4, pp. 647-657
1 2 Train, K (1998)."Recreation Demand Models with Taste Variation," Land Economics, Vol. 74, No. 2, pp. 230-239.

This article is issued from Wikipedia - version of the 11/18/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.