On-line version ISSN 2224-7890
Print version ISSN 1012-277X
S. Afr. J. Ind. Eng. vol.23 n.1 Pretoria Jan. 2012
In this paper, a decision support model for capacity planning of bus rapid transit stations is proposed. Static, deterministic models are used at present for station capacity design, whereas bus passenger flow through stations is dynamic with an associated stochastic element. Passenger arrivals vary over time, and are driven by bus schedules and process variation. The proposed model produces capacity design curves for given input data sets and bus schedules, and assists the station designer without being prescriptive. A case study of model application to the proposed Thibault Station in Cape Town, South Africa, is presented.
'n Besluitsteunmodel vir kapasiteitsbeplanning van sneldiensbusstasies word in hierdie artikel voorgehou. Tans word statiese, deterministiese modelle gebruik vir stasie-kapasiteitsontwerp, terwyl passasiersvloei deur busstasies dinamies is met 'n stochastiese element. Passasiersaankomste varieer oor tyd en word bepaal deur busskedules en proses-variasie. Die voorgestelde model lewer kapasiteitontwerpkurwes vir gegewe insetdata en busskedules en ondersteun die stasieontwerper sonder om voorskriftelik te wees. Die toepassing van die model word geïllustreer met behulp van 'n gevallestudie by die beplande Thibault-stasie in Kaapstad, Suid-Afrika.
Traffic congestion has increased dramatically over the past two decades, threatening the economy of many developing countries and the quality of life of their citizens. Traffic congestion is defined as a phenomenon on transport networks that arises as use increases, and is characterised by slower travelling speeds, longer trip times, and increased queuing. It occurs when traffic demand is greater than the capacity of a road .
In 2003, the Texas Transportation Institute recorded that congestion in the top 85 urban areas in the United States had cost the country $3.7 billion in travel delays, and resulted in 2.3 billion gallons of wasted fuel. Internationally, countries are searching for ways to decrease congestion on roads. This has led to the constant development of new technology, and to different congestion-management strategies being developed and tested worldwide. Some include high-occupancy vehicle (HOV) lanes, congestion pricing, carpooling, vanpooling, ridesharing, bikeways, transit lanes, and various modes of public transport - the last being the prevalent congestion-relieving alternative. Public transport in which city management could invest includes metro rail, light rapid transit (LRT), monorail, suburban rail, standard bus systems, and bus rapid transit (BRT) systems (Wright & Hook ). This paper focuses on the last of those, the BRT system.
A BRT system is a high-capacity public transportation system that carries passengers from one point to another, providing a bus service that is faster and more efficient than an ordinary bus line. At present BRT systems are predominantly planned and designed according to a bus rapid transit planning guide developed by the Institute for Transportation and Development Policy (ITDP), whose head office is in New York (ITDP ). In this guide, the methods (Wright & Hook ) used to calculate the capacity of a BRT system are based on deterministic equations that provide capacity measures for an entire BRT system. We identified the need to investigate the capacity parameters of a single bus station while considering dynamic, stochastic input.
In this paper, we describe the development of a decision support model for BRT system planners, designers, and managers that can assist with the planning and designing of BRT stations, particularly in relation to station capacity. The term 'capacity' refers to the ability of the station to host entering and departing bus passengers - i.e., the floor area required over time. We demonstrate how to estimate percentiles of the capacities, as percentiles afford the decision maker or designer more information than the classic estimated mean.
The importance of this work for the transportation industry is that:
- system planners and designers can use the model to investigate different scenarios and evaluate the impact of each in terms of cost, relieving congestion, and other environmental impacts;
- system planners and designers can use this model to assist with decision-making about the system sizing and capacity requirements of stations;
- the future costs of proposed systems can be estimated based on the projected performance of the simulation models; and
- system managers could use this tool to investigate an existing system to reveal problem areas as well as opportunities for improvement.
We first present some literature related to our study. This is followed by the problem description, a description of the model development using a real-world case study, and the results and conclusions. Finally, we discuss future research on this model.
A literature review showed little information on BRT station capacity planning. The most extensive literature pertains to setting design standards (e.g. Highway Capacity Manual 2000 , Transit Capacity and Quality of Service Manual (Kittelson & Associates )). Pedestrians do not want to wait in dense crowds that cause limited and uncomfortable movement. Therefore standards and procedures were introduced to maintain a desirable pedestrian level of service. These standards and procedures for evaluating pedestrian capacity and levels of service are provided in the Transit Capacity and Quality of Service Manual (Kittelson & Associates ) and in Pedestrian Planning and Design (Fruin ). Furthermore, relationships between pedestrian flow measures, such as speed, space, and delay, are contained in the Highway Capacity Manual 2000 .
Capacity planning of BRT systems focuses more on the corridor's capacity as a whole. The entire network is modelled so that the overall capacity of the system can be evaluated. BRT corridors are usually modelled with simulation, which can be classified according to the level of detail with which they represent the system: microscopic, mesoscopic, and macroscopic simulation (Lieberman & Rathi ). The type of simulation model used depends on the required level of detailed analysis. A microscopic model describes the system entities and their interactions at a high level of detail. A mesoscopic model also describes most of the entities at a high level of detail, but their activities and interactions at a much lower level of detail. Macroscopic models describe the entities and their activities at a low level of detail; lane change manoeuvres, for example, would not be represented (Lieberman & Rathi ). Examples of these traffic simulation models used on corridors include NETSIM (microscopic), NETFLO (mesoscopic), TRANSYT-7F, SYNCHRO, and CORFLO (macroscopic) (Sabra et al. ). These traffic simulation models do not take passenger flow into account, and instead see bus stations as a 'delay' in the entire corridor. This was one of our reasons for investigating the passenger capacity of a single bus station.
We propose estimating percentiles for the capacities, as these are more informative than estimated means. If X has the cumulative distribution function F(x), the p-percentile is defined as F(p) = P(X < p) = q. In most courses in elementary statistics, students are only taught to do inferences based on the mean. This is inadequate, because in practical situations it is often necessary (and sensible) to estimate percentiles of a system performance parameter. A simple example is the required volume of a buffer tank, given a stochastic inflow and outflow of, say, water. Assume the tank is a cylinder with a fixed diameter; then the tank height must be determined to provide some buffer volume. Suppose the nett height due to the flow can be estimated via analytical methods, and the nett height is found to be gamma (α, β) distributed with α=2 and β=3. The mean is αβ = 6, and if the tank is designed for this value, it will overflow 40% of the time. If the designer wants the tank to absorb flow 95% of the time, the tank height should be 14.2 units.
Estimating percentiles is more difficult than estimating means, mainly due to data storage and sorting requirements. Research on this topic dates back to 1976 (Iglehart ; Seila ). Raatikainen  proposed a sequential procedure that requires no storing and sorting, but is difficult to implement. Chen & Kelton  propose two heuristic sequential algorithms and apply these to delay queue processes. We present our percentile estimations in Section 6.
3. PROBLEM OVERVIEW
A primary objective of BRT system operations is to maximise passenger throughput. Factors that affect the speed and ease with which passengers travel throughout the system include rapid and efficient bus operations, adequate facilities, the physical layout of stations, and bus designs. Design parameters of a station that affect the station capacity include the inter-arrival times of passengers, the station layout, the number of bus parking bays or platforms, the number of disembarking and boarding passengers, bus capacities, and bus schedules. Capacity, measured in persons per hour per direction (pphpd), is currently calculated using simple deterministic equations that are accepted worldwide as the standard. Deterministic design parameter values derived from these equations are used to model the entire BRT system. However, the bus station can be considered a complex queuing system, as there are two processes for passenger arrivals and for departures. In the first arrival process, passengers arrive at the station on foot and are removed by buses, while in the second arrival process, random numbers of passengers are delivered to the station by bus and leave the station on foot. This is schematically shown in Figure 1.
During design, models are mainly used to investigate the flow of buses throughout the entire system. However, we identified the need to investigate the effects of these deterministic design parameter values on the capacity of a single bus station. Because these parameters are assigned predetermined, fixed values, any variation in the operations of a bus station is ignored, which is unrealistic. Therefore, it is necessary to investigate what effect these parameters will have on the capacity of a single bus station when introducing elements of variation. The variation is realised through pseudo-randomness and consequent stochastic modelling. An example of a random input is the number of passengers entering the station at some time during the day. This time-dependent event and others can be studied using a specific probability distribution to determine the number of passengers entering the station at that time of day. Considering these, we propose a dynamic, stochastic model to accommodate the time-dependent variation of the real-life system (bus station).
Next, a modelling case study of a real-life station is presented, which served to guide the model development and evaluate the resulting proposed stochastic modelling approach.
4. CASE STUDY: THIBAULT STATION
The study is based on Thibault Station, a bus station that forms part of the BRT system in Cape Town, South Africa. At the time of writing, Thibault Station was under construction. A stochastic model of Thibault Station was developed in which the random arrival of passengers and the scheduled arrival of buses were realised so that the impact on station capacity could be examined. The stochastic model is a time-dependent, discrete-event computer simulation model that allows for what-if questions and the replication of passenger arrivals using known input data. This stochastic model, together with an appropriate user interface, forms the proposed decision support model.
This section provides information on the system operations that were simulated for Thibault Station. A concept model that helps to plan the simulation model is presented, followed by an explanation of the model architecture and model functions.
4.1 Concept model description
Thibault Station is an enclosed bus station with glass sliding doors, staffed ticket booths, pre-board fare verification (by means of turnstiles), and four platforms. The operations of Thibault Station involve:
Passenger arrivals and flow:
- There are two entrances to the station. Entrance 1 is the main entrance where passengers buy tickets from a ticket booth before proceeding to the platforms where outgoing passengers embark and incoming passengers disembark. Entrance 2, which does not have a ticket booth, is situated at the opposite end of the station. This entrance only allows access to passengers who already have a ticket. The majority of these passengers have monthly or season tickets, allowing them quicker access to the platforms. They can also enter through Entrance 1.
- Passengers verify their tickets at the turnstiles and then proceed to their desired platform waiting area.
- The buses operate in a separate dedicated BRT bus lane, allowing little interference from other modes of traffic.
- Two routes depart from Thibault Station, with different types of buses serving each route. The TO1 route uses bi-articulated buses with a passenger capacity of 120. The TO2 route uses smaller buses seating only 40 passengers.
- The buses operate according to a schedule; bus arrivals are thus considered to be deterministic.
- The fleet size is changed throughout the course of the day in accordance with variation in demand at various times during the day.
- The station has four platforms to facilitate passenger arrivals and departures.
- Platform 1 is allocated to TO1 buses, and transfers passengers from various urban bus stations to the CBD of Cape Town. Platform 3 is also a TO1-allocated platform, and routes out to Blaauwberg, approximately 30 km away in a northerly direction.
- TO2 buses stop at Platform 2 and proceed to the CBD. Platform 4, situated opposite Platform 2, is also a TO2-allocated bus stop, and transfers passengers to Cape Town international airport.
- There are two parking berths at each platform, allowing parking space for an additional bus waiting for the platform.
A schematic layout to explain the movement of the buses to and from their dedicated platforms, as well as the position of the platforms at Thibault Station, is shown in Figure 2. TO1 buses follow the cycle represented by the broken lines, while TO2 buses follow the dotted lines. The two drop-off stations (A and B in Figure 2) are dummy stations, and were used in the model to represent all other stations that form part of a cycle. The cycle time of each bus on each route to these dummy stations was incorporated into the model, which ensures model validity.
The entities in the model are the passengers and the buses. These are the 'discrete units of traffic' in the simulation model, as stated by Schriber & Brunner .
Model restrictions include:
- Vehicles operating between 04:00 and 23:00, every day of the week.
- TO1 buses having a capacity of 120 passengers.
- TO2 buses having a capacity of 40 passengers.
- The number of boarding passengers can only be as large as the number of available seats on the particular bus.
Model assumptions include:
- The various bus types always stop at their dedicated platforms. If buses arrive simultaneously at the same platform, or if one arrives while the parking berth is already occupied, that bus will stop behind the bus already occupying the designated parking berth to await its turn.
- Dwell, boarding, and disembarkation times are assumed to be incorporated into the passenger arrival times.
- Buses always arrive according to scheduled times.
- The bus capacities, already defined for the different types of buses, represent the passenger seating capacity, and are taken as the maximum number of passengers that the vehicle can accommodate. Standing passengers are not included in the model.
The concept model is based on the operations mentioned above, and was used to build the simulation model.
4.2 Model architecture
The decision support model has a user interface implemented in Microsoft® Excel, while the simulation package Arena of Rockwell Automation was used to implement the stochastic part of the model. A second interface between MS Excel and Arena made the user input available to Arena in the required format. This can be seen in Figure 3.
The simulation model is built according to the operations specified in the concept model, and is aimed to assist with capacity planning at Thibault Station. Although this model is based on the operations of Thibault Station - i.e., Thibault's specific platform layout, bus schedule, passenger arrival rates, and so on - it is sufficiently generic to be applied and used at similar BRT stations.
The key functions of the simulation model are that:
- users are able to define input data for the model to suit the specifications of their station under investigation, and
- analysis of capacity parameters can be performed, and required station capacities can be predicted, from the model results.
The next section describes the functionality of the developed simulation model by presenting extracts of the 'user interface' spreadsheets.
4.3 Model functionality
The user interface spreadsheets, which are briefly explained below, allow users to define and change certain input data. The user interface spreadsheet shown in Figure 4 accepts passenger arrival data at Thibault Station. The data vary over the day in hourly resolution, and are expressed as a percentage of the daily maximum expected number of passenger arrivals.
The user interface spreadsheet for the input data related to the arrival of buses is shown in Figure (part of the TO1 route specification). An MS Excel spreadsheet was developed for each route, which also serves as a method of assigning the attributes to the arriving buses. The simulation model reads the 'expected' arrival times as constants, thus emulating the bus schedule.
Next, the number of passengers arriving by bus is assigned a random quantity less than bus capacity, which represents the number of passengers disembarking. This yields a secondary source of arrivals at the station, and also affects the station capacity.
The decision support model was built, verified, and validated to ensure the correctness and reasonableness of the model's functions and outcomes. A deterministic version of the model was most useful, since outcomes could be predicted based on known input.
The performance measurements provided by the model are:
- the percentiles of the number of passengers waiting and duration per station platform, and
- the average waiting time of the passengers.
These were all estimated using point estimators and 95% confidence intervals, according to the procedure of Law & Kelton . Some scenarios that were evaluated are described in the next section.
5. EXAMPLE OF SCENARIOS
To demonstrate the purpose and usefulness of the decision support model, two scenarios were constructed that differ in terms of their mean inter-arrival times of passengers entering the station. Scenario 1 is based on the original input data obtained from the station designer, while for Scenario 2 the passenger inter-arrival times were decreased by 25%. This implies a higher arrival rate of passengers for Scenario 2, while the bus schedule and the number of buses allocated are maintained. These scenarios are used to show the various ways in which this model could be used to obtain system response information, as well as the consequences of decisions made.
We shall present two types of output: estimated means, and estimated percentiles of system performance measures. We introduced the idea of percentiles in Section 2, and here we briefly outline how we estimated the percentiles for this problem. The simulation package Arena allows for the specification of histograms as output, and can determine confidence intervals per histogram class, provided a sufficient number of simulation replications were executed. From the histogram counts it is possible to determine the percentage of time a parameter was at a certain level, and the percentiles can then be deduced from these results. One possible drawback of this approach is that the histogram classes must be specified beforehand (Seila ).
To illustrate the approach, we simulated an M|M|1 queue with traffic intensity ρ = 0.9. The number of entities in the system X in steady state follows the geometric distribution f(x) = ρχ(1-ρ), x=0,1,2,... These exact probabilities were calculated for the first 30 values of X (0 to 29), and are shown as 'Exact Pr(x)' in Table 1. The column labelled 'Sum Pr(x)' in the table is the exact cumulative distribution values for X. The simulated values are listed in the columns labelled 'Estim. Pr(x)' and 'Sum Estim. Pr(x)'. The estimations are based on 5,000 pseudo-independent replications of 2,000 time units each, with 500 time units as warm-up. The absolute deviations ('Abs. Dev.') are reasonably small, so the simulation model is an acceptable estimator for the percentiles. It can, for example, be concluded from Table 1 that there will be 10 or fewer entities in the system with a probability of 0.69. Note that the cumulative distribution is discontinuous in this case, due to the discrete nature of X.
Some scenario results for the Thibault Station design are presented next.
6. SCENARIO RESULTS
The estimated expected number of passengers occupying each platform for every hour of the operating day is shown in Figure 6 for current arrival estimations (Scenario 1), and a 25% decrease of inter-arrival times (Scenario 2). We did not test for statistically significant differences, as the numerical differences are large over several time periods.
The estimated passenger waiting times per platform and the daily average number of passengers waiting for the respective scenarios are shown in Table 2. The effect of passenger increase is clear: decreasing the arrival times by 25% results in an average increase of 37% in the waiting times and an 83% increase in the number of passengers waiting at Platform 1. (Note that Platform 3 is located opposite Platform 1, and Platform 2 and 4 form an opposite pair - see Figure 2)
One can also obtain the time frequencies for the daily number of passengers waiting. The cumulative frequency diagram of Platform 1 (Scenario 1) is shown in Figure 5, while the frequency duration of passenger presence on Platform 1 (Scenario 1) is shown in Table 3.
From Table 3, the following can be deduced:
- For 57% of the operating time in a day, the occupation will range between 10 and 20 passengers. See tag (1) in Table 3.
- Platform 1 occupation ranges between 20 and 30 passengers for 3.303 hours of the day. See tag (2) in Table 3.
- For 16.52 hours of the 19-hour day, the occupation will be less than 40 passengers. See tag (3) in Table 3.
A further interpretation is possible using the results in Table 4.
The passenger frequencies show the number of observations (passengers) counted in the specified ranges. Valuable conclusions can be made by using the passenger frequencies together with the time frequencies. The following is an example for the range between 20 and 30 passengers:
- For 74.4% of the operating time (see Table 3), up to 59.9% (Table 4) of the total number of passengers made use of the service that day, where they will wait in a group comprising fewer than 30 passengers.
- Up to 59.9% of the total number of passengers will make use of the service 74.4% (see Table 3) of the operating time, where they will wait in a group comprising fewer than 30 passengers. This means the remainder of the expected passengers (40.1%) will make use of the service in only 25.6% of the operating hours, which is 19 - 14.87 = 4.87 hours of the day.
Results were also generated for the waiting areas of Platform 2, 3, and 4. As mentioned before, Platform 3 is located opposite Platform 1 at Thibault Station (see Figure 2), and Platform 2 and 4 form an opposite pair. Aggregated results for the two pairs of platforms were collected during the simulation runs (Platforms 1 and 3, Platforms 2 and 4), and the results can be used to design the aggregated capacity of the station. The combined observed response for the Platform 1-3 pair (Scenario 1) is shown in Figure 8 and Table 5 (Time frequencies).
Table 6 shows the passenger frequencies for the combined number of passengers on Platform 1 and Platform 3 (Scenario 1).
Similar interpretations of the data given for the passenger count of Platform 1 hold for Figure 6 and Table 6; consider for example the cumulative line in Figure 6. Examples of how this information could be used as a basis for design decisions are as follows:
- Different alternatives can be weighed up against one another: for an area to accommodate passengers 80% of the operating time, the cumulative line shows that it must be designed and built to fit up to 60 passengers. For the other 20% of the time the area will not be able to accommodate the expected demand and these passengers will have to wait at alternative waiting areas. This is known as the 80th percentile at which the station is designed, and is illustrated in Figure 6 by the dotted line indicated as 'A'. Knowing this information, the space requirements for that area can be calculated. Using values of standard practices taken from the Highway Capacity Manual 2000 , the basic minimum space requirement for a single standing passenger is determined at 0.3 m2. For a passenger facility, an area of 0.75 m2 is used as a buffer zone for each pedestrian. Using this figure, the space required to accommodate passengers 80% of the time would be 52 χ 0.75 m2 = 39 m2.
- Conversely, the cumulative line could also be used to obtain information about the percentile for a specific space: if the proposed station is designed with a waiting area of 75 m2, the number of passengers it can hold is calculated as 100 passengers (75 m2 ÷ 0.75 m2/passenger). Referring back to the graph in Figure 6, the cumulative line shows a 95th percentile (operating time) for 100 passengers. The dashed line labelled 'B' shows how this value can be read from the graph. This means that an area of 75 m2 will be sufficient to accommodate passengers 95% of the time.
There is a trade-off between the percentile for which a station is designed and the cost of the station. Is the designer going to design a station at the 80th percentile, which will be less expensive, and accommodate fewer passengers, or at the 95th percentile, which will accommodate more passengers? A larger percentile will require more space, which in turn increases the costs. These graphs are useful for such decisions, and more can be created using the proposed decision support model.
The information presented for Scenario 1 is valuable for decision-making purposes; however, it will also be valuable to compare Scenario 1 with Scenario 2 to see what impact a 25% difference in inter-arrival rates will have on the proposed station design - i.e., a busier process. The decision support model allows for exactly that: testing the sensitivity of the passenger estimates.
The previous analysis for the waiting areas of Platform 1 and Platform 3 combined (Scenario 1) can be repeated for Scenario 2. The results of the time percentiles for Platform 1 and Platform 3 combined (Scenario 2) are presented in Figure 7 and Table 7.
The 80th time percentile for Scenario 1 requires an area that holds up to 52 passengers. The 80th time percentile for Scenario 2 (Figure 7) requires an area that holds between 100 and 120 passengers. There is a vast difference between the two scenarios. If the station was built according to Scenario 1 and the passenger estimates were incorrect, the impact could be very negative for the same number of buses.
The scenarios were constructed to investigate the effect the variation would have on the performance of the system. More importantly, it also showed how the decision support model can be used as a tool to investigate and evaluate many different aspects relating to the capacity of a station.
This method generates easily understandable graphs that display significant information about the performance measurements of a system, the impacts that might be experienced, and the capacity capabilities for different parameters. Outcomes of scenarios can be generated with relative ease, investigated, and compared against each another.
This paper presented a decision support model for capacity estimations of BRT stations. A case study of Thibault Station in Cape Town, South Africa, was presented and showed how this model can be used to investigate capacity parameters for the station. The model delivers useful design measures, such as the average number of passengers per waiting area, the average waiting times, and the frequencies and percentiles of passenger counts.
Through the scenarios presented in the case study, it is possible to affirm the purpose and usefulness of the model. The model offers a variety of options to its users, giving them control over investigations to create concrete evidence for design decision support. The direct implications of the proposed decision support model are:
- users are able to define input data for the model to suit the specifications of the station they are investigating;
- the performance of a station can be predicted for different scenarios; and
- capacity parameters can be analysed.
We have demonstrated the need to estimate percentiles when doing simulation studies, because that type of output parameter affords the decision maker more information than point estimators of the mean.
8. FUTURE RESEARCH
The following suggestions are made for future work:
- Building a corridor by replicating the simulation model:
The model developed in this study can be used as a basis for expanding the simulation model into an entire BRT corridor, consisting of more than one station. The model can be duplicated and linked, on which numerous capacity studies could be done.
- Integrating a bus schedule into the simulation model:
This study did not focus on implementing a bus schedule into the simulation model, as a schedule was provided by the designers. A mathematical programming technique can be used to develop an optimal bus schedule. This could then be linked to the simulation model on which further investigations could be done, such as the costs related to the optimal bus schedule requirements versus its efficiency.
- Developing a fully generic stochastic simulation model:
The current stochastic model is based on the station configuration of Thibault Station, i.e. four platforms and four waiting areas. Model adjustments would need to be made for other station configurations. The same holds for the input data spreadsheets that were developed for the model. Although the spreadsheets are user-friendly and values are easily adjustable, the bus schedule spreadsheets are designed for a cycle consisting of four station stops in a cycle. For a cycle consisting of fewer or more station stops, MS Excel and Arena adjustments are needed. A more generic approach could be valuable in future studies.
 Chen, E.J. & Kelton, W.D. 2006. Quantile and tolerance-interval estimation in simulation, European Journal of Operational Research, 168(2), pp 520-540. [ Links ]
 Fruin, J.J. 1971. Pedestrian planning and design. Metropolitan Association of Urban Designers and Environmental Planners. [ Links ]
 Iglehart, D.L. 1976. Simulating stable stochastic systems; VI. quantile estimation, Journal of the Association for Computing Machinery, 23, pp 347-360. [ Links ]
 Kittelson & Associates, Inc. 1999. TCRP web document 6 - Transit capacity and quality of service manual. Transportation Research Board, National Research Council. [ Links ]
 Law, A.M. & Kelton, W.D. 2000. Simulation modeling and analysis, 3rd edition, McGraw-Hill. [ Links ]
 Lieberman, E. & Rathi, A.K. 1975. TRB Special Report 165 Traffic Flow Theory, Transportation Research Board. [ Links ]
 Raatikainen, K.E.E. 1990. Sequential procedure for simultaneous estimation of several percentiles, Transactions of the Society for Computer Simulation, 7(1), pp 21-44. [ Links ]
 Sabra, Z., Wallace, C.E. & Lin, F. 2000. Traffic analysis software tools, Transportation Research Board, pp 12-16. [ Links ]
 Schriber, T.J. & Brunner, D.T. 2005. Inside discrete-event simulation software: How it works and why it matters, in Proceedings of the 2005 Winter Simulation Conference, (M.E. Kuhl, N.M. Steiger, F.B. Armstrong & J.A. Joines, eds), Simulation Society of INFORMS, www.informs-sim.org, pp 167-177. [ Links ]
 Seila, A.F. 1982. A batching approach to quantile estimation in regenerative simulations, Management Science, 28(5), pp 573-581. [ Links ]
 The American Heritage Dictionary of the English Language. 2003. U.w. 'traffic congestion'. Houghton Mifflin. [ Links ]
 Transportation Research Board. 2000. Highway capacity manual 2000, National Research Council. [ Links ]
 Wright, L. & Hook, W. 2007. Bus rapid transit planning guide. Institute for Transportation and Development Policy, p. 48. [ Links ]
 Yang, F., Ankenman, B.E. & Nelson, B.L. 2008. Estimating cycle time percentile curves for manufacturing systems via simulation, INFORMS Journal on Computing, 20, pp 628-643. [ Links ]