Optimization of an Improved Intermodal Transit Model Equipped with Feeder Bus and Railway Systems Using Metaheuristics Approaches

: One of the serious concerns in network design is creating an efﬁcient and appropriate network capable of efﬁciently migrating the passenger’s mode of transportation from private to public. The main goal of this study is to present an improved model for combining the feeder bus network design system and the railway transit system while minimizing total cost. In this study, the imperialist competitive algorithm (ICA) and the water cycle algorithm (WCA) were employed to optimize feeder bus and railway services. The case study and input data were based on a real transit network in Petaling Jaya, Kuala Lumpur, Malaysia. Numerical results for the proposed model, including the optimal solution, statistical optimization results and the convergence rate, as well as comparisons are discussed in detail.


Introduction
Transportation is a multimodal, multi-problem and multi-spectral system, as it involves different categories and activities, such as policy-making, planning, designing, infrastructure construction and development.Currently, considering the significant developments in technology, economy and society, an efficient transportation system plays a key role in passengers' satisfaction and the reduction of costs.
Many people use public transportation systems to reach their destination; however, others employ personal vehicles.Passengers are more likely to use a transit service that is highly reliable.Travelers may switch to other transportation modes if a transit service does not provide the expected levels of service [1].To prevent the increase of private transport entering city centers, effective alternative travel modes must be provided [2].Nuzzolo and Comi [3] investigated methods of transit network modelling that can be implemented in transit decision-making support system tools to improve their performances, according to the latest innovations in information technology and telematics.
In addition, a good public transportation system has been recognized as a potential means of reducing air pollution, decreasing energy consumption, increasing mobility and improving traffic congestion.In order to improve complicated public transportation systems, a well-integrated transit system in urban areas can play a crucial role in passengers' satisfaction and reduce operating costs.This system usually consists of integrated rail lines and a number of feeder routes connecting transfer stations.
In general, previous approaches to tackle transit network problems can be divided into two major groups: analytic and network approaches.These approaches differ in their purposes and have different advantages and disadvantages.They should be considered as complementary rather than alternative [4].Numerous studies have attempted to implement analytic models [5][6][7][8][9][10].On the other hand, some of the researchers adopted network approaches instead of analytical methods [2,[11][12][13][14].Many studies have been carried out to identify solutions using the aforementioned approaches.These can be categorized into four groups, namely mathematical, heuristic, metaheuristic and hybrid techniques [15].Kuah and Perl [6] presented a mathematical method for designing an optimal feeder bus network to access an existing rail line.Furthermore, Chang and Hsu [10] developed a mathematical model to analyze the passenger waiting time in an intermodal station in which the intercity transit system was served by feeder buses.They presented the analytic model for quantifying the relationships of passenger waiting time to the reliability of feeder bus services and the capacity of intercity transit.
A large number of research papers has been published in recent years utilizing heuristic methods due to their flexibility.Shrivastav and Dhingra [12] developed a heuristic algorithm to integrate suburban stations and bus services, along with the optimization of coordinated schedules for feeder bus services using existing schedules for suburban trains.Sumalee et al. [16] proposed a stochastic network model for a multimodal transport network that considers auto, bus, underground and walking modes.Chowdhury [9] proposed a model for better coordination of the intermodal transit system.Furthermore, Steven and Chien [17] suggested the use of a specific feeder bus service to provide shuttle service between a recreation center and a major public transportation facility.They suggested an integrated methodology (i.e., analytical and numerical techniques) for the development and optimization of decision variables, including bus headway, vehicle size and route choice.
In terms of metaheuristic methods, Kuan [13] applied genetic algorithms (GAs), ant colony optimization (ACO), simulated annealing (SA) and tabu search (TS) to resolve a feeder network design problem (FNDP) for a similar work conducted by Kuah and Perl [11], which improved previously-proposed solutions.
In another study carried out by Shrivastava and O'Mahony [18], optimum feeder routes and schedules for a suburban area were determined using GAs.Mohaymany and Gholami [19] suggested an approach for solving multi-modal FNDP (MFNDP), whose objective was to minimize the total operator, user and society costs.They used ACO for constructing routes and modifying the optimization procedure in order to identify the best mode and route in the service area.
Hybrid methods are categorized as another type of solution method that combines the abilities of different computational techniques to solve complex problems.Shrivastava and O'Mahony [20] developed the Shrivastava-O'Mahony hybrid feeder route generation algorithm (SOHFRGA).The idea was to develop public bus routes and to coordinate schedules in a suburban area.
Numerous researchers have attempted to design a more efficient feeder network and to provide feeder services connecting major transportation systems and welfare facilities.The main target of this paper is to represent an improved model and to present an efficient transit system to increase the efficiency of feeder network designs in order to minimize costs.
The structure of this paper is organized as follows: Sections 2 and 3 present a brief description and definition of the problem and explain the details of the mathematical model, respectively.Section 4 provides succinct representations of the applied ICA and WCA in order to optimize the transit system problem considered.The computational optimization results obtained from the methods used along with a discussion and comparisons are presented in Section 5. Finally, the concluding results and suggestions for future research are given in Section 6.

Problem Definition and Assumptions
In large metropolitan areas, particularly those with high transit demand, an integrated transit system consisting of rail lines and a number of feeder routes connected at different transfer stations is essential.Consequently, designing a proper feeder network that can provide access to an existing rail system and coordinate the schedule of transit service can be a significant issue.
The development of improved integrated intermodal systems can result in a higher quality of service and passenger satisfaction by providing better coverage, reduced access time, minimal delay and shorter travel times.From the transit operators' point of view, the operating costs may be reduced by an overall coordination between different public transport modes.Profit can also be increased by shorter route maintenance and eliminating the duplication of routes by trains and buses.
This study is focused on designing a set of feeder bus routes and determining operating frequency on each route, such that the objective function of the sum of the operator, user and social cost is minimized.The mathematical formulation of the improved model and the details of the constraints are presented in the following sections.
An intermodal transit network consisting of a rail line and feeder bus routes connecting the transfer stations is assumed to serve the examined area.The optimal transit system will be determined based on an assumed route structure (i.e., one rail line and feeder bus routes are linked with straight lines between nodes) and the peak hour demand situations in the entire service area.To formulate the mathematical model for an intermodal transit system and its application in the case study, the following assumptions are made: (1) The transit network was designed with feeder buses and a fixed rail line.
(2) Transit demand is assumed to be independent of the quality of transit service (i.e., fixed demand).
The demand pattern for feeder bus routes is many-to-one.(3) The location of nodes (i.e., bus stops and rail stations) is given.Some of the model parameters (e.g., vehicle sizes, operating speed, cost) are specified.(4) All feeder routes can be used in both directions for the transit service.

Model Formulation
To propose a mathematical formulation for the model based on the problem statement of this study, the total cost function is expressed in Equation (1).The total cost function is the sum of the user, operator and social costs, which can be formulated as follows: where C T , C u , C o and C s represent the total cost, user costs, operation costs and social costs, respectively.The well-structured cost classification for the proposed model is shown in Figure 1.Table 1 tabulates each cost mentioned more comprehensively later in this section.
For nomenclature and convenience purposes, all variables and parameters used for the modified objective function are defined in Table 2.The user cost is the expense imposed on passengers using the transit system (contains feeder and train services).This cost is comprised of access, waiting and in-vehicle traveling costs, denoted by C a , C w and C ui , respectively, in the following equation: In light of the user cost, which is the summation of feeder bus and train cost, Equation ( 2) can be re-written as follows: C u " pC aF `CaT q `pC wF `CwT q `pC uiF `CuiT q (3) Generally, all elements of the user cost can be formulated as the product of an hourly demand, average time spent in each travel time category (i.e., access time, wait time and in-vehicle time) and the users' value of time, which are explained in the following subsections.

Access Costs (C a )
Feeder and train passengers who have access to stops and stations mainly incur the access cost.The access cost is generally experienced by local and train passengers accessing the transfer station.
The access cost for feeder bus passengers is the product of local demand, q i , with average access time t aF and the value of time µ a , where t aF can be estimated from the average distance divided by the average access speed.The average access time for train passengers (t aT ) can be formulated similarly.t aT is dependent on the distance between the platforms of bus and train services and access speed.Assume that access speed and the value of time for feeder bus and train passengers are identical.Thus, the access cost for feeder route k can be formulated as follows: C a " µ a pQ k ˆtaF `Qk ˆtaT q (4) The users' value of time (µ a ) is an important parameter in determining the user cost and is usually dependent on the economic situation (e.g., annual income).

Waiting Cost (C w )
The waiting time includes passengers waiting for the buses and trains.Additionally, it starts counting when a passenger arrives at the bus stop or rail station and stops when the person boards the vehicle [21].The waiting cost is the product of average wait time, demand and the value of users' wait time (µ w ).Average wait time can be estimated by a fraction of the headway.In this model, the average wait times for the feeder bus at the stops and for trains at the stations are assumed to be one half of the headway.Hence, the user waiting cost can be represented using Equation (5) as follows: Similarly, the product of demand, in-vehicle time and the value of time can define the user in-vehicle cost (C ui ).The C ui is formulated based on the average journey time and is calculated in two main parts: the run time and the dwell time.Running costs for all passengers (C rui ) are equal to the link travel distance from stop i to station j in route k (L ijk ) divided by the average bus real speed (V k ).
The dwell time is the boarding and alighting time at the feeder bus stops (t dF ) and rail stations (t dT ).The observation of feeder bus stops and rail stations revealed that the dwell time is an important part of in-vehicle travel time.This time will increase the user, operation and social costs for both feeder bus and train travel and consequently has a significant effect on the total cost of the transit network.Dwell time will increase user costs by increasing the in-vehicle time for a boarding passenger.In addition, this time cost will increase operation costs by increasing fuel consumption, maintenance and personnel costs.Accordingly, with the increase in pollution, noise, greenhouse gases, etc., the social costs will also be higher.
Since the time spent on boarding and alighting has an important role in user in-vehicle time, we tried to present a new concept for determining such costs.Moreover, because of the variation in the time spent on boarding and alighting, which is dependent on the dwell time at each of the bus stops, the geometric series equation was adopted to develop a more accurate model for distributing the dwell cost of the bus stops along the routes.
The average cost of dwell time is determined by demand multiplied by the rate of passenger boarding and alighting.The derivation of the dwell cost for feeder buses and trains is discussed in detail in Appendix A. Therefore, the in-vehicle cost, including in-bus and in-train cost, for route k, is given as follows: where C rui is the running cost for all passengers, given as follows: and C dui is the average cost of dwell time as described in Appendix A, given in the following equation: The first and second terms in Equations ( 7) and ( 8), respectively, denote the feeder bus and train user costs.In Equation ( 6), C rui represents the running cost for all passengers, which is equal to the link travel distance from stop i to station j in route k (L ijk ) divided by the average bus speed on route k (V k ).t Tj denotes the riding time between station j and the destination of the train regardless of boarding and alighting times.

Operating Cost (C o )
The operating cost (C o ) is the summation of railway and feeder bus operation costs.It can be described by the unit time or distance cost (in hours or km) in connection with the transit service provided.Thus, C o can be formulated as the sum of C oi , C m , C p and C f .These costs include the cost of trains and buses; therefore, it can be formulated as follows: C o " pC oiF `CoiT q `pC mF `CmT q `pC pF `CpT q `pC fF `CfT q (9)

Feeder Bus Maintenance Cost (C mF )
The feeder bus maintenance cost (C mF ) consists of maintenance, repair and tire costs.This cost depends on the fleet size and round trip distance, formulated as follows: The feeder bus personnel cost (C pF ), including the drivers and administrative costs, is dependent on the fleet size, hourly pay and insurance rate.Since the time for boarding and alighting (i.e., dwell time), as well as bus slack time have important roles in spending time for personnel, in this study, an effort has been made to represent the improved concept for the determination of these costs.Hence, in order to increase the accuracy of the cost function (objective function), adding slack time (S kj ) into the schedule of bus route k at station j and average rest time were considered for each bus at stations.Moreover, on mathematical formulation, the dwell times were added into the calculation of personnel costs with respect to the interrelationship among cost terms.The derivation of the C pF is given in Appendix A. Therefore, C pF for feeder bus route k can be formulated as follows: The first and second terms in Equation ( 11) rely on the feeder running time and the dwell time in route k, respectively.Accordingly, the third term denotes the personnel cost when drivers are in the rest time or queue.

Feeder Bus Fixed Costs (C fF )
The feeder bus fixed cost contains initial fleet costs, such as vehicle ownership costs, license, insurance and so forth.It is formulated according to the fleet size and hourly fixed cost for the vehicle given for route k as follows: The feeder bus operating in-vehicle cost (C oiF ) is dependent on the travel time and round trip distance.C oi (for bus or train) is formulated based on the running cost (C roi ) and the dwell cost (C doi ).The running cost for the bus is formulated according to the round trip distance against the rail, which is the round trip time.It is assumed that the stop delay time incurred at bus stops and intersections should be taken into consideration.
As explained in Section 3.1.3., the average cost of dwell time was defined by demand multiplied by the passenger boarding and alighting rate.Furthermore, these costs were also determined similarly in Section 3.1.3.Thus, the C oiF for feeder bus route k can be formulated as given in the following equations: C roiF " λ l p2F k ˆLK q (14) The derivation of C doiF is discussed later in Appendix A.

Train Operating Cost (C oT )
The operating cost for a rail system can be obtained through multiplying the fleet size by the value of the train operating cost (λ IT ).The fleet size can be obtained from the trip time multiplied by the train frequency (F T ), and the rail trip time consists of running and dwell time.The train running time is trip distance divided by average running speed (V T ).In addition, the rail dwell time is the product of the number of inflow or outflow passengers on the route and the average service time for passengers boarding and alighting from a vehicle.Thus, the train operation cost can be formulated as given in the following equation: where the first term in Equation ( 16) corresponds to the train running time and the second term denotes the train dwell time.As a fixed rail line is assumed and operation cost depends on route station distance and demand, one operating value for all operating costs is considered in order to simplify the model in this study.λ IT represents all elements of operating cost, including fixed, maintenance, personnel and in-vehicle costs ($/veh-h).The derivation of this cost is represented in Appendix A.

Social Costs (C s )
Social costs consist of many parameters that non-users pay indirectly.For instance, accident costs, pollution costs, infrastructure costs, noise, greenhouse gases, and so forth.This cost is assumed to be dependent on in-vehicle operating costs for feeder services and formulated as follows: Each cost term consists of several parameters and items, which, consequently, have different effects on total cost.Furthermore, determining some of the parameters and items requires the cooperation of other organizations.Therefore, the interrelationship between some of the cost terms and other related costs is considered.It is assumed that there is an interrelationship among some of the costs, such as "social cost" and "feeder operating in-vehicle cost".Thus, based on previous studies [19,22,23], in order to simplify the proposed model, the social cost is assumed to be 20% of the "feeder operating in-vehicle cost" in this study.

Total Cost for a Route (C Tk )
After calculating all cost components for route k, the total cost function C Tk for route k is expressed as given in the following equation: ff λl p2F k ˆLK q `λIT rpQ k ˆtdT q `pF T ˆTT qs `λI pQ k ˆtdF q `λm p2F k ˆLK q λp "´2

Objective Function and Constraints of the Model
The total system cost of the intermodal transit model consists of user parameters (i.e., the value of the time for user's access, wait and in-vehicle cost, etc.), operation parameters, social parameters and the number of decision variables (i.e., number of routes).
This transit network model must satisfy users, operators and social terms.Thus, the objective function is defined as the sum of the user, operator and social costs, which is given in the following equation: Therefore, the objective function can be formulated after substitution of all cost terms as given follows: subject to: where decision variables contain two binary variables, called Y ij and X ihk , which stand for the transit network, and a continuous variable for the feeder bus frequency (F k ).The determination of F k , as one of the decision variables, depends on the transit network configuration.Thus, the optimal feeder bus frequency using the analytical solution can be determined by setting the first derivative of the total cost function (C TK ) with respect to the feeder bus frequency, equating it to zero and solving it.Therefore, the optimal bus frequency can be taken as: Furthermore, the minimum required frequency for route k is given as follows: Thus, the given frequency for route k is obtained by selecting the maximum value for the optimum frequency (F opt,K ) and required frequency (F req,K ).Some limitations are considered for the proposed improved model to represent an effective transit network model satisfying route feasibility, frequency, and so forth.Equations ( 21)-( 25) correspond to the route feasibility in the network design.Several researchers used these constraints in their studies [2,11,13].Equation (21) explains that each bus stop should be placed in a single route (many-to-one pattern).Furthermore, Equation (22) ensures that each generated route must be connected to only one railway station.Accordingly, in Equation ( 23), each bus is assumed to pass all of the stops in its route node.Equation (24) describes how each feeder bus route should be linked to only one railway station.The constraint given in Equation (25) specifies that a bus stop can be assigned to a station in which the corresponding route terminates at one of the rail stations.Constraints on the minimum and maximum length of feeder routes are given in Equation (26).Similarly, limitations for the minimum and maximum frequencies are indicated in Equation (27).Equation (28) shows the allowable maximum number of vehicles in the fleet, and Equation ( 29) represents the restriction for the minimum frequency in order to satisfy the demand.

Applied Optimization Methods
The transit network design problems are categorized as NP-hard problems with a nonlinear objective function and constraints.Searching for the best feasible routes in order to minimize the cost function is crucial in solving the feeder network design and scheduling problem (FNDSP).Therefore, optimization approaches, which are mostly metaheuristics, are of great importance.There are many methods being used to solve transit network design problems.Based on the literature, there are pros and cons for all of these optimization methods [22].
WCA and ICA have shown great potential for solving optimization problems, as they have been used for global stochastic searches [24][25][26].These two metaheuristic algorithms were employed to optimize the model for the case study considered in this paper.Brief explanations of each optimizer are provided in the following subsections.

Imperialist Competitive Algorithm
The imperialist competitive algorithm (ICA) is inspired by the social-political process of imperialism and imperialistic competition.Similar to many optimization algorithms, the ICA starts with an initial population.Each individual in the population is called a 'country'.Some of the best countries with minimal cost are considered imperialist states, and the rest are colonies of those imperialist states.All of the colonies are distributed among the imperialist countries based on their power.
To define the algorithm, first, initial countries of size N Country are produced.Then, some of the best countries (with the size of N imp ) in the population are selected to be the imperialist states.Therefore, the rest of the countries with size N col will form the colonies that belong to the imperialists.Then, the colonies are divided among the imperialists [24] in such a way that the initial number of each empire's colonies has to be proportional to its power.Hence, the initial number of colonies for the n-th empire will be [27]: where NC n is the initial number of colonies for the n-th empire and N col is the total number of initial colonies.To divide the colonies, NC n of the colonies are randomly chosen and assigned to the n-th imperialist [27].After dividing all colonies among the imperialists and creating the initial empires, these colonies start moving toward their relevant imperialist country.This movement is a simple model for assimilation policy.Furthermore, the total power of an empire is defined by the sum of the cost of the imperialist, and some percentage of the mean cost of its colonies, as given below [27]: where TC n is the total power of the n-th empire and ξ is a positive small number.After computing the total power of empires, the weakest colony (or colonies) of the weakest empire is usually targeted by other empires and competition begins on possessing this colony.Each imperialist participating in this competition, according to its power, has a probable chance of possessing the cited colony.
To start the competition, at first, the weakest empire is chosen, and then, the possession probability of each empire is estimated.The possession probability P p is related to the total power of the empire (TC).In order to evaluate the normalized total cost of an empire (NTC), the following equation is used [27]: During the imperialistic competition, the weak empires will slowly lose their power and become weaker over time.At the end of the process, only one empire will remain that governs all colonies [27].

Water Cycle Algorithm
The water cycle algorithm (WCA) is inspired by nature and is based on the observation of the water cycle and how rivers and streams flow downhill towards the sea in the real world.Similar to other metaheuristic algorithms, the WCA begins with an initial population called the population of streams.First, we assume that we have rain or precipitation.The best individual (i.e., best stream) is chosen as a sea [25].
Then, a number of good streams (N sr ) are chosen as rivers.Depending on their magnitude of flow (i.e., cost/fitness function), rivers and the sea absorb water from streams.Indeed, streams flow into rivers and rivers flow to the sea.Furthermore, it is possible that some streams directly flow to the sea.Therefore, new positions for streams and rivers may be given as follows [26]: where rand is a uniformly-distributed random number between 0 and 1 (1 < C < 2).If the solution given by a stream is better than its connecting river, the positions of the river and stream are exchanged.Such an exchange can similarly occur for rivers and the sea, as well as for the sea and streams.
For the exploration phase, if the normal distances among rivers, streams and the sea are smaller than a predefined value (d max ), new streams are generated flowing to the rivers and sea (i.e., evaporation condition).A schematic view of the WCA is illustrated in Figure 2, where circles, stars and a diamond correspond to the streams, rivers and the sea, respectively [25].Detailed comparisons concerning similarities and differences between the PSO and WCA, as well as other optimizers have been discussed in the literature [28].by other empires and competition begins on possessing this colony.Each imperialist participating in this competition, according to its power, has a probable chance of possessing the cited colony.
To start the competition, at first, the weakest empire is chosen, and then, the possession probability of each empire is estimated.The possession probability Pp is related to the total power of the empire (TC).In order to evaluate the normalized total cost of an empire (NTC), the following equation is used [27]: During the imperialistic competition, the weak empires will slowly lose their power and become weaker over time.At the end of the process, only one empire will remain that governs all colonies [27].

Water Cycle Algorithm
The water cycle algorithm (WCA) is inspired by nature and is based on the observation of the water cycle and how rivers and streams flow downhill towards the sea in the real world.Similar to other metaheuristic algorithms, the WCA begins with an initial population called the population of streams.First, we assume that we have rain or precipitation.The best individual (i.e., best stream) is chosen as a sea [25].
Then, a number of good streams (Nsr) are chosen as rivers.Depending on their magnitude of flow (i.e., cost/fitness function), rivers and the sea absorb water from streams.Indeed, streams flow into rivers and rivers flow to the sea.Furthermore, it is possible that some streams directly flow to the sea.Therefore, new positions for streams and rivers may be given as follows [26]: ) ) where rand is a uniformly-distributed random number between 0 and 1 (1 < C < 2).If the solution given by a stream is better than its connecting river, the positions of the river and stream are exchanged.Such an exchange can similarly occur for rivers and the sea, as well as for the sea and streams.
For the exploration phase, if the normal distances among rivers, streams and the sea are smaller than a predefined value (dmax), new streams are generated flowing to the rivers and sea (i.e., evaporation condition).A schematic view of the WCA is illustrated in Figure 2, where circles, stars and a diamond correspond to the streams, rivers and the sea, respectively [25].Detailed comparisons concerning similarities and differences between the PSO and WCA, as well as other optimizers have been discussed in the literature [28].

Description of the Study Area
The mathematical formulation and optimization algorithms are applied to a real study area in Petaling Jaya (PJ), Kuala Lumpur, Malaysia.PJ is a major Malaysian city originally developed as a satellite township in Kuala Lumpur.The objective was to minimize the total cost of the feeder bus network in the area.The case study region, shown in Figure 3, is an area of 5.5 km by 6.5 km in the south of PJ in Malaysia and includes the Kelana Jaya Line of the Kuala Lumpur LRT.There are four stations in the study region.The existing bus stops in each traffic zone are considered feeder bus stops covered by a feeder line.A certain amount of demand corresponding to the traffic zone is assigned to each of the bus stops in that particular zone.

Description of the Study Area
The mathematical formulation and optimization algorithms are applied to a real study area in Petaling Jaya (PJ), Kuala Lumpur, Malaysia.PJ is a major Malaysian city originally developed as a satellite township in Kuala Lumpur.The objective was to minimize the total cost of the feeder bus network in the area.The case study region, shown in Figure 3, is an area of 5.5 km by 6.5 km in the south of PJ in Malaysia and includes the Kelana Jaya Line of the Kuala Lumpur LRT.There are four stations in the study region.The existing bus stops in each traffic zone are considered feeder bus stops covered by a feeder line.A certain amount of demand corresponding to the traffic zone is assigned to each of the bus stops in that particular zone.

Data Collection
In order to execute the transit network problem, four important datasets must be available, namely the list of all nodes (i.e., bus stops and rail stations), the network available connectivity list, the transit demand matrix and the cost parameters.A total of 54 nodes is defined to describe the service area and associated network connectivity.The list of locations associated with these 54 nodes is tabulated in Tables 3 and 4. Network connectivity is generated from street links that connect these 54 nodes, and these are suitable for bus operations.The generation of the demand matrix is based on a questionnaire survey data collection.The demand matrix was determined by extracting the abstained results from the survey.
The questionnaire is designed to collect a respondent's origin and destination.Targeted respondents are LRT passengers that queue up at bus stops and LRT stations in different locations in the study area.Generally, larger sample sizes provide more accurate survey results.Nonetheless, due to the constraints of limited resources and time, the sample size for transit service is confined to 20 percent of passengers counted in each LRT station using public buses.
The random sampling technique is employed in this survey to make sure that each member of the population has an equal chance of being selected as a respondent.The locations of LRT stations in the PJ area used for conducting the questionnaires, which widely covers the study area, are shown in Table 5.
The data for normal weekdays were applied in the research.The survey time slot was for the three-hour morning peak period from 6:30 a.m. to 9:30 a.m. and was designed to capture the feeder bus passengers of morning peak times.To determine some of the data and the value of parameters in the proposed mathematical model, the observation method was used in this research.In this section, the design and procedures employed for conducting the observation in train stations and bus stops are presented.Observations are conducted based on the LRT stations and existing feeder bus routes.Observations are categorized into two types of questions.The first question was about the time spent boarding and alighting at the LRT station and bus stops based on time per passenger.The second question was about the average feeder bus speed in existing routes.The cost parameters are based on data collection for the current study, as well as ridership and financial reports publicized by Barton and Valley Metro [29,30].

Results and Discussions
The presented model, explained in detail in Section 3, was applied to the transit services, including bus feeder services connecting the rail stations in the case study (i.e., PJ area).The locations of nodes (i.e., bus stops and rail stations) are given in Tables 3 and 4. Furthermore, the demand of each bus stop is listed in Table 6.The values for the model parameters (e.g., vehicle sizes, operating speed and costs) are specified in Table 7.The WCA and ICA techniques have demonstrated their viability as powerful optimization tools with great potential for solving optimization problems [31][32][33][34].The proposed model and the corresponding optimization methods were coded and run in MATLAB programming software.The optimization procedure for the transit service model involved 50 independent runs, which were performed for each of the considered optimizers.
After performing sensitivity analyses for both optimizers with 50 independent runs, initial parameters for the WCA were a population size of 100, N sr of eight and d max of 1 ˆ10 ´5.Accordingly, for the ICA, the initial parameters consisted of a country population of 100, a number of imperialist country of eight and a revolution rate of 0.4.
The application of different optimization algorithms resulted in solutions with diverse precision values.In fact, the solutions state the accuracy of the applied methods and the method's ability to determine the optimum results.There is a close relationship between the number of function evaluations (NFEs) and the best solutions obtained.
This means that the ideal situation contains the least number of NFEs and is a more accurate solution.With regard to the convergence trend of optimization algorithms and in order to draw a fair comparison between the optimizers, a maximum of 100,000 NFEs was considered the stopping condition for both optimizers.
The number of generated routes is also considered a design variable and varies in a generated population.In fact, the total number of design variables can be changed in this problem in each iteration.Having numerous and changeable design variables can be considered a special feature that can categorize this model as a dynamic optimization model.
The optimization results for the presented optimization algorithms are compared and discussed in this section.Table 8 shows the comparison of the best solutions attained for all cost terms using the applied optimization engines for the improved model.The obtained total cost (C T ) is highlighted in bold in Table 8 for the two reported algorithms.The best solution obtained in the PJ area is provided by the WCA, as shown in Table 8.Additionally, the main costs are illustrated graphically in Figure 4.A total cost of RM 23,494.8 per hour is achieved, with an average service frequency of 6.2 trips per hour (buses on average arriving at intervals of 9.67 min).Accordingly, Table 9 demonstrates the comparison of the statistical optimization results obtained for the two reported optimizers for the FNDSP.The best solution obtained in the PJ area is provided by the WCA, as shown in Table 8.Additionally, the main costs are illustrated graphically in Figure 4.A total cost of RM 23,494.8 per hour is achieved, with an average service frequency of 6.2 trips per hour (buses on average arriving at intervals of 9.67 min).Accordingly, Table 9 demonstrates the comparison of the statistical optimization results obtained for the two reported optimizers for the FNDSP.Table 9 shows that the WCA obtained the best cost (minimum cost) for the FNDSP.The WCA performed better compared to the ICA, having better solution stability.The detailed statistical optimization results associated with each term in the modified cost function using the applied algorithms are presented in Table 10.
As shown in Table 10, it can be concluded that the WCA is superior over the ICA optimizer for finding all cost terms (except the CP) with minimum statistical optimization results.Figure 5 shows the deviation percentage associated with the corresponding optimization algorithms for 50 independent runs.Acceptable stability can be seen in the results among different runs.
In fact, Figure 5 confirms the reliability of the presented optimization methods.The stochastic nature of the method applied to produce the initial values for each iteration makes it natural not to obtain the same results through different independent runs.However, the final results (CT) are similar for various runs, with about 0.92% and 1% for the WCA and ICA, respectively.The ICA is only about 0.08% worse than the WCA.Both methods are comparable.The convergence rate and cost history (i.e., cost reduction) of the applied optimization algorithms have been compared and are illustrated in Figure 6.Table 9 shows that the WCA obtained the best cost (minimum cost) for the FNDSP.The WCA performed better compared to the ICA, having better solution stability.The detailed statistical optimization results associated with each term in the modified cost function using the applied algorithms are presented in Table 10.
As shown in Table 10, it can be concluded that the WCA is superior over the ICA optimizer for finding all cost terms (except the C P ) with minimum statistical optimization results.Figure 5 shows the deviation percentage associated with the corresponding optimization algorithms for 50 independent runs.Acceptable stability can be seen in the results among different runs.
In fact, Figure 5 confirms the reliability of the presented optimization methods.The stochastic nature of the method applied to produce the initial values for each iteration makes it natural not to obtain the same results through different independent runs.However, the final results (C T ) are similar for various runs, with about 0.92% and 1% for the WCA and ICA, respectively.The ICA is only about 0.08% worse than the WCA.Both methods are comparable.The convergence rate and cost history (i.e., cost reduction) of the applied optimization algorithms have been compared and are illustrated in Figure 6.  1 "SD" stands for standard de viation, and the value s are in RM.     1 "SD" stands for standard de viation, and the value s are in RM.Considering the trend of convergence for each method, the WCA is capable of determining faster optimum solutions with a higher level of precision in comparison with the ICA, as can be seen in Figure 6.It can be observed that the convergence rate for the WCA is faster than the ICA at earlier Considering the trend of convergence for each method, the WCA is capable of determining faster optimum solutions with a higher level of precision in comparison with the ICA, as can be seen in Figure 6.It can be observed that the convergence rate for the WCA is faster than the ICA at earlier iterations.Figure 7 illustrates the location and variation of the total cost (CT) and its main components, namely Cu, Co and Cs, or 50 independent runs of each optimization algorithm.
It can be observed that the lowest levels for the average cost terms are 24,023.97,16,917.10,6993.46 and 113.45, respectively, for CT, Cu, Co and Cs with the WCA.The differences between the cost terms of both algorithms for the average cost values are 121, 110.6, 8.9 and 1.5, respectively, for CT, Cu, Co and Cs.It can be highlighted that the lowest level of the average cost in terms of both algorithms becomes nearly the same.However, the ICA shows minimal variation in levels between the first and third quartiles compared to the WCA.In terms of CT for the ICA, the second and third quartile boxes are approximately the same size.The box plot for that dataset would look like one for a normal distribution, however, with a number of outliers beyond one whisker.Table 11 provides the best solution obtained by the WCA among all of the runs, which is illustrated in Figures 8-10.It can be highlighted that the lowest level of the average cost in terms of both algorithms becomes nearly the same.However, the ICA shows minimal variation in levels between the first and third quartiles compared to the WCA.In terms of C T for the ICA, the second and third quartile boxes are approximately the same size.The box plot for that dataset would look like one for a normal distribution, however, with a number of outliers beyond one whisker.Table 11 provides the best solution obtained by the WCA among all of the runs, which is illustrated in Figures 8-10.As shown in Table 11, the transit network consists of 17 feeder bus routes with an average service frequency of 6.24 trips per hour (the average headway is 9.61 min).The case solution includes 3062.4 passenger-km of travel.The provided total cost consists of 71.05% user, 28.53% operator and 0.42% social costs, which is illustrated graphically in Figure 8.It can be observed that the main costs are ranked as user, operator and social costs, respectively, from the maximum to minimum.As shown in Table 11, the transit network consists of 17 feeder bus routes with an average service frequency of 6.24 trips per hour (the average headway is 9.61 min).The case solution includes 3062.4 passenger-km of travel.The provided total cost consists of 71.05% user, 28.53% operator and 0.42% social costs, which is illustrated graphically in Figure 8.It can be observed that the main costs are ranked as user, operator and social costs, respectively, from the maximum to minimum.Furthermore, all costs based on RM for each of the cost terms are shown in Figure 9.It is obvious from Figure 9 that the maximum cost of RM 16,692.6 belongs to CU, while Cs has the minimum cost of RM 98.8.Furthermore, the best transit network obtained by the WCA is given graphically in Figure 10.Furthermore, all costs based on RM for each of the cost terms are shown in Figure 9.It is obvious from Figure 9 that the maximum cost of RM 16,692.6 belongs to C U , while C s has the minimum cost of RM 98.8.Furthermore, the best transit network obtained by the WCA is given graphically in Figure 10.As shown in Table 11, the transit network consists of 17 feeder bus routes with an average service frequency of 6.24 trips per hour (the average headway is 9.61 min).The case solution includes 3062.4 passenger-km of travel.The provided total cost consists of 71.05% user, 28.53% operator and 0.42% social costs, which is illustrated graphically in Figure 8.It can be observed that the main costs are ranked as user, operator and social costs, respectively, from the maximum to minimum.Furthermore, all costs based on RM for each of the cost terms are shown in Figure 9.It is obvious from Figure 9 that the maximum cost of RM 16,692.6 belongs to CU, while Cs has the minimum cost of RM 98.8.Furthermore, the best transit network obtained by the WCA is given graphically in Figure 10.The application and optimization of the proposed model to the PJ transit network provided more accurate and efficient solutions for various conditions in the transit systems by employing additional terms and constraints in the objective function.In other words, an effort was made to widen the scope of the research by considering all aspects of satisfaction (i.e., user cost satisfaction, operation cost satisfaction and social cost satisfaction).
Taking into consideration the proposed objective function and imposed constraints, which have already been explained in Section 3, deriving these levels of cost terms shows that the proposed model can be considered a potentially feasible model to overcome current difficulties in the public transit system.This model may lead to the creation of a more realistic model for simulating real-life problems by providing fresh empirical data for future works.The application and optimization of the proposed model to the PJ transit network provided more accurate and efficient solutions for various conditions in the transit systems by employing additional terms and constraints in the objective function.In other words, an effort was made to widen the scope of the research by considering all aspects of satisfaction (i.e., user cost satisfaction, operation cost satisfaction and social cost satisfaction).
Taking into consideration the proposed objective function and imposed constraints, which have already been explained in Section 3, deriving these levels of cost terms shows that the proposed model can be considered a potentially feasible model to overcome current difficulties in the public transit system.This model may lead to the creation of a more realistic model for simulating real-life problems by providing fresh empirical data for future works.
The results of the application and optimization of the transit network design problem in the study area provided more accurate and efficient solutions for various conditions in transit systems.The outputs of these solutions have demonstrated that the presented model has been verified, and the applied optimizers (i.e., WCA and ICA) are considered suitable for obtaining moderate quality solutions under certain computational cost.This confirms the reliability of the presented methodology.Therefore, this model could be considered an alternative model for real transit networks.

Conclusions
In this paper, an improved model was suggested for transit network problems, including rail system and feeder bus network designs, as well as frequency setting problems.The main purpose of this paper was to develop a real-life model (actualizing the cost function and adding additional constraints) for handling the feeder bus design and frequency setting problems.
The case study for the research was based on the actual transit network in the Petaling Jaya area, in Kuala Lumpur, Malaysia.Finding the optimum feasible routes in order to reduce the cost function is a vital and difficult task for solving the transit network design problem, which is classified as an The results of the application and optimization of the transit network design problem in the study area provided more accurate and efficient solutions for various conditions in transit systems.The outputs of these solutions have demonstrated that the presented model has been verified, and the applied optimizers (i.e., WCA and ICA) are considered suitable for obtaining moderate quality solutions under certain computational cost.This confirms the reliability of the presented methodology.Therefore, this model could be considered an alternative model for real transit networks.

Conclusions
In this paper, an improved model was suggested for transit network problems, including rail system and feeder bus network designs, as well as frequency setting problems.The main purpose of this paper was to develop a real-life model (actualizing the cost function and adding additional constraints) for handling the feeder bus design and frequency setting problems.
The case study for the research was based on the actual transit network in the Petaling Jaya area, in Kuala Lumpur, Malaysia.Finding the optimum feasible routes in order to reduce the cost function is a vital and difficult task for solving the transit network design problem, which is classified as an NP-hard problem.For this reason, the importance of optimization techniques, particularly metaheuristics, is understood.
Therefore, two recent optimization algorithms, namely the water cycle algorithm (WCA) and the imperialist competitive algorithm (ICA), were considered.The analysis of the objectives shows that the obtained statistical optimization results acquired by the WCA were superior to those attained by the ICA.In terms of solution stability, the WCA also slightly outperformed the ICA.
The optimum number of routes obtained using the WCA was 17, with an average frequency of 6.2 feeder buses per hour in the network.Applying the optimum network resulted in the lowest level of total cost at RM 23,494.8using the WCA, whereas the corresponding costs obtained by the ICA were about 0.8 percentage points greater than found by the WCA.This confirms the trustworthiness of the presented system.Therefore, the proposed improved model could be considered an alternative model for real transit networks.

Figure 1 .
Figure 1.The cost structure of the propose d improve d mode l.

Figure 1 .
Figure 1.The cost structure of the proposed improved model.

Figure 3 .
Figure 3.The case study regions for the Petaling Jaya (PJ) (e xtracte d from Google Maps, 2014).Figure 3. The case study regions for the Petaling Jaya (PJ) (extracted from Google Maps, 2014).

Figure 3 .
Figure 3.The case study regions for the Petaling Jaya (PJ) (e xtracte d from Google Maps, 2014).Figure 3. The case study regions for the Petaling Jaya (PJ) (extracted from Google Maps, 2014).

Figure 4 .
Figure 4. Comparison of obtained best results and main cost terms (RM).

Figure 4 .
Figure 4. Comparison of obtained best results and main cost terms (RM).

Figure 5 .
Figure 5.The deviation percentage associated with the corre sponding optimization algorithms.

Figure 6 .
Figure 6.Comparison of the convergence rate and cost history (in RM) with re spect to the number of ite rations for the considered optimize rs.

Figure 5 .
Figure 5.The deviation percentage associated with the corresponding optimization algorithms.

Figure 5 .
Figure 5.The deviation percentage associated with the corre sponding optimization algorithms.

Figure 6 .
Figure 6.Comparison of the convergence rate and cost history (in RM) with re spect to the number of ite rations for the considered optimize rs.Figure 6.Comparison of the convergence rate and cost history (in RM) with respect to the number of iterations for the considered optimizers.

Figure 6 .
Comparison of the convergence rate and cost history (in RM) with respect to the number of iterations for the considered optimizers.
7 illustrates the location and variation of the total cost (C T ) and its main components, namely C u , C o and C s , or 50 independent runs of each optimization algorithm.It can be observed that the lowest levels for the average cost terms are 24,023.97,16,917.10,6993.46 and 113.45, respectively, for C T , C u , C o and C s with the WCA.The differences between the cost terms of both algorithms for the average cost values are 121, 110.6, 8.9 and 1.5, respectively, for C T , C u , C o and Cs.Susta ina bility 2016, 8, 537 19 of 27

Figure 7 .
Figure 7. Location and variation of the total cost (CT) and its main compone nts, name ly, Cu, Co and Cs.

Figure 7 .
Figure 7. Location and variation of the total cost (C T ) and its main components, namely, C u , C o and C s .

Figure 8 .
Figure 8. Summary of the main costs (in RM) obtaine d by the WCA.

Figure 9 .
Figure 9. Obtaine d best results and each cost term (in RM) using the WCA.

Figure 8 .
Figure 8. Summary of the main costs (in RM) obtained by the WCA.

Figure 8 .
Figure 8. Summary of the main costs (in RM) obtaine d by the WCA.

Figure 9 .
Figure 9. Obtaine d best results and each cost term (in RM) using the WCA.Figure 9. Obtained best results and each cost term (in RM) using the WCA.

Figure 9 .
Figure 9. Obtaine d best results and each cost term (in RM) using the WCA.Figure 9. Obtained best results and each cost term (in RM) using the WCA.

Figure 10 .
Figure 10.Be st solution obtaine d by the WCA (the base map source is extracted from Google Earth, 2014).

Figure 10 .
Figure 10.Best solution obtained by the WCA (the base map source is extracted from Google Earth, 2014).

Figure A2 .
Figure A2.Passenger demand at each rail station in the transit syste m.

Figure A2 .
Figure A2.Passenger demand at each rail station in the transit system.

Table 1 .
Illustration of the total cost with all terms in the proposed improved model.

Table 2 .
De scription of the parameters used in the propose d improve d mode l.CdoiF Feeder dwell operating cost ($/h) CroiT Train running ope rating cost ($/h) CdoiT Train dwe ll ope rating cost ($/h) Cf Fixe d costs ($/h) Cm Mainte nance cost ($/h) Cs Social cost ($/h) CoF Feeder bus operation cost ($/h)

Table 2 .
Description of the parameters used in the proposed improved model.
doiT Train dwell operating cost ($/h) C f Fixed costs ($/h) C m Maintenance cost ($/h) C s Social cost ($/h) C oF Feeder bus operation cost ($/h)

Table 3 .
Location of rail stations and bus stops in the PJ study area.

Table 5 .
Selected LRT stations for public transit passenger questionnaires.

Table 6 .
Passenger demand at bus stops in the Petaling Jaya study area.

Table 7 .
Selected values for the parameters used in the Petaling Jaya study area.

Table 8 .
Comparison of the best solutions obtained for the transit service model using the reported methods.All cost values are based on Malaysian Ringgit (RM).

Table 8 .
Susta ina bility 2016, 8, 537 17 of 27 the applied optimization engines for the improved model.The obtained total cost (CT) is highlighted in bold in Table 8 for the two reported algorithms.Comparison of the best solutions obtained for the transit service model using the reported me thods.All cost value s are based on Malaysian Ringgit (RM).

Table 9 .
Comparison of statistical re sults gaine d by the optimize rs unde r conside ration.

Table 9 .
Comparison of statistical results gained by the optimizers under consideration.

Table 10 .
Statistical optimization results for each cost term.All cost values are in RM. "SD" stands for standard deviation, and the values are in RM. 1

Table 10 .
Statistical optimization re sults for e ach cost te rm.All cost value s are in RM.

Table 11 .
Be st solution obtaine d by the WCA.All cost value s are in RM.

Table 11 .
Best solution obtained by the WCA.All cost values are in RM.