Original Research ARTICLE
Consumer Behavior in the Choice of Mode of Transport: A Case Study in the Toledo-Madrid Corridor
- 1Econometrics Area, Department of Spanish and International Economy, Econometrics, History and Economic Institutions, University of Castilla-La Mancha, Cobertizo San Pedro Mártir, Toledo, Spain
- 2Business Administration Department, University of Castilla-La Mancha, Talavera de la Reina, Toledo, Spain
Within the context of the consumption of goods or services the decisions made by individuals involve the choice between a set of discrete alternatives, such as the choice of mode of transport. The methodology for analyzing the consumer behavior are the models of discrete choice based on the Theory of Random Utility. These models are based on the definition of preferences through a utility function that is maximized. These models also denominated of disaggregated demand derived from the decision of a set of individuals, who are formalized by the application of probabilistic models. The objective of this study is to determine the behavior of the consumer in the choice of a service, namely of transport services and in a short-distance corridor, such as Toledo-Madrid. The Toledo-Madrid corridor is characterized by being short distance, with high speed train available within the choice options to get the airport, along with the bus and the car. And where offers of HST and aircraft services can be proposed as complementary modes. By applying disaggregated transport models with revealed preference survey data and declared preferences, one can determine the most important variables involved in the choice and determine the arrangements for payment of individuals. These payment provisions may condition the use of certain transport policies to promote the use of efficient transportation.
Transport is one of the most important services of a developed society. The growing need for mobility of people motivated mainly by the spatial differences of the locations to which people need access to makes transport modeling a matter of spatial importance for any developed country. One of the main problems to be analyzed is how individuals move or what their mobility patterns are. This issue is fundamental to the proper planning of the transport system. An efficient transport system must serve the mobility needs of individuals, for this it must use the necessary tools to be able to plan this mobility.
The individual faces, decisions daily between different alternatives of choice, whether goods or services, conditioned by the qualities or attributes of the different options available (McFadden, 2001b; Ye et al., 2007; Chowdhury and Ceder, 2016; Houdek, 2016). In order to determine the variables that affect their choice and to determine the probability of choosing between different available options, disaggregated demand models are used (McFadden, 1981; Schakenbos et al., 2016).
Disaggregated demand models are also called discrete choice models, due to the characteristics variable of choice that can be the mode of transport; this variable is discrete, due to the qualitative nature of the individual's response and also a finite number of responses are shown (Ben-Akiva and Lerman, 1985; Train, 2003; Navarette and Ortuzar, 2013). These models have the ability to predict individual decisions and joint decisions and thus serve as a basis for policy planning (St-Louis et al., 2014).
In order to try to predict the behavior of users regarding modal choice, we will focus on the disaggregated demand models based on the Theory of Choice, in which the traveler maximizes its utility (McFadden, 2001a; Train, 2003). These disaggregated demand models are based on the theory of discrete choice, to determine the probability of choosing the different alternatives which the individual counts with (Martín et al., 2011).
Most of the models used for travel behavior applications are based on utility theory (McFadden, 1974; Domencich and McFadden, 1975; Manski, 1977; Williams, 1977; de Dios Ortúza and Willumsen, 2001), which assumes that the preference of choice of an alternative is captured by a value, called utility, and decision-making selects the alternative in the set of choices most satisfactory (Taniguchi et al., 2014).
This concept, used by the microeconomic theory of the consumer, presents strong limitations for practical applications, since the complexity of human behavior suggests that the decision rule must include a probabilistic dimension (Simon, 1957; Sandoval, 1994; McFadden, 2001b).
Some models assume that the decision rule is intrinsically probabilistic, and even a complete understanding of the problem does not surpass uncertainty (Thurstone, 1927; Quandt, 1956; Luce, 1959; Tversky, 1969; McFadden, 1981). Others consider that the decision rules of individuals are deterministic and motivate uncertainty to the analyst's limited ability to observe and grasp all the dimensions of the election process because of its complexity (Aguado Franco, 2012).
A fundamental assumption in the choice process is that decision making is assumed to have “rational behavior” and implies that a decision maker is a “maximize” of utility and, in the same circumstances, will repeat the same choice, in addition to the transitivity of preferences (Domencich and McFadden, 1975; McFadden, 2001b; de Dios Ortúza and Willumsen, 2001; Train, 2003).
For the estimation of disaggregated demand models, survey data are required. Two types of surveys, revealed preferences, and declared preference surveys can be differentiated. The revealed preference surveys, which try to reflect the current behavior of individuals in their travel decisions and the declared preference surveys, provide us with data that try to reflect how travelers share a certain hypothetical situation (Espino, 2003).
The Toledo-Madrid corridor has HST, bus and car, as alternative modes of transport. In addition, the completion of the Atocha-Chamartín-Airport railway corridor will allow the connection between the Spanish high-speed network and Madrid's international airport. This will favor the cooperation of operators to offer better transport services to travelers (Muro, 2012): new figures of transport tickets, shorter connection times, integrated baggage management, information visualization, etc. (Eurocontrol, 2005b).
In this context, the research work on the “Toledo case” analyzes the current state of mobility in the Toledo-Madrid corridor and, on the other hand, changes the situation in a scenario of new transport infrastructures. This makes it possible to formulate policies at the level of passenger mobility to mitigate the most harmful effects of transport, such as congestion, environmental damage (INFRAS IWW, 1995, 2000, 2004; Intergovernmental Panel on Climate Change (IPCC), 1999; Aviation Environment Federation (AEF), 2000; Whitelegg and Williams, 2001; IATA, 2003; Steer Davies Gleave, 2006; CE Delft, 2008, 2011), by promoting complementarity between modes (Watkiss et al., 2001). Specifically with the promotion of intermodality between aircraft and train, the European Commission aims to address some of the negative effects of transport (Eurocontrol, 2004b).
The Individual Behavior in the Choice of a Service
One of the main objectives of social research is to determine the behavior of the consumer in the decision-making process and this is especially relevant when we are talking about services (McFadden, 2001b), whose supply is neither storable nor cumulative such as transport (Domencich and McFadden, 1975; de Dios Ortúza and Willumsen, 2001).
According to Ben-Akiva and Lerman (1985), an election can be considered as the result of a sequential decision-making process, which includes the following steps: definition of the problem of choice, generation of alternatives, evaluation of the attributes of each alternative, choice, and application of choice.
A specific theory of choice is therefore a collection of procedures and elements based on some general hypotheses. The elements that define the process of decision making are the decision maker and its characteristics, the alternatives or the determination of the options available for decision making, attributes, and decision rules, which describes the process used to choose a decision alternative (McFadden, 2001b; Ye et al., 2007; Martín et al., 2011; Chowdhury and Ceder, 2016; Houdek, 2016).
The decision-making or decision-making unit can be a person or a group of people, such as a family or an organization. If one considers a group of people as a decision maker we can ignore all the internal interactions within the group, and consider only the decisions of the group as a whole. In transport modeling, it has been customary to use aggregate models, which are calibrated with data that has been grouped or aggregated in some form (e.g., using average income per zone), although in the last decades a more disaggregated focus has been considered taking into account the decision unit, the individual (Domencich and McFadden, 1975). In this case, to explain the heterogeneity of preferences among decision makers at the individual level, a disaggregated model should include the socioeconomic characteristics of each decision maker, such as age, sex, and income (Collins and Chambers, 2005).
The analysis of individual decision-making requires not only the knowledge of what has been chosen, but also of what has not been chosen. Therefore, it is necessary to make hypotheses about the available options, or alternatives, that an individual considers during the election process. The set of considered alternatives is called the set of choice. We can differentiate two types of set of choice: continuous and discrete choice sets (Ben-Akiva and Lerman, 1985). The first case defines the set of choice taking into account all possible combinations of goods and services available to the individual1. On the contrary, a discrete choice set contains a finite number of mutually exclusive alternatives that can be specifically enumerated. The choice of a means of transport is a typical example of a choice of a discrete choice set2 (Taniguchi et al., 2014).
The identification of the list of alternatives is a complex process, since not all alternatives will be available all the time for all decision makers. Therefore, it is necessary to differentiate between the universal set of choice, which contains all possible alternatives and the set of choice. The latter is a subset of the universal set of choice available to a particular individual, since alternatives to the universal choice set not available to the individual are excluded.
Normally, to determine these alternatives, deterministic criteria of availability of alternatives are used (for example, having a driving license determines the availability of the alternative of car). In addition to availability, knowledge of the existence of an alternative is a very important factor for the choice made. In this sense, Ben-Akiva and Lerman define that the causes of an alternative is feasible depending on availability in the market, budgetary constraints, available time or other informal factors such as knowledge of the service of a given alternative (Ben-Akiva and Lerman, 1985; Chowdhury et al., 2015; Schakenbos et al., 2016).
Behavior aspects sow uncertainty in the production of alternatives and motivate the use of probabilistic models of choice generation that allow us to obtain the probability of each alternative in the universal set (Swait, 1984).
An attribute is a trait that characterizes an alternative, with a certain value for individuals. The neoclassical economic theory considers that the individuals choose an alternative, from the different amounts of the goods that are included in them. In contrast, the discrete choice theory considers each alternative as a set of attributes and in function of these the attraction of an alternative is expressed by an attribute vector, called the utility function (Ben-Akiva and Lerman, 1985). Moreover, each alternative has its own attributes, and these can be generic for all alternatives or specific for a specific alternative. An attribute is not necessarily an objective quantity that can be measured, but can derive from a subjective measurement from perceptions. Just as the characteristics of the decision maker are taken into account, the analyst must include the attributes of each alternative (Sandoval, 1994).
In decision making, the last of the elements is the decision rule, defined as the process used in decision making, to evaluate the information available on the attributes of each alternative, in the set of choice and to determine a single choice. There are numerous decision rules in the literature that can be grouped into the following four categories (Ben-Akiva and Lerman, 1985):
(1) Dominance: one alternative is dominant over another, if it is better in at least one attribute and not worse in all the others. This rule in many situations can lead to the choice of several alternatives. To avoid this you can add to each attribute values that consider some better than others (minimum levels are defined for the election).
(2) Satisfaction: for each attribute a level is given that serves as a criterion of satisfaction for the individual. If an alternative does not reach the minimum criterion in some of its attributes it can be eliminated as a choice.
(3) Lexicographical rule: it is given that the attributes are ordered according to the level of importance. The decision maker chooses the alternative that is more attractive based on the most important attributes. You can also use this rule to eliminate the worst alternatives in each attribute in order of importance. The combination of the decision rule, in terms of satisfaction and lexicographic rules, is called “elimination by aspects”. This process begins with the most important attribute and eliminates the alternatives that do not reach the level of satisfaction. If two or more alternatives are available, the second most important attribute is continued and so on (Tversky, 1972).
(4) Utility: This type of decision rule assumes that the attributes are measurable. This means that the attractiveness of an alternative is expressed by a vector of attribute values by means of a scalar, called utility. The individual will seek to maximize this utility (or minimize its costs).
These four decision rules, at the same time can be grouped into two types of behavior: compensatory or non-compensatory behavior. The compensatory type takes into account the set of all the attributes, so that changes in one or more attributes can be compensated by others, and in case of a decrease in one attribute can be compensated for an improvement in another, as is the case of the utility rule. In non-compensatory behaviors, rules or levels are defined to restrict the choice of some alternatives. Among these are the first three decision rules described (McFadden, 2001b).
The classical (neoclassical) and discrete-choice theories analyze individual consumer behavior (Ben-Akiva and Lerman, 1985), from a utility function that individuals must maximize, based on the analysis of revealed preferences, although with important differences between them, in the definition of decision rules and alternatives.
Most transport choice models have been based on this behavior (utility maximization), although other decision rules have been tried (Cantillo and Ortúzar, 2005), especially for the case of the declared preferences, where the choices that do not follow the principles of maximum utility are measured (Sælensminde, 2002; Rouwendal and Blaeij, 2004).
The mobility of people is a complex phenomenon, due to the large number of factors that influence the decisions of individuals. For this reason, it is necessary to analyze the displacements according to three interacting aspects a spatial approach, a social approach and a perspective approach (Rodríguez, 1991).
The spatial approach focuses on the different land use and the distribution of activities in space. This creates the necessity of the individual of the individual to move, depending on the location of the activities you want to perform.
The social approach refers to the fact that an individual's movements are exclusive of their characteristics and are therefore the result of the socioeconomic characteristics of the individual who performs it. Because of this it is necessary to have very present variables like the sex, the age, cultural level, income, health, etc.; for the analysis of the mobility.
Finally, we can speak of a perceptive approach, which refers to the image that the individual has formed. This image varies in each individual or group of people, leading to different assessments of the same data in decision making.
In the classic model a sequence of decisions of the individuals is done in four stages which are the generation of trips, zonal distribution, modal distribution, and assignment, but another sub-model has recently been included that is the last one mentioned, the election of the hour that gives place to the models of time distribution.
At present, this method is recognized to be too strict, since the decisions of the individuals are not taken following this sequence, but each sub-model depends on the type of the function of utility assumed to explain all these choices of trip. There are some current approaches that differ from this four-step sequential methodology and simultaneously address the stages of choice of frequency, destination, and mode of travel, but they are still at a level of research and have not been implemented. In addition, we must highlight models based on the activities of families or decision-making centers, which take into account constraints on budget and time choice (de Dios Ortúza and Willumsen, 2001) stand out.
Methodological Framework for Modeling Consumer Choices
From an economic perspective, consumer theory is the economic modeling of the behavior of an economic agent that consumes goods and services. This theory relates preferences, indifference curves and budget constraints to consumer demand curves.
The fundamentals of the individual choice theory of the neoclassical model of consumer behavior are based on the fact that the individual chooses an amount of goods, which form their basket of goods, in order to maximize their utility which translates into their level of satisfaction, subject to their income restriction (Torres, 2008). Therefore, given a consumer whose preferences fulfill these assumptions, there is a deterministic utility function, U, which represents these preferences that will be ordered so that the individual will choose the alternative that gives him greater satisfaction3.
Given that consumer economic theory has been developed without taking into account the nature of the alternatives, some extensions have been made such as Strotz (1957, 1959), which analyzes consumer behavior by introducing the concept of “utility tree” (Strotz, 1957, p. 270). Muth (1966) proposes that the goods and services of the market are considered by the consumer as inputs of a domestic production function, whose production is the satisfaction of the consumer needs derived from consuming certain goods or Becker (1965) that extends the traditional theoretical formulation by adding the time constraint, in addition to the budget, to the production function.
Finally, one of the most important extensions is that one made by Lancaster (1966). This author considers three assumptions that break with the traditional approach and states. That (1) goods per se do not contribute to usefulness, but are the characteristics they possess, which provide usefulness to the consumer, (2) a good has more than one characteristic, and many of these characteristics can be shared by more than one good; (3) and through the combination of different goods different characteristics can be obtained from those corresponding to the separate goods. These assumptions represent the definition of utility in terms of the attributes of the goods.
Taking into account all these assumptions, consumer economic theory has important limitations to explain the consumer behavior, which justifies models based on Discrete Choice Theory. This reformulation of the behavior of the consumer will not occur until 1981, with the inclusion of consumer goods of a discrete nature (McFadden, 1981), with which the Theory of Discrete Choice begins, whose characteristics are that:
• The individual consumes continuous (or divisible) goods.
• The individual chooses between a set of mutually exclusive discrete alternatives (non-divisible goods).
• The discrete alternatives are represented by a feature vector Qj.
The problem that arises when using Discrete Choice Theory is that it represents a mechanism of deterministic choice and these deterministic election mechanisms do not fit the analysis of choice problems in real situations (Block and Marschak, 1960). In this line of study, this approach has been widely criticized, both in the field of psychology by Thurstone (1927), Luce (1959), and Tversky (1969), as in the field of economics by Quandt (1956) and McFadden (1981). These authors justify the inclusion of uncertainty in modeling, being the source of randomness determinant when specifying the different models.
Therefore probabilistic mechanisms of choice are required in order to analyze individual options. In order to do this, we use Probabilistic Election Theory, specifically to the Theory of Random Utility formalized by McFadden (1974, 1981), Domencich and McFadden (1975), and Manski (1977) and whose most recent developments are by (de Dios Ortúza and Willumsen, 2001).
The Random Utility Approach assumes that the individual always selects the most useful alternative, although the utilities are not known to the analyst and are treated as random variables. In mathematical terms, this is expressed by separating the total utility Uin into a deterministic component Vin, called Systematic Utility and a random component εi that captures uncertainty, called Random Perturbation:
Vjq is the systematic component that represents an appreciable part of the utility of the individual q to choose alternative j. Systematic utility depends on the attributes of the alternative and the socioeconomic characteristics of the person.
εjq is the random term that accounts for unobserved factors.
The individual chooses the most useful alternative. Alternative i is chosen, if and only if:
The analyst can only obtain the probability of choosing alternative i as:
This is the probability that each random term εnj – εni, is less than an observed quantity Vni − Vnj, so it is a cumulative distribution function. From the density function f(εn) we can calculate this probability as:
Being I (·) a function that has a value of one if the term given in brackets is true (if the individual has chosen alternative i) and zero in another case. According to the different hypotheses formulated about the distribution, εnj will have a closed value of this integral (Simple Logit or Hierarchical Logit) or it will have to be evaluated numerically by simulation (Probit or Mixed Logit), and therefore the different Models of discrete choice (Train, 2003).
Defined the set of choice, we can proceed to define the utility function that will give rise to different types of models. The typology of this type of model is broad and can be classified (Medina, 2003 according to the number of alternatives of the endogenous variable (dichotomous response models and multiple response models) or according to the type of function to estimate the probability (Linear Probability Model, Logit Model, and Probit Model): Whether the alternatives are exclusive or incorporate ordinal information (models with non-ordered data and models with ordered data) and whether the regressors refer to aspects of individuals or to alternatives in non-ordered models (multinomial and conditional models).
For the theoretical justification of the models of discrete choice based on the Theory of Random Utility, we can find three types of models of discrete choice depending on the different hypotheses that are taken for the distribution of the random term. First the linear probability model, assuming a uniform distribution, then the Probit Model, assuming a normal distribution and in third place the Logit Model, assuming a logistic distribution.
The Linear Probability Model is the linear fit regression model that is applied to a binary dependent variable. This model is estimated by ordinary minimums to the square and is easy to estimate and interpret. The estimated parameters measure the predicted change in probability of success, vs. a unit increment of Xi. Although there are problems in the estimation of the regression model when the endogenous variable is binary. The specific problems encountered with this model are heteroskedasticity of the perturbation term, the predicted probabilities are inconsistent, since it cannot be guaranteed that they are bounded between zero and one, the non-normality of the perturbation and the coefficient of Determination is not appropriate. Due to these problems, what is interesting is a model that reproduces properly the behavior of a probability function and the alternatives we find the Logit model and the Probit model, similar numerically. When the normal distribution is used as a probability function, the so-called Probit model is obtained, while the use of the logistic distribution provides the Logit model.
The main advantage of the Probit model is its ability to capture all correlations between alternatives, but its main problem is the complexity of its formulation, so there are very few applications that have been developed (Daganzo, 1979) and more than three alternatives are used for the calculation Simulation procedures (Train, 2003).
The logistic regression model, also called Logit models, is much more popular thanks to its analytical flexibility. Thus, the possible hypothesis for deriving some Logit models from the operation of individual choice is analyzed here: the most general case of interest, the Logit Multinomial model (McFadden, 1974), in its particular case, beyond binary situations; The Logit Hierarchical or Nested model (Williams, 1977) and the Logit Mixed model (McFadden and Train, 2000).
Modeling the Demand for Services: Application to the Choice of Transport Services
Transport contributes significantly to the economic development and allows the market to function in a global way. It should be noted that most modes of transport do not affect society only in a positive ways, there are negative side effects such as congestion, noise and air pollution.
The promotion of intermodality using airplanes and trains is intended to solve some negative effects of transport, such as the impact of congestion on the environment, economy, safety and passengers (Eurocontrol, 2004b). For example, regarding to the access to the airport with the promotion of intermodality, beneficial effects are expected on the economy, especially in the regional economy, and on the environment are expected (Watkiss et al., 2001; Li and Hensher, 2012).
The estimation and internalization of external costs of transport have been important issues for research in the transport sector in the European Union. A Concern for the environment has grown over the last decade and the European Commission has raised the issue of internalization in several of its directives, such as the 1995 Green Paper (European Commission, 1995), the White Paper on 2001 (European Commission, 2001), and the Mid-Term Review in 2006 (European Commission, 2006). The latter two, underline the need for equitable and efficient pricing in terms of external costs. On the other hand, although the impact on the environment or congestion are indisputable, in the literature, there are some reports that try to quantify them, some of them are: (Intergovernmental Panel on Climate Change (IPCC), 1999; Aviation Environment Federation (AEF), 2000; INFRAS IWW, 2000, 2004; Whitelegg and Williams, 2001; IATA, 2003; Steer Davies Gleave, 2006; CE Delft, 2008) and (CE Delft, 2011).
In general, all of these investigations seek to establish pricing systems that capture the external costs associated with transport, so as to reduce negative impacts, in order to improve the efficiency of the transport system, to ensure fairness between modes and improve safety and sustainability.
In recent years, considerable progress has been made in the construction of transport demand models based on the theoretical principles of choice. Within random utility theory, it can be shown that the structure of the models depends on the perceived similarity between discrete choice alternatives. Moreover, this aspect can be interpreted mathematically in terms of the correlation between the components of random utility functions. The use of Logit-type disaggregated demand models allows us to determine the likelihood of choosing one way or another depending on the previously defined attributes, as well as quantifying the willingness to pay of these subjects based on the different characteristics previously given.
To do this, we must first determine which are the parameters that measure the utility of the services or that define the attributes of the service (Martín et al., 2011; Navarette and Ortuzar, 2013; Chowdhury et al., 2015; Schakenbos et al., 2016; Ettema et al., 2017). In order to select the mode of transport, the factors influencing the choice of mode can be classified into three groups (Eurocontrol, 2004a) by the characteristics of the mode of transport (vehicle availability and/or ownership, driving license, home structure, income), the trip characteristics (the mode of choice is heavily influenced by the purpose of the trip, time of travel), and the characteristics of transport infrastructures and services (relative travel time, relative monetary cost, comfort, and convenience, reliability and regularity, protection, and safety).
The information necessary for the estimation of these models is obtained through surveys of the mobility of the different individuals at a given moment (Martínez and Muro, 2011). Two types of surveys can be differentiated into revealed preferences and declared preferences.
Revealed preferences (RP) show data on the current behavior of individuals and give us information about their travel decisions. These provide us with descriptive information about the characteristics of the traveler on a particular route. Until the mid-1980s this type of data was the most used to model the transport demand. The main drawbacks to the use of this data in the modeling are the costs of the sample and its limited capacity to understand the behavior of the traveler, the observations of the actual elections may not have sufficient variability for the construction of good models. Because of this it is difficult to be able to detect the relative importance of certain factors on behavior; (de Dios Ortúza and Willumsen, 2001; Espino, 2003; Espino et al., 2006). For the Toledo-Madrid case, information is obtained on the type of mobility of travelers in the current situation with a socio-economic profile and mobility parameters (Martínez and Muro, 2014).
The stated preferences (SP) give information on the behavior of the individual before certain hypothetical situations raised by the researcher (Navarette and Ortuzar, 2013). Unlike RP data, which give information about the trips that an individual habitually makes, they inform about the trips that the individual would realize in certain conditions. The SPs started in the field of market studies and began to be applied to the field of transport demand analysis in the late 1970s. The possibility of designing SP experiments allows, theoretically, to solve the problems they present The RP (de Dios Ortúza and Willumsen, 2001). One of the problems presented by these models is that the analyst cannot assure that the individual performs what he has answered in an SP survey, so it is important to construct realistic alternatives and design comprehensible exercises to present the individual.
In an SP exercise, three elements can be mainly distinguished (Espino, 2003). The first are the situation in which the individual finds himself to declare his preferences which can be a real situation (a journey that is carried out at that moment) or hypothetical (a journey that would take place in the future given a series of conditions), and constitutes the context of the decision. Second, the alternatives that are usually hypothetical are selected, although some of them might exist today, they are presented in the exercise as a function of a set of attributes. Third, it defines how individuals can state their preferences. The most frequent techniques are the hierarchy of the answer, the punctuation of the alternatives and the choice between alternatives.
In the case of Toledo-Madrid, the different alternatives that the user has to carry out this journey are traveling by public transport, either by High Speed Train (HST) or by regular bus (BUS), or by private transport, Private vehicle (CAR).
The basic methodology for the presentation of the options of the declared preferences (SP) is the realization of an efficient design of the scenarios. Specifically, the declared preference questionnaire is designed following a time and price reduction scheme, in three blocks or scenarios, to determine the variation of the individuals' choice. The SP surveys are structured with two approaches: a first comparative approach between the bus and the HST and a second approach comparing the car and the HST. The scenarios are conformed with the variables of trip price, travel time, comfort and cost of parking, this last variable referred only to the case of the car.
The SP survey proposes a simulation exercise is exposed of alternative scenarios that layout the mobility in the future. For the elaboration of these scenarios the combination of price and travel time variables are used. Prior to the establishment of the options, a characterization of the services and infrastructures currently available for each mode of transport is made. This current quote will serve as a basis for calculating the different alternatives for the arrival to Madrid Airport from Toledo and choose the best.
Thus, the methodology followed was to calculate the price and average travel time of each mode of transport from the origin, which in this case is Toledo, to the destination in the current situation and compare it with the different hypothetical scenarios. Therefore, the scenarios that are presented in the SP questionnaires are calculated based on the average prices and times calculated in that route, and compared in each pair of transport modes (Muro, 2012).
Consequently, three route scenarios are established:
• Toledo-Airport without trans-shipment with the of the new infrastructure, making the trip in HST from Toledo to Madrid Airport making two stops, one at Atocha station and another at Chamartín station but without any transfer.
• Toledo-Airport with a transfer: with the new infrastructure, making the trip in HST from Toledo to Madrid Airport by making a stop at Atocha station and another stop with transshipment at Chamartín station.
• Toledo-Airport with two transfers with the use of the new infrastructure, making the trip in HST from Toledo to the Madrid Airport carrying out a stop with transshipment in the Atocha station and another stop with transshipment in the station of Chamartín.
The chosen bus and bus (CAR) optimum alternatives are fixed and compared in blocks of three, according to the three scenarios raised for the HST. The prices and travel times for the mode of transport bus and private vehicle are constant for all scenarios and are based on current values. The prices of the bus are calculated according to ticket prices and the prices of the private vehicle are calculated by applying a cost per kilometer journey that incorporates the maintenance of the vehicle, insurance, etc., adding to the cost of the toll. This price does not include parking costs since it is presented as a different variable measured in euros per day of parking.
In summary, for the calculation of these values of time and price, the following hypotheses have been taken into account:
1. Existence of a single ticket for the case of the HST, with information on possible transfer for all scenarios. Thus, the time of transfer between trains is minimized if necessary.
2. The travel time is composed of the sum of the time from the origin, the waiting and access times, and the travel time.
3. The prices related to the different stages of the trip in HST are calculated as variations of the base price, which is calculated based on the current trip.
From the econometric point of view, the main difference of both types of data (RP and SP) are the types of error that each present. The RP data suffer from errors in the measurement of independent variables and those of SP present errors in the dependent variable, motivated by the doubt about whether the individual will actually perform what he is declaring. The joint estimation of both types of data is based on the different errors mentioned and can be specified if we consider error terms with different variance. Ben-Akiva and Morikawa formulate mixed data estimates that combine revealed preferences and stated preferences with the aim of avoiding their disadvantages and taking advantage of different sources of information (Ben-Akiva and Morikawa, 1990).
Specification and Estimation
In order to explain the mathematical approach of the models of discrete choice applied to the case of transport choice, one can start with the assumption based on the choice of the individual between two available alternatives (binary response models), and then extend this assumption to the j alternatives available (multiple response models). The dichotomous response models have two categories and these usually indicate that an event has occurred, that some feature is present or that an option is chosen.
In the modeling of the choice of transport, the objective will be to analyze the choice of the different alternatives that the user has to make this journey, and these are for example to travel by public transport, either in High Speed Train (HST) or by regular bus (BUS), or by private transport, in a private vehicle (CAR).
As this study focuses on the probability of choosing a mode of transport, the available alternatives will be the high-speed train against that the bus and the car. This allows specifying and estimating binary logistic regression models that will serve as the basis for a first analysis of the goodness of the database. Thus, two models are specified that will have a functional form of binomial logistic regression model, whose alternatives will be HST-BUS and HST-CAR, respectively.
The functional form of the binomial logistic regression models with the HST-BUS / HST-CAR alternatives will be represented by:
Where U is the total utilities of the high-speed train (a) and bus/car (r), respectively.
If the deterministic utility component, V, is assumed to be linear with regard to the parameters, as in most disaggregated model approaches, then:
, is the vector of coefficients.
, is the vector of independent variables.
So if we have two alternatives, for example the binomial Logit model will be represented
If we suppose that α = 1, since in the case linear utilities with regard to the parameters, the parameter α cannot be distinguished from the general scale of the β.
Taking the logarithm on both sides, from the equation above, we obtain a formulation that is suitable for applying linear regression techniques:
Since the probability of choosing the alternative r is given by: P(r) = 1 − P(a)
The final functional form of the model with an additive linear utility function is as follows:
The main explanatory variables that are considered in the transport demand analysis are, differential between the attributes and the socioeconomic variables (de Dios Ortúza and Willumsen, 2001; Cokasova, 2003; Espino, 2003; Eurocontrol, 2005a; Steer Davies Gleave, 2006; Muro, 2012) which are:
(1) Attribute variables:
- Travel cost:
The price variable price measured in euros will be the price paid for each of the means of transport used as public transport, and for the private vehicle, the cost of gasoline based on the kilometers traveled, including the cost of the amortization of the vehicle (insurance, maintenance, etc.). The travel cost variable represent the price paid, in euros, for the route according to the form of transport chosen. In addition to the private car, the price of the car is different, so the total cost of the car will be the sum of the cost of the trip plus the cost of parking.
- Travel time:
The total travel time is one of the most important factors that the traveler takes into account for choosing a mode of transport (Steer Davies Gleave, 2006). It is therefore a major factor to increase the demand, as travel time is reduced (IATA, 2003). The variable of travel time measured in minutes can be represented as the total time of the trip or by the time dedicated to each of the stages of the journey, such as waiting time, transfer, etc. The variable time considered in this study is the total time spent on the trip done by the chosen mode of transport.
This variable measures the level of service and is represented by dummy variables, which include values whose coding contains the different levels of service perceived by the individual. In this case, comfort tries to measure how the individual perceives the transshipments that are presented as hypothetical scenarios of the future services of the HST to Madrid airport (for the rest of the alternatives they would be the current travel options). Comfort is the dummy variable defined as the number of transfers chosen for the alternative train. So that in the first scenario, with no transfers, the comfort is 1; in the second scenario, with one trans for, the comfort is 2; and the third scenario, with two transfers the comfort is 3. If you do not choose the HST option the comfort level will be 0.
(2) Socioeconomic variables:
The variable sex is included as a dichotomous dummy variable, to analyze the influence of them in the decisions
- Level of income was specified by the individual corresponding to the monthly salary range was divided into five categories being the first of an income of up to 600 € then income in the range of 601 to 1,000 €, following that in the range of 1,001 to 2,000 next in the range of 2,001 to 3,000 € and last above 3,000.
This variable is expressed in four ranges: the value 1 being from 18 to 33 years; 2 from 34 to 49 years; 3 from 50 to 64 years, and 4 over 65 years. So we are in the specification of three dummy variables that represent value one when it meets that range and zero otherwise.
The utility specification is expressed as a Logit Binomial model, to measure the choice of mode of transport between the HST/BUS (MNL1) alternatives, in the Toledo-Madrid corridor (Madrid Airport):
Utility HST1: UHST1 = ASCHST1 + θPPHST1 + θTTHST1
Utility BUS: UBUS = ASCbus + θPPbus + θTTbus + ∑θCjC
ASCHST1 and ASCbus, the specific constants of each alternative, in which ASCHST1 is fixed.
θP, θT, and θCj are the parameters associated to each of the explanatory variables. In this specification, the two most important variables (cost, θP and time, θT), where C, refers to all socioeconomic variables and characteristics of travel such as income, sex, age, motive travel, and comfort.
Cj, is a generic variable that represents each and every one of the variables not explicitly collected but mentioned above, where j refers to S, sex; ING to income; and COM, comfort depending on the transfers.
PHST, Pbus, is the price of alternatives HST1 and Bus, respectively.
THST, Tbus, is the travel time of the alternatives HST1 and Bus, respectively.
The specification of the utilities of the two alternatives of the HST-CAR Model (MNL2) will be:
Utility HST2: UHST2 = ASCHST2 + θPPHST2 + θTTHST2
Utility CAR: UCAR = ASCcar + θPPcar + θTTcar + θPAPAcar + ∑θCjCj
Being, ASCHST2; ASCcar, the specific constants of each alternative.
θP; θT; θPA, and θCj are the parameters associated to each of the explanatory variables.
Cj, is a generic variable that represents each and every one of the variables not explicitly collected but mentioned previously.
PHST, Pcar, is the price of alternatives HST2 and CAR, respectively.
THST, Tcar, is the travel time of alternatives HST2 and CAR, respectively.
The MNL1 models Table 1 are the result of the estimation of the data of the alternatives HST1 and BUS. The MNL1.1, MNL1.2, MNL1.3, MNL1.4, MNL1.5, MNL1.6, and MNL1.7 models are variants of specified models with different explanatory variables. In general, these variables are statistically significant, as it shows in Table 1, and with correct signs (negative values for its inverse relationship).
The results of the first binomial models show low determination coefficients, which are a reflection of the inadequate fit of the Logit Binomial model at around 1% in the case of the HST-BUS alternatives (MNL1.1: 0.057; MNL1.2: 0.080; MNL1.6: 0.81). This can be caused by the lower amount of people who chose the BUS option vs. the HST option.
The MNL2 models Table 2 are the results of the estimation of the data of the alternatives HST1 and BUS. The MNL1.1, MNL1.2, MNL1.3, MNL1.4, MNL1.5, MNL1.6, and MNL1.7 models are variants of specified models with different explanatory variables. In general, these variables are statistically significant, as it is shown in Table 2, and with correct signs (negative values for its inverse relationship), except in the MNL2.2 where we can find the travel cost with positive values. For this reason this model is nullified and in the final models only the cost of the total trip is included.
The MNL2 models (Table 2) are the results of the estimation with the data of the alternatives HST2 and CAR. The MNL2.1, MNL2.2, MNL2.3, MNL2.4, MNL2.5, MNL2.6, and MNL2.7 models are the same as the previous variants of the original. These models show determination coefficients around 53% (MNL2.1: 0.522, MNL2.3: 0.522; MNL2.7: 0.531).
From the results of the estimated models we proceed to calculate the payment arrangements of the users. These payment arrangements are calculated as the quotient of the parameters estimated for time and price, so they are interpreted as the willingness to pay or waiting to save in travel time (Espino et al., 2006).
The provisions for the payment of the binomial HST-BUS, fluctuate between 0.397 and 0.431, which means that individuals would be willing to pay from 23 to 25 € to save an hour of travel.
In the case of Logit Binomial HST-CAR models, Table 2, the results are very high being 600 € per hour what is saved in travel time. This may be motivated by individuals who answered in a lexicographic form, in which regardless the options, the users have chosen car, although this was but in all the options and by the difficulty of capturing the usual users of the car.
Next, if we assume that there are three available alternatives of choice for the individual in this path, a model is proposed with three alternatives, which are HST, bus and car, whose functional form will be that of a multinomial logistic regression model.
The specification is based on the approach of four alternatives, two from each of the sub-databases:
Utility HST1: UHST1 = ASCHST1 + θPPHST1 + θTTHST1
Utility BUS: UBUS = ASCBUS + θPPBUS + θTTBUS
Utility HST2: UHST = ASCHST2 + θPPHST2 + θTTHST2
Utility CAR: UCAR = ASCCAR + θPPCAR + θTTCAR + θPAPACAR + ∑θCjCj
The joint treatment of both databases requires, in the specification of the multinomial model, a specific treatment and equal to the problematic applicable to Mixed data (RP and SP).
Thus, the Logit Multinomial (MNL3) model with the HST-BUS-CAR alternatives will have the form:
The probability of choosing each of the alternatives will be:
Therefore, a multinomial Logit model (MNL3) is specified, with sub-bases 1 and 2, so that a specification is made with the mixed data processing with different errors and therefore variances, so that ε is the stochastic error Of the HST-BUS data (Base 1) and η the data of the base HST-CAR (Base 2). It can be expressed as: , where μ is an unknown parameter.
So the utilities of both databases are expressed
θ, α, and w are the parameters to be estimated.
and are common attributes of alternative j for databases 1 and 2, respectively, while and are non-common attributes of alternative j for each data set.
By multiplying the utility function of the Base 2 data by the unknown parameter μ, what is obtained is that the stochastic error of this type of data has the same variance as the data of Base 1.
Taking into account these utility functions we are allowed to homogenize the type of error homogenized, due to the multiplication of the SP parameters where by multiplying the utility function of the sub-base data 2 by the unknown parameter μ, we ensure that the stochastic error of this type of data has the same variance as the data of sub-base 1.
In mixed data, it is assumed that SP data should have more noise than RP data. If this is the case, the value of μ, that is known as the scale coefficient of the model, will be between 0 and 1. If the value is >1 it would indicate that the data with the highest noise level are those of RP. In our case and due to the errors detected, we assume that base 2 (HST-CAR) will have more noise (Ben-Akiva and Morikawa, 1990)5.
Taking into account the mixed data (MNL3.1) the specification of the model would be: UHST1 = OP * (ASCHST1 + BETA1 *PHST1 + BETA2 *THST1) + (1-OP) *μ* (ASCHST2 + BETA1 *PHST2 + BETA2 *THST2) UBUS = OP * (ASCBUS + BETA1 * PBUS + BETA2 *TBUS) + (1-OP)* μ* (ASCBUS + BETA1 *PBUS + BETA2 *TBUS) UHST2 = OP * (ASCHST1 + BETA1 *PHST1 BETA2 *THST1) + (1-OP) *μ* (ASCHST2 + BETA1 *PHST2 BETA2 *THST2) UCAR =OP * (ASCCAR + BETA1 * CTCAR + BETA2 *TCAR) + (1-OP) *μ* (ASCCAR + BETA1 * PTCAR + BETA2 *TCAR)
We proceeded to estimate multinomial models with three and four alternatives. The results of the estimation of the models are shown in Tables 3, 4 (Bierlaire, 2003), for the multinomial simple models (MNL) with four and three alternatives, respectively.
Subsequently, models that combine the two databases have been estimated applying the specific problem of mixed data (Ben-Akiva and Morikawa, 1990), estimating Multinomial Logit models, with three and four alternatives. Once detected that the scale value of mixed data is significant, which means that there are differences in errors in the two databases, we start scaling the data, verifying that the results improve with the scaled data. From this moment on when estimating the Multinomial Logit models, the data of the sub-base 2 are scaled by this calculated coefficient in the model MNL3.1.
In Table 3 it can be verified that the variable Theta (μ) can be verified to be statistically significant. For this reason, the database 1 must be scaled in order to be able to establish alternatives HST1 and HST2 as a single HST alternative. Thus, the three estimated alternatives models will have the scaled data. (Tables 4, 6).
In Table 4 the estimated model is shown with three alternatives. From the MNL4.1 model with three alternatives without scaled data we can verify how in subsequent models with scaled data the results are improved in terms of Pseudo-R2 (ρ2). Therefore, from the scaled data the proposed specifications are estimated with three alternatives of Multinomial Logit and Mixed Logit models. From each of these results shown the provisions to be paid are calculated as the quotient of the estimated parameters for time and price.
Finally, the specification of a Logit Mixed model (ML), of fixed parameters, is carried out. The Logit Mixed model with error components is used when analyzing arbitrary substitutability or correlation patterns. In this case the utility of the alternative i is specified as:
α, is a vector of fixed coefficients and represents the parameter vector of the explanatory variables of alternative i.
μq and εiq, are independent, being μq, is a vector of random coefficients with mean 0 and covariance W y εiq, are random variables, iid Gumbel(0, β).
With this formulation the random term is:
It can be verified that with this formulation there is correlation between any pair of alternatives:
The probability of choice can be obtained in a similar way to the previous case, conditioned by the value of μq. Due to the fact that there is no guarantee of being able to solve the integral, the probability is obtained by simulation.
Both models (of random parameters and error components) are equivalent when the vector parameter β breaks clown into its mean α plus the deviations μq (βq = α + μq). Reciprocally, if ziq = xiq the error component model is equivalent to a random parameter model: and if ziq = xiq, we would obtain a random parameter model with fixed coefficients for xiq and random coefficients with zero mean for ziq.
The specification of a Logit Mixed (ML) model, with a panel effect to analyze the correlation between variables, using a fixed effects model, where four alternatives are specified first and further on three are specified.
First, a Logit Mixed model (ML5.1) is proposed to determine the panel effect, with mixed data and with the same error component for the two databases, so that the model will be:
Being: OP = 1 with values of sub-base 1 and 0 for values of sub-base 2 (OP-1 = Base 2)
And later it is specified with different errors for each base (ZERO_SIGMA1 and ZERO_SIGMA2)
The specifications of the three alternatives are performed, first without scaling the data, ML 6.2:
And with the Base 2 data scaled, with the scale factor obtained in the MNL3.1 estimate (μ = 0.0830), the model is specified with three alternatives such that ML.6.4:
Table 6. Results of ML models estimated with three alternatives, with panel effect and with scaled data.
The results of the Multinomial Logit models with three choice alternatives (HST-CAR and BUS) improve in terms of the coefficient of determination, showing more realistic payment arrangements. Even in the case of when we include the comfort variable (measured according to the number of transfers), together with the total price of the trip (which separately includes the cost of the trip and the cost of parking for the car, is significant. As well as that obtained from models whose variables are not significant and with signs of time not expected, the results improve in a substantial way.
Conclusions and Recommendations
This article shows an application of the demand models or discreet choice of Logit type to obtain the provisions for the payment of the consumers of a transport service, specifically in the corridor Toledo-Madrid/Madrid Airport. In this case it is analyzed what will be the choice of the travelers that make this journey, considering that they have three alternatives at present, the high speed train, car and the bus.
The theoretical basis of these emodels starts from the theory of random utility specifically in the theory of the maximization of utility. This model is based on the existence of a rational consumer that represents the average behavior of the set of consumers, based on a series of parameters that measure that utility. In the case of transport, these variables are usually the travel time, the cost or price of the service and the convenience of the trip, as well as different socioeconomic variables of the individual.
The model data are obtained through a survey of declared preferences. In it an exercise of simulation of alternative scenarios is proposed that configure mobility in the future. This mobility is based on the new infrastructure of connection with the airport. So there are different service options on the route which are two transfers, one transfer or no transfer.
In general, of the results presented above, we highlight the variables of the equation and the measures of goodness of the model.
Logit models have been estimated that measure the probability of choosing the HST in front of the bus or the car independently and according to the socioeconomic characteristics of the subject and the attributes of the different modes of transport. The variables included in each of the models are: time, price, sex and income. Thus, the attributes for each of the transport modes are (1) HST, price and travel time, (2) bus, price and travel time and (3) car taking into account price, parking cost and travel time.
Specifically, the results show that the willingness to pay HST-Car binomial models is too large (600 € per hour saved on the trip), motivated by individuals who have responded in a lexicographic way and motivated by the smaller number of people who have chosen option bus in front of the HST option. The provisions for the payment of the binomial HST-Bus, are more reasonable and oscillate between 23 and 25 € per hour.
In multinomial models the provisions for payment are more realistic, except when we include the comfort variable, together with the total price of the trip, with which models are obtained whose variables are not significant and even with the signs of the time not expected, estimating The model without independent terms with which you get a result of 26.3 € per hour.
From the results of the modeling we can analyze how the consumer is willing to pay 26.3 € for each hour of travel saved which gives us information about the importance of travel time over the price paid.
In conclusion, based on this information, a sustainable transport policy focusing on HST-airplane intermodal transport is recommended in order to facilitate the reduction of travel times, by encouraging the use of efficient public transport as the high-speed train to reach the airport.
A transport policy focused on replacing private transport with public transport (HST), benefits, on the one hand, users who will have shorter travel times and in general, the transport system to reduce the congestion of access to Madrid and to the Madrid Airport (Li and Hensher, 2012).
The main variables to encourage the integration of air and high-speed modes of transport are the visibility of the offer, integrated management of the reservation and boarding passes, check-in and baggage control, full travel responsibility, the management of passenger loyalty programs, ease of connection in intermodal mode, travel time and price (Muro, 2012). As for this offer, and above all regarding the connections between Toledo and Madrid airport, a new line of direct buses to the airport has emerged (Muro and Pérez, 2016) and recently, in February 2017, the offer of combined train and flying tickets (Train and Fly)6 from the city of Toledo, which indicates improvements in air-rail intermodality (air-rail intermodality).
AM: modeling and coordination. IP: methodological framework. SG: introduction and conceptual framework.
This paper has been funded by Aid Research Groups at the University of Castilla-La Mancha 2017. Research Group “Observatorio de la Innovación en Distribución Comercial”. This paper has been funded by The Ministry of Economy and Competitivity (Spain), Research Project with reference: ECO2014-59688-R, Programa Estatal de Investigación, Desarrollo e Innovación Orientada a los Retos de la Sociedad, Plan Estatal de Investigación Científica y Técnica y de Innovación 2013–2016. And the revision of the translation has been financed with the collaboration of the contract-program of the Faculty of Juridical and Social Sciences of Toledo.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are grateful to our teacher, Professor Timoteo Martínez Aguado, for the collaboration and help he has always given us and, in particular, for his contributions in this work.
1. ^Based on the neoclassical theory of consumer behavior.
2. ^From this set of decision we derive discrete choice models, which we will see below.
5. ^Assuming that the errors distribute as Gumbel with zero mean and different variance.
Aguado Franco, J. (2012). Análisis del Comportamiento del Consumidor: Egoísmo, Altruismo, Cooperación Y Otras Posibles Motivaciones Sociales, Contribuciones a la Economía. Available online at: http://EconPapers.repec.org/RePEc:erv:contri:y:2012:i:2012-11:11
Block, H. D., and Marschak, J. (1960). “Random orderings and stochastic theories of responses,” in Contributions to Probability and Statistics, eds I. Olkin, S. G. Ghurye, W. Hoeffding, W. G. Madow, and H. B. Mann (Stanford University Press), 97–132.
Chowdhury, S., Ceder, A., and Schwalger, B. (2015). The effects of travel time and cost savings on commuters' decision to travel on public transport routes involving transfers. J. Trans. Geogr. 43, 151–159. doi: 10.1016/j.jtrangeo.2015.01.009
Espino, R. (2003). Análisis y Predicción de la Demanda de Transportes de Pasajeros: Una Aplicación al Estudio de dos Corredores de Transporte en Gran Canaria. Tesis doctoral, Universidad de las Palmas de Gran Canaria.
European Commission (1995). Towards Fair and Efficient Pricing in Transport: Policy Options for Internalising the External Costs of Transport in the European Union. Green Paper, Directorate-General for Transport-DG VII. Bruselas, COM (1995) 691.
European Commission (2006). Por una Europa en Movimiento. Movilidad Sostenible Para Nuestro Continente. Revisión Intermedia del Libro Blanco del Transporte de la Comisión Europea de 2001, 22 Junio de 2006. COM (2006) 314 Final.
CE and Delft, Infras Fraunhofer ISI. (2008). Handbook on Estimation of External Costs in the Transport Sector, Produced Within the Study Internalisation Measures and Policies for All external Cost of Transport (IMPACT). European Commission DG TREN.
Li, Z., and Hensher, D. A. (2012). Congestion charging and car use: a review of stated preference and opinion studies and market monitoring evidence. Trans. Policy 20, 47–61. doi: 10.1016/j.tranpol.2011.12.004
Martínez, T., and Muro, A. I. (2011). Caso de Estudio3: Corredor Toledo-Madrid, Tercer Informe de Seguimiento, Proyecto AEROAVE: La integración del transporte aéreo con las redes ferroviarias de larga distancia. Recuperado de: http://www.vialibre-ffe.com/pdf/AEROAVE.pdf
McFadden, D. (1981). “Econometric models of probabilistic choice,” in Structural Analysis of Discrete Choice Data with Econometric Applications, eds C. Manski and D. McFadden (Cambridge: The MIT Press), 198–272.
Watkiss, P., Jones, R., Rhodes, D., Hardy, A., Handley, C., and Walker, C. (2001). A Comparative Study of the Environmental Effects of Rail and Shoer-Haul Air Travel. Commission for Integrated Transport.
Keywords: consumer behavior, choice of service, transportation, modeling, discrete choice, logit, willingness to pay
Citation: Muro-Rodríguez AI, Perez-Jiménez IR and Gutiérrez-Broncano S (2017) Consumer Behavior in the Choice of Mode of Transport: A Case Study in the Toledo-Madrid Corridor. Front. Psychol. 8:1011. doi: 10.3389/fpsyg.2017.01011
Received: 25 March 2017; Accepted: 01 June 2017;
Published: 20 June 2017.
Edited by:Alicia Izquierdo-Yusta, University of Burgos, Spain
Reviewed by:Cristina Peñasco, Consejo Superior de Investigaciones Científicas (CSIC), Spain
Victor Martin, Complutense University of Madrid, Spain
Copyright © 2017 Muro-Rodríguez, Perez-Jiménez and Gutiérrez-Broncano. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana I. Muro-Rodríguez, firstname.lastname@example.org