Next Article in Journal
Rainfall Runoff Mitigation by Retrofitted Permeable Pavement in an Urban Area
Previous Article in Journal
Inter-Basin Water Transfer Green Supply Chain Equilibrium and Coordination under Social Welfare Maximization
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Transport Choice Modeling for the Evaluation of New Transport Policies

1
DeustoTech-Fundación Deusto, Avda. Universidades, 24, 48007 Bilbao, Spain
2
Facultad Ingeniería, Universidad de Deusto, Avda. Universidades, 24, 48007 Bilbao, Spain
*
Author to whom correspondence should be addressed.
Sustainability 2018, 10(4), 1230; https://doi.org/10.3390/su10041230
Submission received: 1 March 2018 / Revised: 6 April 2018 / Accepted: 12 April 2018 / Published: 17 April 2018
(This article belongs to the Section Sustainable Transportation)

Abstract

:
Quantifying the impact of the application of sustainable transport policies is essential in order to mitigate effects of greenhouse gas emissions produced by the transport sector. One of the most common approaches used for this purpose is that of traffic modelling and simulation, which consists of emulating the operation of an entire road network. This article presents the results of fitting 8 well known data science methods for transport choice modelling, the area in which more research is needed. The models have been trained with information from Biscay province in Spain in order to match as many of its commuters as possible. Results show that the best models correctly forecast more than 51% of the trips recorded. Finally, the results have been validated with a second data set from the Silesian Voivodeship in Poland, showing that all models indeed maintain their forecasting ability.

1. Introduction

Transport systems performance evaluation and forecasting is a difficult task currently under consideration by many research centers. This complexity emerges mainly due to its stochastic nature, whereby a large number of variables and personal preferences come into play [1]. Transport infrastructures are highly expensive and if the solution does not fit the necessities of the citizens, the infrastructure might need to be upgraded, increasing the cost or, even worse, becoming a worthless investment [2]. Therefore, a cheap and flexible method for testing future transport developments is to use computer simulations supported by advanced mathematical models [3].
To calculate the impact of the transport sector, traditional approaches are relying on simulations, a paradigm that has become a popular and effective technique to analyze a wide range of dynamically changing systems. The essence of simulating is to carry out a series of computational experiments using a model that describes the operations of the real system (how drivers, vehicles, signals, pedestrians, etc. interact) and to generate output data that characterizes the subject system [4]. Traffic simulation allows one to check, monitor, and evaluate the behavior of the real systems under different realistic conditions in an artificial-computer-based environment [5]. Most common traffic simulations focus on simulating everyday traffic activity in search of congestions, assessing different ways of enabling traffic in intersections (crossings, roundabouts, traffic lights, etc.) or applying dynamic speed limits to avoid jams.
Over the last decades, many universal approaches have been widely used for modeling stochastic processes, such as Monte Carlo techniques [6], Markov chain models [7], or queuing theory methods [8]. Today’s state-of-the-art simulation packages (listed later) add features such as graphical interfaces and object-oriented programming and even allow one to simulate both continuous and discrete events. Regarding transport simulation, two mutually interrelated pieces of information are needed [9]: the demand that characterizes the need for movement (passengers or freight) and a description of the transport network, traffic zones, and vehicles. By combining these data sets, state-of-the-art models follow the standard concept of the four-step model (FSM) [10]. FSM is a primary tool for the analysis of transport systems, assessment of their performance, and prediction of their future behavior. This model is composed of the following phases:
  • Trip generation: Generate a specific, whether fixed or approximated, amount of vehicles that represents the population.
  • Trip distribution: Spread the population along the scenario in a realistic distribution.
  • Transport choice: Schedule their transport mode and itineraries throughout the city to reach their destination.
  • Traffic simulation: Estimate the global impact derived from the configuration of a single scenario and compare the results with those obtained from others.
In this sense, multi-agent systems (MASs) have been an active area of research in the last three decades [11]. What started as isolated systems for simple problem-solving has turned into full-scale simulation platforms that are able to integrate several state-of-the-art technologies. Specially, for transport systems evaluation, agent-based simulation models have been proven to be very useful [12]. In these models, each vehicle is represented as a software object, a so-called intelligent agent, that has its own set of parameters and attributes that define how the agent will behave, interact, make decisions, and communicate with the other agents in the system [13]. This approach allows one to build a thoroughly detailed description of the model population. According to their accuracy and scope, simulation models can be classified as follows:
  • Micro-simulation models: These models describe traffic with a high level of detail and distinguish separate elements in the traffic flow, such as types of vehicles and pedestrians. Though the use of this high level of detail entails very precise analysis, it must be limited to a small area or intersection due to its complexity. The most popular micro-simulation tools are PTV Vissim [14], Aimsun [15], and Corsim [16].
  • Meso-simulation models: These models describe traffic with an intermediate level of detail, distinguishing separate elements in the traffic flow but not taking into account the interactions between them. They are less precise and can be applied to cover larger areas than those of the micro-simulation. The most popular meso-simulation tools are Dynasmart [17] and Transims [18].
  • Macro-simulation models: These models describe traffic with a high level of aggregation, as a uniform traffic flow. They are based on deterministic relationships between the parameters characterizing the traffic flow, such as volume, speed, or density. Macroscopic simulation has been developed to model an entire transport network and/or system. The most popular macro-simulation tools are Emme/2 [19], PTV Visum [20], and Transcad [21].
A great amount of research effort has been put into understanding transport mode choice modeling and how the travelers’ background, trips and transport facilities influence the attractiveness of transport [22]. The evaluation of transport systems depends completely on how people adapt to it. This implies modeling social preferences that have particular importance when traveling, such as time, cost, comfort, or environment-friendly awareness. Even as the aforementioned simulators are relatively easy to run and have many proven success stories, most focus primarily on the calculation and estimation of routes and congestion avoidance, leaving out the social behavior that influences the traveler’s choice of using one transport method or another.
For those simulators that model transport selection, discrete choice regression models (probit and logit models) [23] are the most widespread techniques [24]. In fact, the three simulators that include transport choice models (Aissum, TRANSIMS, and TransCAD) implement these types of models. However, the addition of new transport modes, the alteration of the transport network, or the evaluation of the application of new incentives and restrictions requires calibrating (and sometimes building) the entire model again.
Within this scope, this article compares eight methodologies to model human decision making when choosing between transport modes. The objective is to emphasize the models that perform better and assess the difficulty of extending them to assess the introduction of new policies or transport methods. All models will be trained with a real data set, and their results will be validated with another data set. This process is of great importance for analyzing the impacts that sustainable transport policies will have when applied over the specific geographical and demographic characteristics of a certain area.
The rest of the paper is organized as follows. Section 2 describes the materials used and details of the proposed transport choice methodologies. Section 3 presents and analyzes the results derived from the application of these methodologies to two different scenarios. The final section presents the conclusions and discussion that emerge from the analysis of the results.

2. Materials and Methods

In order to better understand the methodology presented, the next subsections will detail its implementation according to the stages defined in the FSM.

2.1. Trip Generation and Distribution

In order to generate the information to train the models, agent-based simulations were extensively used. The main characters of the simulation are the citizens, where each citizen agent has a set of personal preferences that will directly influence the transport choice for traveling. For the generation and distribution of the trips for the experimentation, two data sets were used. On the one hand, Biscay’s commute data (depicted in Figure 1a and Figure 2a) was extracted from the census performed every 10 years by the Spanish National Statistics Institute [25]. This data set covers the region of Biscay, in the Basque Country, with an area of 2217 km2. The census also details the points of origin and destination of traveling for the most populated 11 municipalities of the region, with a total population of about 1,100,000 citizens. This data set will be used to train the different models. On the other hand, Silesia commutes (depicted in Figure 1b and Figure 2b) were extracted from the green traveling project’s surveys [26]. The data set covers 19 municipalities of the central part of Silesian Voivodeship with a population of 4,710,000 citizens. This data set will be used to validate the results.
When the data set is brought to life by the simulation, each citizen agent registers its point of origin, point of destination, departure time, and the itinerary for the five transport modes. The platform preprocesses the input data to geolocate each citizen agent on the environment through a weighted distribution. This process assigns the origin of each citizen agent to household buildings in the municipality and the destination to available amenities, commercial, or industrial buildings. The weighted distributions are based on the buildings’ size and levels, i.e., the bigger and the taller the building, the higher the probability of assigning a citizen agent to it. Once located, citizens will calculate the five itineraries for CAR, MOTORCYCLE, TRANSIT, BICYCLE, and WALK. For each itinerary, the main characteristics of every choice (see Section 2.2 for details) will be stored in order to build their decision process. Some itineraries might be null due to a lack of public transport stop nearby or an inability to calculate the route.

2.2. Transport Mode Choice

There are many features that determine the modal split of a territory. More than 45 years ago, the authors of [27] identified the main features that influence commuters’ choices, and since then other important characteristics have been added to the list [28]. In spite of using different techniques, researchers agree that itinerary dependent variables can be grouped into four main categories that represent the physical and factual features of a trip [29]: duration, price, length, and environmental impact. Table 1 describes these categories and shows how they are usually further divided in other second-level features. Nonetheless, these four variables do not cover some other features listed in the literature, which are directly linked to the traveler’s background. These background variables, presented in Table 2 do not appear in the census nor in the surveys and are therefore beyond the scope of this study (see Section 4 for details).
Missing values have been the main source of problem. In these data sets, missing values appear as NaN (not a number), representing that the itinerary is undefined. There are valid reasons to have missing values in the data set (for example, it is not possible to travel from one place to other using public transport), but a closer look reveals that the missing values are randomly distributed and do not follow any logical pattern (for example, itineraries that can be performed by MOTORCYCLE and not by CAR). This suggests that missing values originated from the routing algorithm time-out errors due to the inability to calculate complex routes. For this reason, all NaN values were replaced with an arbitrary number that was greater than any of the real measured values (in this case, all NaNs were replaced by the maximum value in the data set).
The rest of the variables do not present any obvious problems. In fact, they clearly show an expected distribution. For example, Figure 3 shows the length of the trips made with BICYCLE. As expected, there is a point where the trips using BICYCLE are too costly and fall sharply. Another expected result is, for example, the distribution of time expended in CAR trips. As seen in the left panel of Figure 3, the distribution almost follows a χ 2 distribution. As citizens tend to be normally geographically distributed [30], the distance traveled follows a χ 2 distribution, as does the time spent.
The influence of these variables on transport choice was modeled through eight modeling techniques: k-nearest neighbors (KNNs) [31], multinomial logit (M) [32], support vector machines (SVMs) [33], neural networks (NNs) [34], naïve Bayes (B) [35], fuzzy logic (CE) [36], expert knowledge (EK) [37], and random search (RA).

2.2.1. Multinomial Logit Models

Following the so-called latent construct interpretation, multinomial models work by adjusting a helper function p that takes as input a vector of external factors that affect the choice and gives as output the selected one. For a two class case, the function p takes the form:
p ( X ) : = 1 β X + ε > 0 0 e l s e
where X is a vector of the variables considered to make the regression, β is the vector of regressors to be estimated, and ε is a random variable following the logistic (logit) or normal (probit) distribution. Details can be consulted in [32].

2.2.2. Support Vector Machines

This technique builds a black-box model that can be tailored for classification and regression problems [33]. A more detailed explanation is beyond the scope of this document and readers interested can consult [33] for details. Since the transport choice is a classification problem, the default approach ε -classification was followed. This technique features three hyper-parameters: C controls the penalization to the errors in the training sample, ε controls the margin of tolerance to errors, and γ controls the bandwidth of the radial function. These hyper-parameters have been optimized using a grid search following the advice given in [38].

2.2.3. Neural Networks

This technique builds a black-box model that tries to mimic the behavior of the natural neurons. This way, a neural network is composed of a set of artificial networks arranged into three layers: an input layer, a hidden layer, and an output layer (see Figure 4). The outputs of one layer are connected to the inputs of the subsequent layers, building in this way a network of neurons. Finally, each neuron computes the following function:
y i : = φ b i + k = 0 m w k i x k
where x is the input vector w, b are the parameters of the model (wights and bias), and φ : = 1 1 + e x denotes the sigmoid function [39]. Further details can be consulted in [34]. The hidden layer features several numbers of neurons, and the parameters of the model have been adjusted through backpropagation.
The main hyper-parameter of this model is the number of neurons m in the hidden layer. In this case, a grid search has been followed to optimize this value such that m = 100 is the best value.

2.2.4. The Naïve Bayes (B) model

This model extends Bayes’ theorem to perform classifications. This way, this model estimates a probability p c ( X ) of every transport mode given the vector of variables considered to make the prediction X. For this end, the algorithm just needs to calculate an estimation of the probability of X conditional of every transport mode c and then use Bayes’ Theorem. The output class will be the one with the greatest probability. Details of this model can be found in [35].

2.2.5. k-Nearest Neighbours

This method uses the notion of the closeness of the input to perform classifications. This way, given a new element y to classify, the algorithm searches the k-nearest elements known in the training set with respect to a given metric (usually the Euclidean) and classifies y as the weighted mean value of the nearest neighbours (in a regression setting) or the most common class among the k-nearest elements (in a classification setting as this one). Using a grid search procedure, it was found that k = 5 was the best value for this project.

2.2.6. Fuzzy Logic

All methods presented above are black-box models or are not straightforward in considering new transport modes or new input variables. Fuzzy logic [36], however, has been shown to achieve high success in the emulation of individual preferences for transport mode choice [40,41], therefore enabling a wider level of parametrization for each citizen’s profile. A well trained fuzzy logic engine and accurately defined rule sets are expected to provide consistent results even in the case of variation in the characteristics of the transport network.
Fuzzy logic engines build discrete choice models through a simple natural language rule-based approach [42]. Classical logic only allows conclusions to be true or false. In contrast, fuzzy logic defines a set of rules that map numeric data into linguistic terms and creates fuzzy sets, indicating the extent to which each term is part of. That is, a variable can have several values (called terms) that can be overlapped and shaped. Thereby, it is possible to model how each feature increases or decreases the likelihood of choosing a given transport mode. Fuzzy logic engines are parametrized though their inputs, outputs, and rules. For this research, the fuzzy logic was built as follows:
  • Inputs, terms, and membership functions: For the engine to translate from linguistic terms to numbers, inputs are composed of a name that identifies the input, the possible terms or values that the input can take, and a membership function that limits the extent to which each term can be classified. The inputs of the model presented here are built from the four previously mentioned itinerary features: DURATION, PRICE, LENGTH, and ENVIRONMENTAL IMPACT. Each of these inputs feature a fixed set of three terms: LOW, MEDIUM, and HIGH, shaped as triangles. The anchors of these triangles ( α , β , γ , and δ ), shown in Figure 5, are key values that describe how citizens understand these qualitative variables.
  • Outputs, terms, and membership functions: Similarly to the inputs, the outputs are composed of a name that identifies the output, the possible terms or values the output can take, and a membership function for each term. The model presented here features a single output with five fixed, non-overlapping terms that represent the transport modes among which to choose: CAR, MOTORCYCLE, PUBLIC TRANSPORT, BICYCLE, and WALK, as represented in Figure 6.
  • Rules: Rules follow the following structure: if INPUT is INPUT_TERM and … then OUTPUT is OUTPUT_TERM. The engine needs a set of rules, each with one output but with one or more terms, as depicted in Table 3. Additionally, the engine provides fuzzy value modifiers or hedges through natural language keywords to help describe border cases without the need to include additional terms: ANY, NOT, EXTREMELY, SELDOM, SOMEWHAT, or VERY. The defuzzifier will be the maximum value between the three terms, that is, the term that receives the maximum value and thus the one chosen by the citizen. Please note that only AND conjunctions will be used between rule terms, since disjunctions can be represented by splitting a rule in two or using the NOT hedge.

Expert Fuzzy Knowledge (EK)

A first approach of the definition of the fuzzy logic inputs and rules was created according to both literature [37] and internal transport surveys carried out for the green traveling project [43,44]. From this information, a simple set of rules was extracted. These rules can be consulted in Table 4.

A Co-Evolutionary Fuzzy Logic Algorithm (CE)

While expert knowledge can lead to accurate solutions and is crucial for a better understanding of the problem, it is not flexible enough to variations in the inputs used to build that knowledge. Therefore, it is always advisable to further search whether a better solution can be reached. Trying to manually improve the expert knowledge rules to match a wider number of citizens from the census involves applying many changes to interrelated inputs and rule sets in a very large search space. One of the approaches to automate this is to use an evolutionary algorithm and to divide the search space in subsets of problems.
Genetic programming techniques have already been used to train fuzzy logic models [45]. Here, a cooperative co-evolution programming method (CE) is presented. CE is a form of evolutionary computation method that divides a large problem into sub-components and solves them independently in order to solve the large problem [46]. The sub-components are implemented as sub-populations and the only interaction between sub-populations is in the cooperative evaluation of each individual of the sub-populations. The cooperative evaluation of each individual in a sub-population is done by concatenating the current individual with the best individuals from the rest of the sub-populations. The algorithm presented hereby tackles both the adjustment of the the input linguistic variables membership functions and the search of the rule sets that most faithfully fit the census data. These two variables are strongly related since any change in the input terms’ extent requires adapting the rules to the new meaning of the linguistic variable. Thus, the CE algorithm suits perfectly the problem of the optimizing two interrelated variables in two sub-components:
  • Input-evolving sub-component: This is composed of a population of 200 input sets. Each input set features the four input variables already mentioned—DURATION, PRICE, LENGTH, and ENVIRONMENT. Each input has three fixed and ordered terms—LOW, MEDIUM, and HIGH—whose triangle anchors α , β , γ , and δ (see Figure 5 for details) can be modified in search of limits that better fit the citizens’ appreciation for these qualitative terms.
  • Rule evolving sub-component: This is composed of a population of 200 rule sets. Each rule set features a list with a variable number of rules, and, for each rule, its terms and output are modified to extract the combination of rules that more faithfully represents citizens’ transport choice from the census.
Each combination of input membership functions and rules defines what is called a fuzzy experiment. On evaluation, each fuzzy experiment configures the fuzzy logic engine to test how accurate the outputs of the engine are compared to the census. On the evolution of the input sets, whenever a change is applied in any of the triangle anchors, a fuzzy experiment will be created with corresponding input set, the fixed output, and the best found rule set to calculate the input set’s fitness. Likewise, on the evolution of the rule sets, whenever a change is applied to the rules (whether it is an addition, a deletion, or a change in the terms), a fuzzy experiment will also be created with a corresponding rule set, and the best input set for calculating its fitness will be found.
Although evolutionary algorithm operators traditionally use mutation and crossover, rule set mutations produce only slight changes in the population. Therefore, a more abrupt approach has been designed by relying only on the crossover operator. The crossover operators for evolving rule sets and evolving input sets are detailed below.
  • Operators used for Rules:
    • rule set slice crossover Given two rule sets, this operator creates a new rule set, taking the first to n rules from the first set and the n to last rules from the second set.
    • rule set combine crossover Given two rule sets, this operator creates a new rule set, appending for each position a randomly selected rule from the same position of the two provided rule sets.
    • rule set slice terms crossover Given two rule sets, this operator clones the first rule set and selects a p rule in the set. The p rule’s terms are replaced by the first to n terms from the p rule of the first set and n to last from the p rule of the second provided set.
  • Fuzzy Logic Input Membership Functions:
    • input set combine crossover Given two input sets, this function creates a new input set, appending for each position a randomly selected element from the same position of the two provided input sets.
    • input set shape mutation Given an input set, this function modifies one of the α , β , γ , or δ values from the shapes in Figure 5. It will displace at the same time the first term’s triangle’s third point, the second term’s triangle’s mid-point, and the third term’s triangle’s first point.
The fitness of a fuzzy experiment is given by the amount of citizens that match the transport choice originally registered in the census (that is, the accuracy of the forecasting model). The goal of the CE algorithm is, therefore, to maximize the fitness in order to estimate as many correct transport choices as possible. For this, a confusion matrix is calculated within each fitness evaluation. Confusion matrices are square matrices with a number of rows and columns equal to the number of output terms. Columns denote the amount of real elements, and rows denote the predicted elements of that class, that is, a column contains a split of 100% of the elements. Namely, cell m i j denotes the percentage of elements classified as i that are of class j. The diagonal of this matrix (namely the elements m i i ) describes the correctly classified elements. For instance, in Table A1, 19.8004% of C (CAR) instances were forecast correctly, but 14.4995% of C users were forecast as M (MOTORCYCLE) users. The global fitness is composed of the best element of the input and rules populations as seen in Figure 7.
Since the number of elements for each vehicle type is not balanced, a fuzzy experiment that matched only CAR would obtain a better fitness than other operating on a more distributed modal split. Therefore, the previously defined correctly classified elements has been improved to include balanced correctly classified elements. In this sense, let t k denote the relative partial matches of class k, namely,
t k : = i k m i k i = 1 n m i k .
Then, the fitness function f is defined as the sum of the logarithm of the partial matches t k . Namely,
f : = k = 1 n ln ( 1 + t k ) .
Please note that we add 1 to the partial matches to avoid undefined values of the fitness function when a model does not correctly classify any element. Thus, with this approximation, not only are all vehicle types balanced, but rules that match the census in a more distributive way also weigh more.
Additionally, a gradually growing subset of the census is used to calculate the fitness. This subset is extended as the matched percentage of that subset increases. That is, the first iteration’s fitness is calculated with a 10% subset of the entire census data set. Following iterations, the subset will be kept equal or increased parallel to the matched percentage. In this way, as matched elements increase, more elements are introduced in the census, against which the next iteration’s fitness can be evaluated. Moreover, in order to go through all the elements of the census, at every 15 iterations, the subset to extract from the census is moved by 10%. This changes in the reference data with which to evaluate the fitness and generates the repeating pattern visible in Figure 7.
As a summary, the pseudo code for the CE algorithm is described in Figure 8. The algorithm is divided into two main procedures:
  • initialize the population of rule sets and input sets that will take part in the co-evolutionary process, and
  • initialize the co-evolutionary process itself where the rule sets and input sets evolve together through the application of the evolutionary operators, the calculation of the fuzzy experiments’ fitness on each iteration, and the update of the percentage of census being used for evaluation.

Random Search (RA)

In order to assess whether the results had been produced by chance or, indeed, due to the evolutionary process, a simple yet powerful method, random search, was also trained and validated. For this purpose, the same amount of random fuzzy models as that of fitness evaluations for the CE were produced. As 200 rules are in the population and 200 generations have been made, 400,000 random search rules have been generated. All models, in general, and the CE algorithm, in particular, are expected to produce far better results than this method.

2.3. Transport Simulation

Since the current research focuses primarily on transport choice modeling, the description and evaluation of the transport simulation stage will be featured in future works. Section 4 describes how this phase will be tackled in future works.

3. Experimental Results

This section presents the results for the execution of the eight data science methods described in Section 2. The experiments were carried out using the following:
  • The GeoWorldSim platform [47] with the open-source fuzzy logic library Fuzzylite [48]. GeoWorldSim is a powerful software that eases the integration of multi-agent systems with reference simulations tools such as Matlab or EPA-NET.
  • The Caret package [49] of the R Statistical Suite [50]. This package uniforms the access to several libraries and stages needed to build data-based models. For example, it provides a uniform way to tune the hyper-parameters, train the models, and produce a forecast among other facilities.
Simulations were executed on two computers: large experiments were run on a system with an AMD Opteron 6168 CPU and 32 GB of RAM, and short experiments were run on a system with an Intel Core i7-2600 CPU and 16 GB of RAM.
Table 5 shows the global accuracy results, that is, the number of trips correctly classified. The first row (5-TM) shows the results with the five transport modes, while the second row (3-TM) shows the results after grouping CAR and MOTORCYCLE into PRIVATE VEHICLE as well as WALK and BICYCLE into WALK.
As seen, all techniques are able to transfer the results almost seamlessly since the differences between the accuracy of the training (T) and validation (V) sets are rather small. As explained previously, each method was trained with data from the census of Biscay (Basque Country, Spain) and validated over the data set of Silesian Voivodeship (Poland). The models can be classified in three groups: the cluster of the best models, the contrast models, and the rest of the models. As expected, the worst model is the random search (RA). The best models are the most complex ones: SVM, NN, and M. Surprisingly, the KNN method, while being very simple, achieves the best results. The rest of the models—EK, CE, and B—complete the clusters of the other models. In all cases, the best models can correctly predict the transport choice of almost half of the trips registered in the data set.
For more detail, Table A1 contains the confusion matrices of all models included in this article. As previously stated, within these matrices, columns show the observed modes, and rows show the forecasted real modes. In theses tables, it is clear that some of the models have difficulties differentiating between CAR - MOTORCYCLE and WALK - BICYCLE. This may be caused by the existence of other additional variables, such as socio-demographics, car/motorcycle ownership, or climatology factors that are not taken into consideration in this article. Therefore, since the census lacked information about these features, the five transport modes were grouped into three categories. This way, CAR and MOTORCYCLE were gathered into PRIVATE VEHICLE, WALK - BICYCLE into WALK, and TRANSIT is maintained as its own class.
With these three transport modes, all techniques are able to correctly transfer the results from the training to the validation set (Table 5). As before, three clusters of techniques appear. On the one hand, RA continues to be the worst method, but it is close to B. On the other hand, EK has greatly improved its performance and produces the best cluster. In any case, the best models have not improved their performance, yet they are still able to correctly predict the transport choice of half of the trips made.
In order to validate the previous claims, a bootstrapping validation was used [51]. Thus, 100 samples were randomly built following the original distribution for the Silesia data sets. Figure 9 shows a box plot with the results of the different models bringing to light the significant differences among them. In fact, a Friedman test confirms this hypothesis (p-value < 10 16 ). In order to assess the differences among the different models, a post-hoc analysis has been made following the procedure described in [52]. This procedure clusters all methods with similar results and gives them the same letter. It should be taken into consideration that these groups are not sorted, namely, the best algorithms are not labeled with an a. The results are shown in Figure 10. In this figure, the mean value of the different algorithms is plotted and groups are identified with a letter. As can be seen, the test confirm our claims. The three groups are clearly visible.
Again, when analyzing the results in detail, problems arise for some of the models. Looking at the confusion matrices (Table A2) highlights the fact that some models continue to have problems with some of the classes. For example, the SVM cannot correctly predict any WALK trip (even if it has correctly predicted half of the trips). The same problem occurs, when looking at the M and NN models, though it is not as acute. In the end, it seems that only the KNN, EK, and CE models do not suffer these problems in the setting.
Finally, Table 6 shows the forecasts that the different models share. As before, 100 bootstrapping samples were built from the Silesia data set. Column Real shows the results of the survey carried out in 2015 [26], and Columns L.C.I. and U.C.I. contain the lower and upper confidence interval at the 0.05 significant level ( α = 0.5 ) for the forecast value made. As expected, given the low variance of the bootstrapped samples (see Figure 9) the modal split forecasts of all models are quite similar (the confidence intervals are quite narrow). Nevertheless, the real value is not included in the confidence interval.

4. Conclusions and Future Work

This article presents eight methodologies for building transport choice models according to FSM. The methodologies were trained for one scenario, Biscay (Basque Country, Spain), and later evaluated with a second scenario, Silesian Voivodeship (Poland). As explained in Section 3, the results suggest the existence of three groups of models: On the one hand, as expected, the most complex models tested (SVM, CE, and M) perform quite well, hitting almost half of the trips. Next, there is another set of complex algorithms but with slightly worse results (NN, EK, and KNN). The results of the KNN algorithm was in fact unexpected, as it is a very simple algorithm. On the other hand, the control methods (B and RA) do not achieve good results, as expected. In order to draw meaningful conclusions, a qualitative analysis should also be made:
  • Some of the models produce unbalanced predictions, like M, SVM and NN, which completely ignore some of the transport modes. Therefore, given the similar prediction capabilities, the models with more balanced predictions should be used, namely, the EK, CE, and KNN models.
  • The number of parameters and complexity of the models are not the same either. In this sense, simpler models should be preferred to more complex models. However, EK, CE, and KNN are non-parametric, and their complexity is difficult to assess. On the one hand, the number of parameters of KNN is indeed the amount of points in the training data set, but the complexity of the model is quite low (understanding here the complexity in terms of its VC dimension or similar measures [53]). On the other hand, the amount of parameters in the EK and CE models is quite low (compared to the number of parameters of KNN), but the complexity of the model is higher [54].
  • Comparing the EK and CE fuzzy rule sets, both have similar amounts of rules (EK has 33 and CE has 30). However, the EK rules are only composed of one or two terms, which makes it easy to follow and modify, while the CE evolved rules usually have about 10 or more terms with a relation between them that is not very clear.
  • It is always better to use a model that is easy to be understood than a data-driven model. Following this advice, EK and M should be preferred to the rest of the models.
Based on the above considerations, it is clear that the best model is EK, as it is possible to be understood, is easy to be extended so as to cover new transport policies, has a complexity that is similar to other models that produce similar numerical results, and its predictions are well balanced among the classes.
Even though the results are good, the models need to be improved. At the moment, the methodology does not take into account socio-economic variables, which are not present on the tested data sets (see Table 2). Recent works ([55,56]) have achieved successful results by using demographic and socio-economic features. As shown in [56] (where they have almost reached a 90% chance of success), socio-economic information from commuters enables the clustering and creation of custom utility functions for each group. Furthermore, climatic variables have been identified as relevant [28]. Future refinements of the model may introduce both sets of variables in search of reducing error. A first approach should consist of adding these non-trip related features to the global census models and evaluate their improvement. In case this is not satisfactory, the next step will be to try clustering the commuters and adjusting the models for each group.
Additionally, the ownership of private vehicles is a relevant parameter not detailed in the used data sets. In order to introduce this information, the model will use Monte Carlo simulations so that several distributions of these variables are simulated and the results assessed.
Finally, once the transport choice model is fitted, the next step is to determine the extent to which different non-physical incentives, such as changes in the price of the journey, discounts on public transport, additional taxes on the use of private vehicles, reductions in the duration of the journey, and congestion charges, affect commuters’ choices. These incentives alter the features of the itineraries, resulting in a different modal split. In order to check the suitability of these incentives, the citizen agents within the simulation will be given a set of preferences that will determine whether they are prone to accept or ignore a certain policy, based on the itinerary features already modeled in the transport choice model. This approach will involve defining the thresholds at which the applications of transport policies will be effective, all the while seeking an equilibrium between the objectives of the policies and the citizens’ preferences in order to avoid social rejection or oversizing.

Acknowledgments

This work was partially supported by (a) the GREENTRAVELLING project (Era NET Transport 6/12/IN/2014/195 label) partially funded by Diputación Foral de Bizkaia through the Plan de Promoción de la Innovación (6/12/-IN/2014/195) and the Basque Government through the GAITEK program IG-2014/0000133), (b) the S-MILE project (Era NET Transport ENTIII Flagship 2015 label) funded by Basque Government through the HAZITEK program (ZL-2017/00081), and (c) the Industrial Ph.D. grant given by the University of Deusto (2015–2018).

Author Contributions

Ander Pijoan, Oihane Kamara-Esteban, Ainhoa Alonso-Vicario and Cruz E. Borges. conceived and designed the experiments, performed the analysis of the results and contributed to the writing of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Result tables

Table A1. Confusion matrices for all the models for the five transport modes (in %). Source: own research.
(a) Expert Knowledge Confusion Matrix
(a) Expert Knowledge Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB20.930225.958813.548316.904222.913839.81806.323500.00050.0051
C14.341019.800411.612918.56396.85613.59723.98768.53395.1276453.5860
M18.992214.499523.870921.089623.049113.30288.169920.84801.72674112.15994
T12.790617.570914.838714.08984.239916.402912.159960.130724.025933.56175
W32.945722.170236.129029.352342.940926.879249.358910.487469.12920.6872
(b) Co-evolutionary Fuzzy Algorithm Confusion Matrix
(b) Co-evolutionary Fuzzy Algorithm Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB25.581418.4526628.7096822.9061420.758124.34213.65663.22584.06525.2098
C5.0387626.158173.8709685.0180513.74548743.289455.660920.967750.004214.4617
M0000046.315735.688864.516140.378248.2227
T24.0310131.1495924.5161329.079425.6859210.39470.141500.13850.5517
W45.3488424.2395942.9032342.9963969.810475.65784.851911.29035.413731.5539
(c) Support Vector Machine Confusion Matrix
(c) Support Vector Machine Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB0000000000
C34.108563.657636.774138.228311.141161.769368.958343.790865.711777.1642
M0000000000
T65.891436.342363.225861.771688.858838.230631.041656.209134.288222.8357
W0000000000
(d) Multinomial Confusion Matrix
(d) Multinomial Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB3.0076000000000
C34.210564.5774118230624863.138471.468344.11768.336315.8941
M000000.18430.649100.53310.0051
T62.781935.34451923236196936.571827.808355.882331.058183.6567
W00.07790100.10530.074200.07230.4439
(e) Neural Network Confusion Matrix
(e) Neural Network Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB00.01550000.02630.023300.03160.0487
C33.082763.5640115225029842.733052.207826.143748.880217.4520
M000000.10530.062300.05930.2540
T62.406036.07731933257188356.766747.257472.875850.517081.1723
W1.50370.3429236360.36860.44900.98030.51161.0727
(f) Naïve Bayesian Confusion Matrix
(f) Naïve Bayesian Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB3.38350.28060.32250.216400.92152.40441.63391.08440.8469
C9.398429.70067.74197.46884.194826.513931.634117.973830.548117.5213
M000000.76350.562200.60380.6980
T18.796927.284020.322526.35755.142018.667724.34707.843121.92180.1154
W68.421042.734671.612965.957090.663053.133241.052172.549045.841780.8181
(g) k-Nearest Neighbours Confusion Matrix
(g) k-Nearest Neighbours Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB0000000000
C33.458664.483937.096733.104821.244952.948961.764238.235258.532226.5013
M0000000000
T54.135331.602754.838760.779360.983344.707736.109058.823539.170561.8109
W12.4063.91338.06456.115817.77172.34332.12672.94112.297211.6877
(h) Random Search Confusion Matrix
(h) Random Search Confusion Matrix
BiscaySilesia
ObservedObserved
BCMTWBCMTW
ForecastB0000000000
C10010010010010000.14341.96085.174223.0056
M0000010085.71798.039291.150776.2244
T000000003.67510
W00000014.1396000.7699
Table A2. Confusion matrices for all the models for the three transport modes (in %). Source: own research.
(a) Expert Knowledge Confusion Matrix
(a) Expert Knowledge Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP27.23068.577110.861764.305260.493232.9044
T47.686751.663327.55830.05280.05870.4739
W25.082639.759561.579835.641939.447966.6216
(b) Co-evolutionary Confusion Matrix
(b) Co-evolutionary Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP27.88516.96375.292923.339726.487415.1919
T44.824546.112221.050540.558355.267359.8538
W27.290346.924073.656536.10218.245324.9543
(c) Suport Vector Machine Confusion Matrix
(c) Suport Vector Machine Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP63.979740.014414.141463.977160.465820.4924
T36.020259.985585.737336.022839.534179.5075
W000.1212000
(d) Multinomial Confusion Matrix
(d) Multinomial Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP65.913144.109615.191976.545674.021422.5083
T34.012455.872284.808023.332925.845177.0076
W0.07430.018000.12140.13340.4840
(e) Neural Network Confusion Matrix
(e) Neural Network Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP58.001136.388214.909047.296147.037024.9403
T41.612163.196883.959551.628351.941870.7169
W0.38660.41491.13131.07541.02114.3426
(f) Naïve Bayes Confusion Matrix
(f) Naïve Bayes Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP28.97087.68534.808034.606232.243819.8774
T26.933326.35756.626224.302921.91271.7585
W44.095765.957088.565641.090845.843478.3639
(g) Random Search Confusion Matrix
(g) Random Search Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP99.936699.943110022.3468100100
T00077.653200
W0.06330.05680000
(h) K-Nearest Neighbours Confusion Matrix
(h) K-Nearest Neighbours Confusion Matrix
BiscaySilesia
ObservedObserved
PTWPTW
ForecastP65.987536.460425.373764.200661.209931.7244
T29.089855.998553.494933.038635.761754.8664
W4.92267.541021.13132.76073.028213.4091

References

  1. De Donnea, F.X. The Determinants of Transport Mode Choice in Dutch Cities: Some Disaggregate Stochastic Models; Universitaire Pers Rotterdam: Rotterdam, The Netherlands, 1971. [Google Scholar]
  2. Holtsmark, B.; Skonhoft, A. The Norwegian support and subsidy policy of electric cars. Should it be adopted by other countries? Environ. Sci. Policy 2014, 42, 160–168. [Google Scholar] [CrossRef]
  3. Karlaftis, M.G.; Vlahogianni, E.I. Statistical methods versus neural networks in transportation research: Differences, similarities and some insights. Transp. Res. Part C Emerg. Technol. 2011, 19, 387–399. [Google Scholar] [CrossRef]
  4. Fernández-Isabel, A.; Fuentes-Fernández, R. An agent-based platform for traffic simulation. In Soft Computing Models in Industrial and Environmental Applications, 6th International Conference SOCO 2011; Springer: Berlin/Heidelberg, Germany, 2011; pp. 505–514. [Google Scholar]
  5. Bazzan, A.L.; Klügl, F. A review on agent-based technology for traffic and transportation. Knowl. Eng. Rev. 2014, 29, 375–403. [Google Scholar] [CrossRef]
  6. Rubinstein, R.Y.; Kroese, D.P. Simulation and the Monte Carlo Method; John Wiley & Sons: Hoboken, NJ, USA, 2016; Volume 10. [Google Scholar]
  7. Keilson, J. Markov Chain Models—Rarity and Exponentiality; Springer Science & Business Media: Berlin, Germany, 2012; Volume 28. [Google Scholar]
  8. Allen, A.O. Probability, Statistics, and Queueing Theory; Academic Press: Cambridge, MA, USA, 2014. [Google Scholar]
  9. Meyer, M.D.; Miller, E.J. Urban Transportation Planning: A Decision-Oriented Approach; McGraw-Hill: New York, NY, USA, 1984. [Google Scholar]
  10. McNally, M.G. The four-step model. In Handbook of Transport Modelling, 2nd ed.; Emerald Group Publishing Limited: Bingley, UK, 2007; pp. 35–53. [Google Scholar]
  11. Wooldridge, M.; Jennings, N.R. Intelligent Agents: Theory and Practice; Cambridge University Press: Cambridge, UK, 1995. [Google Scholar]
  12. Owen, L.E.; Zhang, Y.; Rao, L.; McHale, G. Street and traffic simulation: Traffic flow simulation using CORSIM. In Proceedings of the 32nd conference on Winter simulation, Orlando, FL, USA, 10–13 December 2000; pp. 1143–1147. [Google Scholar]
  13. Saad, F. In-depth analysis of interactions between drivers and the road environment: Contribution of on-board observations and subsequent verbal report. In Proceedings of the 4th Workshop of ICTCT, University of Lund, Vienna, Austria, November 1992. [Google Scholar]
  14. PTV AG. PTV Vissim 7—USER MANUAL; PTV AG: Karlsruhe, Njemačka, 2015. [Google Scholar]
  15. Barceló, J.; Casas, J. Dynamic network simulation with AIMSUN. In Simulation Approaches in Transportation Analysis; Springer: Boston, MA, USA, 2005; pp. 57–98. [Google Scholar]
  16. Halati, A.; Lieu, H.; Walker, S. CORSIM-corridor traffic simulation model. In Proceedings of the Traffic Congestion and Traffic Safety in the 21st Century: Challenges, Innovations, and Opportunities, Chicago, IL, USA, 8–11 June 1997. [Google Scholar]
  17. Mahmassani, H. Dynamic traffic assignment and simulation for advanced network informatics (DYNASMART). In Proceedings of the 2nd International Seminar on Urban Traffic Networks, Capri, Italy, 5–8 July 1992. [Google Scholar]
  18. Smith, L.; Beckman, R.; Anson, D.; Nagel, K.; Williams, M. TRANSIMS: Transportation Analysis and Simulation System; Technical Report; Los Alamos National Lab.: Los Alamos, NM, USA, 1995.
  19. Babin, A.; Florian, M.; James-Lefebvre, L.; Spiess, H. Emme/2 an interactive graphic method for road and transit planning. In Proceedings of the RTAC Annual Conference Preprints, Washington, DC, USA, 18–22 January 1981; Volume 1. [Google Scholar]
  20. PTV Group Ptv Visum. 2016. Available online: http://your.visum.ptvgroup.com/FreeDemoVersion (accessed on 4 November 2016).
  21. Ji, M.; Yu-long, P. Development and application of TransCAD for urban traffic planning. J. Harbin Univ. Civ. Eng. Archit. 2002, 5, 028. [Google Scholar]
  22. Balcombe, R.; Mackett, R.; Paulley, N.; Preston, J.; Shires, J.; Titheridge, H.; Wardman, M.; White, P. The Demand for Public Transport: A Practical Guide; Transportation Research Laboratory: Wokingham, UK, 2004. [Google Scholar]
  23. Train, K. Qualitative Choice Analysis: Theory, Econometrics, and an Application to Automobile Demand, 1st ed.; The MIT Press: Cambridge, MA, USA, 1985; Volume 1. [Google Scholar]
  24. Rich, J.; Holmblad, P.; Hansen, C. A weighted logit freight mode-choice model. Transp. Res. Part E Logist. Transp. Rev. 2009, 45, 1006–1019. [Google Scholar] [CrossRef]
  25. Instituto Nacional de Estadística. Censos 2011. 2011. Available online: http://www.ine.es/censos2011_datos/cen11_datos_resultados.html (accessed on 13 April 2018).
  26. Green Travelling Consortium. A Platform to Analyse and Foster the Use of Green Travelling Options; Technical Report, Project Proposal, The ERA-NET Transport III: Future Travelling; Avda. Universidades: Bilbao, Spain, 2013. [Google Scholar]
  27. Domencich, T.A.; McFadden, D. Urban Travel Demand—A Behavioral Analysis; North-Holland Publishing Co.: Amsterdam, The Netherlands, 1975. [Google Scholar]
  28. Chiu Chuen, O.; Karim, M.R.; Yusoff, S. Mode Choice between Private and Public Transport in Klang Valley, Malaysia. Sci. World J. 2014, 2014, 394587. [Google Scholar] [CrossRef] [PubMed]
  29. Clauss, T.; Döppe, S. Why do urban travelers select multimodal travel options: A repertory grid analysis. Transp. Res. Part A Policy Pract. 2016, 93, 93–116. [Google Scholar] [CrossRef]
  30. Pijoan, A.; Borges, C.E.; Oribe-Garcia, I.; Martin, C.; Alonso-Vicario, A. Agent Based Simulations for the Estimation of Sustainability Indicators. Procedia Comput. Sci. 2015, 51, 2943–2947. [Google Scholar] [CrossRef]
  31. Altman, N.S. An introduction to kernel and nearest-neighbor nonparametric regression. Am. Stat. 1992, 46, 175–185. [Google Scholar]
  32. Lemeshow, S.; Sturdivant, R.X.; Hosmer, D.W. Applied Logistic Regression; Wiley: Hoboken, NJ, USA, 2013. [Google Scholar]
  33. Cristianini, N.; Shawe-Taylor, J. An Introduction to Support Vector Machines and Other Kernel-Based Learning Methods; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  34. Schalkoff, R.J. Artificial Neural Networks; McGraw-Hill: New York, NY, USA, 1997; Volume 1. [Google Scholar]
  35. Lewis, D.D. Naive (Bayes) at forty: The independence assumption in information retrieval. In Machine Learning: ECML-98; Springer: Berlin/Heidelberg, Germany, 1998; pp. 4–15. [Google Scholar]
  36. Mendel, J.M. Uncertain Rule-Based Fuzzy Logic Systems: Introduction and New Directions; Prentice Hall PTR: Upper Saddle River, NJ, USA, 2001. [Google Scholar]
  37. Souche, S. Measuring the structural determinants of urban travel demand. Transp. Policy 2010, 17, 127–134. [Google Scholar] [CrossRef] [Green Version]
  38. Chang, C.C.; Lin, C.J. LIBSVM: A Library for Support Vector Machines. ACM Trans. Intell. Syst. Technol. 2011, 2, 27:1–27:27. [Google Scholar] [CrossRef]
  39. Von Seggern, D.H. CRC Standard Curves and Surfaces with Mathematica; CRC Press: Boca Raton, FL, USA, 2006. [Google Scholar]
  40. Kumar, M.; Sarkar, P.; Madhu, E. Development of fuzzy logic based mode choice model considering various public transport policy options. Int. J. Traffic Transp. Eng. 2013, 3, 408–425. [Google Scholar] [CrossRef]
  41. Tuzkaya, U.R.; Önüt, S. A fuzzy analytic network process based approach to transportation-mode selection between Turkey and Germany: A case study. Inf. Sci. 2008, 178, 3133–3146. [Google Scholar] [CrossRef]
  42. Berkan, R.C.; Trubatch, S. Fuzzy System Design Principles; Wiley-IEEE Press: Hoboken, NJ, USA, 1997. [Google Scholar]
  43. Staniek, M.; Sierpiński, G. Cost Criteria as a Means to Support Travelling Mode Related Decisions—A Case Study for the Central Part of Silesian Voivodeship (Poland). In Sustainable Transport Development, Innovation and Technology; Springer: Cham, Switzerland, 2016; pp. 137–150. [Google Scholar]
  44. Sierpiński, G.; Staniek, M.; Celiński, I. Travel behavior profiling using a trip planner. Transp. Res. Procedia 2016, 14, 1743–1752. [Google Scholar] [CrossRef]
  45. Herrera, F. Genetic fuzzy systems: Taxonomy, current research trends and prospects. Evolut. Intell. 2008, 1, 27–46. [Google Scholar] [CrossRef]
  46. Potter, M.A.; De Jong, K.A. A cooperative coevolutionary approach to function optimization. In Parallel Problem Solving from Nature; Springer: Berlin/Heidelberg, Germany, 1994; pp. 249–257. [Google Scholar]
  47. Pijoan, A.; Kamara-Esteban, O.; Borges, C.E. Environment modelling for spatial load forecasting. In Agent Environments for Multi-Agent Systems IV; Springer: Verlag, New York, 2014; pp. 188–206. [Google Scholar]
  48. Rada-Vilela, J. Fuzzylite: A Fuzzy Logic Control Library. 2017. Available online: https://www.fuzzylite.com/ (accessed on 13 April 2018).
  49. Kuhn, M. caret: Classification and Regression Training; R Package Version 6.0-78; ASCL: Leicester, UK, 2017. [Google Scholar]
  50. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2017. [Google Scholar]
  51. Efron, B. Bootstrap Methods: Another Look at the Jackknife. Ann. Stat. 1979, 7, 1–26. [Google Scholar] [CrossRef]
  52. Piepho, H.P. An Algorithm for a Letter-Based Representation of All-Pairwise Comparisons. J. Comput. Graph. Stat. 2004, 13, 456–466. [Google Scholar] [CrossRef]
  53. Vapnik, V.N.; Chervonenkis, A.Y. On the Uniform Convergence of Relative Frequencies of Events to Their Probabilities. Theory Probab. Appl. 1971, 16, 264–280. [Google Scholar] [CrossRef]
  54. Chen, Y.; Wang, J.Z. Support vector learning for fuzzy rule-based classification systems. IEEE Trans. Fuzzy Syst. 2003, 11, 716–728. [Google Scholar] [CrossRef]
  55. Khattak, Z.H.; Magalotti, M.J.; Miller, J.S.; Fontaine, M.D. Using New Mode Choice Model Nesting Structures to Address Emerging Policy Questions: A Case Study of the Pittsburgh Central Business District. Sustainability 2017, 9, 2120. [Google Scholar] [CrossRef]
  56. Ding, L.; Zhang, N. A travel mode choice model using individual grouping based on cluster analysis. Procedia Eng. 2016, 137, 786–795. [Google Scholar] [CrossRef]
Figure 1. Area and transport links of Biscay (Spain) and Silesia (Poland) data sets. Source: own research generated from real boundaries and transport networks.
Figure 1. Area and transport links of Biscay (Spain) and Silesia (Poland) data sets. Source: own research generated from real boundaries and transport networks.
Sustainability 10 01230 g001aSustainability 10 01230 g001b
Figure 2. Modal splits for the Biscay and Silesia data sets. Source: own research from the datasources of [25,26].
Figure 2. Modal splits for the Biscay and Silesia data sets. Source: own research from the datasources of [25,26].
Sustainability 10 01230 g002
Figure 3. Distribution of time (left) and length (right) for CAR and BICYCLE transport modes in Silesian Voivodeship.
Figure 3. Distribution of time (left) and length (right) for CAR and BICYCLE transport modes in Silesian Voivodeship.
Sustainability 10 01230 g003
Figure 4. Structure of a neural network (NN).
Figure 4. Structure of a neural network (NN).
Sustainability 10 01230 g004
Figure 5. Fuzzy engine input membership functions, with the configurable anchors that modify the shape of the terms. Source: own research for better understanding the input membership functions.
Figure 5. Fuzzy engine input membership functions, with the configurable anchors that modify the shape of the terms. Source: own research for better understanding the input membership functions.
Sustainability 10 01230 g005
Figure 6. Fuzzy engine output. Each term represents a mode of transport one may choose. Source: own research for better understanding the fuzzy logic output membership function.
Figure 6. Fuzzy engine output. Each term represents a mode of transport one may choose. Source: own research for better understanding the fuzzy logic output membership function.
Sustainability 10 01230 g006
Figure 7. Iterations of the co-evolutive fuzzy algorithm. Source: own research by plotting the evolution in time of the algorithm.
Figure 7. Iterations of the co-evolutive fuzzy algorithm. Source: own research by plotting the evolution in time of the algorithm.
Sustainability 10 01230 g007
Figure 8. Pseudocode of the co-evolutionary fuzzy logic algorithm. Source: own research simplification of the real code programmed in the algorithm.
Figure 8. Pseudocode of the co-evolutionary fuzzy logic algorithm. Source: own research simplification of the real code programmed in the algorithm.
Sustainability 10 01230 g008
Figure 9. Box plot of the distribution of the models’ accuracy of the three transport modes on 100 bootstrapping samples. In Biscay, the probability of using BICYCLE, CAR, MOTORCYCLE, TRANSIT, and WALK is 1.75%, 43.51%, 2.1%, 37.6%, and 15.04%, respectively. In Silesia, the probability of using BICYCLE, CAR, MOTORCYCLE, TRANSIT, and WALK is 0.82%, 53.96%, 0.12%, 36.33%, and 8.75%, respectively. Therefore, if a model can produce such a forecast, it will achieve an accuracy of 100 %. Source: own research generated from R Statistical Suite [50].
Figure 9. Box plot of the distribution of the models’ accuracy of the three transport modes on 100 bootstrapping samples. In Biscay, the probability of using BICYCLE, CAR, MOTORCYCLE, TRANSIT, and WALK is 1.75%, 43.51%, 2.1%, 37.6%, and 15.04%, respectively. In Silesia, the probability of using BICYCLE, CAR, MOTORCYCLE, TRANSIT, and WALK is 0.82%, 53.96%, 0.12%, 36.33%, and 8.75%, respectively. Therefore, if a model can produce such a forecast, it will achieve an accuracy of 100 %. Source: own research generated from R Statistical Suite [50].
Sustainability 10 01230 g009
Figure 10. Results of an ALL vs. ALL post-hoc analysis of the accuracy of the models. Source: own research by plotting the results from R Statistical Suite [50].
Figure 10. Results of an ALL vs. ALL post-hoc analysis of the accuracy of the models. Source: own research by plotting the results from R Statistical Suite [50].
Sustainability 10 01230 g010
Table 1. Itinerary features used for transport choice. Source: Adapted from [29].
Table 1. Itinerary features used for transport choice. Source: Adapted from [29].
FeatureDescriptionEncompasses
DurationTime necessary (in seconds) to complete the itinerary from the departure point until arriving at the destination.Closeness to public transport stop and parking spaces.
Frequency and reliability of public transport.
In-journey traveling time.
Traffic congestion.
Transfer waiting time between two means of transport.
Waiting environment.
Poor waiting time.
PriceMonetary cost (in euros) of the itinerary including public transport tickets, parking fares and private vehicle costs.Public transport fares.
Private vehicle purchase, maintenance, and fuel costs.
Tolls and congestion charges.
Economic incentives.
LengthDistance of the itinerary (in kilometers). For long- and medium-length travels, this is not a determinant factor. This is only when walking or cycling are an option.Walking distance.
Cycling distance.
Environmental impactContribution to climate change (measured in CO2 emissions). Regarding short distance and itineraries, where transports compete in similar conditions, environmental awareness can be the trigger that determines user decisions.Environmental impact of the trip.
Table 2. Itinerary features are not present in the data sets and are beyond the scope of the study. Source: own research, conducted by analyzing the features mentioned in [27], which are not present in the data sets.
Table 2. Itinerary features are not present in the data sets and are beyond the scope of the study. Source: own research, conducted by analyzing the features mentioned in [27], which are not present in the data sets.
FeatureDescriptionEncompasses
ProfilePersonal demographic, socio-economical, and cultural attributes of the user.Citizen age
Citizen gender
Socio-economic profile
Cultural profile
Family size
Trip type
OwnershipVehicle(s) under the property of the user.Car ownership
Motorcycle ownership
Bicycle ownership
ClimatologyVariables that may have a direct impact on the physical variables of the trip, i.e., if it rains, the duration of a trip may take longer than usual.Environment climatology
Table 3. Example of fuzzy engine rule sets. Source: own research for explaining the structure and examples of fuzzy logic rules.
Table 3. Example of fuzzy engine rule sets. Source: own research for explaining the structure and examples of fuzzy logic rules.
if TRANSIT_DURATION is HIGH and TRANSIT_PRICE is HIGH then TRANSPORT is PRIVATE
if TRANSIT_DURATION is HIGH and PRIVATE_DURATION is LOW then TRANSPORT is PRIVATE
if PRIVATE_PRICE is HIGH and TRANSIT_PRICE is LOW and TRANSIT_DURATION is MEDIUM then
TRANSPORT is TRANSIT
if WALK_LENGTH is LOW and TRANSIT_PRICE is HIGH then TRANSPORT is WALK
Table 4. Fuzzy engine rule sets extracted from surveys and expert knowledge. Source: own research generated from the surveys in [43,44].
Table 4. Fuzzy engine rule sets extracted from surveys and expert knowledge. Source: own research generated from the surveys in [43,44].
if CAR_PRICE is LOW then TRANSPORT is CAR
if CAR_TIME is LOW and TRANSIT_TIME is HIGH then TRANSPORT is CAR
if CAR_TIME is MEDIUM and TRANSIT_TIME is HIGH then TRANSPORT is CAR
if CAR_TIME is LOW and MOTORCYCLE_TIME is MEDIUM then TRANSPORT is CAR
if CAR_TIME is LOW and MOTORCYCLE_TIME is HIGH then TRANSPORT is CAR
if CAR_PRICE is LOW and TRANSIT_PRICE is LOW then TRANSPORT is CAR
if CAR_PRICE is LOW and TRANSIT_PRICE is MEDIUM then TRANSPORT is CAR
if CAR_PRICE is LOW and TRANSIT_PRICE is HIGH then TRANSPORT is CAR
if CAR_PRICE is MEDIUM and TRANSIT_PRICE is MEDIUM then TRANSPORT is CAR
if CAR_PRICE is MEDIUM and TRANSIT_PRICE is HIGH then TRANSPORT is CAR
if CAR_PRICE is HIGH and TRANSIT_PRICE is HIGH then TRANSPORT is CAR
if CAR_DISTANCE is MEDIUM and BICYCLE_DISTANCE is HIGH then TRANSPORT is CAR
 
if TRANSIT_PRICE is LOW and CAR_PRICE is MEDIUM then TRANSPORT is TRANSIT
if TRANSIT_PRICE is LOW and CAR_PRICE is HIGH then TRANSPORT is TRANSIT
if TRANSIT_PRICE is MEDIUM and CAR_PRICE is HIGH then TRANSPORT is TRANSIT
if TRANSIT_TIME is MEDIUM and CAR_TIME is MEDIUM then TRANSPORT is TRANSIT
if TRANSIT_TIME is MEDIUM and CAR_TIME is HIGH then TRANSPORT is TRANSIT
if TRANSIT_TIME is HIGH and CAR_TIME is HIGH then TRANSPORT is TRANSIT
 
if MOTORCYCLE_DISTANCE is LOW and CAR_PRICE is HIGH then TRANSPORT is MOTORCYCLE
if MOTORCYCLE_PRICE is MEDIUM and CAR_PRICE is HIGH then TRANSPORT is MOTORCYCLE
 
if WALK_DISTANCE is LOW then TRANSPORT is WALK
 
if BICYCLE_DISTANCE is LOW then TRANSPORT is BICYCLE
if BICYCLE_TIME is LOW then TRANSPORT is BICYCLE
Table 5. Global accuracy (%) of the algorithm. Column T (testing) shows the results for Biscay’s data set and Column V (validation) for Silesia. Source: own research by extracting the forecast accuracy of the different models in Table A1. EK: expert knowledge; CE: fuzzy logic; SVM: support vector machine; M: multinomial logit; NN: neural network; B: naive Bayes; KNN: knearest neighbor; RA: random search.
Table 5. Global accuracy (%) of the algorithm. Column T (testing) shows the results for Biscay’s data set and Column V (validation) for Silesia. Source: own research by extracting the forecast accuracy of the different models in Table A1. EK: expert knowledge; CE: fuzzy logic; SVM: support vector machine; M: multinomial logit; NN: neural network; B: naive Bayes; KNN: knearest neighbor; RA: random search.
EKCESVMMNNBKNNRA
TVTVTVTVTVTVTVTV
5-TM24.326522.028937.068939.913950.660050.0550.170049.4749.995047.6236.030028.1362.680047.2214.12542.700535
3-TM46.834344.8786749.21845.7252150.880049.3251.070050.6150.950045.7337.990033.6864.250047.5031.27626.26114
Table 6. Bootstraping values (in %) for the modal split forecast by the different models. Source: own research.
Table 6. Bootstraping values (in %) for the modal split forecast by the different models. Source: own research.
ModelTransport MeanRealL.C.IU.C.I.
EKWALK9.5731.1031.12
EKTRANSIT36.330.080.08
EKPRIVATE54.1068.8068.85
CEWALK9.576.896.90
CETRANSIT36.3319.2019.23
CEPRIVATE54.1073.8773.91
SVMWALK9.570.000.00
SVMTRANSIT36.3341.1541.17
SVMPRIVATE54.1058.8158.85
MWALK9.570.110.11
MTRANSIT36.3333.1133.16
MPRIVATE54.1066.7566.77
NNWALK9.574.364.38
NNTRANSIT36.3349.6349.65
NNPRIVATE54.1045.9946.01
BWALK9.5746.1446.16
BTRANSIT36.3321.7921.81
BPRIVATE54.1032.0532.08
KNNWALK9.577.477.48
KNNTRANSIT36.3335.6135.65
KNNPRIVATE54.1056.8856.92
RAWALK9.5744.4644.48
RATRANSIT36.3335.5735.59
RAPRIVATE54.1019.9419.96

Share and Cite

MDPI and ACS Style

Pijoan, A.; Kamara-Esteban, O.; Alonso-Vicario, A.; Borges, C.E. Transport Choice Modeling for the Evaluation of New Transport Policies. Sustainability 2018, 10, 1230. https://doi.org/10.3390/su10041230

AMA Style

Pijoan A, Kamara-Esteban O, Alonso-Vicario A, Borges CE. Transport Choice Modeling for the Evaluation of New Transport Policies. Sustainability. 2018; 10(4):1230. https://doi.org/10.3390/su10041230

Chicago/Turabian Style

Pijoan, Ander, Oihane Kamara-Esteban, Ainhoa Alonso-Vicario, and Cruz E. Borges. 2018. "Transport Choice Modeling for the Evaluation of New Transport Policies" Sustainability 10, no. 4: 1230. https://doi.org/10.3390/su10041230

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop