A Maximum Entropy Fixed-Point Route Choice Model for Route Correlation

Abstract

In this paper we present a stochastic route choice model for transit networks that explicitly addresses route correlation due to overlapping alternatives. The model is based on a multi-objective mathematical programming problem, the optimality conditions of which generate an extension to the Multinomial Logit models. The proposed model considers a fixed point problem for treating correlations between routes, which can be solved iteratively. We estimated the new model on the Santiago (Chile) Metro network and compared the results with other route choice models that can be found in the literature. The new model has better explanatory and predictive power that many other alternative models, correctly capturing the correlation factor. Our methodology can be extended to private transport networks.

1. Introduction

The present study formulates a new route choice model for public transport networks that features significant innovations compared to existing models. The main enhancement in the proposed model is the ability to simultaneously and explicitly integrate the traveler’s lack of information (randomness or uncertainty) and the correlation between route alternatives (due to overlapping). The proposed model is based on a multi-objective mathematical programming problem and its respective scalarized single-objective problem. The multi-objective problem considers the exogenous cost functions of the network, entropy of the route choice, and the covariance matrix for the route flows. An extension to networks with congestion (e.g., private transport with endogenous costs) is also proposed.

A traditional approach to defining the route choice process, and the subsequent traffic equilibrium, is to assume a deterministic behavior. This deterministic equilibrium usually states optimality conditions, such as minimizing transport costs or satisfying Wardrop’s first principle of traffic equilibrium [1]. These models assume that travelers have perfect information and seek to unilaterally minimize their travel costs [2–5]. Typically, a mathematical programming model is formulated and solved by an iterative algorithm. If applied with care and understanding, a deterministic user-equilibrium model provides a simple but effective method of traffic assignment [6–12].

Alternatively, stochastic/probabilistic route choice models differ from deterministic formulations in that they incorporate the uncertainty, randomness and/or the heterogeneity of travelers and alternative routes, and passenger’s imperfect knowledge. This is the approach followed in this study. Reviews of this class of models are found in Daganzo and Sheffi [13], Hazelton [14], Ramming [15], and Prashker and Bekhor [16]. Among these models, there are some that explicitly consider correlations between alternative routes, such as Cascetta et al. [17], Ben-Akiva and Bierlaire [18], Bekhor and Prashker [19,20], and Bovy et al. [21]. A complete review of these studies can be found in Prashker and Bekhor [16] and Prato [22]. Several of the models presented by the authors were considered for performing a comparison with our new model.

In the remainder of this paper, Section 2 contains a brief review of the literature that provides context for understanding the proposed model; Section 3 provides an analytic derivation of the new formulation; Section 4 applies the model to a medium-sized network (the Santiago Metro), comparing the results with existing models and proposing a version for private networks with congestion; and lastly, Section 5 summarizes the results and gives the main conclusions.

2. Literature Review

The formulation of the route choice or traffic assignment stage in transportation modeling has long followed an approach in which users minimize their generalized trip cost on the assumption of perfect knowledge of the transport network. Under this approach, travelers are considered to be homogeneous and each one is fully informed of the cost of each link on the network at any level of flow [1,2]. These assumptions are both rather strong even for a small network, and the results obtained are often not satisfactory. However, due to their simplicity and availability, many transport planners continue to apply such models, especially in large networks.

Because the perfect information assumption is usually not correct, there is a clear need for models that represent users who have incomplete or imperfect information on the transport system in regard to existing routes and their levels of congestion. Various route choice models in the specialized literature are based on system attributes perceived by travelers and their socioeconomic and demographic characteristics [13,15,16,23,24]. In these models, users behave in accordance with the costs they perceive. The socioeconomic and demographic variable data are usually obtained through user surveys or from network data records and are easily justified as an integral part of individuals’ rational decision-making processes. However, because the modeler lacks information on these processes, the modeled choices necessarily embody a degree of variability.

Regarding knowledge of routes, it is widely accepted that individual users do not know (or do not perceive/consider) all of the route possibilities between a given origin-destination pair [17,25–28]. This is a reality that should be incorporated into route choice models. Cascetta et al. [27] propose a Logit-type model of route perception and choice similar to those based on random utility theory. Ben-Akiva et al. [25] model interurban route choice as a two-stage process in which a set of routes is defined in the first stage, and the route choice is performed in the second stage.

Correlations between alternative routes can (to a certain extent) be indirectly addressed, according to Bovy and Hoogendoorn-Lanser [29], by hierarchical nested Logit and multi-nested GEV models in cases where it arises from overlapping segments and/or nodes. Prato and Bekhor [30] and Bekhor et al. [31] study the way in which individuals construct their set of route alternatives and the implications of route similarity for user behavior. Bliemer and Bovy [32] study the interdependence of these two aspects.

2.1. Multinomial Logit for Route Choice

The probability of choosing route p to travel between O-D pair w ( $P_w^p$) can be estimated by discrete choice models in which the cost of a given route as perceived by a user ( $P_w^c$) is assumed to be given by Equation (1):

$$C_w^p=c_w^p+{\epsilon }_w^p$$

(1)

The cost is thus perceived as the sum of a deterministic component ( $c_w^p$) and a random error component ( ${\epsilon }_w^p$). The latter can be interpreted in various ways, one of which is that the error reflects the user’s inaccurate perception of route cost due to, among other things, the scarcity of information. In this case, the first term ( $c_w^p$) represents the real mean cost of the route. The expression for the choice probability ( $P_w^p$) is a function of the assumed distribution of the error terms and whether or not they are independent. If we assume the errors are i.i.d. Gumbel-distributed with a scale parameter θ > 0, we have a multinomial logit (MNL) choice model in which $P_w^p$ is given by Equation (2) (see [5,23,33]):

$$P_w^p=\frac{\text{exp}\left(-\theta c_w^p\right)}{\sum _{r\in p^w}\text{exp}\left(-\theta c_w^r\right)}$$

(2)

where p^w is the set of routes uniting the O-D pair w.

An equivalent optimization problem that generates the stochastic assignment model (2) is (see [5,34,35]):

$$\begin{array}{ll}\underset{\{h_w^p\}}{\text{min}}Z=\sum _w\sum _{p\in p^w}{c}_w^p{h}_w^p+\hfill & \frac{1}{\theta }\sum _w\sum _{p\in p^w}{h}_w^p(\text{ln}{h}_w^p-1)\hfill \\ s.t.\hfill & \hfill \\ \sum _{p\in p^w}{h}_w^p=T_w\hfill & \forall w\hfill \end{array}$$

(3)

where T_w is the total number of trips (exogenous demand) between the O-D pair w and $h_w^p$ is the flow along route p ( $P_w^p=h_w^p/T_w$).

In problem (3), the set of routes p^w between each pair w must be predefined and maintained invariant during the assignment process. The optimality conditions lead to the assignment criterion given by Equation (2) in which users divide up among the alternative routes based on a Logit model.

It is worth noting that the objective function in Equation (3) considers two terms: the first term $\sum _w\sum _{p\in p^w}{c}_w^p{h}_w^p$ represents the total cost in the network, while the second term $\sum _w\sum _{p\in p^w}{h}_w^p\left(\text{ln}\hspace{0.17em}{h}_w^p-1\right)$ represents the negative value of the entropy. The latter is weighted by 1/θ to scale both terms in the objective function. In this way, problem (3) can be interpreted as a bi-objective problem that simultaneously minimizes the total cost and maximizes the total entropy [34].

The principal limitation of model (3) is that it does not incorporate a structure for capturing correlations between routes. Possible direct extensions are either extremely simplistic or difficult to apply correctly to real-world-scale networks (for example, using a hierarchical Logit model) given that they cannot properly capture the various types of correlations between routes that share links with each other in different ways. In urban networks, the routes linking a given O-D pair will typically have many overlaps due to common links, with the result that the independent error assumption implicit in Logit-type models such as Equation (2) is unrealistic.

Different extensions of the multinomial logit model have been proposed to explicitly capture correlations between alternative routes, some of which are presented below. These models are used as a comparative basis for our proposed model.

2.2. C-Logit Model

In Cascetta et al. [17] and Cascetta et al. [27], the authors propose a joint implicit availability/perception (IAP) and route choice (C-Logit) model that also explicitly addresses the correlation issue (i.e., the lack of route independence due to common links) and is analytically tractable even for large-scale networks.

The basic idea behind the model is to address route interdependence via a cost attribute called the “similarity factor”, which is added to the cost of the route in a conventional Logit model, instead of addressing it in terms of error non-independence as in the Probit model. The probability of choosing route p is then given by Equation (4):

$$P_w^p=\frac{\text{exp}\left(-\theta c_w^p+{\mathit{CF}}_w^p\right)}{\sum _{r\in p^w}\text{exp}\left(-\theta c_w^r+{\mathit{CF}}_w^r\right)}$$

(4)

where ${\mathit{CF}}_w^p$ is the “similarity factor” for route p joining pair w and is constructed as follows:

$${\mathit{CF}}_w^p=\beta \cdot \text{ln}\sum _{a\in p\hspace{0.17em}}\left(\frac{l_a}{L_p}\cdot \sum _{r\in p^w}{\delta }_{ar}\right)$$

(5)

where β is a parameter to be calibrated that must be negative, l_a is the length of link a, L_p is the length of route p, and δ_ar is equal to 1 if link a belongs to some route r joining w but is 0 otherwise. Other specifications for ${\mathit{CF}}_w^p$ may be found in Prato [22].

2.3. Path-Size Logit Model

Ben-Akiva and Ramming [36] develop the path-size logit (PSL) model, which also aims to correct for routes that have overlaps and are therefore correlated. It attempts to incorporate behavioral theory in Cascetta’s C-logit model. In this case, $P_w^p$is given by Equation (6):

$$P_w^p=\frac{\text{exp}\left(-\theta c_w^p+\beta \cdot \text{ln}\hspace{0.17em}{\mathit{PS}}_w^p\right)}{\sum _{r\in p^w}\text{exp}\left(-\theta c_w^r+\beta \cdot \text{ln}\hspace{0.17em}{\mathit{PS}}_w^r\right)}$$

(6)

where ${\mathit{PS}}_w^p$ is the correction for the route size.

The correction principle applied in this model is as follows: a route with no links overlapping another route needs no correction and is therefore assigned a size of 1. At the other extreme, if there are J duplicate routes (i.e., total overlap), each one has a size of 1/J. Lastly, the length of routes with partial overlap is based on the sizes of the links, which are appropriately weighted on some criterion, such as the link’s contribution to the total length of the route. Thus, ${\mathit{PS}}_w^p$ can take the following form:

$${\mathit{PS}}_w^p=\sum _{a\in p}\left(\frac{l_a}{L_p}\cdot \frac{1}{\sum _{r\in p^{*}}{\delta }_{ar}}\right)$$

(7)

where the variables are the same as those employed in Equation (5) to define the similarity factor. Further specifications for ${\mathit{PS}}_w^p$ may be determined in Bovy et al. [21].

2.4. Paired Combinatorial Logit Model

The paired combinatorial logit model (PCL) belongs to a family of models that derives from the generalized extreme value (GEV) model [33]. The general PCL model can be derived from the following generation function:

$$G(y_1,y_2,\dots ,y_n)=\sum _{k=1}^{n-1}{\sum _{j=k+1}^n(1-{\sigma }_{kj})\left(y_k^{1/1-{\sigma }_{kj}}+y_j^{1/1-{\sigma }_{kj}}\right)}^{1-{\sigma }_{kj}}$$

(8)

where:

• y_k characterizes each alternative (in our case, a route between an O-D pair),
• σ_kj is a similarity index between alternatives k, and
• n is the number of alternatives.

The model was adapted for route choice by Bekhor and Prashker [19], and, in our case, we use the expression used by Chen et al. [37]: $y_k=y_w^p=e^{-\theta c_w^p}$. Next, the probability of choosing route p when traveling between O-D pair w is Equation (9):

$$P_w^p=\frac{\sum _{q\in p^w,q\ne p}{e}^{\frac{-\theta c_w^p}{{}^{1-{\sigma }_w^{pq}}}}(1-{\sigma }_w^{pq}){\left(e^{\frac{-\theta c_w^p}{{}^{1-{\sigma }_w^{pq}}}}+e^{\frac{-\theta c_w^q}{{}^{1-{\sigma }_w^{pq}}}}\right)}^{-{\sigma }_w^{pq}}}{{\sum _{r=1}^{p^w-1}\sum _{m=r+1}^{p^w}(1-{\sigma }_w^{rm})\left(e^{\frac{-\theta c_w^r}{{}^{1-{\sigma }_w^{rm}}}}+e^{\frac{-\theta c_w^m}{{}^{1-{\sigma }_w^{rm}}}}\right)}^{1-{\sigma }_w^{rm}}}$$

(9)

In this case, for a choice set of p^w alternative routes, there are p^w(p^w − 1)/2 pairs of alternatives.If ${\sigma }_w^{pq}$ equals 0 for every pair (p,q) in w routes, then the PCL model reduces to a MNL model such as Equation (2). We use the similarity index proposed by Chen et al. [37]:

$${\sigma }_w^{pq}=\frac{L_w^{pq}}{\sqrt{L_w^p\cdot L_w^q}}$$

(10)

where $L_w^{pq}$ is the length of overlap between routes p and q and $L_w^p$ and $L_w^q$ are the respective lengths of routes p and q, respectively. It is possible (and convenient) to include an additional degree of freedom in this model by including a parameter β in Equation (10) as shown in Equation (11). This new parameter must be estimated together with the rest of the model’s parameters:

$${\sigma }_w^{pq}=\beta \cdot \frac{L_w^{pq}}{\sqrt{L_w^p\cdot L_w^q}}$$

(11)

2.5. Cross-Nested Logit Model

The cross-nested model (CNL) also belongs to the GEV family of models and was adapted for the route choice context by Prashker and Bekhor [38] and Vovsha and Bekhor [39]. In that adaptation, the model considers a two-level hierarchical structure. The upper level includes all of the links in the network, and the lower level includes all of the routes that belong to p^w. Next, every route is assigned to the nests that represent the links that belong to it. In this case, the probability of choosing route p is given by Equation (12):

$$P_w^p=\frac{e^{-\theta c_w^p}\sum _{m\in M_w^p}{\alpha }_w^{mp}{\left(\sum _{q\in p^w}{\alpha }_w^{mq}{e}^{-\theta c_w^q}\right)}^{\mu -1}}{\sum _{m\in M_w^p}{\left(\sum _{q\in p^w}{\alpha }_w^{mq}{e}^{-\theta c_w^q}\right)}^{\mu }}$$

(12)

where:

• m characterizes the links and therefore the nests,
• $M_w^p$ is the set of links that belongs to route p in the O-D pair w,
• ${\alpha }_w^{mp}$ are parameters that represent the degree of inclusion of alternative p in nest m, and
• µ is the nesting coefficient. If µ = 1, the model reduces to a MNL model.

The CNL model is adapted to a route choice context, defining the parameters ${\alpha }_w^{mp}$ depending on the topology of the network. Prashker and Bekhor [38] proposed the following expression:

$${\alpha }_w^{mp}={\left(\frac{L_w^m}{L_w^p}\right)}^{\gamma }{\delta }_w^{mp}$$

(13)

where $L_w^m$ is the length of link m, $L_w^p$is the length of route p, ${\delta }_w^{mp}$ equals 1 if link m belongs to route p, and γ is a parameter that must be calibrated, which reflects the perception of the travelers regarding the similarity of the alternative routes. For the estimation of the CNL model, we use γ = 1 following Bekhor et al. [31].

3. Route Choice Model with Correlated Routes

In this section, we develop a new route choice model that simultaneously incorporates: (i) users with imperfect knowledge of the network; (ii) correlations between alternatives (in this case, routes); and (iii) the effect of demand or flow levels on network costs (congestion). As noted earlier, the proposed formulation is entropy-based with quadratic constraints.

3.1. Mathematical Formulation of Model

The proposed stochastic equilibrium assignment model that captures route correlation is based on the following multi-objective optimization problem (other applications of multi-objective models in transportation can be found in [34,35,40])

$$\begin{array}{l}\text{min}\hspace{0.17em}{F}_1=\sum _w\sum _{p\in p^w}{c}_w^p{h}_w^p\hfill \\ \text{min}\hspace{0.17em}{F}_2=\sum _w\sum _{p\in p^w}{h}_w^p\left(\text{ln}{h}_w^p-1\right)\hfill \\ \text{min}\hspace{0.17em}{F}_3=\sum _w\sum _{p\in p^w}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^p-t_w\right)\left(h_w^q-t_w\right)\right]\hfill \\ s.t.:\hspace{0.17em}\sum _{p\in p^w}{h}_w^p=T_w,\forall w\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\left({\gamma }_w\right)\hfill \end{array}$$

(14)

Objective F₁ relates to the total system cost. Objective F₂ attempts to maximize the entropy to determine the most likely routes; in probabilistic terms, it finds the most feasible route combinations for travelers in equilibrium. Combining F₁ and F₂ with the flow conservation constraints gives the stochastic assignment model expressed in Equation (3).

Objective F₃, the novel element in the proposed model, explicitly incorporates correlations of flows between different routes, whether or not they join the same O-D pair w. The objective is constructed as a weighted sum of the divergences of the flows on the individual defined routes joining w from the average flow on those routes, where $t_w=\frac{T_w}{N_w}$ is the average flow, N_w is the number of defined routes (i.e., the cardinality of p^w) and T_w is the (fixed) total number of trips between w. The parameters ${\eta }_w^{pq}$ are exogenous and determine the degree of correlation (0 to 1) between the routes p and q of w. The values of the parameters can be defined in a number of ways (see [15,17,41]).

As with F₂, objective F₃ is an information criterion [42]. If, for example, the flows on all traveled routes were uniform (that is, if $h_w^p=t_w,\forall p$), the value of F₃ would be 0 and thus contain no information. F₂ also takes its lowest possible value if $h_w^p=t_w,\forall p$.

An alternative optimization problem [34,43] for Equation (14) that generates the proposed stochastic equilibrium model is the following:

$$\begin{array}{l}\text{min}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}F=\sum _w\sum _{p\in p^w}{c}_w^p{h}_w^p+\frac{1}{\theta }\sum _w\sum _{p\in p^w}{h}_w^p\left(\text{ln}{h}_w^p-1\right)+\frac{1}{\rho }\sum _w\sum _{p\in p^w}\sum _{\begin{array}{l}p\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^p-t_w\right)\left(h_w^q-t_w\right)\right]\hfill \\ s.t.:\sum _{p\in p^w}{h}_w^p=T_w,\forall w\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\left({\gamma }_w\right)\hfill \end{array}$$

(15)

The terms 1/θ and 1/ρ are the respective relative weights of the two information criteria F₂ and F₃ with respect to the reference objective F₁. θ and ρ are parameters to be estimated.

The first-order conditions for Equation (15) are:

$$c_w^p+\frac{1}{\theta }\text{ln}{h}_w^p+\frac{1}{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]+{\gamma }_w=0$$

(16)

where $c_w^p$ could be, for a public transportation application, replaced by $L_w^p=\sum _k{\beta }_k{X}_{w,k}^p$. From Equation (16), we obtain:

$$h_w^p=\text{exp}\left(-\theta c_w^p-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]-\theta {\gamma }_w\right)$$

(17)

$$\sum _{p\in p^w}{h}_w^p=T_w=\text{exp}\left(-\theta {\gamma }_w\right)\cdot \sum _{p\in p^w}\text{exp}\left(-\theta c_w^p-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]\right)$$

(18)

Dividing Equations (17) and (18) we have:

$$h_p^w=T_w\frac{\text{exp}\left(-\theta c_w^p-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]\right)}{\sum _{r\in p^w}\text{exp}\left(-\theta c_w^r-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne r\end{array}}\left[{\eta }_w^{rq}\left(h_w^q-t_w\right)\right]\right)}$$

(19)

Because ${\eta }_w^{pq}$ is an exogenous model parameter (defined by the modeler) and t_w is assumed to be constant (for calibration and modeling purposes), the term $\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in q^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{t}_w\right]$ is also constant. Letting $\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in q^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{t}_w\right]={\alpha }_w^q$, which is an intercept or specific constant, Equation (19) can then be rewritten as:

$$h_w^p=T_w\frac{\text{exp}\left({\alpha }_w^p-\theta c_w^p-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{h}_w^q\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^r-\theta c_w^r-\frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{rq}{h}_w^q\right]\right)}$$

(20)

This non-linear expression is similar in structure to the model specified by Ben-Akiva and Ramming [36], given here as Equation (6). The main difference is that the right-hand side of Equation (20) includes the endogenous variable $h_w^q$ and is therefore a fixed-point function [44], whereas in Equation (6), the right-hand side contains only the model’s exogenous variables.

An alternative can be specified for F₁ that can model service levels in public transportation networks. An example of such a replacement is ${\tilde{F}}_1=\sum _w\sum _{p\in p^w}\left[L_w^p{h}_w^p\right]$, where $L_w^p=\sum _k{\beta }_k,X_{w,k}^p$ is a generalized cost given by the weighted sum of attributes or explanatory variables $X_w^p$ that could represent, for example, the trip time, wait time, or the cost of route p between pair w.

We can define $V_w^p={\alpha }_w^p-\theta c_w^p-\frac{\theta }{\rho }\sum _{\underset{q\ne p}{q\in p^w}}\left[{\eta }_w^{pq}{h}_w^p\right]$ as the “utility level” of the proposed Logit-based model. By the route choice probability function in Equation (20), the marginal utility of the attribute $X_{w,k}^p$ can be written as:

$$\frac{\partial V_w^p}{\partial X_{w,k}^p}={\beta }_k+\frac{\partial V_w^p}{\partial X_{w,k}^p}\cdot \frac{1}{T_w}\cdot \frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}{\eta }_w^{pq}\cdot h_w^p\cdot h_w^q$$

(21)

Clearing, we obtain:

$$\frac{\partial V_w^p}{\partial X_{w,k}^p}=\frac{{\beta }_k}{1-\frac{1}{T_w}\cdot \frac{\theta }{\rho }\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}{\eta }_w^{pq}\cdot h_w^p\cdot h_w^q}$$

(22)

Because the denominator of Equation (21) is independent of attribute k, the marginal rates of substitution are generic and have the same functional form as traditional discrete choice models:

$$\frac{\partial V_r/\partial X_{w,k_1}^p}{\partial V_r/\partial X_{w,k_2}^p}=\frac{{\beta }_{k_1}}{{\beta }_{k_2}}$$

(23)

Thus, in the proposed fixed-point model with spatial correlation, the marginal rates of substitution between attributes can be obtained without any additional complexity.

To implement model (20), the parameters ( ${\alpha }_w^p$, θ, ρ^*), where ${\rho }^{*}=\frac{\theta }{\rho }$ must first be estimated and the fixed-point Equation (20) and then solved given that the variable $h_w^p$ appears on both sides of the equation.

However, because the endogenous variable ( $h_w^p$) appears on the right-hand side of the model, the parameters ( ${\alpha }_w^p$, θ, ρ^*) cannot be estimated using the maximum likelihood because the presence of the endogenous variable violates the assumption of independent marginal probability functions necessary for defining the likelihood function. To bypass this problem, we resort to the use of an instrumental variable to replace $h_w^p$ as explained below.

3.2. Estimation of Model Parameters

An instrument or instrumental variable [45] is an exogenous variable that is highly correlated with an explanatory variable exhibiting endogeneity and can therefore be used as a replacement for the latter without a loss of asymptotic properties in the estimated parameters. In our case, a suitable instrument $h_w^{p0}$ to replace $h_w^p$ is:

$$h_w^{p0}=T_w\frac{\text{exp}\left(-\lambda c_w^p\right)}{\sum _{r\in p^w}\text{exp}\left(-\lambda c_w^r\right)}$$

(24)

This formula is a classic MNL model and can be easily estimated. Once this is done, the value of $h_w^{p0}$ is substituted into the right-hand side of Equation (20), as follows:

$$h_w^p=T_w\frac{\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{h}_w^{q0}\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^r-\theta c_w^r-{\rho }^{*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{rq}{h}_w^{q0}\right]\right)}$$

(25)

Unlike Equation (20), the parameters ( ${\alpha }_w^p$, θ, ρ^*) in Equation (25) can be directly estimated by the maximum likelihood [46–48] because $h_w^{q0}$ does not exhibit endogeneity. Using these estimates, which we now denote ( ${\alpha }_w^{p0}$, θ⁰, ρ⁰*), we can estimate $h_w^{p1}$ by Equation (26):

$$h_w^{p1}=T_w\frac{\text{exp}\left({\alpha }_w^{p0}-{\theta }^0{c}_w^p-{\rho }^{0*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{h}_w^{q0}\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^{r0}-{\theta }^0{c}_w^r-{\rho }^{0*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}{\eta }_w^{rq}{h}_w^{q0}\right)}$$

(26)

The original model (26) is thus estimated iteratively from the following recursive relation (which represents the equilibrium of the fixed-point function):

$$h_w^{p(n+1)}=T_w\frac{\text{exp}\left({\alpha }_w^{p(n)}-{\theta }^{(n)}{c}_w^p-{\rho }^{(n)*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}{h}_w^{q(n)}\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^{r(n)}-{\theta }^{(n)}{c}_w^r-{\rho }^{(n)*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{rq}{h}_w^{q(n)}\right]\right)}$$

(27)

The iterative estimation process concludes when convergence in obtained; that is, $\left({\alpha }_w^{p\left(n\right)},{\theta }^{\left(n\right)},{\rho }^{\left(n\right)*}\right)\approx \left({\alpha }_w^{p\left(n-1\right)},{\theta }^{\left(n-1\right)},{\rho }^{\left(n-1\right)*}\right)$. If the estimator of parameter ρ^(n)* in Equation (27) is significantly different from 0, the null hypothesis of no correlation between route p and any of the other routes joining O-D pair w is rejected.

In an analogous fashion to the derivation of Equation (21), upon iteratively solving model (20), the marginal utility of attribute $X_{w,k}^p$ is then:

$$\frac{\partial V_w^{p(n+1)}}{\partial X_{w,k}^p}={\beta }_k^{(n)}+\frac{\partial V_w^{p(n)}}{\partial X_{w,k}^p}\cdot \frac{1}{T}\cdot {\rho }^{(n)*}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}{\eta }_w^{pq}\cdot h_w^{p(n)}\cdot h_w^{q(n)}$$

(28)

This expression generates the marginal utilities recursively in successive iterations. The model converges to the equilibrium of the fixed-point function as follows:

$$\begin{array}{ccc}\frac{\partial V_w^{p(n+1)}}{\partial X_{w,k}^p}\approx \frac{\partial V_w^{p(n)}}{\partial X_{w,k}^p}& \to & \frac{\partial V_r/\partial V_{w,k_1}^{p(n+1)}}{\partial V_r/\partial X_{w,k_2}^{p(n+1)}}=\frac{{\beta }_{k_1}^{(n)}}{{\beta }_{k_2}^{(n)}}\end{array}$$

(29)

The marginal rates of substitution thus continue to be generic and equal in their functional form to those of the traditional discrete choice models.

However, because the endogenous variable ( $h_w^p$) appears on the right-hand side of the model, the parameters ( ${\alpha }_w^p$, θ, ρ^*) cannot be estimated using the maximum likelihood because the presence of the endogenous variable violates the assumption of independent marginal probability functions necessary for defining the likelihood function. To bypass this problem, we resort to the use of an instrumental variable to replace $h_w^p$ as explained below.

3.3. Existence and Uniqueness of the Fixed-Point Problem Solution

To demonstrate the existence of a solution to the fixed-point problem expressed by Equation (27), we apply the Brouwer fixed-point theorem [49]. This theorem may be stated as follows:

Let H be a non-empty compact convex set of a finite-dimensional Euclidean space k, and let f : H → H be a continuous function. Then, f has a fixed point, i.e., ∃h ∈ H : h = f(h).

Below, it is shown that the system of equations describing the equilibrium in model (27) satisfies the theorem’s existence and completeness hypotheses.

Let $h={\left(h_w^p\right)}_{p\in P^w,w\in W}$, h ∈ H, with $H=\underset{w\in W}{\otimes }{T}_w{\Delta }^{N_w}$ where ${\Delta }^n\left\{x\in X_{+}^n:{\mathit{1}}^{T}x=1\right\}$ is the probability simplex. H is a nonempty convex compact set because it is a finite cross product of nonempty convex compact sets.

The function f: H_w → H_w is defined as:

$$f\left({\mathit{h}}_{\mathit{w}}\right)={\left(T_w\frac{\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^q\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^r-\theta c_w^r-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{rq}{h}_w^q\right]\right)}\right)}_p$$

(30)

This expression is a construction of two continuous functions and is therefore itself continuous. The parameters ( ${\alpha }_w^p$, θ, ρ^*) are estimated using the maximum likelihood given an h_w and thus are continuous functions with respect to the latter, thereby satisfying the Brouwer’s theorem hypotheses and proving the existence of at least one fixed point.

The function f is continuous because it is a composition of continuous functions. In this case, the parameters ( ${\alpha }_w^p$, θ, ρ) are continuous functions of h, which is given by the continuity and strict concavity ( ${\alpha }_w^p$, θ, ρ) of the log-likelihood function to be used in the multinomial logit model. The estimation of our model, which corresponds to the point where succession described in the previous section converges, belongs to the set H and will be a fixed point of the function f.

Using Brouwer’s Theorem, we can establish the existence of at least one fixed point, i.e., ∃h ∈ H : h = f(h), which is exactly the point of equilibrium that we want:

In formal terms:

$$\exists {\mathit{h}}^{\mathbf{*}}\in H:f\left({\mathit{h}}^{\mathbf{*}}\right)=\left(T_w\frac{\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{q^{*}}\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^r-\theta c_w^r-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{rq}{h}_w^{q*}\right]\right)}\right)$$

(31)

We now prove uniqueness for the case of ρ^* ≥ 0. Assume that there are 2 different equilibrium points h^*, g^* ∈ H. It is therefore the case that:

$$\left({\mathit{h}}^{*}\ne {\mathit{g}}^{*}\wedge \sum _p{h}_w^{p*}=\sum _p{g}_w^{p*}=T_w\right)$$

(32)

$$\Rightarrow \exists p^{+},p^{-}\in \mathbb{N}:\left(p^{+}\ne p^{-};1\le p^{+},p^{-}\le n;h_w^{p*-}<g_w^{p*-};h_w^{p*-}>g_w^{p*+}\right)$$

(33)

$$\frac{h_w^{p*}}{g_w^{p*}}=\frac{\sum _{p\in p^w}\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{g}_w^{q*}\right]\right)}{\sum _{p\in p^w}\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{p*}\right]\right)}\times \text{exp}\left\{{\rho }^{*}\left(\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{p*}\right]-\sum _{q\in p^w}\left[{\eta }_w^{pq}{g}_w^{p*}\right]\right)\right\}$$

(34)

Without a loss of generality, we can assume that:

$$\exists {\mathit{h}}^{*}\in H:f\left({\mathit{h}}^{*}\right)=\left(T_w\frac{\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{q*}\right]\right)}{\sum _{r\in p^w}\text{exp}\left({\alpha }_w^r-\theta c_w^r-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{rq}{h}_w^{q*}\right]\right)}\right)$$

(35)

For the case of p⁺, we have:

$$h_w^{p*+}>g_w^{p*+}\Rightarrow \frac{h_w^{p*+}}{g_w^{p*+}}>1$$

(36)

$$\Rightarrow \frac{\sum _{p\in p^w}\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{g}_w^{p*+}\right]\right)}{\underset{\le 1}{\underbrace{\sum _{p\in p^w}\text{exp}\left({\alpha }_w^p-\theta c_w^p-{\rho }^{*}\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{p*+}\right]\right)}}}\times \underset{<1}{\underbrace{e^{-\rho *\left(\sum _{q\in p^w}\left[{\eta }_w^{pq}{h}_w^{p*+}\right]-\sum _{q\in p^w}\left[{\eta }_w^{pq}{g}_w^{p*+}\right]\right)}}}<1$$

(37)

This is true given assumption (35) and provided that ρ^* ≥ 0. From this, we deduce that $-{\rho }^{*}\left(h_w^{p^{*}}-g_w^{p^{*}}\right)<0\Rightarrow e^{-{\rho }^{*}\left(h_w^{p^{*}}-g_w^{p^{*}}\right)}<1.$. There is thus a contradiction, and we conclude that a fixed point not only exists but is unique.

4. Numerical Results and Extensions

In Section 4.1, we present a route choice case study comparing the results of a real application of the proposed model to those generated by the classic MNL model in Equation (2), C-Logit model in Equation (4), path-size logit model in Equation (6), PCL model in Equation (9) and the CNL model in Equation (12). Next, in Section 4.2, we develop an extension of the proposed formulation with route correlations to traffic assignments on congested networks (private transportation).

4.1. Application to a Medium-Sized Network: The Santiago Metro

This case study applies the various models to route choice on the Metro rail transport system in the Chilean city of Santiago. Successive transfer points are chosen between origin and destination pairs that in many instances are joined by more than one feasible route (see Figure 1). The analysis focuses on the morning and evening peak periods (7 am to 9 am and 6 pm to 8 pm) when approximately 790,000 trips are taken across the system, 44% of which include transfers.

Trip data were obtained through an O-D survey conducted on the Metro in which 92,800 system users, or approximately 12% of the total, participated. Because only those trips for which there existed more than one alternative route were retained in the data set, the number of individuals or observations eventually employed was 16,029, or approximately 40% of users who transferred.

When the survey was performed in October 2008, the Santiago Metro consisted of five lines and 85 stations, seven of which were transfer points. Of the 7140 O-D pairs in the network, 4985 (70%) of them required transfers. The reasonability criterion applied in including a route as a possible alternative for a given pair was that at least one surveyed traveler was observed to have used it. Data on the alternative routes for O-D pairs across the system that had more than one route on this criterion are given in

Table 1. Origin-destination pairs and alternative routes, 2008.

No. of Observed Alternative Routes	% of All O-D Pairs	% of All Trips
2	97%	93%
3	3%	7%
4	<1%	<1%

. Although in the majority of cases there were only two alternatives, some pairs had as many as four. Denser networks than this one would undoubtedly have a greater proportion of pairs with alternative routes.

For the proposed model (20), the generalized cost of traveling between pair w on route p is defined as $L_w^p={\beta }_{\mathit{time}}{X}_{\mathit{time},w}^p+{\beta }_{\mathit{access}}{X}_{\mathit{access},w}^p+{\beta }_{\mathit{trans}}{X}_{\mathit{trans},w}^p+{\beta }_{rr}{D}_{rr,w}^p+{\beta }_{\mathit{old}}{D}_{\mathit{old},w}^p$, where $X_{\mathit{time},w}^p$ is total trip time, $X_{\mathit{access},w}^p$ is the access time (walking and waiting), $X_{\mathit{trans},w}^p$ is the number of transfers along the route, $D_{rr,w}^p$ is a binary variable that takes the value 1 if route p is reasonable for travelling in w (based onDial [23]) and takes the value 0 if not, $D_{\mathit{old},w}^p$ is a binary variable that takes the value 1 if the route is less than 10 years old (a proxy of how well known the route is) and takes the value 0 if it is older. The number of transfers on a given route alternative is considered as an indicator of the disutility of transferring. Because the Metro uses a flat fare system, the fare is the same regardless of the route and can therefore be omitted from the cost function.

The ${\eta }_w^{pq}$ term, the spatial correlation factor due to route overlap between routes p and q joining pair w, is defined as suggested by Yai et al. [41]:

$${\eta }_w^{pq}=\frac{D_w^{pq}}{\sqrt{D_w^p\cdot D_w^q}}$$

(38)

where $D_w^{pq}$ is the length of overlap between routes p and q and $D_w^p$ and $D_w^q$ are the respective lengths of routes p and q. This expression stems from the traditional definition of the similarity factor in the C-logit model. As with the parameter β in that model (see Equation (5) above), the parameter –ρ^* must be negative, ensuring that ρ^* is positive.

The proposed model was compared to a conventional MNL model that included only network service level variables. That is, ρ^* was set to 0. The comparison was performed using statistical tests. Thus, the six route choice models estimated for the Santiago Metro were the following:

(a)

Multinomial logit (MNL). This is the base model because it does not account for correlations and is also used to construct the route choice proxy variables.

(b)

C-logit.

(c)

Path-size logit (PSL).

(d)

PCL.

(e)

CNL.

(f)

Fixed-point model (FPM) with spatial correlations, which is our proposed model.

The estimation results are set out in

Table 2. Route Choice Model Estimation Results.

Variable	MNL		C-Logit		PSL		PCL		CNL		FPM
	Parameter	Test-t	Parameter	Test-t	Parameter	Test-t	Parameter	Test-t	Parameter	Test-t	Parameter	Test-t
Trip time	−0.124	−38.1	−0.126	−37.8	−0.125	−37.5	−0.102	−35.2	−0.407	−11.9	−0.097	−29.2
Walk and waiting time	−0.192	−7.6	−0.212	−8.0	−0.203	−7.7	−0.190	−9.6	−0.518	−7.5	−0.136	−5.5
No. of transfers	−0.853	−13.4	−0.811	−12.3	−0.828	−12.6	−0.588	−11.4	−2.031	−8.9	−0.610	−9.4
Reasonable route	−0.445	−11.1	−0.458	−11.3	−0.454	−11.2	−0.268	−8.2	−0.968	−8.1	−0.239	−5.7
Old route	0.458	10.3	0.448	10.1	0.452	10.2	0.319	9.6	0.979	7.7	0.336	7.2
Spatial correlation	-	-	−0.351	−2.4	0.222	1.5	0.571	19.2	0.306	12.4	−1.980	−17.9
Log-likelihood	−7225.35		−7222.51		−7224.31		−7117.90		−7060.65		−7062.12
Adjusted rho squared	0.375		0.375		0.375		0.384		0.390		0.389
Sample Size	16,029		16,029		16,029		16,029		16,029		16,029

. The proposed FPM converged in only four iterations at the 0.01% tolerance level. As can be seen, all explanatory variables had the correct sign and were statistically significant. In the C-logic, PCL, CNL and FPL models, the correlation parameter was significant; in the PSL model, the statistical significance was lower.

The models with better goodness-of-fit (log-likelihood value) were the CNL and our proposed FPL model, and both were very similar. In addition to the log-likelihood function, the following indicators were used to compare the MNL and fixed-point models based on their goodness of fit [50,51]:

• percent of correct predictions (PCP);
• residual sum of squares (RSS): $S=\sum _i{\left(Y_i-P_i\right)}^2$, where Y_i is 1 if an alternative i is chosen and 0 otherwise and P_i is the probability predicted by the model of choosing alternative I; and
• weighted residual sum of squares (WRSS): $S^{\prime }=\frac{\sum _i{\left(Y_i-P_i\right)}^2}{P_i\left(1-P_i\right)}$.

The results of the three indicators for the estimated models are given in

Table 3. Goodness-of-fit indicators for route choice models.

Indicator	MNL	C-Logit	PSL	PCL	CNL	FPM
PCP	81.4%	81.4%	81.4%	81.4%	81.6%	81.6%
RSS	2273.1	2272.0	2272.5	2253.9	2227.5	2236.5
WRSS	28,072.2	28,339.0	28,320.3	23,147.8	19,745.0	19,757.9

. The PCP indicator is practically the same for all models; therefore, it does not allow any analysis regarding forecasting capability. In regard to RSS and WRSS, once again the best models are the CNL and FPL models.

4.2. Extension of Proposed Fixed-Point Model to Traffic Assignment with Route Correlation

The proposed fixed-point model can be extended to the traffic assignment problem in any congested transportation networks. The multi-objective optimization problem would then take the following form:

$$\begin{array}{ll}\text{min}\hfill & F_1=\sum _a\underset{0}{\overset{f_a}{\int }}{c}_a(x)dx\hfill \\ \text{min}\hfill & F_2=\sum _w\sum _{p\in p^w}{h}_w^p(\text{ln}{h}_w^p-1)\hfill \\ \text{min}\hfill & F_3=\sum _w\sum _{p\in p^w}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^p-t_w\right)\left(h_w^q-t_w\right)\right]\hfill \\ s.t.:\hfill & \sum _{p\in p^w}{h}_w^p=T_w,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\forall w\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}({\gamma }_w)\hfill \\ \hfill & \sum _{\begin{array}{l}p\in p^w\\ a\in p\end{array}}{h}_w^p=f_a,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\forall a\hfill \end{array}$$

(39)

Objective F₁ is the classic Beckmann transformation of the traffic assignment problem [2]. Objectives F₂ and F₃ were described above in Section 3.1 as information criteria. Combining F₂ and F₃ with flow conservation constraints gives the Fisk stochastic assignment model [52].

A substitute optimization problem for Equation (14) that defines a new stochastic traffic assignment model with route correlation is the following:

$$\begin{array}{ll}\text{min}\hfill & Z=\sum _a\underset{0}{\overset{f_a}{\int }}{c}_a(x)dx+\frac{1}{\theta }\sum _{p\in p^w}\sum _w{h}_w^p(\text{ln}{h}_w^p-1)+\frac{1}{\rho }\sum _w\sum _{p\in p^w}\sum _{\begin{array}{l}q\in p^w\\ q\ne p\end{array}}\left[{\eta }_w^{pq}\left(h_w^p-t_w\right)\left(h_w^q-t_w\right)\right]\hfill \\ s.t.:\hfill & \sum _{p\in p^w}{h}_w^p=T_w,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\forall w\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}({\gamma }_w)\hfill \\ \hfill & \sum _{\begin{array}{l}p\in p^w\\ a\in p\end{array}}{h}_w^p=f_a,\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\hspace{0.17em}\forall a\hspace{0.17em}\hspace{0.17em}\hfill \end{array}$$

(40)

The first-order conditions for Equation (40) are:

$$h_w^p=T_w\frac{\text{exp}\left(-\theta C_w^p-\frac{\theta }{\rho }\sum _{q\in P_w}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]\right)}{\sum _{p\in P_w}\text{exp}\left(-\theta C_w^p-\frac{\theta }{\rho }\sum _{q\in P_w}\left[{\eta }_w^{pq}\left(h_w^q-t_w\right)\right]\right)}$$

(41)

where $C_w^P=\sum _{\begin{array}{l}p\in p^w\\ a\in p\end{array}}{c}_a\left(f_a^{*}\right)$.

5. Conclusions

A new transportation network route choice model was developed that explicitly incorporates the phenomenon of correlations between routes. The model is applicable to public transportation networks and extensible to the traffic assignment problem. The presence of correlations between route alternatives was contrasted empirically by means of classical econometric techniques.

A multi-objective problem was stated and a substitute problem then formulated whose optimality conditions yielded a logit specification with an endogenous variable constituting a fixed-point model that is estimated and solved iteratively. The estimation was performed by the maximum likelihood and by the use of instrumental variables due to the presence of endogeneity in the model’s explanatory variables. The functional form of the proposed fixed-point model combined with the utilization of instrumental variables guaranteed both the existence and uniqueness of the solution.

The proposed model was compared with other route choice models reported in the literature in a case study of the Santiago Metro. The results obtained were both satisfactory and superior to many other existing formulations, and statistically equivalent to CNL model. Unlike the latter, the proposed fixed-point model was able to capture the correlation between routes and provided a better goodness of fit. A future research should be related to compare FPM with CNL in large networks. However, FPM model is able to include congestion in traffic equilibrium conditions; this is an advantage of our new model.

Although the proposed formulation could be more complex to estimate due to the iterative process that must be used to obtain the parameters, in practice the iterations converge rapidly. Furthermore, the additional estimation complexity does not complicate the derivation of project evaluation indicators such as the marginal rates of substitution, which retain the simple functional form of the MNL models.

Lastly, future research should address the apparent advantages of the proposed model on larger networks and traffic assignment problems.