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 
will choose alternative 
from the set 
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 \]](images/etsug_mdc0092.png){kind=small}
where the choice set has 
elements and
![\[ \Gamma (x) = \exp (-\exp (-x)) \]](images/etsug_mdc0093.png){kind=small}
![\[ \gamma (x) = \exp (-x)\Gamma (x) \]](images/etsug_mdc0094.png){kind=small}
are the cumulative distribution function and probability density function of the Type I extreme-value distribution. The variance
of the error term for the th alternative is 
. If the scale parameters, 
, 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)] \]](images/etsug_mdc0095.png){kind=small}
where
![\[ d_{ij} = \left\{ \begin{array}{cl} 1 & \mr {if \; individual} \; i \; \mr {chooses \; alternative} \; j \\ 0 & \mr {otherwise} \end{array} \right. \]](images/etsug_mdc0084.png){kind=small}
Since the log-likelihood function contains an improper integral function, it is computationally difficult to get a stable
estimate. With the transformation 
, 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*}](images/etsug_mdc0097.png)
Using the Gauss-Laguerre weight function, 
, 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}) \]](images/etsug_mdc0099.png){kind=small}