Heteroscedastic Extreme-Value Model

The heteroscedastic extreme-value (HEV) model (Bhat, 1995) allows the random components of the utility function to be nonidentical. Specifically, the HEV model assumes independent but nonidentical error terms distributed with the Type I extreme-value distribution. The HEV model allows the variances of the random components of utility to differ across alternatives. Bhat (1995) argues that the HEV model does not have the IIA property. The HEV model contains the conditional logit model as a special case. The probability that an individual $i$ will choose alternative $j$ from the set $C_{i}$ of available alternatives is

\[  P_{i}(j) = \int _{-\infty }^{\infty }\prod _{k\in C_{i},k\ne j} \Gamma \left[\frac{\mathbf{x}_{ij}\bbeta - \mathbf{x}_{ik}\bbeta + \theta _{j}w}{\theta _{k}}\right]\gamma (w)dw  \]

where the choice set $C_{i}$ has $n_{i}$ elements and

\[  \Gamma (x) = \exp (-\exp (-x))  \]
\[  \gamma (x) = \exp (-x)\Gamma (x)  \]

are the cumulative distribution function and probability density function of the Type I extreme-value distribution. The variance of the error term for the $j$th alternative is $\frac{1}{6}\pi ^{2}\theta _{j}^{2}$. If the scale parameters, $\theta _{j}$, of the random components of utility of all alternatives are equal, then this choice probability is the same as that of the conditional logit model. The log-likelihood function of the HEV model can be written as

\[  \mathcal{L} = \sum _{i=1}^{N}\sum _{j \in C_ i}d_{ij}\ln [P_{i}(j)]  \]

where

\[  d_{ij} = \left\{  \begin{array}{cl} 1 &  \mr {if \;  individual} \;  i \;  \mr {chooses \;  alternative} \;  j \\ 0 &  \mr {otherwise} \end{array} \right.  \]

Since the log-likelihood function contains an improper integral function, it is computationally difficult to get a stable estimate. With the transformation $u=\exp (-w)$, the probability can be written

\begin{eqnarray*}  P_{i}(j) &  = &  \int _{0}^{\infty }\Pi _{k\in C_{i},k\ne j} \Gamma \left[\frac{\mathbf{x}_{ij}\bbeta - \mathbf{x}_{ik}\bbeta - \theta _{j}\ln (u)}{\theta _{k}}\right]\exp (-u)du \\ &  = &  \int _{0}^{\infty }G_{ij}(u)\exp (-u)du \end{eqnarray*}

Using the Gauss-Laguerre weight function, $W(x)=\exp (-u)$, the integration of the log-likelihood function can be replaced with a summation as follows:

\[  \int _{0}^{\infty }G_{ij}(u)\exp (-u)du = \sum _{k=1}^{K}w_{k}G_{ij}(x_{k})  \]