#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Unravelling travellers’ route choice behaviour at full-scale urban network by focusing on representative OD pairs in computer experiments


Authors: Humberto González Ramírez aff001;  Ludovic Leclercq aff001;  Nicolas Chiabaut aff001;  Cécile Becarie aff001;  Jean Krug aff001
Authors place of work: Univ. Lyon, Univ. Gustave Eiffel, IFSTTAR, ENTPE, LICIT, Lyon, France aff001
Published in the journal: PLoS ONE 14(11)
Category: Research Article
doi: https://doi.org/10.1371/journal.pone.0225069

Summary

In a city-scale network, trips are made in thousands of origin-destination (OD) pairs connected by multiple routes, resulting in a large number of alternatives with diverse characteristics that influence the route choice behaviour of the travellers. As a consequence, to accurately predict user choices at full network scale, a route choice model should be scalable to suit all possible configurations that may be encountered. In this article, a new methodology to obtain such a model is proposed. The main idea is to use clustering analysis to obtain a small set of representative OD pairs and routes that can be investigated in detail through computer route choice experiments to collect observations on travellers behaviour. The results are then scaled-up to all other OD pairs in the network. It was found that 9 OD pair configurations are sufficient to represent the network of Lyon, France, composed of 96,096 OD pairs and 559,423 routes. The observations, collected over these nine representative OD pair configurations, were used to estimate three mixed logit models. The predictive accuracy of the three models was tested against the predictive accuracy of the same models (with the same specification), but estimated over randomly selected OD pair configurations. The obtained results show that the models estimated with the representative OD pairs are superior in predictive accuracy, thus suggesting the scaling-up to the entire network of the choices of the participants over the representative OD pair configurations, and validating the methodology in this study.

Keywords:

Network analysis – Behavior – Roads – Transportation – Computer networks – Decision making – Statistical distributions

Introduction

Urban congestion occurs when traffic demand locally exceeds the network capacity. The local demand is the combination of the global travel demand between the different origins and destinations and the travellers’ route choices, which define how many trips are made at the same place in a given time window. Thus, at a city-scale level, i.e., considering all the OD pairs and links in the network, route choice is a key determinant of urban transportation network performance. Route choice behaviour has been extensively studied in the transportation literature from two main but different angles. The first, related to human factors and mainly founded in discrete choice models [13], is focused on the identification of the determinants of travellers’ individual choices. This line of research is based on investigating travellers’ behaviour through experiments that consist in either observing their choices in the field (revealed preference) or asking them what would be their choices in hypothetical scenarios (stated preference). The second line of research, tackles the problem at full-scale and aims to solve the network loading problem to determine static or dynamic traffic states over all the network links. In this case, the interactions between the demand and the route choices on all the OD pairs in the network are considered altogether to define general principles that determine the network equilibrium. This is, for example, the case of the deterministic network equilibrium principle [4], that states that the travellers are selfish optimisers who only try to minimise their travel costs when choosing a route amongst all the alternatives; at the equilibrium, all the used routes that connect an OD pair have the same minimal cost.

Theoretically, the study of route choice at an individual and network level are consistent, however, in practice, there is a lack of connection between the two [5]. The reason is that, on the one hand, studies of route choice behaviour are focused in specific determinants of travellers’ route choice and, therefore, are based on simple scenarios (two or three routes in few OD pair configurations) that do not cover the multiplicity of situations that are found in a city-scale transportation network. In these experiments, particular attention has been paid to the study of how travellers learn from experience [68], the impact of advanced travel information systems (ATIS) [915], and the effect of travel time variability and risk attitudes in the travellers choices [1618]. On the other hand, in the network loading problem, representations have been designed as a simplified mathematical abstractions that permit to calculate the network loading under different behavioural principles, such as the deterministic user equilibrium [4], stochastic user equilibrium [19] or bounded rational user equilibrium [20]. These representations often assume that the only variable influencing travellers’ route choice behaviour is the travel time, ignoring other local factors, related to the network OD configuration, that have been recognised to influence route choice behaviour [2125]. One of the main reasons for the gap between research in individual route choice behaviour and network loading is the lack of observations at large scale over a sufficient number of OD configurations, that would allow discrete choice models to scale-up at the network level and thus enable the design of network equilibrium founded in a more user-oriented approach. The ambition of this study is to fulfil this gap by the selection of OD pairs that are representative of the OD configurations that are found in a transportation network, and then use these OD pairs in computer experiments to collect data on travellers’ route choice behaviour.

In a city-scale network, trips are made in thousands of OD pairs connected by several routes (in the case concerning this study, the city of Lyon in France, the network has 96,096 OD pairs and 559,423 routes), resulting in a large number of diverse routes, consequence of the topology of the network. For example, the route alternatives connecting an OD pair located in the central part of a city are likely to have short length, a high number of intersections and turns, but are unlikely to include segments of freeways. In contrast, the routes connecting an OD pair that traverses the city are longer and are more likely to include routes with fewer number of intersections and segments of freeways. From the point of view of the design of experiments, this implies that the number of scenarios must be reduced to a small but representative set of scenarios, such that the choices of travellers in any scenario found in the network can be approximated by a choice model estimated with this small set. More specifically, a representative set of OD pairs and routes is such that, for any randomly sampled OD pair in the network it is possible to find an OD pair in the representative set with similar attributes. Thus, assuming that the choices of travellers are similar for similar situations, an estimated model on the representative OD pairs could adequately reproduce the choices in the rest of the OD pairs. The question addressed in this work is: how to find a set of OD pairs and routes, such that it is representative of the OD configurations and route attributes found in the network, while being small enough so that a sufficient number of observations on route choices can be collected through computer route choice experiments?

The solution proposed in this work is based on k-means clustering [26] of the full set of OD pairs and routes in the network. In cluster analysis, the observations, in this case OD pairs and routes, are grouped in clusters characterised for having elements that are similar among themselves, but dissimilar to the elements in the other clusters. In the problem pertaining this article, the elements in a cluster will show similar orientation, length, % of freeway, directness and number of turns, and thus a cluster Ci will be, for example, composed mainly of OD pairs of short length in the central part of the city, with direct routes and low % of freeway composition, whereas another cluster Cj will be composed of OD pairs representing long trips traversing the city, with some non-direct routes composed mainly of freeways. Assuming that there are k clusters, the elements of a cluster can be regarded as belonging to a same class of OD pair configurations, and the whole network as being composed of elements of k different classes. Therefore, the OD pairs and routes in the network can be represented by elements in the k clusters. A natural choice to represent the elements in a cluster is the mean element in the cluster (cluster centroid), as it is the point with minimum euclidean distance to all the elements in the cluster. Thus, the cluster centroids are chosen as representative of the clusters’ elements, and the k clusters’ centroids as representative of the OD pairs and routes in the whole network. These OD pairs and routes are then used in computer experiments to collect data on travellers route choice behaviour. Note that the set representative OD pairs and routes found with k-means is representative of the attributes of the network, so the question that arises here is if a model estimated over this representative set can adequately reproduce the choices in the rest of the OD pairs in the network. To answer this question, three discrete choice models are estimated with the observations over the representative set. The discrete choice model used in this work is a joint mixed logit model (MXL), which under certain conditions, as is the case in this study, is equivalent to to the panel data formulation of MXL models [2, 2729]. The predictive accuracy of these models are compared against the predictive accuracy of the same models, but estimated with randomly chosen sets of OD pair configurations in a sort of cross validation procedure.

The results of the above methodology are that the models estimated with the observations over the representative OD pair configurations are better in predicting the route choices on unseen OD pairs, i.e., on OD pairs not used for the estimation process. On the one hand, these results demonstrate how a careful selection of OD pairs for experiments on route choice behaviour can improve the results of a choice model in a broader set of OD pairs and, on the other hand, that cluster analysis can be used to find these OD pair configurations. These findings have direct implications for urban traffic simulators, which solve the network loading problem to determine the time-evolving traffic states in the network. The scalable route choice model proposed in this paper can be implemented in such simulators without adding significantly computational complexity, compared to the usual simple equilibrium rule, e.g., user equilibrium. Furthermore, the use of clustering techniques to find the most relevant OD pairs and routes in the network, provides an efficient method to calibrate route choice models that can be easily replicated in any urban transportation network.

Materials and methods

Participants in the route choice experiments

The data on route choice behaviour in this article comes from six route choice experiments carried out between February 2018 and February 2019. The participants in the experiments were students at the University of Lyon taking part in the courses of traffic theory (66%), staff from the IFSTTAR (French Institute of Science and Technology for Transport, Development and Networks) and other universities, who received an invitation by e-mail to remotely join the experiments via a web browser (34%). All participants have signed, before the experiments begin, an informed consent form describing the objectives of the study, the data collection and processing, and the confidentially rules. Participants could opt out of the experiment at any time. No personal data were mandatory to participate to the experiments as people had the opportunity to identify themselves by a login of their choice. Finally, all data were fully anonymised and processed as such. At the beginning of the experiments, the participants were briefed about the objective of the experiment and the interface of the experimental platform; for the participants that joined the experiments via web, a document with the instructions was shared. The participants were instructed to choose the route that they consider the best to complete a trip on time.

Three of the six experiments were specifically implemented for the purpose of this work, so they were configured to obtain observations in the 9 representative OD pairs. The rest of the experiments were implemented for previous studies, so they were configured over 21 OD pairs different from the 9 representative set; data coming from these experiments was used to validate the methodology in this work. Throughout the six experiments, 3,334 choices of 483 participants were recorded, from which 802 choices of 73 participants were made over the nine representative OD pairs. In the experiments, the participants were confronted to several route choice problems in the different OD pairs, the task of the participants was to choose one of the three alternative routes to complete the trip before a given time.

Obtaining representative OD pairs and routes

Data: The Lyon network

The origins and destinations in the network of the city of Lyon, France, come from the zoning by the National Institute of Statistics and Economic Studies (INSEE) [30], and the major entry/exit points to the network. The zones are the geostatistic units used for the trip demand estimations and represent the origins and destinations of the trips generated or terminated inside the network. The entry and exit points represent the origins and destinations of the demand coming or going outside the network. In total, there are 285 zones, 29 entry points and 28 exit points in the Lyon network. The total number of origins is 313 and of destinations 310 (this quantity does not correspond exactly to the sum of zones plus entries/exits as there are zones that may have no outgoing or incoming trips), giving a total number of 96,096 OD pairs (see Fig 1). The most likely routes joining the origins to the destinations are derived with the A* algorithm looking for the k-shortest paths in free-flow travel time. This is roughly equivalent to minimising the travel distance, but accounting for the influence of spacial limits. The total number of routes in the network, obtained with this algorithm, is 559,423, with an average number of 5.82 routes per OD pair.

Lyon road network.
Fig. 1. Lyon road network.
The zones are depicted in different colours with their centroids in blue. The entry/exit points to the network are depicted with yellow points. There is an average of 5.82 routes connecting each origin and destination. The geodata used to render the plot is from ©OpenStreetMap contributors, licensed under the Open Data Commons Open Database License (ODbL). The zones are from Contours…Iris® licensed under the Open Licence v1.0 from Etalab. The maps were rendered using The R Project for Statistical Computing which is distributed as Free Software under the terms of the Free Software Foundation’s GNU General Public License.

Data representation

The selected route features in this study are the informed travel time, the length, directness, number of turns per kilometre and the percentage of freeway in the route composition. This features were selected as they are variables relevant in travellers’ route choice behaviour [2125], and because they are the attributes that participants can observe in the computer route choice experiments. In the experiments, the number of routes connecting each OD pair is limited to three. This limitation, however, does not restrict the scope of the experiments or diminishes the quality of the results for two reasons. First, choice sets with many alternatives may be burdensome for participants as they may have trouble identifying the differences between the routes. Second, the low variability between routes attributes due to the small number of alternatives in the choice set is compensated by the presence of many OD pairs, that are considered by jointly estimating a random utility model.

An OD pair and three routes connecting the origin and destination, defined as OD-routes, are characterised by the variables describing the origin and the destination (latitude, longitude and the euclidean distance between them), and the variables describing the three routes connecting them (the length of the route, the number of turns per kilometre, the directness of the route, and the percentage of freeway in the route). An OD-routes is then defined by 17 variables: 5 OD pair specific and 12 describing the routes (4 for each route). An OD-routes is represented as a vector in which the attributes of the three routes appear ordered by length, from shortest to longest. A depiction of the OD-routes objects is shown in Fig 2.

OD-routes vector.
Fig. 2. OD-routes vector.
The vector is composed of the attributes of the OD pair and the three routes connecting the origin and destination, with length(Route_1) ≤ length(Route_2) ≤ length(Route_3).

For the clustering of the OD pairs and routes, and thus the route choice experiments, the short routes (less than 1.5 km) and highly overlapping routes belonging to the same OD pair (sharing more than 70% of their links) were not considered. The reason is that very short trips lack of real alternatives: usually there is an unique route to travel from origin to destination. The highly overlapping routes are removed from the analysis because, from a route choice experiment perspective, the similarity between the routes may cause participants not to consider some routes as real alternatives and, furthermore, highly overlapping routes lack of the variability required for a choice model to capture the impact of each route attribute in the choices. After removal of the very short trips and the high overlapping routes, the OD-routes are obtained by considering all the possible combinations of three routes from the set of routes joining that particular OD pair. For example, if there are 5 routes joining an OD pair, then the total number of OD-routes that are obtained is ( 5 3 ) = 10. The total number of OD-routes in the network is 624,490.

Clustering of the OD-routes

Before clustering, the data was normalised so that all the variables describing the OD-routes have the same weight in determining the dissimilarity between observations; this step is necessary when the range of the variables are not comparable, as is the case in the OD-routes where the directness of the routes takes values in the interval (0, 1), but the length of the routes takes values in the interval (0, 35). The OD-routes are clustered using k-means with euclidean distance, determining the optimal number of clusters, k*, using the elbow method [26]. The idea behind the elbow method is to select the optimal number of clusters k*, such that the mean dissimilarity of the elements in the clusters does not decrease significantly with the k* + 1 clustering. The measure of dissimilarity of the elements in a cluster is the within-cluster sum of squares (WCSS), i.e., the sum of the square distance between the elements in a cluster. One of the OD-routes among the 1% nearest to the theoretical cluster centroid is selected as the cluster centroid. This is done because the theoretical centroid, i.e., the mean of the variables of the elements in the clusters, may not be part of the data.

Route choice experiments

The route choice experiments were carried out using a computer platform, the mobility decision game (MDG), that has been developed in the LICIT laboratory to investigate travellers’ decisions in transportation networks at large scale. The network description in the MDG is based on the full map of a real road network: the city of Lyon in France in our experiment. To generate the scenarios in which the choices are made, the MDG interacts with a single dynamic microscopic simulator, based on the LWR traffic model [31], which generates and handles all the trips that populate the network. To produce the simulated scenarios, the microscopic traffic simulator takes as an input the OD pairs in the network, their corresponding trip demands and the most important alternatives connecting the origins to the destinations. In the MDG, the participants access simultaneously to the experiment through a web interface, consisting of the map of the city of Lyon, France, road network. Each participant receives periodically new specific missions, that consist in travelling from an origin to a destination by choosing one of the three alternative routes that are proposed. The choices of the participants are considered as updated trip specifications by the microscopic simulator that runs in a central server. This alters the traffic conditions in the network. During a MDG session, multiple OD pairs are assigned to the participants, allowing to observe the choices of the same participants in different OD pairs. Furthermore, some of the participants receive traffic information as travel time estimates for the different route options in a given mission, allowing to assess the impact of travel time information in the decisions of travellers. Thus, the MDG permits to investigate the determinants of the participants’ decisions under different conditions (traffic conditions, traffic information and route characteristics).

Route choice model

Random utility models (RUMs) have been broadly used to understand and predict the route choice of travellers [7, 11, 12, 17, 18]. Joint RUMs arise in situations in which decisions of the same individuals are observed in several related choice problems, and correlation among their decisions is suspected. This is the case of surveys, where the answers of individuals to different questions may be correlated; or in route choice, with decisions of travellers in different OD pairs. A special case in joint RUMs is when the choice problems share part of their variables. In this situation, the coefficients of the shared variables in the model can be assumed to be equal across the choice problems. This problem is encountered when combining different data sources, as in [32], where the authors developed the techniques to jointly estimate a multinomial logit model (MNL) combining reveal preference and stated preference data to study the switching of mode of transportation of travellers. Other related works can be found in [3337]. When the choice problems share all of their variables, then the joint RUMs consists of an unique representation of the utility, given by the variables and their respective coefficients, which is equivalent to a panel data RUM. This is the case here, as several models are estimated, one for each OD pair. However, as the OD pairs are described by the same variables, the coefficients can be assumed to be equal for all the OD pairs. Therefore, the utility of the joint model is reduced to a single representation, and the model can be estimated as a panel data model.

The joint model for route choice, used in this study, is based on the mixed multinomial logit model (MXL) for panel data [2, 2729]. The MXL is a generalisation of the MNL in which the coefficients are assumed to be random, accounting for the heterogeneity in individuals’ preferences. Furthermore, estimation of MXLs are easily extended to account for observations of repeated choices by the same individuals, i.e., panel data. Formally, when decisions of individuals are observed in several choice situations, the utility that an individual i gets from alternative j in the MXL model is written as


where s = 1, …, Si indexes the choice situation in which the observation is made, εijs are i.i.d. Gumbel distributed, and βiFβ(b, Σ) is a vector of random coefficients. Estimating the MXL means obtaining the parameters b (mean) and Σ (covariance matrix) of the distribution of the β’s. The independence of the βi’s implies that the individuals are heterogeneous: their tastes vary following the distribution Fβ(μ, Σ). The correlation between the responses of the same individual to different choice situations, s and s′, is given by C o v ( U i j s , U i k s ′ ) = x i j s ′ Σ x i k s ′.

In Eq (1), the coefficients βi are indexed only by individual, i, and not by choice situation, s, which means that the tastes of the same individual do not vary between choice situations. This is equivalent to the panel data formulation of MXLs [2]. Specifically, in the route choice problem, travellers’ preferences towards the route attributes are equal for all routes regardless of the OD pair, but the preferences of two different travellers may be different. Also, in the panel data formulation in Eq (1), the alternatives j are the same in all the choice situations s, i.e., the choice sets are the same across choice situations. This is not true for route choices in several OD pairs, where alternatives are OD pair specific, i.e., j s ∈ C ( s ), where C ( s ) is the choice set in situation s. However, since the variables describing the alternatives are the same for all the choice sets, then the joint mixed logit model can be treated as a MXL for panel data. Conditioning on βi, the joint probability of individual i choosing the sequence of alternatives j 1 , j 2 , … , j S i is given by


where yijs = 1 when alternative j s ∈ C ( s ) is chosen and yijs = 0 otherwise; and Y i = ( y i j 1 , y i j 2 , … , y i j S i ) T. The first equality in Eq (2) is guaranteed by the independence of the εijs’s, and the second because they are identically Gumbel distributed.

Bayesian estimation of panel MXL models

In this work, Bayesian inference is used to estimate L(Yi|βi). The classical inference method to estimate the MXL requires integrating the expression in Eq (2) to obtain the unconditional probability L(Yi). Nevertheless, since the integral has no closed form, it needs to be numerically approximated, which could present convergence problems and be computationally expensive. An alternative approach to estimate MXLs, is to regard them as Bayesian hierarchical models, which have the advantage of avoiding the numerical multiple integration [2, 38, 39].

In Bayesian methods, the parameters of the model are assumed to be random variables rather than fixed values. Inference, in this context, refers to obtaining the joint distribution of the parameters that best fits the data. To estimate the joint distribution of the parameters, first, a prior distribution, h, representing the researchers’ beliefs over the values of the parameter, is defined. Then, when data becomes available, the prior is updated through the likelihood function to obtain the posterior distribution, H. As a result of the Bayes’ theorem, the posterior distribution is proportional to the prior multiplied by the likelihood. In the general case, the posterior distribution H of the parameters is


where X represents the alternative and individuals’ attributes; Y the observed choices and ϕn is the multivariate normal density function of the random coefficients parametrised by b (mean) and Σ (covariance matrix). The expression in brackets is the likelihood of the observed choices and h is the joint prior distribution of the model’s parameters. The joint priors, h, for the three MXL models estimated in this work, M1, M2 and M3 (see the Results section for the specification of the models), are, respectively,

where ϕ is the density function of the normal distribution, ϕn of the n-variate normal distribution, fIG the density of the inverse-Gamma distribution and fIW of the inverse-Wishart distribution. The inverse-Gamma is the conjugate prior for the variance of the normal distribution, and the inverse-Wishart its generalisation for the multivariate case.

The right hand side in Eq (3) has no closed form, however samples from the joint posterior distribution H can be obtained using the Gibbs sampling method [40]. In this study, the Gibbs sampler software JAGS [41] and the R [42] package rjags were used to obtain 10000 samples of the posterior distribution H after a burn-in period of 20000 samples. The values of the hyperparameters μ0, σ 0 2, r0, λ0, Σ0, I0 and k0, which define the priors, were chosen to be weakly-informative (very high variances). In other words, it is assumed high uncertainty on the real values of the parameters that are being estimated. They are shown in Table 1.

Tab. 1. Hyperparameters of the prior distribution h.
Hyperparameters of the prior distribution <i>h</i>.

Results

Clustering

To determine the optimal number of clusters, k-means algorithm, with k = 1, …, 30, was performed over the 624,490 OD-routes in the network. The mean within-cluster sum of squares (WCSS) is plotted against the number of clusters k in Fig 3. In the results, the optimal number of clusters is not clear, according to the elbow method: big improvements happen for the first values of k (k ≤ 4); for values 5 ≤ k ≤ 9 the improvement is mediocre; and for k ≥ 10 the improvements are rather small. In terms of the purpose of this article, choosing a small number of k has the risk of sub-representing the OD-routes in the network and, more important, a small number of OD-routes in the route choice experiments implies that the variability in the route attributes is also small, posing a problem in estimating a route choice model (overfitting). In this sense, choosing high values of k is preferable, even if some of the clusters are similar. However, the needed number of observations in the route choice experiment increases with the number of OD-routes, implying higher costs in the organisation of the experiments, not to mention the difficulties to recruit participants. In view of these limitations, the number of clusters is set to k = 9. The clustering results are presented in Table 2 and the centroids are depicted in the map shown in Fig 4.

Determining the number of clusters <i>k</i>.
Fig. 3. Determining the number of clusters k.
The sum of squared errors for k-means clustering of the 624,490 OD-routes with k = 1, …, 30. After k = 9 the decrease in the mean WCSS is marginal.
Selected OD-routes for the route choice experiments.
Fig. 4. Selected OD-routes for the route choice experiments.
The geodata used to render the plot is from ©OpenStreetMap contributors, licensed under the Open Data Commons Open Database License (ODbL). The maps were rendered using The R Project for Statistical Computing which is distributed as Free Software under the terms of the Free Software Foundation’s GNU General Public License.
Tab. 2. Cluster analyis results.
Cluster analyis results.

With k = 9, the variability of the full set of OD-routes is reduced in 46.4%. If well, this reduction may not be big in terms of clustering analysis, it can be seen (Fig 5) that the road attributes of the cluster centroids cover likely values to be observed in the network. To be more specific, 83% of the values of the attributes of the OD-routes in the network lie in the range of the centroids: 83% for the euclidean distance and the directness, 89% for the freeway composition, 85% for the number of turns per kilometre and 90% for the route length. Furthermore, the resulting p-values of the two-sample Kolmogorov-Smirnov test (Fig 5) are high: pvalue > 0.1 for the five variables, suggesting that there is not enough statistical evidence (with a significance level of α = 0.1) to reject the null hypothesis that the values of the attributes of the centroids and the full network come from the same distribution. This implies that a random selected OD pair or route in the network is likely to have attributes similar to one of the nine OD-routes used in the route choice experiments. In this sense, the nine cluster centroids can be regarded as representative of the network.

Cluster centroids as representative OD-routes.
Fig. 5. Cluster centroids as representative OD-routes.
The distribution of the attributes of the selected OD pairs are similar to that of the whole network. The p-values of the Kolmogorov-Smirnov, presented in red in the top of each panel, indicate the lack of statistical evidence (with a confidence level of 0.90) to reject the hypothesis that the two distributions are the same.

A further characterisation of the nine clusters, based on their elements’ attributes (see S1 and S2 Figs), is proposed as follows:

  • Clust. C1: Medium-range direct trips going from south to north, with routes having small number of turns per kilometre and some freeway segments.

  • Clust. C2: Short non direct trips mainly in the central part of the network, with routes having a lot of turns per kilometre and no freeway segments.

  • Clust. C3: Medium-range direct trips going from north to south, with routes having small number of turns per kilometre and some freeway segments.

  • Clust. C4: Medium-range direct trips mostly in the central part of the network, with routes having average number of turns per kilometre and with longest route highly composed of freeway segments.

  • Clust. C5: Long trips going from east to west, with routes having a small number of turns and with large portions of freeway.

  • Clust. C6: Long trips going from west to east, with routes having a small number of turns and with large portions of freeway.

  • Clust. C7: Medium-range non direct trips in the central part of the network, with routes with average number of turns per kilometre and high portions of freeway.

  • Clust. C8: Short direct trips mainly in the central part of the network with routes with low number of turns per kilometre (among short trips) and no freeway segments.

  • Clust. C9: Short non direct trips mainly in the central part of the network, with routes with low number of turns per kilometre (among short trips) and some freeway segments.

Experiment results

Three route choice experiment sessions were carried out using the nine OD pairs and routes obtained from the clustering analysis of the network. In total, 73 individuals participated in the three sessions, from these participants, 56 (77%) received estimates of the travel times in each route. Participants recorded a total number of 802 choices in the nine defined OD-routes, with an average number of 11 choices per participant, and an average number of 89 choices in each OD-routes. The choices of the participants are presented in Fig 6, where it can be immediately noticed that travel time information changes the behaviour of the participants.

Route choice distribution in the nine cluster centroids.
Fig. 6. Route choice distribution in the nine cluster centroids.
The choices of the informed participants are different from those of the not informed participants.

Route choice model estimation

Three panel data mixed logit models are estimated using the observations collected in the route choice experiments. Five variables are used in the specification of the models. Four of these variables correspond to the variables used in the selection of the representative OD-routes for the route choice experiments, which will help to test if the choices in the representative OD-routes (cluster centroids) can approximate the choices in other OD-routes. The fifth variable is the estimated travel time that the participants received during the experiments. These variables are known to influence the route choice behaviour of travellers and that can be observed by the participants in the computer route choice experiments.

Let the individuals and alternatives be indexed by i and j, respectively. Since participants were allowed to repeat decisions in the same OD pair, the choice situation, indexed by s, represents the pair (od, t), where od is the OD pair in which the decision was made and t indexes the moment of the choice. The explanatory variables considered in the model are

  • FRWj, the percentage of freeway that composes the route j;

  • DIRj, the directness of the route j, defined as the length of j divided by the euclidean distance between origin and destination;

  • TNRj, the number of turns per kilometre in the route j;

  • ITTjs, the informed travel time in the route j in OD pair and moment s, the variable is normalised by OD pair by dividing the informed travel time by the free flow travel time in the fastest of the three routes;

  • INFi, binary variable indicating if participant i received information; and

  • LENj, the length (in km) of the route j.

The specifications of the three models M1, M2 and M3 are




In models M1 and M2, the coefficients βip for p = 1, ‥, 6 are independent and normally distributed, i.e., β i p ∼ N ( b p , σ p 2 ). In model M3, the coefficient β6 is fixed for all individuals (not random), and the coefficients βip are correlated for p = 1, ‥, 4, i.e., βi.N4(b, Σ), but independent from β i 5 ∼ N ( b 5 , σ 5 2 ). Model M1 is the simplest MXL model considering the five variables. In models M2 and M3 the interactions between the route length and the travel time information are taken into account, allowing for the preference towards the length of the route to change depending on the informed travel time. In model M3 the correlations between the coefficients βip for p = 1, …, 4 are also estimated. In MXL models, the parameters that are estimated are the means and variances (covariances) of the coefficients’ distributions, b p ^, σ ^ and Σ ^. The estimated parameters for the three models are shown in Table 3; more detailed result of the posterior distribution of the parameters can be found in S1 Table, and the details of the computational effort for the estimation process in S2 Table.

Tab. 3. MXL models estimation results.
MXL models estimation results.

The estimated parameters b p ^ represent the mean preferences in the population. The positive sign of the estimates b 1 ^ and b 2 ^ in the three models is interpreted as the average traveller prefers routes with high composition of freeways, and direct routes. On the contrary, the negative signs of b 3 ^ and b 5 ^ mean that the average traveller avoids routes with many turns and higher travel times. These results are in line with the findings in [24], and provide more evidence in favour of travel time as the most important variable in route choice. Note that in the three models b 3 ^ ≈ 0, but with large standard deviations σ 3 ^, meaning that (i) the sign is positive for a large number of participants (near half), and that (ii) even when the mean of the coefficient is close to zero, this variable is still important for a large percentage of the participants, specially in model M3, where Pr(|βi3| > 1) = 0.42. The case for the length of the route is different, as the standard deviations are smaller: for models M1 and M2 this implies that the length of the route is not important for the majority of the participants, Pr(|βi4| < 0.2) = 0.93 and Pr(|βi4| < 0.2) = 0.81, respectively; but for model M3 it is, Pr(|βi4| < 0.2) = 0.29. Finally, note that in models M2 and M3 the mean preference for the length of the routes can be written as (b4 + b6ITTjs INFi), with b4 < 0 and b6 > 0, meaning that the informed travel time diminishes the preference for shorter routes.

Choices on representative OD-routes

Until now, the discussion on the representativeness of the nine selected OD-routes (the cluster centroids) has been in terms of the route attributes. In this section, the representativeness of the OD-routes is assessed in terms of how well a choice model, estimated using the cluster centroids, can be generalised to the entire road network or, in other words, how well it scales-up the travellers’ choices to other OD pairs in the network. The hypothesis is that if the choices in the nine cluster centroids are representative of the choices in the entire network, then the predictive accuracy of a model, estimated with observations in the nine cluster centroids, should be higher than the predictive accuracy of models (with the same specification) estimated with observations in random sets of OD-routes. To this end, data collected in other route choice experiments carried out with the MDG platform is used. The data consists of route choice observations in 21 OD-routes, defined with a different methodology for previous experiments, and not comprising the representative OD pairs.

The methodology to validate the representative OD-routes is based on bootstrapping for model validation: at each step, a random part of the data is left-out of the estimation process, and then used to measure the predictive accuracy of the model. However, in this case, the predictive accuracy of the models obtained at each iteration are compared to the predictive accuracy of the model estimated with the nine cluster centroids. Let C be the set of choice observations in the 9 cluster centroids and T the set of observations in the 21 test OD-routes. Denote by M* the model (it can be either M1, M2 or M3) estimated with observations on the nine cluster centroids, C. At iteration r, r = 1, …, 40,

  • obtain Tr ⊂ (CT), composed of all the observations from nine randomly sampled OD-routes;

  • estimate the model Mr (M1, M2 or M3) with the observations in Tr;

  • for each OD-routes od ∈ (CT) − Tr, compute the prediction error of models M* and Mr, i.e., er(M*, od) and er(Mr, od), where


  • obtain the mean prediction error for iteration r, defined as


    where the weight wod is the percentage of OD-routes in the cluster to which od belongs, multiplied by the inverse of the number of OD-routes in (CT) that belong to that cluster. The weighting is done to adjust for the probability of observing an OD-routes in the network like od. This follows since some clusters are over-represented in T, as the OD-routes in T were not randomly selected from the network, but they were selected following a different methodology in previous studies.

Note that the error measure, er(M*, od) has a direct interpretation in terms of traffic assignment: the percentage of trips that are wrongly distributed amongst the three alternative routes.

The MPEr(M*) is compared against the MPEr(Mr), r = 1, …, 40 for the three model specifications. In Fig 7, MPEr(M*) is plotted against MPEr(Mr), with blue dots when MPEr(M*) ≤ MPEr(Mr), and red otherwise. The models estimated with the clusters’ centroids performed better in predicting the choices of travellers than most of the models estimated with randomly selected OD-routes. To be more specific, MPEr(M*) ≤ MPEr(Mr) in 35 out of 40 cases (87.5%) for models M1 and M2, and in 31 cases (77.5%) for model M3. Furthermore, in the cases when the model estimated with the centroids performed worst, i.e., MPEr(M*) ≤ MPEr(Mr), the errors were close to those of the models estimated with randomly selected OD-routes. Define the improvement of M* with respect to Mr as αr = (MPEr(Mr) − MPEr(M*))/MPEr(Mr). Then, the mean improvements, α ¯, are 14% for model M1, 14.5% for model M2 and 9.9% for model M3. In 20% of the test cases, αr is at least 26%, 25% and 22% for models M1, M2 and M3, respectively; and αr reaches 52%, 48% and 43% in the worst case scenarios. This results highlights the importance of a careful selection of OD pairs in route choice model estimation. As the MPE represents the percentage of trips that are not assigned to the right route, and since the total number of trips at a city level can be very high, about 1 million in the Lyon Metropolis during one day, even low αr values may have an impact on how the traffic is distributed on the network.

Mean predictive errors.
Fig. 7. Mean predictive errors.
The MPEs of model M* are smaller than the MPEs of models Mr in the majority of the cases (blue dots). Furthermore, in the cases where the MPEs of models M* are bigger (red dots), the differences are small (close to identity line).

If the MPE is analysed by whether or not the participants received travel time information (Fig 8), it can be seen that the models M* are better than the models Mr for the not informed participants than for the informed ones. For models M1 and M3, MPEr(M*) ≤ MPEr(Mr) in 97.5% of the cases, and for model M2 in 95%; and when the participants were informed, MPEr(M*) ≤ MPEr(Mr) in 65% for model M1, 77.5% for model M2 and 67.5% for M3. In the case of the not informed participants the values of α ¯ are 20%, 14% and 15%, respectively for models M1, M2 and M3. The high performance of the models M* for the not informed participants implies that the models are capable of approximating the choices of this group in a variety of scenarios, i.e., the models estimated with the nine centroids generalise well to other OD-routes for this group. Moreover, considering that the informed travel time was not part of the variables used in the clustering of the OD-routes (% of freeway, directness, no. of turn per kilometre, distance), this result suggests that the choices in the cluster centroids are representative of the choices in the entire network, thus validating the methodology proposed in this article.

Mean predictive errors by information group.
Fig. 8. Mean predictive errors by information group.
The models estimated with the cluster centroids are clearly better in predicting the choices for the not informed participants.

The predictive errors of the representative models, M*, and the test models, Mr, can be disaggregated by OD-routes. In the results, shown in Fig 9, it is clear that the magnitude and the variance of the predictive errors depend on the OD-routes where the choices are being predicted. The choices in some OD-routes are difficult to predict, regardless of the training set used to estimate the models. Furthermore, there is no clear pattern indicating that these errors are associated with the road characteristics of the OD-routes: two OD-routes belonging to the same cluster, i.e., having similar route attributes, may have a low and a high prediction error. Such is the case of OD-routes c3_od2 and c3_od3, both belonging to cluster C3, but with errors below 0.1 for the former and above 0.2 for the later. Similar cases can be found in cluster C5 and C8. An important observation is that the models estimated with the cluster centroids M* are not as accurate in predicting the choices in individual OD-routes as the models Mr, for some values of r. In fact, their prediction errors are amongst the lowest 25% in only 8 out of 21 test OD-routes for models M1 and M3, and in 6 for model M2. However, at the same time, the individual errors are almost never amongst the highest 75%: in 0 OD-routes for model M1, in 2 for model M2, and in 1 for M3. Moreover, when the individual errors are averaged to obtain the MPE (as in the previous analysis), the models M* outperform the models Mr for the majority of values of r. This result implies that a model M r 0 having low prediction errors for some OD-routes has also high prediction errors in other OD-routes, and therefore its mean predictive accuracy is reduced. In this sense, the models estimated with the cluster centroids, M*, are preferred, as they will show a relative better global prediction accuracy without incurring in large errors in individual OD-routes.

Distributions of the predictive errors the 21 validation OD-routes.
Fig. 9. Distributions of the predictive errors the 21 validation OD-routes.
The level and the variability of the errors amongst the different OD-routes imply that the choices in some OD-routes are difficult to predict, regardless of the training set used to estimate the models.

The models estimated with the representative OD-routes, M*, are compared in terms of their prediction errors over the 21 validation OD-routes. The error distribution of the three models, depicted in Fig 10, show that, practically, there is no difference in the predictive accuracy. This means that the interaction between the informed travel time and length of the route in models M2 and M3 does not improve the predictive accuracy; nor considering the correlations in model M3 does.

Distributions of the prediction errors of the models <i>M</i>* on the 21 validation OD-routes.
Fig. 10. Distributions of the prediction errors of the models M* on the 21 validation OD-routes.
There are no significant differences between the error distributions.

Discussion

In this study, it was demonstrated that the choices of participants in a route choice experiment over a small but representative set of OD configurations can be scaled-up to the entire network. To obtain the set of representative OD configurations, a new methodology based on k-means cluster analysis is proposed. First, the OD configurations in the network, i.e., the OD pairs and three connecting routes, are represented in vector form according to the attributes of the OD pairs and routes. Then, these OD configurations are clustered in order to obtain a partition of the road network and the cluster centroids selected as representative of the entire network. The main hypothesis is that the choices of travellers over the entire network can be approximated with route choice models estimated using data collected for the representative set. The obtained results point in this direction.

In the current study, for the city of Lyon in France, 9 OD pairs and their connecting routes were used as representative of 624,490 OD configurations. These nine representative OD configurations cover around 83% of the values of the attributes of the OD-routes in the network. The predictions of the models estimated with the representative set were superior in most of the test cases (87.5% and 77.5% in the general case). For the not informed participants, whose decisions were based on the same attributes used in the clustering, the predictions are better in at least 95% of the test cases. By estimating the route choice model with the cluster centroids, the mean prediction errors are reduced by up to 14.5% for model M1 (similar results are observed for models M2 and M3). The reduction of the prediction error is more than 22% for the 20% of the test cases, and it goes up to 51% in the worst case. This demonstrates that a careful selection of the OD configurations significantly improves the prediction accuracy, independently of the model specification. Another significant finding, is that the models estimated with the representative OD configurations are more robust than the ones obtained from the models with random OD configurations. The models estimated with the representative set never show extreme errors for individual OD pairs, contrary to the models estimated with random sets of OD configurations. This implies that the models estimated with the representative set will show a relative better global prediction accuracy without incurring in large errors on individual OD-routes. This result is important when predicting the trip distribution over the network, as high errors in individual OD pairs may have significant impact in local traffic conditions, causing spreading.

The last finding is that estimating the models with the representative OD pairs leads to an average prediction error of 12.7%. This value can be considered quite low when considering the scale of the city, the heterogeneity of OD configurations, and the actual performance of user equilibrium approaches.

From the clustering analysis in this study, it is clear that there are OD pairs in the network that are not well represented by the representative set of nine OD configurations. Therefore, it cannot be claimed that the choices in these non-represented OD pairs can be well approximated by the set of nine OD pairs found in this study. However, these non-represented OD pairs are those with attributes not covered by the representative set, which are no more than 17% of the OD configurations in the network. Note that this result does not hinder the usefulness of the proposed methodology, as it can be extended by either using other clustering techniques that allow taking into account for these atypical OD configurations or by including more clusters in the representative set.

Supporting information

S1 Fig [tif]
Distribution of the route attributes in the nine clusters.

S2 Fig [tif]
Origins and destinations in the network.

S1 Table [pdf]
Complete estimation results of the MXL models.

S2 Table [pdf]
Computational effort in the estimation of the MXL models.


Zdroje

1. Manski C, McFadden D. Structural Analysis of Discrete Data with Econometric Applications. The MIT Press; 1981.

2. Train KE. Discrete Choice Methods with Simulation. Cambridge University Press; 2003.

3. Walker J, Ben-Akiva M. Generalized random utility model. Mathematical Social Sciences. 2002;43(3):303–343. doi: 10.1016/S0165-4896(02)00023-9

4. Wardrop JG. Road paper. Some theoretical aspects of road traffic research. Proceedings of the Institution of Civil Engineers. 1952;1(3):325–362. doi: 10.1680/ipeds.1952.11259

5. Yildirimoglu M, Kahraman O. Searching for empirical evidence on traffic equilibrium. PLoS ONE. 2018;13(5):1–16. doi: 10.1371/journal.pone.0196997

6. Iida Y, Akiyama T, Uchida T. Experimental analysis of dynamic route choice behavior. Transp Res B. 1992;26(1):17–32. doi: 10.1016/0191-2615(92)90017-Q

7. Bogers EaI. Joint modeling of ATIS, habit and learning impacts on route choice by laboratory simulator experiments. Delft University of Technology; 2005.

8. Selten R, Chmura T, Pitz T, Kube S, Schreckenberg M. Commuters route choice behaviour. Games and Economic Behavior. 2007;58(2):394–406. doi: 10.1016/j.geb.2006.03.012

9. Adler JL, McNally MG. In-laboratory experiments to investigate driver behavior under advanced traveler information systems. Transportation Research Part C. 1994;2(3):149–164. doi: 10.1016/0968-090X(94)90006-X

10. Lotan T. Effects of familiarity on route choice behavior in the presence of information. Transportation Research Part C: Emerging Technologies. 1997;5(3-4):225–243. doi: 10.1016/S0968-090X(96)00028-9

11. Mahmassani HS, Liu YH. Dynamics of commuting decision behaviour under Advanced Traveller Information Systems. Transportation Research Part C: Emerging Technologies. 1999;7(2-3):91–107. doi: 10.1016/S0968-090X(99)00014-5

12. Ben-Elia E, Shiftan Y. Which road do I take? A learning-based model of route-choice behavior with real-time information. Transportation Research Part A: Policy and Practice. 2010;44(4):249–264.

13. Ben-Elia E, Avineri E. Response to Travel Information: A Behavioural Review. Transport Reviews. 2015;35(3):352–377. doi: 10.1080/01441647.2015.1015471

14. Abdel-Aty MA, Kitamura R, Jovanis PP. Using stated preference data for studying the effect of advanced traffic information on drivers’ route choice. Transportation Research Part C: Emerging Technologies. 1997;5(1):39–50. doi: 10.1016/S0968-090X(96)00023-X

15. Srinivasan K, Mahmassani H. Modeling Inertia and Compliance Mechanisms in Route Choice Behavior Under Real-Time Information. Transportation Research Record: Journal of the Transportation Research Board. 2000;1725(January):45–53. doi: 10.3141/1725-07

16. De Moraes Ramos G, Daamen W, Hoogendoorn S. Modelling travellers’ heterogeneous route choice behaviour as prospect maximizers. Journal of Choice Modelling. 2013;6:17–33. doi: 10.1016/j.jocm.2013.04.002

17. Avineri E, Prashker JN. Sensitivity to travel time variability: Travelers learning perspective. Transportation Research Part C: Emerging Technologies. 2005;13(2):157–183. doi: 10.1016/j.trc.2005.04.006

18. de Palma A, Picard N. Route choice decision under travel time uncertainty. Transportation Research Part A: Policy and Practice. 2005;39(4 SPEC. ISSS.):295–324.

19. Sheffi Y. Urban transportation networks. Prentice-Hall, Inc.; 1985.

20. Mahmassani HS, Chang GL. On Boundedly Rational User Equilibrium in Transportation Systems. Transportation Science. 1987;21(2):89–99. doi: 10.1287/trsc.21.2.89

21. Bovy PHL, Stern E. Route Choice: Wayfinding in Transport Networks. Kluwer Academic Publishers; 1990.

22. Ramming MS. Network knowledge and route choice. Massachusetts Institute of Technology; 2002.

23. Bekhor S, Ben-Akiva ME, Ramming MS. Evaluation of choice set generation algorithms for route choice models. Annals of Operations Research. 2006;144(1):235–247. doi: 10.1007/s10479-006-0009-8

24. Papinski D, Scott DM, Doherty ST. Exploring the route choice decision-making process: A comparison of planned and observed routes obtained using person-based GPS. Transportation Research Part F: Traffic Psychology and Behaviour. 2009;12(4):347–358. doi: 10.1016/j.trf.2009.04.001

25. Zhu S, Levinson D. Do people use the shortest path? An empirical test of wardrop’s first principle. PLoS ONE. 2015;10(8):1–18. doi: 10.1371/journal.pone.0134322

26. Hastie T, Tibshirani R, Friedman J. The Elements of Statistical Learning. vol. 1 of Springer Series in Statistics. 2nd ed. New York, NY: Springer New York; 2009. Available from: http://www.springerlink.com/index/10.1007/b94608.

27. McFadden D, Train K. Mixed MNL Models for Discrete Response. Journal of Applied Econometrics. 2000;15(5):447–470. doi: 10.1002/1099-1255(200009/10)15:5%3C447::AID-JAE570%3E3.0.CO;2-1

28. Bhat CR, Castelar S. A unified mixed logit framework for modeling revealed and stated preferences: formulation and application to congestion pricing analysis in the San Francisco Bay area. Transportation Research Part B: Methodological. 2002;36(7):593–616. doi: 10.1016/S0191-2615(01)00020-0

29. Brownstone D, Bunch DS, Train K. Joint mixed logit models of stated and revealed preferences for alternative-fuel vehicles. Transportation Research Part B: Methodological. 2000;34(5):315–338. doi: 10.1016/S0191-2615(99)00031-4

30. Institut national de la statistique et des études économiques. Découpage infracommunal: Table d’appartenance géographique des IRIS; 2018. Available from: https://www.insee.fr/fr/information/2017499.

31. Leclercq L. Hybrid approaches to the solutions of the “Lighthill-Whitham-Richards” model. Transportation Research Part B: Methodological. 2007;41(7):701–709. doi: 10.1016/j.trb.2006.11.004

32. Ben-Akiva M, Morikawa T. Estimation of switching models from revealed preferences and stated intentions. Transportation Research Part A: General. 1990;24(6):485–495. doi: 10.1016/0191-2607(90)90037-7

33. Bradley M, Daly A. Estimation of logit choice models using mixed stated preference and revealed preference information. Paper presented to the 6th International Conference on Travel Behavior, Quebec. 1991.

34. Earnhart D. Combining Revealed and Stated Data to Examine Housing Decisions Using Discrete Choice Analysis. Journal of Urban Economics. 2002;51(1):143–169. doi: 10.1006/juec.2001.2241

35. Adamowicz W, Boxall P, Williams M, Louviere J. Stated Preference Approaches for Measuring Passive Use Values: Choice Experiments and Contingent Valuation. American Journal of Agricultural Economics. 1998;80(1):64–75. doi: 10.2307/3180269

36. Adamowicz W, Louviere J, Williams M. Combining Revealed and Stated Preference Methods for Valuing Environmental Amenities. Journal of Environmental Economics and Management. 1994;26(3):271–292. doi: 10.1006/jeem.1994.1017

37. Hensher DA, Bradley M. Using stated response choice data to enrich revealed preference discrete choice models. Marketing Letters. 1993;4(2):139–151. doi: 10.1007/BF00994072

38. Regier DA, Ryan M, Phimister E, Marra CA. Bayesian and classical estimation of mixed logit: An application to genetic testing. Journal of Health Economics. 2009;28(3):598–610. doi: 10.1016/j.jhealeco.2008.11.003 19345433

39. Balcombe K, Chalak A, Fraser I. Model selection for the mixed logit with Bayesian estimation. Journal of Environmental Economics and Management. 2009;57(2):226–237. doi: 10.1016/j.jeem.2008.06.001

40. Levin DA, Peres Y. Markov Chains and Mixing Times. 2nd ed. American Mathematical Society; 2017.

41. Plummer M. JAGS: A program for analysis of Bayesian graphical models using Gibbs sampling JAGS: Just Another Gibbs Sampler. DSC 2003 Working Papers (Draft Versions). 2003; p. 1–8.

42. R Core Team. R: A Language and Environment for Statistical Computing; 2018. Available from: https://www.R-project.org/.


Článek vyšel v časopise

PLOS One


2019 Číslo 11
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#