Discrete choice

In economics, discrete choice models, or qualitative choice models, describe, explain, and predict choices between two or more discrete alternatives, such as entering or not entering the labor market, or choosing between modes of transport. Such choices contrast with standard consumption models in which the quantity of each good consumed is assumed to be a continuous variable. In the continuous case, calculus methods (e.g. first-order conditions) can be used to determine the optimum amount chosen, and demand can be modeled empirically using regression analysis. On the other hand, discrete choice analysis examines situations in which the potential outcomes are discrete, such that the optimum is not characterized by standard first-order conditions. Thus, instead of examining "how much" as in problems with continuous choice variables, discrete choice analysis examines "which one". However, discrete choice analysis can also be used to examine the chosen quantity when only a few distinct quantities must be chosen from, such as the number of vehicles a household chooses to own and the number of minutes of telecommunications service a customer decides to purchase. Techniques such as logistic regression and probit regression can be used for empirical analysis of discrete choice.

Discrete choice models theoretically or empirically model choices made by people among a finite set of alternatives. The models have been used to examine, e.g., the choice of which car to buy, where to go to college, which mode of transport (car, bus, rail) to take to work among numerous other applications. Discrete choice models are also used to examine choices by organizations, such as firms or government agencies. In the discussion below, the decision-making unit is assumed to be a person, though the concepts are applicable more generally. Daniel McFadden won the Nobel prize in 2000 for his pioneering work in developing the theoretical basis for discrete choice.

Discrete choice models statistically relate the choice made by each person to the attributes of the person and the attributes of the alternatives available to the person. For example, the choice of which car a person buys is statistically related to the person's income and age as well as to price, fuel efficiency, size, and other attributes of each available car. The models estimate the probability that a person chooses a particular alternative. The models are often used to forecast how people's choices will change under changes in demographics and/or attributes of the alternatives.

Discrete choice models specify the probability that an individual chooses an option among a set of alternatives. The probabilistic description of discrete choice behavior is used not to reflect individual behavior that is viewed as intrinsically probabilistic. Rather, it is the lack of information that leads us to describe choice in a probabilistic fashion. In practice, we cannot know all factors affecting individual choice decisions as their determinants are partially observed or imperfectly measured. Therefore, discrete choice models rely on stochastic assumptions and specifications to account for unobserved factors related to a) choice alternatives, b) taste variation over people (interpersonal heterogeneity) and over time (intra-individual choice dynamics), and c) heterogeneous choice sets. The different formulations have been summarized and classified into groups of models. When discrete choice model are combined with structural equation models to integrate psychological (latent) variables, they are referred as hybrid choice models.

Applications

 * Marketing researchers use discrete choice models to study consumer demand and to predict competitive business responses, enabling choice modelers to solve a range of business problems, such as pricing, product development, and demand estimation problems. In market research, this is commonly called conjoint analysis.
 * Transportation planners use discrete choice models to predict demand for planned transportation systems, such as which route a driver will take and whether someone will take rapid transit systems. The first applications of discrete choice models were in transportation planning, and much of the most advanced research in discrete choice models is conducted by transportation researchers.
 * Disaster planners and engineers rely on discrete choice models to predict decision take by householders or building occupants in small-scale and large-scales evacuations, such as building fires, wildfires, hurricanes among others.   These models help in the development of reliable disaster managing plans and safer design for the built environment.
 * Energy forecasters and policymakers use discrete choice models for households' and firms' choice of heating system, appliance efficiency levels, and fuel efficiency level of vehicles.
 * Environmental studies utilize discrete choice models to examine the recreators' choice of, e.g., fishing or skiing site and to infer the value of amenities, such as campgrounds, fish stock, and warming huts, and to estimate the value of water quality improvements.
 * Labor economists use discrete choice models to examine participation in the work force, occupation choice, and choice of college and training programs.
 * Ecological studies employ discrete choice models to investigate parameters that drive habitat selection in animals.

Common features of discrete choice models
Discrete choice models take many forms, including: Binary Logit, Binary Probit, Multinomial Logit, Conditional Logit, Multinomial Probit, Nested Logit, Generalized Extreme Value Models, Mixed Logit, and Exploded Logit. All of these models have the features described below in common.

Choice set
The choice set is the set of alternatives that are available to the person. For a discrete choice model, the choice set must meet three requirements:


 * 1) The set of alternatives must be collectively exhaustive, meaning that the set includes all possible alternatives. This requirement implies that the person necessarily does choose an alternative from the set.
 * 2) The alternatives must be mutually exclusive, meaning that choosing one alternative means not choosing any other alternatives. This requirement implies that the person chooses only one alternative from the set.
 * 3) The set must contain a finite number of alternatives. This third requirement distinguishes discrete choice analysis from forms of regression analysis in which the dependent variable can (theoretically) take an infinite number of values.

As an example, the choice set for a person deciding which mode of transport to take to work includes driving alone, carpooling, taking bus, etc. The choice set is complicated by the fact that a person can use multiple modes for a given trip, such as driving a car to a train station and then taking train to work. In this case, the choice set can include each possible combination of modes. Alternatively, the choice can be defined as the choice of "primary" mode, with the set consisting of car, bus, rail, and other (e.g. walking, bicycles, etc.). Note that the alternative "other" is included in order to make the choice set exhaustive.

Different people may have different choice sets, depending on their circumstances. For instance, the Scion automobile was not sold in Canada as of 2009, so new car buyers in Canada faced different choice sets from those of American consumers. Such considerations are taken into account in the formulation of discrete choice models.

Defining choice probabilities
A discrete choice model specifies the probability that a person chooses a particular alternative, with the probability expressed as a function of observed variables that relate to the alternatives and the person. In its general form, the probability that person n chooses alternative i is expressed as:


 * $$ P_{ni} \equiv \Pr( \text{Person } n \text{ chooses alternative } i)  = G(x_{ni}, \;x_{nj, j \neq i}, \; s_n, \;\beta), $$

where


 * $$x_{ni} $$ is a vector of attributes of alternative i faced by person n,


 * $$x_{nj, j \neq i} $$ is a vector of attributes of the other alternatives (other than i) faced by person n,


 * $$ s_n $$ is a vector of characteristics of person n, and


 * $$ \beta $$ is a set of parameters giving the effects of variables on probabilities, which are estimated statistically.

In the mode of transport example above, the attributes of modes (xni), such as travel time and cost, and the characteristics of consumer (sn), such as annual income, age, and gender, can be used to calculate choice probabilities. The attributes of the alternatives can differ over people; e.g., cost and time for travel to work by car, bus, and rail are different for each person depending on the location of home and work of that person.

Properties:
 * Pni is between 0 and 1
 * $$\forall n:\; \sum_{j=1}^J P_{nj} = 1, $$  where J is the total number of alternatives.
 * (Expected fraction of people choosing i ) $$ = {1 \over N} {\sum_{n=1}^N P_{ni}}, $$ where N is the number of people making the choice.

Different models (i.e., models using a different function G) have different properties. Prominent models are introduced below.

Consumer utility
Discrete choice models can be derived from utility theory. This derivation is useful for three reasons:


 * 1) It gives a precise meaning to the probabilities Pni
 * 2) It motivates and distinguishes alternative model specifications, e.g., the choice of a functional form for G.
 * 3) It provides the theoretical basis for calculation of changes in consumer surplus (compensating variation) from changes in the attributes of the alternatives.

Uni is the utility (or net benefit or well-being) that person n obtains from choosing alternative i. The behavior of the person is utility-maximizing: person n chooses the alternative that provides the highest utility. The choice of the person is designated by dummy variables, yni, for each alternative:


 * $$ y_{ni} = \begin{cases}

1 & U_{ni} > U_{nj} \quad  \forall j \neq i\\ 0 & \text{otherwise}\end{cases}$$

Consider now the researcher who is examining the choice. The person's choice depends on many factors, some of which the researcher observes and some of which the researcher does not. The utility that the person obtains from choosing an alternative is decomposed into a part that depends on variables that the researcher observes and a part that depends on variables that the researcher does not observe. In a linear form, this decomposition is expressed as
 * $$ U_{ni}= \beta z_{ni} + \varepsilon_{ni} $$

where


 * $$z_{ni} $$ is a vector of observed variables relating to alternative i for person n that depends on attributes of the alternative, xni,  interacted perhaps with attributes of the person, sn, such that it can be expressed as $$z_{ni}=z(x_{ni}, s_n) $$  for some numerical function z,
 * $$\beta $$ is a corresponding vector of coefficients of the observed variables, and
 * $$ \varepsilon_{ni} $$ captures the impact of all unobserved factors that affect the person's choice.

The choice probability is then


 * $$\begin{align}

P_{ni}& = \Pr(y_{ni} = 1) \\ &= \Pr \left ( \bigcap_{j \neq i}U_{ni} > U_{nj}, \right ) \\ & = \Pr \left ( \bigcap_{j \neq i} \beta z_{ni} + \varepsilon_{ni} >  \beta z_{nj} + \varepsilon_{nj}, \right ) \\ & = \Pr \left ( \bigcap_{j \neq i} \varepsilon_{nj}- \varepsilon_{ni} < \beta z_{ni}- \beta z_{nj}, \right ) \end{align}$$

Given β, the choice probability is the probability that the random terms, εnj − εni (which are random from the researcher's perspective, since the researcher does not observe them) are below the respective quantities $$\forall j \neq i: \beta z_{ni} - \beta z_{nj}.$$ Different choice models (i.e. different specifications of G) arise from different distributions of εni for all i and different treatments of β.

Only differences matter
The probability that a person chooses a particular alternative is determined by comparing the utility of choosing that alternative to the utility of choosing other alternatives:


 * $$ P_{ni}= \Pr(y_{ni} = 1) = \Pr \left ( \bigcap_{j \neq i} U_{ni} > U_{nj} \right ) = \Pr  \left ( \bigcap_{j \neq i} U_{ni} - U_{nj} > 0 \right )$$

As the last term indicates, the choice probability depends only on the difference in utilities between alternatives, not on the absolute level of utilities. Equivalently, adding a constant to the utilities of all the alternatives does not change the choice probabilities.

Scale must be normalized
Since utility has no units, it is necessary to normalize the scale of utilities. The scale of utility is often defined by the variance of the error term in discrete choice models. This variance may differ depending on the characteristics of the dataset, such as when or where the data are collected. Normalization of the variance therefore affects the interpretation of parameters estimated across diverse datasets.

Prominent types of discrete choice models
Discrete choice models can first be classified according to the number of available alternatives.
 * * Binomial choice models (dichotomous): 2 available alternatives
 * * Multinomial choice models (polytomous): 3 or more available alternatives

Multinomial choice models can further be classified according to the model specification:
 * * Models, such as standard logit, that assume no correlation in unobserved factors over alternatives
 * * Models that allow correlation in unobserved factors among alternatives

In addition, specific forms of the models are available for examining rankings of alternatives (i.e., first choice, second choice, third choice, etc.) and for ratings data.

Details for each model are provided in the following sections.

A. Logit with attributes of the person but no attributes of the alternatives
Un is the utility (or net benefit) that person n obtains from taking an action (as opposed to not taking the action). The utility the person obtains from taking the action depends on the characteristics of the person, some of which are observed by the researcher and some are not. The person takes the action,, if Un > 0. The unobserved term, εn, is assumed to have a logistic distribution. The specification is written succinctly as:


 * $$\begin{cases} U_n = \beta s_n + \varepsilon_n  \\ y_n = \begin{cases} 1 & U_n > 0 \\ 0 & U_n \leqslant 0 \end{cases} \\ \varepsilon \sim \text{Logistic} \end{cases} \quad \Rightarrow \quad  P_{n1} = \frac{1}{1+\exp(-\beta s_n)} $$

B. Probit with attributes of the person but no attributes of the alternatives
The description of the model is the same as model A, except the unobserved terms are distributed standard normal instead of logistic.


 * $$\begin{cases} U_n = \beta s_n + \varepsilon_n \\  y_n = \begin{cases} 1 & U_n > 0 \\ 0 & U_n \leqslant 0 \end{cases} \\ \varepsilon \sim \text{Standard normal} \end{cases} \quad \Rightarrow \quad P_{n1}  = \Phi(\beta s_n),$$

where $$\Phi$$ is cumulative distribution function of standard normal.

C. Logit with variables that vary over alternatives
Uni is the utility person n obtains from choosing alternative i. The utility of each alternative depends on the attributes of the alternatives interacted perhaps with the attributes of the person. The unobserved terms are assumed to have an extreme value distribution.


 * $$ \begin{cases} U_{n1} = \beta z_{n1} + \varepsilon_{n1} \\ U_{n2} = \beta z_{n2} + \varepsilon_{n2} \\ \varepsilon_{n1}, \varepsilon_{n2} \sim \text{iid extreme value} \end{cases} \quad \Rightarrow \quad P_{n1}= \frac{\exp(\beta z_{n1})} {\exp(\beta z_{n1})+\exp(\beta z_{n2})} $$

We can relate this specification to model  A  above, which is also binary logit. In particular, Pn1 can also be expressed as
 * $$P_{n1} =\frac{1}{1+\exp(-\beta (z_{n1}-z_{n2}))}$$

Note that if two error terms are iid extreme value, their difference is distributed logistic, which is the basis for the equivalence of the two specifications.

D. Probit with variables that vary over alternatives
The description of the model is the same as model C, except the difference of the two unobserved terms are distributed standard normal instead of logistic.

Then the probability of taking the action is
 * $$P_{n1} = \Phi(\beta (z_{n1}-z_{n2})),$$

where Φ is the cumulative distribution function of standard normal.

E. Logit with attributes of the person but no attributes of the alternatives
The utility for all alternatives depends on the same variables, sn, but the coefficients are different for different alternatives:

The choice probability takes the form
 * Since only differences in utility matter, it is necessary to normalize $$\beta_i =0 $$ for one alternative. Assuming $$ \beta_1=0 $$,
 * εni are iid extreme value
 * εni are iid extreme value
 * $$P_{ni}= {\exp(\beta_i s_n) \over \sum_{j=1}^J \exp(\beta_j s_n)},$$

where J is the total number of alternatives.

F. Logit with variables that vary over alternatives (also called conditional logit)
The utility for each alternative depends on attributes of that alternative, interacted perhaps with attributes of the person:


 * $$\begin{cases} U_{ni} = \beta z_{ni} + \varepsilon_{ni} \\ \varepsilon_{ni} \sim \text{iid extreme value} \end{cases} \quad \Rightarrow \quad P_{ni} = {\exp(\beta z_{ni}) \over \sum_{j=1}^J \exp(\beta z_{nj})},$$

where J is the total number of alternatives.

Note that model E can be expressed in the same form as model F by appropriate respecification of variables. Define $$w_{nj}^k=s_n \delta_{jk} $$ where $$\delta_{jk}$$ is the Kronecker delta and sn are from model E. Then, model F is obtained by using


 * $$z_{nj} = \left \{ w^1_{nj}, \cdots, w^J_{nj} \right \}  \quad \text{and} \quad \beta = \left \{ \beta_1, \cdots,\beta_J \right \},$$

where J is the total number of alternatives.

Multinomial choice with correlation among alternatives
A standard logit model is not always suitable, since it assumes that there is no correlation in unobserved factors over alternatives. This lack of correlation translates into a particular pattern of substitution among alternatives that might not always be realistic in a given situation. This pattern of substitution is often called the Independence of Irrelevant Alternatives (IIA) property of standard logit models. A number of models have been proposed to allow correlation over alternatives and more general substitution patterns:


 * Nested Logit Model - Captures correlations between alternatives by partitioning the choice set into 'nests'
 * Cross-nested Logit model (CNL) - Alternatives may belong to more than one nest
 * C-logit Model - Captures correlations between alternatives using 'commonality factor'
 * Paired Combinatorial Logit Model - Suitable for route choice problems.
 * Generalized Extreme Value Model - General class of model, derived from the random utility model to which multinomial logit and nested logit belong
 * Conditional probit - Allows full covariance among alternatives using a joint normal distribution.
 * Mixed logit  - Allows any form of correlation and substitution patterns. When a mixed logit is with jointly normal random terms, the models is sometimes called "multinomial probit model with logit kernel". Can be applied to route choice.

The following sections describe Nested Logit, GEV, Probit, and Mixed Logit models in detail.

G. Nested Logit and Generalized Extreme Value (GEV) models
The model is the same as model F except that the unobserved component of utility is correlated over alternatives rather than being independent over alternatives.


 * The marginal distribution of each εni is extreme value, but their joint distribution allows correlation among them.
 * The probability takes many forms depending on the pattern of correlation that is specified. See Generalized Extreme Value.
 * The probability takes many forms depending on the pattern of correlation that is specified. See Generalized Extreme Value.

H. Multinomial probit
The model is the same as model G except that the unobserved terms are distributed jointly normal, which allows any pattern of correlation and heteroscedasticity:


 * $$\begin{cases} U_{ni} = \beta z_{ni} +\varepsilon_{ni} \\ \varepsilon_n \equiv (\varepsilon_{n1},\cdots,\varepsilon_{nJ}) \sim N(0,\Omega) \end{cases} \quad \Rightarrow \quad P_{ni} = \Pr \left ( \bigcap_{j \neq i}\beta z_{ni}+\varepsilon_{ni} > \beta z_{nj} + \varepsilon_{nj} \right )  = \int I\left ( \bigcap_{j \neq i}\beta z_{ni}+\varepsilon_{ni} > \beta z_{nj} + \varepsilon_{nj} \right ) \phi(\varepsilon_n | \Omega) \;d \varepsilon_n,$$

where $$ \phi(\varepsilon_n | \Omega) $$ is the joint normal density with mean zero and covariance $$\Omega $$.

The integral for this choice probability does not have a closed form, and so the probability is approximated by quadrature or simulation.

When $$\Omega $$ is the identity matrix (such that there is no correlation or heteroscedasticity), the model is called independent probit.

I. Mixed logit
Mixed Logit models have become increasingly popular in recent years for several reasons. First, the model allows $$\beta$$ to be random in addition to $$\varepsilon$$. The randomness in $$\beta$$ accommodates random taste variation over people and correlation across alternatives that generates flexible substitution patterns. Second, advances in simulation have made approximation of the model fairly easy. In addition, McFadden and Train have shown that any true choice model can be approximated, to any degree of accuracy by a mixed logit with appropriate specification of explanatory variables and distribution of coefficients.

The choice probability is
 * $$\beta \sim f(\beta | \theta) $$ for any distribution $$ \it f $$, where $$\theta$$ is the set of distribution parameters (e.g. mean and variance) to be estimated,
 * εni ~ iid extreme value,
 * εni ~ iid extreme value,
 * $$P_{ni}= \int_\beta L_{ni} (\beta) f(\beta | \theta) \, d\beta,$$

where
 * $$ L_{ni} (\beta) = {\exp(\beta z_{ni}) \over {\sum_{j=1}^J \exp(\beta z_{nj})}}$$

is logit probability evaluated at $$ \beta, $$ with $$ J $$ the total number of alternatives.

The integral for this choice probability does not have a closed form, so the probability is approximated by simulation.

Estimation from choices
Discrete choice models are often estimated using maximum likelihood estimation. Logit models can be estimated by logistic regression, and probit models can be estimated by probit regression. Nonparametric methods, such as the maximum score estimator, have been proposed. Estimation of such models is usually done via parametric, semi-parametric and non-parametric maximum likelihood methods, but can also be done with the Partial least squares path modeling approach.

Estimation from rankings
In many situations, a person's ranking of alternatives is observed, rather than just their chosen alternative. For example, a person who has bought a new car might be asked what he/she would have bought if that car was not offered, which provides information on the person's second choice in addition to their first choice. Or, in a survey, a respondent might be asked:


 * Example : Rank the following cell phone calling plans from your most preferred to your least preferred.
 * * $60 per month for unlimited anytime minutes, two-year contract with $100 early termination fee
 * * $30 per month for 400 anytime minutes, 3 cents per minute after 400 minutes, one-year contract with $125 early termination fee
 * * $35 per month for 500 anytime minutes, 3 cents per minute after 500 minutes, no contract or early termination fee
 * * $50 per month for 1000 anytime minutes, 5 cents per minute after 1000 minutes, two-year contract with $75 early termination fee

The models described above can be adapted to account for rankings beyond the first choice. The most prominent model for rankings data is the exploded logit and its mixed version.

J. Exploded logit
Under the same assumptions as for a standard logit (model F), the probability for a ranking of the alternatives is a product of standard logits. The model is called "exploded logit" because the choice situation that is usually represented as one logit formula for the chosen alternative is expanded ("exploded") to have a separate logit formula for each ranked alternative. The exploded logit model is the product of standard logit models with the choice set decreasing as each alternative is ranked and leaves the set of available choices in the subsequent choice.

Without loss of generality, the alternatives can be relabeled to represent the person's ranking, such that alternative 1 is the first choice, 2 the second choice, etc. The choice probability of ranking J alternatives as 1, 2, ..., J is then


 * $$\Pr(\text{ranking } 1, 2, \ldots, J) = {\exp(\beta z_1) \over \sum_{j=1}^J \exp(\beta z_{nj})} {\exp(\beta z_2) \over \sum_{j=2}^J \exp(\beta z_{nj})} \ldots {\exp(\beta z_{J-1}) \over \sum_{j=J-1}^J \exp(\beta z_{nj})}$$

As with standard logit, the exploded logit model assumes no correlation in unobserved factors over alternatives. The exploded logit can be generalized, in the same way as the standard logit is generalized, to accommodate correlations among alternatives and random taste variation. The "mixed exploded logit" model is obtained by probability of the ranking, given above, for Lni in the mixed logit model (model I).

This model is also known in econometrics as the rank ordered logit model and it was introduced in that field by Beggs, Cardell and Hausman in 1981. One application is the Combes et al. paper explaining the ranking of candidates to become professor. It is also known as Plackett–Luce model in biomedical literature.

Ordered models
In surveys, respondents are often asked to give ratings, such as:


 * Example : Please give your rating of how well the President is doing.
 * 1: Very badly
 * 2: Badly
 * 3: Okay
 * 4: Well
 * 5: Very well

Or,
 * Example : On a 1-5 scale where 1 means disagree completely and 5 means agree completely, how much do you agree with the following statement. "The Federal government should do more to help people facing foreclosure on their homes."

A multinomial discrete-choice model can examine the responses to these questions (model G, model H, model I). However, these models are derived under the concept that the respondent obtains some utility for each possible answer and gives the answer that provides the greatest utility. It might be more natural to think that the respondent has some latent measure or index associated with the question and answers in response to how high this measure is. Ordered logit and ordered probit models are derived under this concept.

K. Ordered logit
Let Un represent the strength of survey respondent n's feelings or opinion on the survey subject. Assume that there are cutoffs of the level of the opinion in choosing particular response. For instance, in the example of the helping people facing foreclosure, the person chooses for some real numbers a, b, c, d.
 * 1, if Un < a
 * 2, if a < Un < b
 * 3, if b < Un < c
 * 4, if c < Un < d
 * 5, if Un > d,

Defining $$ U_n = \beta z_n + \varepsilon, \; \varepsilon \sim$$ Logistic, then the probability of each possible response is:


 * $$\begin{align}

\Pr(\text{choosing } 1) & = \Pr(U_n   d) = \Pr(\varepsilon  >  d - \beta z_n) = 1 - {1 \over 1+\exp(-(d - \beta z_n))} \end{align}$$

The parameters of the model are the coefficients β and the cut-off points a − d, one of which must be normalized for identification. When there are only two possible responses, the ordered logit is the same a binary logit (model A), with one cut-off point normalized to zero.

L. Ordered probit
The description of the model is the same as model K, except the unobserved terms have normal distribution instead of logistic.

The choice probabilities are ($$\Phi$$ is the cumulative distribution function of the standard normal distribution):


 * $$\begin{align}

\Pr(\text{choosing } 1) &= \Phi (a - \beta z_n) \\ \Pr(\text{choosing } 2) &= \Phi (b - \beta z_n) - \Phi (a - \beta z_n) \\ &\cdots \end{align}$$