In many applications of discrete choice models, modelers are interested in estimating consumer’s marginal “willingness-to-pay” (WTP) for different attributes. WTP can be estimated in two ways:

  1. Estimate a discrete choice model in the “preference space” where parameters have units of utility and then compute the WTP by dividing the parameters by the price parameter.
  2. Estimate a discrete choice model in the “WTP space” where parameters have units of WTP.

While the two procedures generally produce the same estimates of WTP for homogenous models, the same is not true for heterogeneous models where model parameters are assumed to follow a specific distribution, such as normal or log-normal (Train and Weeks 2005). For example, in a preference space specification, a normally distributed attribute parameter divided by a log-normally distributed price parameter produces a strange WTP distribution with large tails. In contrast, a WTP space specification allows the modeler to directly assume WTP is normally distributed. The package was developed to enable modelers to choose between these two utility spaces when estimating multinomial logit models.

The random utility model in two spaces

The random utility model is a well-established framework in many fields for estimating consumer preferences from observed consumer choices (Louviere, Hensher, and Swait 2000; Train 2009). Random utility models assume that consumers choose the alternative jj a set of alternatives that has the greatest utility uju_{j}. Utility is a random variable that is modeled as uj=vj+εju_{j} = v_{j} + \varepsilon_{j}, where vjv_{j} is the “observed utility” (a function of the observed attributes such that vj=f(𝐱j)v_{j} = f(\mathrm{\mathbf{x}}_{j})) and εj\varepsilon_{j} is a random variable representing the portion of utility unobservable to the modeler.

Adopting the same notation as in Helveston et al. (2018), consider the following utility model:

uj*=𝛃*𝐱j+α*pj+εj*,εj*Gumbel(0,σ2π26)\begin{equation} u^*_{j} = \mathbf{β}^{*'} \mathrm{\mathbf{x}}_{j} + \alpha^{*} p_{j} + \varepsilon^{*}_{j}, \quad\quad \varepsilon^{*}_{j} \sim \textrm{Gumbel}\left(0, \sigma^2\frac{\pi^2}{6}\right) \label{eq:utility} \end{equation}

where 𝛃*\mathbf{β}^{*} is the vector of coefficients for non-price attributes 𝐱j\mathrm{\mathbf{x}}_{j}, α*\alpha^{*} is the coefficient for price pjp_{j}, and the error term, εj*\varepsilon^{*}_{j}, is an IID random variable with a Gumbel extreme value distribution of mean zero and variance σ2(π2/6)\sigma^2(\pi^2/6). This model is not identified since there exists an infinite set of combinations of values for 𝛃*\mathbf{β}^{*}, α*\alpha^{*}, and σ\sigma that produce the same choice probabilities. In order to specify an identifiable model, the modeler must normalize equation (). One approach is to normalize the scale of the error term by dividing equation () by σ\sigma, producing the “preference space” utility specification:

(uj*σ)=(𝛃*σ)𝐱j+(α*σ)pj+(εj*σ),(εj*σ)Gumbel(0,π26)\begin{equation} \left(\frac{u^*_{j}}{\sigma}\right) = \left( \frac{\mathbf{β}^{*}}{\sigma} \right)' \mathrm{\mathbf{x}}_{j} + \left( \frac{\alpha^{*}}{\sigma} \right) p_{j} + \left( \frac{\varepsilon^{*}_{j}}{\sigma} \right), \quad\quad \left( \frac{\varepsilon^{*}_{j}}{\sigma} \right) \sim \textrm{Gumbel}\left(0, \frac{\pi^2}{6}\right) \label{eq:utilityPreferenceScaled} \end{equation}

The typical preference space parameterization of the multinomial logit (MNL) model can then be written by rewriting equation () with uj=(uj*/σ)u_j = (u^*_j / \sigma), 𝛃=(𝛃*/σ)\mathbf{β}= (\mathbf{β}^{*} / \sigma), α=(α*/σ)\alpha = (\alpha^{*} / \sigma), and εj=(εj*/σ)\varepsilon_{j} = (\varepsilon^{*}_{j} / \sigma):

$$\begin{equation} u_{j} = \mathbf{β}' \mathrm{\mathbf{x}}_{j} + \alpha p_{j} + \varepsilon_{j} \hspace{0.5in} \varepsilon_{j} \sim \textrm{Gumbel}\left(0,\frac{\pi^2}{6}\right) \label{eq:utilityPreference} \end{equation}$$

The vector 𝛃\mathbf{β} represents the marginal utility for changes in each non-price attribute, and α\alpha represents the marginal utility obtained from price reductions. In addition, the coefficients 𝛃\mathbf{β} and α\alpha are measured in units of utility, which only has relative rather than absolute meaning.

The alternative normalization approach is to normalize equation () by α*- \alpha^* instead of σ\sigma, producing the “willingness-to-pay (WTP) space” utility specification:

(uj*α*)=(𝛃*α*)𝐱j+(α*α*)pj+(εj*α*),(εj*α*)Gumbel(0,σ2(α*)2π26)\begin{equation} \left(\frac{u^*_{j}}{- \alpha^*}\right) = \left(\frac{\mathbf{β}^{*}}{- \alpha^{*}}\right)' \mathrm{\mathbf{x}}_{j} + \left(\frac{\alpha^{*}}{- \alpha^{*}}\right) p_{j} + \left(\frac{\varepsilon^{*}_{j}}{- \alpha^{*}}\right), \quad\quad \left(\frac{\varepsilon^{*}_{j}}{- \alpha^{*}}\right) \sim \textrm{Gumbel} \left(0, \frac{\sigma^2}{(- \alpha^{*})^2}\frac{\pi^2}{6} \right) \label{eq:utilityWtpScaled} \end{equation}

Since the error term in equation is scaled by λ2=σ2/(α*)2\lambda^2 = \sigma^2/(- \alpha^{*})^2, we can rewrite equation () by multiplying both sides by λ=(α*/σ\lambda= (- \alpha^{*} / \sigma) and renaming uj=(λuj*/α*)u_j = (\lambda u^*_j /- \alpha^*), 𝛚=(𝛃*/α*\mathbf{ω}= (\mathbf{β}^{*} /- \alpha^{*}), and εj=(λεj*/α*)\varepsilon_j = (\lambda \varepsilon^*_j / - \alpha^*):

$$\begin{equation} u_{j} = \lambda \left( \mathbf{ω}' \mathrm{\mathbf{x}}_{j} - p_{j} \right) + \varepsilon_{j} \hspace{0.5in} \varepsilon_{j} \sim \textrm{Gumbel}\left(0, \frac{\pi^2}{6}\right) \label{eq:utilityWtp} \end{equation}$$

Here 𝛚\mathbf{ω} represents the marginal WTP for changes in each non-price attribute, and λ\lambda represents the scale of the deterministic portion of utility relative to the standardized scale of the random error term.

The utility models in equations and represent the preference space and WTP space utility specifications, respectively. In equation , WTP is estimated as 𝛃̂/α̂\hat{\mathbf{β}} / \hat{- \alpha}; in equation , WTP is simply 𝛚̂\hat{\mathbf{ω}}.

References

Helveston, John Paul, Elea McDonnell Feit, and Jeremy J. Michalek. 2018. Pooling stated and revealed preference data in the presence of RP endogeneity.” Transportation Research Part B: Methodological 109: 70–89.
Louviere, Jordan J., David A. Hensher, and Joffre Swait. 2000. Stated Choice Methods: Analysis and Applications. Edited by J. J. et al. Louviere. Cambridge, Massachusetts: Cambridge University Press.
Train, Kenneth E. 2009. Discrete Choice Methods with Simulation. 2nd ed. Cambridge University Press.
Train, Kenneth E., and Melvyn Weeks. 2005. Discrete Choice Models in Preference and Willingness-to-Pay Space.” In Appl. Simul. Methods Environ. Resour. Econ., 1–16.