Evaluation of Micro EV’s Spreading to Local Community by Multinomial Logit Model
 Author: Seki Yoichi, Manrique Luis C., Amagai Kenji, Takarada Takayuki
 Organization: Seki Yoichi; Manrique Luis C.; Amagai Kenji; Takarada Takayuki
 Publish: Industrial Engineering and Management Systems Volume 11, Issue2, p148~154, 30 June 2012

ABSTRACT
Micro Electric Vehicles are considered as a solution for reducing CO_{2} emissions, however, it is difficult to evaluate its impact in a local community when it has been introduced. In this study, we evaluated how to spread the Micro EV within the community, using the utility derived from a multinomial logit model, and analyze the effect on CO_{2} emissions. The householder’s utility model is based on an investigation about Kiryu citizen's activities of shopping, transportation methods, etc. Using the geographic information system, we get the distances of each householder and the stores, and estimate a multinomial logit model about the combination choices of shopping stores and transportation method.

KEYWORD
Multinomial Logit Model , Direct Product Alternatives , Utility , CO2 Emission , Geographic Information System , Data Mining

1. INTRODUCTION
Japanese society has two issues for sustainability. One is the rapid aging of population, especially, in local cities where this issue is severe as shown in Figure 1, because young persons tend to leave for larger cities. It has resulted in the decay of urban area of the local small cities, along with the exodus of young residents from the inner cities to the suburbs of the cities during the period of rapid economic growth in the postwar era of Japan. The spread population gives rise to operational difficulties of public transportations in local cities. Under these circumstances, most of the residents in local city use his/her private automobiles for daily short trips. The other issue is the energysaving for antiglobalwarming. To resolve this issue, it is needed to reduce CO_{2} emissions, that is, to save fossil fuel. However, automobiles consume much more fossil fuel than public transportation methods.
One of the solutions is a multimodal urban region (Bertolini and Clercq, 2003). Combined use of public transformation methods for longer trips and personal transportation methods for shorter trips produces a lowcarbon society. Micro electric vehicles (micro EV’s) are considered as one method for personal transformation in this scenario. Micro EV is an EV for one personal use. It has small and light body, low speed, and has characteristics of an automobile and a bicycle. It is easy even for elderly persons to drive, and its energy can be supplied from solar energy (Toshio, 2011). From December 2011, use of micro EV was demonstrated in Kiryu, a local city in Japan, founded by the Minister of Land, Infrastructure, Transportation and Tourism (MLIT, 2010; Hada, 2009).
However, it is difficult to evaluate its impact in a local community when it has actually been introduced. In this study, we evaluated how to spread the micro EV within the community, using the utility derived from the multinomial logit model and analyzed the effect on CO_{2} emissions.
First, we propose a generalization of multinomial logit model (BenAkiva and Lerman, 1995; Train, 2009) for direct product alternatives of shopping stores and transportation methods to the stores. Using a result of questionnaire survey, we estimate this model, and clarify the utility structure for the combination of shopping store and transportation method. This model explains discrete choices of combinations according with other variables such as distance to the store, his/her age, etc.
Next, we assume some conditions of utility of micro EV’s, and estimate the probability of citizens’ choices of the combinations. This new method was added as a new choice of transportation. From this probability, we estimate CO_{2} emissions for shopping trips in Kiryu city.
Our approach in this paper may be used for evaluating how widely a new technology can be spread to a society, and be able to forecast the demand level in the saturated case.
2. MULTINOMIAL LOGIT MODEL FOR DIRECT PRODUCT ALTERNATIVES
2.1 Multinomial Logit Model
The multinomial logit model is one of discrete choice models which explain the probabilities to select alternatives in a finite choose set C. These models assume that alternatives in the choice set are mutually exclusive, exhaustive Based on the utilitymaximizing theory that states that decision maker chooses the alternative which maximizes his/ her utility, the probability that decision maker
i choose an alternativec iswhere
U _{ic} is the utility of the alternativec for the decision makeri .V _{ic} is a representative part ofU _{ic}, which is the function of attributes of the alternativec and the decision makeri. ε _{ic} is a random part ofU _{ic}, which summarizes the contribution of unobserved variables and distributes independently and identically extreme value distribution. We can find the closed form of the probability as the Eq. (2).If we observe
K+K ' attributes from the alternatives and the decision makers, and among them,K attributes are categorical factor attributes andK ’ attributes are continuous variablesx_{k’ic} . We usually formulate the functionV _{ic} of attributes as Eq. (3), for example, in the regression model.where
β _{0} is the intercept constant,β _{kl}’s andβ _{k’l}’s are the coefficients for the attributes of the categorical factor and the continuous variables, respectively.l_{kic} is the level of thek ^{th} factor attribute of choiceic andx _{k’ ic} is the value ofk ’^{th} continuous attribute of choiceic . However, in the multinomial logit model,β _{0} and one of theβ _{kl}’s for eachk must be set 0, to avoid the indetermination of coefficients because of equivalent ratio of Eq. (3) between the choices. In this paper, we setβ _{k1} =0 , i.e., the coefficients of the first categories are set zero.2.2 The Model for Direct Product Alternatives
The choice set of the multinomial logit model is a direct product of two sets, for example, set of shopping stores and the set of transportation methods. Let the two sets be S and T. The choice set is C = {(
s, t ) s ∈S,t ∈ T}. In this example, a decision maker chooses a combination of shores and transportation methodt to the selected store. Utilities are defined for the combinations.Estimation method of the model is the same for the most part to ordinary logit model. However, we need some attention to parameterization of coefficients. Because utility is a latent index, and its values are relative in each decision maker, we must set coefficients of base categories to zero, as mentioned above. So attributes of the two sets are parameterized, as setting their coefficients of first categories to zero, respectively. Furthermore, if we need to estimate any interaction term that includes two attributes, the coefficient of the interaction whose level is the first category, at least in one attribute, must be set zero. An example of this formulation is given in section 4.
3. SURVEY IN KIRYU CITY
The information needed for estimation of the model is based on a survey conducted in Kiryu in October 2009. Kiryu is a local city in Gunma prefecture, and is located about 100 km north of Tokyo.
This study is part of the main project “Construction of the town of Kiryu for the future with antiglobalwarming through the regional power” supported by Japan Science and Technology Agency, and the main goal is to reduce 80% of CO_{2} emitted in the city. We research to understand the effect of modal shift of shopping from the information given on the questionnaires.
This questionnaire was sent to random sampled 10,000 householders from a total of 50,163 household
ers of Kiryu and returned 2,269 questionnaires; it means the 4.52% of the total or the 22.69% of the sample. We divided Kiryu city into 4 areas in order to analyze the relationships between them. These areas are represented as: old town, new town, across river, and detached area in Figure 2.
In this survey, we asked about shopping to the main person who usually goes to shopping in his/her home. The main part of the question is about the frequency per month. We must specify the visited store to obtain the distance from his/her home. However, it is too difficult to ask all of the stores located in Kiryu city. For that reason, we limited alternatives in the questionnaire to the representative shopping places in Kiryu. We have selected 41 representative stores and shopping areas in/near the city, and we have selected some shopping activities done in the store, such as: foods, clothes, daily goods, windowshopping, etc. The frequencies are asked for every store and for every kind of shopping activity. In addition with these frequencies, we asked about the means of transportation used to go to every store. The categories of transportation are bicycle, car, bus, train, and non use (walk).
In addition, we surveyed the demographic information of the respondents; sex, age, location, whether he/she had a car license.
We used a geographic information system (ArcGIS ver. 9.5; Esri, Redlands, CA, USA) to get the location of each householder and the stores from their addresses, and also to get the Euclidean distance between them.
4. ESTIMATION OF SHOPPING TRIP UTILITY
In this paper, we restrict the kind of shopping activity to food shopping, which is the most frequent activity, and estimate a multinomial logit model for alternatives which are direct product of set of shopping stores and the set of transportation methods. The set of stores are 12 shopping places as listed in Table 1. These places are selected from 41 stores as the stores selling foods, and all branches having the same name are summarized to one store.
Also, in the original data we have distinguished transportation methods: car, bus, train, and bicycle. Nevertheless, in order to reduce the number of categories that are few checked, we have created a new transportation category called public transportation, it takes into account the information about bus and train. We consider that the respondent who checked no transformation method walks to the store, and we merge this case to the bicycle. Finally, we distinguish three transportation methods; walk/bicycle, public transportation, and car.
Our dataset consists of 4,871 selected alternatives that are selected by 1,618 householders. It means that each householder selects 3.01 pairs of stores and methods from 12 × 3 = 36 pairs on average. For each selected alternative, numbers of trips in the month are observed; the total of these numbers is 26,092 times.
Using the above data set, we developed different kind of models, which explain the representative utility by 5 attributes and the interactions between them and choices. The used attributes are distance to the store and age, sex, place area, car license status of respondents. Estimated models are evaluated by Akaike information criteria (AIC). From these results, we select the next model.
where
i : Householder.s : Store choice.s ∈{old town shopping area, supermarket 1, …, 9, huge shopping mall 1, 2}t : Transportation method choice.V _{ist} : Utility of visit to stores by the methodt for the individuali .d _{is} : Distance from the homei to the stores .A _{i} : Age of householderi .r _{i} : Place area where home ofi is located.r _{i} ∈{old town, new town, across river, detached area}.O _{s}: Is stores in old town.O _{i} ∈{0, 1}l _{i}: Car license status of householderi .l _{i} ∈{hold, not hold}β _{k} : Coefficient ofk ^{th} attribute .The estimated coefficients for the model are shown in the Table 2.
First, effects related with the kind of transportation method are shown in the Figure 3. This figure indicates the effect to utility of distance from home
i to stores and possession of driving license:β_{1t} + (β_{3} +β_{4t} )log(d_{is} ) +β_{5tl} for eacht andl . For the trip less than 1.37km, for the average person who does not have a car driving license, walk/bicycle, car and public transformation have much utility in this order, and the opposite for the longer trip. However, checking the interaction between transportation method and car license, if he/she has a car license, car’s utility become higher (+1.4), and public transportation’s drop 1.85. Accordingly, car drivers always have maximum utility for car in practical distance, and they don’t tend to use public transportations. As for the effect of age to utility:β_{6t} , elderly persons tend to avoid car driving.Next, we consider the effect of areas. Householders residing in the new town or across river area tend to use some transportation method; public one or car, compared with old town and detached area. The reason is considered that four railways have come to Kiryu city, and they cover across river area and some part of the new city area. Furthermore, bus network covers these areas, though it is not so frequent. Householders in these areas, who need to trip longer distance than old town area, use the some transportation method. On the other hand, in the detached area, access to these transportation methods is inconvenient and the services are less frequent, so they depend on car selfsufficiently. We represent this effect by term
β_{7tr} .With respect to the location of stores, we distinguish whether the store is in the old town. The result of estimated
β_{8t} indicated that public transportation, not car, tends to be selected for the trips to the store in old town.Finally, we confirm the effect of store:
β_{2s} andβ_{9s} . Every store have higher utility than the old shopping street, and have higher utility to younger person, except for oldst and LS, that are located in the old town area.5. ESTIMATION OF CO_{2} EMISSION
5.1 How to Estimate the Spread of New Techno logy
Once we have estimated utilities of alternatives for each householder, using the multinomial logit model, we can estimate any probability that a householder choices a alternative; a combination of a store and a transportation method. Adding some assumption about the utility of new technology, we can estimate the probability that the new technology is selected by each person.
In our case, we estimate the probabilities of micro EV in foodbuying trip. To calculate the probability, we assume below assumptions:
The utility of the micro EV is the mean of car and bicycle in every situation of householder i and store s, as below equation. The utilities of the other methods are invariant under the introduction of the micro EV.
The citizens can buy this transportation method, if he/she likes it.
The citizens don’t change the frequency of foodbuying trip.
The population volume and its demographic situation of the citizen will remain constant.
Once we could calculate the utility value of micro EV, we can estimate the probability of each combination of store and method using Eq. (2). This calculation can be sound, if we could assume an appropriate condition to the utility of the new technology to the original state.
We have already one month frequency of foodbuying for each householder, so we can estimate expected frequency of the alternatives by multiplying the selection probability by the frequency. Each alternative; a trip to a store, has distance of trip. Therefore, we can estimate the one month accumulated trip distance of each householder for each transportation method. The results are listed in Table 3.
5.2 Effect the Spread of Micro EV
First, we see the trip distance per person in Table 3. In surveyed data, car distance is longest among all methods, and at old town it is shortest among the other areas. This feature is the same in the estimated result added micro EV.
After the introduction of micro EV, micro EV became the second method, and it is used in 24.8% of all trip in distance on the average weighted by number of persons. Car remains the first method, however, its use is reduced 20.9%. Micro EV is used in old town most frequently, because of short mean distance of trips.
5.3 Assessment of the Emission of CO_{2}
We assess the rough amount of CO_{2} emissions in Kiryu for foodbuying in a year using Eq. (6).
where,
i : Householder.j : Store.m : Transportation method.d_{ij} : Distance betweeni andj .n_{ij} : Frequency of visits in a month.r_{ijm} : Ratio of visiting methodm .ρ _{m} : Methodsm ’s average CO_{2} emission (g/km).0.04523: Rate of respondents among all citizens.
12: Number of months in a year.
It is difficult to assess precisely Eq. (6), because we do not have precise information about the type of transportation method, for example, fuel consumption of householders’ cars. So, we use the representative value shown in Table 4. In this table the emissions related with bicycle for CO_{2} are 0, because it does not use any fuel. For the public transportation, we use the value of bus as a leading method. As to this value, we assume an average number of passengers equal to 10, and divide the value 403.1 in Table 4 by 10. As to the micro EV, it is also difficult to assess its CO_{2} emission (De Boncourt, 2011), we assume that its emission rate is estimated at 1/20 of car’s one.
Using above assumptions, we assess the amount of CO_{2} emission in Kiryu city at 188.3 (t/year) after micro EV was introduced, against 142.3 (t/year) in the present situation. This is 46.0 (t/year) (24.4%) reductions as indicated in Table 5. Emission per person is smallest in old town, and is largest in detached area. Furthermore, in detached area, where population is sparse and residents need long trips for living, the emission does not reduce so much after the introduction of EV. Therefore, it is important for antiglobalwarming to construct compact city, where almost residents live in a compact area in the city, and move by bicycle or energysaving EV.
6. CONCLUSION
The multinomial logit model is useful, because it helps us to understand the process of choice. This model allows us to calculate the utility using different variables of alternatives’ and decision makers’ attributes. We generalized this model for the combinations of direct product alternatives such as stores and transportation methods.
Using the proposed model, we clarified the structure of CO_{2} emission from foodbuying shopping in a local city, and assessed the impact of introduction of micro EV. We estimated these CO_{2} emissions from personal behavior. This approach gives us only a gross estimate of CO_{2} emission, compared with detailed estimate from consumption levels of different fuels in the city. However, when we think how we can reduce the emission, this approach may be rewarding, because it is informative about why the emission occurs.

[Figure 1.] Population pyramid of Kiryu city.

[Figure 2.] Areas of Kiryu city and sampled householders.

[Table 1.] Shopping places investigated

[Table 2.] Estimated model coefficients

[Figure 3.] Utility function of moving distance of each transformation method.

[Table 3.] Impact of introduction micro electric vehicles (EV) to trip distance

[Table 4.] CO2 emissions rate

[Table 5.] Amount of CO2 emission and reduction