arXiv:1506.07424v1 [] 24 Jun 2015

Simulating the Effects of
Various Road Infrastructure Improvements to
Vehicular Traffic in a Busy Three-road Fork
Marian G. Arada

Merly F. Tataro

Jaderick P. Pabico

Polytechnic University of the
Philippines-Taguig
Taguig City 1632, Metro
Manila

Polytechnic University of the
Philippines-Taguig
Taguig City 1632, Metro
Manila

Institute of Computer Science
University of the Philippines
Los Baños
College 4031, Laguna

ABSTRACT
Using microsimulations of vehicular dynamics, we
studied the effects of several proposed infrastructure
developments to the mean travel delay time ∆ and mean
speed Σ of vehicles passing a busy three-road fork, particularly in the non-signalized roundabout junction
of Lower Bicutan, Taguig City, Metro Manila. We
designed and implemented multi-agent-based microsimulation models to mimic the autonomous driving
behavior of heterogeneous individuals and measured the
effect of various proposed infrastructure developments
on ∆ and Σ. Our aim is to find out the best infrastructure development from among three choices being
considered by the local government for the purpose of
solving the traffic problems in the area. We created
simulation models of the current vehicular traffic situation in the area using the mean travel times τ of
statistically sampled vehicles to show that our model
can simulate the real-world at a significance level of
α = 0.05. Based on these models, we then simulated
the effect of the proposed infrastructure developments
on ∆ and Σ and used these metrics as our basis of
comparison. We found out that the proposed widening
of one fork from two lanes to three lanes has the most
improved metrics at the same α = 0.05 compared to
the metrics we observed in the current situation. Under
this infrastructure development, the ∆ increases linearly (R2 = 0.98) at the rate of 1.03 s, while the Σ
decreases linearly (R2 > 0.99) at the rate of 0.14 km/h
per percent increase in the total vehicle volume V.

Categories and Subject Descriptors
K.3.7.1 [Computing Methodologies]:
Artificial
Intelligence—Multi-agent systems; K.5.1.2 [Computing
Methodologies]: Modeling and Simulation—Model

verification and validation; L.9.5 [Applied Computing]: Operations Research—Transportation

1. INTRODUCTION
The increasing number of vehicles that use the road
at the same time is causing congestions on road networks especially at road sections where nominal widths
change such as intersections, three-way forks, and
roundabouts. One real-world example of such width
change is the Bicutan Roundabout (BR) in Taguig
City, Metro Manila. Commuters passing through this
roundabout experience delay in travel everyday, especially during rush hours on week days, because of
vehicular congestion in the area. BR is located at the
Upper Bicutan area and has an inscribed diameter of
34 m. It is bounded to the East by the Department
of Science and Technology (DOST) campus, to the
North by the Philippine National Construction Corporation (PNCC) campus, to the West by the Philippine
National Railways (PNR) crossing, and to the South
by North Daanghari. BR serves as the T-intersection
between General Paulino Santos Avenue (GPSA) which
lies along the ENE-WSW line and SLEX’s East Service Road (ESR) which lies along the NEN-SWS line
of the area (Figure 1). The nominal center of BR is
approximately 100 m WSW from DOST’s main gate
along GPSA, 150 m SWS of the ESR and San Martin
de Porres fork, and 106 m ENE of PNR crossing.
Jeepneys, buses, trucks, taxis, AUVs, motorcycles, tricycles, and bicycles are the usual type of vehicles that
pass the area. Upper Bicutan, a barangay under the
urbanized City of Taguig, is situated on the western
shore of Laguna de Bay and is bordered in the south
by Muntinlupa City, in the southwest by Parañaque
City, in the west by Cainta, Rizal, in the northeast by
Taytay, Rizal and in the north by Makati City, Pasig
City and the Municipality of Pateros. According to the
2010 census, 62,570 out of the city’s total population of
613,343 are residents of Upper Bicutan [13].
The thrust of the city government administration is to
make Taguig City a “Premier City” in the Philippines.
However, the physical and financial development of a
city is greatly dependent, among others, on the efficiency of its transportation support system, particularly
how the city manages the vehicular flow along its road

Figure 1: The area map of Upper Bicutan, Taguig City showing the relative location of BR [9].
networks. The efficient flow of vehicles carrying human
and non-human capitals results in the efficient flow of
resources. This translates to a productive citizenry and
successful businesses. Because of this, the city government must make its move in pursuing efficient and practical steps to at least mitigate, if not totally eradicate,
traffic congestion which must result into improved vehicular flow. Vehicular flow may be quantified by the following metrics and we say that improvement has been
achieved if there is a statistically significant difference
the measurement at a particular significance level α as
compared to the current situation:
1. Average vehicular delay time while in BR (∆) –
This metric measures the total amount of time
that the vehicle is in the stopped position (i.e.,
when the vehicle has a speed of 0 km/h); and
2. Average vehicular speed while in BR (Σ) – This
metric measures the mean speed of a vehicle while
traveling within the BR area.
Ideally, we wanted ∆ to approximate zero while we
wanted Σ to be near but not more than the legal maximum speed. In realistic situations, however, ∆ > 0 and
Σ is even way below the legal minimum speed. These
two metrics complement each other because for some

agents such as those who drive big trucks, an almost
zero ∆ and a very low Σ is much more desirable than
a ∆ > 0 and a relatively high Σ. This is because it
requires much more fuel to move the heavy vehicles from
a full stop than to keep them moving at low Σ. For
agent drivers of light- and medium-weight vehicles, a
∆ > 0 but a very high Σ is desirable because fuel burnt
at low Σ is greater than when idle (for at least within
some reasonable time).
There are two ways the city can make its move:
(1) Improve the infrastructure of the road network
that will best optimize the ∆ and the Σ, or (2) implement responsive but efficient time-dependent traffic
schemes that will significantly improve both metrics.
Both moves are very costly if the government is to
try all possible schemes and infrastructure just for the
simple reason of finding the optimal one.
In terms of finding the optimal infrastructure improvement for BR, one should not do it via a trial-and-error
way because it means the following: (1) building the
infrastructure, (2) letting it tried by the commuters for
a reasonable amount of time just to statistically compile
the ∆ and the Σ metrics, and then (3) removing it to
build the next one. By the time the government has
compiled the data it needed, it has already spent sig-

nificant amount of money, time, and other resources.
Besides, the data collected might have already been
biased due to what we call “data collection dependency”
brought about by the memory of the agents in their
experiences. This means that for a trial-and-error with
n iterations, the data collected on the (i + 1)th iterate
might have been interactively affected by the improvement made during the ith iterate, unless of course, after
every iterate we revert back to the current infrastructure for the same amount of time just so both the drivers
and the commuters will forget their ith experiences.
In terms of implementing responsive yet efficient traffic
schemes, most local governments have the habitual
liking of doing it in the trial-and-error way, as experienced by all of this paper’s authors. If for example
the local Traffic Management Office (TMO) is considering n traffic schemes, they will: (1) implement the
ith scheme for a t period of time (where t is empirically selected to be long enough for the drivers and
the commuters to adapt to it), (2) compile the ∆
and the Σ metrics, and (3) repeat the process for the
(i + 1)th scheme. The local TMO will then compare
the schemes and decide based on the metrics collected
without considering (maybe even knowing) that there
is an interactive effect on the (i + 1)th scheme the
experiences and adaptation of the driver and commuter
agents during the implementation of the ith scheme.
One efficient way that may be used to try out ideas
and schemes in solving traffic congestions is to use
computer-based simulation, a model created to reflect
real-world situations. Once the model is constructed,
solutions can be tried out to determine its impact on
traffic congestion [3], without the need for a costly
real-world trial-and-error implementation.
In this
study, we used a multi-agent framework to simulate the
object-following behavior of heterogeneous agents that
share a road network. The road network reflects the
scaled-down dimension of a real-world network, such
as that of the BR, servicing scaled-down dimensions
of real-world vehicles. Our model can be characterized as microscopic, stochastic, discrete time-step, and
behavior-based. The model is microscopic because it
simulates the behavior of the vehicles themselves while
on the road in the presence of, and interacting with,
other vehicles. The model is stochastic because it is
based on random events that follow some distributions that we observed in the real-world. The model
is discrete time-step because at each unit time, the
agents concurrently change their respective positions,
directions, and speeds driven only by their respective
simulated behaviors in response to their respective
interactions with the current state of the road network.
The model is behavior-based because the traffic flow
algorithms that we used are based on a psycho-physical
car following model for vehicles moving along the length
of the road [6, 8]. We used a rule-based algorithm for
lateral movement that simulates the lane-changing
behaviors of the agents [1, 5].
We present in this paper the result of our experiments
in searching for the optimal infrastructure improvement
for BR from among three choices that are being consid-

ered by the local government. We present our results in
searching for the optimal traffic schemes for the current
infrastructure in another paper [12]. Our study had
three objectives: (1) To show that our model reflects
the current real-world scenario at a statistical significance of α = 0.05 using the mean travel time τ of the
sampled vehicles as metric, (2) to use the model to simulate the effect on ∆ and on Σ of three infrastructure
developments (ID) in BR currently being considered by
the local government, and (3) simulate the effect on ∆
and on Σ the 10%, 50%, and 100% increase in vehicular volume under the best ID. Based on our results,
we found out that: (1) The average travel time of the
simulated vehicles while at the BR is not significantly
different at α = 0.05 from the observed average travel
time of the sampled real-world vehicles, (2) the best
ID is when the South bound lane of the ESR is widen
up to three lanes continuing to the west bound lane of
GPSA from BR to PNR, and (3) the ∆ increases linearly
(R2 = 0.98) at the rate of 1.03 s per percent increase
in the total vehicle volume V, while the Σ decreases linearly (R2 > 0.99) at the rate of 0.14 km/h per percent
increase in V.

2. MODEL
A computer model that controls an agent’s behavior
(i.e., the driver) in terms of its reaction to the vehicle in
front of it while at the same lane is called a car-following
model [10]. An agent-driven vehicle is said to be following when it is constrained by the speed of the moving
vehicle in front of it, and that driving at the agent’s
desired speed will lead to a collision. When a driver
agent is not constrained by another vehicle, it is said
to be free and travels, in general, at its desired speed.
The actions of the following agent is defined by the
agent’s acceleration, although some models, for example
the car-following model developed by Gipps [8], define
the agent’s actions through the agent’s speed. Some
car-following models only describe the agents’ behavior
when they are currently following another vehicle, while
other models determine the agents’ behavior in all situations. We believe, however, that to model the real-world
decision-making capabilities of human drivers, a carfollowing model should identify both of the following:
1. The current state S the vehicle is in; and
2. What actions A are desirable at this state.
Most car-following models use several states to describe
the following agents’ behavior. Most models use:
1. Sf : A state for free driving, where the vehicles are unconstrained and the respective driver
agents try to achieve their desired speeds (subject
of course to pertinent legal speed limits of the road
network);
2. Sn : A state for normal following, where the following agents adjust their respective speeds with
respect to the speeds of the vehicles in front of
them; and

3. Se : A state for an emergency deceleration,
where agents try to decelerate to avoid a collision
with the vehicle in front.
Throughout the paper, we used the notations summarized in Table 1 to describe the kinematic quantities and
model outputs.

2.1

Sf Free Driving State
The agents at the Sf state try to accelerate or decelerate to achieve its current desired speed. If the current
speed vi of the ith vehicle is higher than its desired
speed videsired , the agents uses the normal deceleration
rate (−anormal
) to slow down to the desired speed. If
i
vi < videsired , the agents use its maximum acceleration
rate amax
to reach videsired at the shortest time possible.
i
normal
The −ai
and amax
are parameters of the simulai
tion and they are dependent on the type of vehicle the
agent is driving. The acceleration rate of the ith agent
is expressed as
 max
vi < videsired
 ai
ai =
0
vi = videsired
 normal
ai
vi > videsired

The time ti it will take for the ith vehicle to achieve its
desired speed is
( vi
vi < videsired
amax
i
ti =
vi
vi > videsired
anormal

3.1 Interview with Stakeholders
We conducted interviews in order for us to understand
the current plans of the two local governments that have
jurisdictions over the BR: the Taguig City Traffic Management Office (TMO) and the Metro Manila Development Authority (MMDA). The interview served as some
sort of leveling off with the intended stakeholders. We
also wanted the result of this research endeavour to be
put to use, particularly because two of the three authors
of this paper pass this area on a daily basis. Thus,
if either or both local governments will implement the
result of our research, we will greatly benefit from such.
We also asked permission from the two governments to
allow us to conduct scientific observations so that we
can quantify the current vehicular flow in BR.

3.2 Observation of Current Situation
We employed the help of twelve enumerators whom we
assigned to different identified points within the BR
junction. Each point of entrance to and exit from BR
has designated persons to record the time of entrance
and exit of the vehicle, and the type of vehicle that
passed through. Each vehicle is identified by their tag
or plate number. We recorded a series of actual observations during peak and non-peak periods for us to be
able to get a statistically accurate data and determine
the highest and lowest traffic volume. From these observations, we were able to obtain the mean travel time τ of
sampled vehicles, as well as the respective distributions
of each vehicle type. To identify the specific points in
the BR, we divided the area into six routes as follows:

i

2.2 Sn Normal Following State
While at the normal following state Sn , the acceleration rate of the ith vehicle is given by an asymmetrical
Gazis-Herman-Rothery (GHR) model [4, 7, 14]. The
acceleration is computed as
±

ai = r ±

vis
(vi−1 − vi ),
(xi−1 − Li−1 − xi )t±

where r ± , s± , and t± are model parameters. The
parameters r + , s+ , and t+ are used if vi ≤ vi−1 , while
the parameters r − , s− , and t− are used if vi > vi−1 .
Notice here that the final acceleration of the ith vehicle
is given as max{ai , ai−1 }.

2.3

Se Emergency State
Under this state, the agents use a deceleration rate that
prevents collision and extends δx. This deceleration rate
is given by
(
0.5(vi −vi−1 )2
min{anormal
, ai−1 − xi−1 −L
}
vi > vi−1
i
i−1 −xi
ai =
normal
normal
min{ai
, ai−1 + 0.25ai
vi ≤ vi−1

3. METHODOLOGY
We present in this section the activities that we performed in this study. We also present here the metrics
that we measured, as well as the statistical analysis that
we employed to analyze the results of the study.

1. Route 1 is the route from PNR to DOST Campus
Eastbound along the GPSA passing through the
BR;
2. Route 2 is the route from PNR to PNCC
Campus, Eastbound along GPSA from the PNR
to BR, and then turning left through the BR, and
up to the Northbound lane of ESR going towards
the PNCC Campus;
3. Route 3 is the route from DOST campus to PNR
Westbound along the GPSA passing through the
BR;
4. Route 4 is the route from the DOST Campus to
the PNR Crossing, Westbound along the GPSA
up to the BR, and then turning right through the
BR, and up to the Northbound lane of ESR going
towards the PNCC Campus;
5. Route 5 is the route from PNCC Campus to
PNR, Southbound along the ESR, turning right
through BR, and then Westbound along the
GPSA to the PNR; and
6. Route 6 is the route from PNCC Campus to
the DOST Campus, Southbound along the ESR,
turning left through BR, and then Eastbound
along the GPSA to the DOST Campus.

Table 1: Alphabetical listing of notations, mathematical variables, and abbreviations used in this
paper including their respective descriptions.
Notation Description
δx
Space headway between xn−1 and xn in m
δv
difference in speed between xn−1 and xn in m/s
an
Acceleration of the nth vehicle in m/s2
Ln
length of the nth vehicle in m
sn
Effective length of the nth vehicle in m
T
Reaction time in s
vn
Speed of the nth vehicle in m/s
vndesired
Desired speed of the nth vehicle in m/s
xn
Longitudinal position of the nth vehicle in m
α
Statistical level of significance
∆
Mean delay time per vehicle in s
Σ
Mean vehicular speed in m/s
τ
Mean vehicular travel time s
τo
τ of the sampled vehicles from observation
τs
Mean τ of the simulated vehicles
V
Total volume of vehicles
ANOVA Analysis of Variance using the F-Statistics
BR
Bicutan Roundabout
DOST
Department of Science and Technology
ESR
East Service Road
GPSA
General Paulino Santos Avenue
ID
Infrastructure Development
MMDA Metro Manila Development Authority
PNCC
Philippine National Construction Corporation
PNR
Philippine National Railways
PUP-T
Polytechnic University of the Philippines-Taguig
SLEX
Southern Luzon Expressway
TUP-T
Technological University of the Philippines-Taguig
TMO
Traffic Management Office

3.3 Simulation of Current Situation
We conducted a replicated microsimulation study of the
vehicular flow under the current BR using the data on
respective distribution of vehicles by type. We replicated the study to n = 10 and computed the mean τ .
We wanted to know if the model can statistically reproduce the observed vehicular traffic data. Statistically,
the respective differences of the τ between the observed
(τo ) and the mean simulated (τs ) runs must not be different from zero at α = 0.05. We used analysis of variance statistics (ANOVA) to evaluate two hypotheses,
the null hypothesis H0 and the alternative hypothesis
Ha as follows:
H0 : There is no significant difference between τo and τs
at α = 0.05.
Ha : There is significant difference between τo and τs at
the same α.

3.4 Microsimulations of Proposed Developments
We created the respective replicated n = 10 microsimulation studies of the vehicular flow using the observed
current vehicular distribution when under the different
planned infrastructure developments (ID) of the BR.
These IDs are based on the C6 Project blueprint, the
medium-term development plan by the City Mayor of

Taguig Lani Cayetano, and DOST’s Monorail Project
which is currently underway. The proposed developments are:
1. Infrastructure Development 1 (ID1): Widening of
the GPSA into three lanes but only those lanes
towards the BR. That is the east bound lane from
PNR to BR and the west bound lane from DOST
to BR.
2. Infrastructure Development 2 (ID2): Widening
of the GPSA into three lanes but only the eastbound lane from PNR to BR continuing from BR
to DOST gate.
3. Infrastructure Development 3 (ID3): Widening of
the south bound lane of the ESR into three lanes
up to the BR continuing to the west bound lane
of GPSA from BR to PNR.
We computed the respective means of ∆ and Σ for all
infrastructure improvements, including that of the current infrastructure we termed here as ID0 (i.e., from
ID0 through ID3) and conducted an analysis of variance (ANOVA) statistics to find out whether the respective means are significantly different from each other at
α = 0.05. That is, we want to test two hypotheses, the
null hypothesis H0 and the alternative hypothesis Ha as
follows:

H0 : The absolute pairwise differences |∆ID0 − ∆ID1 | =
· · · = |∆ID1 − ∆ID2 | = |∆ID1 − ∆ID3 | = |∆ID2 − ∆ID3 | =
0, where the = sign in this sentence would mean “not
significantly different from” or “statistically equal to.”
Ha : Any of the following respective pairwise absolute
differences is true: |∆ID0 − ∆ID1 | > 0, . . . , |∆ID1 −
∆ID2 | > 0, |∆ID1 − ∆ID3 | > 0, or |∆ID2 − ∆ID3 | > 0.
We have the same set of hypotheses for Σ.

3.5 Increased Traffic Volume at Best ID
We rerun the respective replicated microsimulation
studies under ID0, ID1, ID2, and ID3 but with respective increases of 10%, 50%, and 100% in V, making sure
that the distribution by vehicular type is preserved.
Here, we wanted to know if the benefits of the proposed
improvements will be carried with the increase in V and
if so, up to how much increase. The assumed increase
in V is just a natural reaction of the driving agents
when there is a perceive improvement in the current
situation. The improvement of the vehicular flow due to
infrastructure improvements will attract more agents,
increasing V in the area, and thereby might degrade
the expected designed benefits of the improvement in
the long run.

4. RESULTS AND DISCUSSION
4.1 Actual Current Vehicular Flow
Tables 2 and 3 summarize the distribution of vehicle
per route and per vehicle type, respectively, that we
compiled during the actual observation. Route 4 has
the most number of total vehicles that passed during
the observation period with 23.79%. This is followed
by Routes 5, 1, and 2, in non-increasing order. Route
4, undeniably has the most number of vehicles since
Route 4 is the only route that may be used by residents
from the Lower Bicutan area who are trying to go to
the Metropolitan Manila Area via the C-5 Diversion.
Aside from the Lower Bicutan residents, the route is
also extensively used by students and employees of several Higher Education Institutions (HEIs) and government agencies such as the Technological University of
the Philippines-Taguig (TUP-T), the Polytechnic University of the Philippines-Taguig (PUP-T), the LGUfunded Taguig City University, and DOST’s national
offices and laboratories. Route 5 is extensively used
by trucking companies with medium-length trucks carrying heavy cargoes. These vehicles are coming off from
the C-5 Diversion going to the Southern Luzon Area
via the SLEX. Route 1 is obviously used by commuters
coming from the SLEX area going to the residential
Lower Bicutan or to any of the HEIs and government
agencies. Route 2 is the reverse of Route 5, where
mostly medium-length trucks pass through carrying cargoes from the Southern Luzon area via SLEX going to
Metropolitan Manila Area through the C-5 Diversion.
From Table 3, motorcycles and cars combined account
for about 70% of the vehicles that use BR. This is followed by bicycles, jeepneys and vans. The heavier vehicles, such as the many-wheelers and buses only account
for about 4% combined. This distribution by type is
actually expected because the area is used extensively

by residents and commuters rather than by non-human
cargoes. Table 4 shows the mean vehicular dimensions
by vehicle type. We used these dimensions to scale the
respective agents in our simulations.
Table 2: Distribution of the total number of
vehicles that passed through the BR during the
observation period.
Route
Number of Vehicles Percentage (%)
Route 1
778
20
Route 2
719
18
Route 3
524
13
Route 4
934
24
Route 5
815
21
Route 6
156
4
Total
3,926

Table 3: Distribution of the vehicles per type.
Type
Percentage (%)
Motorcycle
38.30
4 × 8 wheeler
0.64
4 × 6 wheeler
2.75
Van
5.88
Jeepney
9.60
Car
31.73
Bus
0.38
Bicycle
10.72
Total
100.00

Table 4: The average dimension of observed
vehicles per type.
Vehicle
Length Width
Type
(m)
(m)
Motorcycle
2.00
1.5
4 × 8 wheeler
6.59
1.5
4 × 6 wheeler
5.41
1.5
Van
5.50
1.5
Jeepney
4.00
1.5
Car
4.50
1.5
Bus
11.54
2.5
Bicycle
1.45
0.5

4.2 Evaluation of the Model against the
Observed
Table 5 summarizes the statistics we computed after
conducting ANOVA on τo and τs . Since we see that the
F statistics αF = 0.6873 > α = 0.05, then we accept
the null hypothesis H0 and we say that the difference
between the mean sampled observed data τo and the
mean simulated data τs is not significantly different from
zero at α ≈ 0.05. Thus, our microsimulation model was
able to mimic the driver’s behavior in the real world
with a guarantee of being correct 1 − α = 0.95 of the
time.

4.3 Comparison of the IDs
Tables 6 and 7 show the respective ANOVA tables of
the mean ∆ and mean Σ for ID0, ID1, ID2, and ID3
Σ
where we find that α∆
F < α and αF < α. Based on the

Table 5: The ANOVA table comparing τo and τs using the F statistics. SOV means Source of Variation
and DF means Degree of Freedom
SOV
DF Sum of Squares Mean Square
F
αF
Replication
9
1,633.05
181.45 39.21 < 0.0001
τo vs. τs
1
0.80
0.80
0.17
0.6873
Error
9
41.65
4.63
Total
19
1,675.50

Σ
computed α∆
F and αF , we accept the alternate hypothesis Ha saying with confidence 1 − αF = 0.9999 that
at least one of the absolute pairwise difference between
∆IDx and ∆IDy , and between ΣIDx and ΣIDy , ∀x 6= y, is
significantly greater than zero. Figures 2 and 3 show the
respective mean ∆ and mean Σ of IDx, ∀x = 0, . . . , 3.
Based on the DMRT, we see that ID1 will not improve ∆
nor Σ. In fact, both ∆ and Σ under ID1 is no different,
statistically, than the ∆ and Σ of the current infrastructure, and thus, implementing ID1 will just be a waste of
taxpayer’s money. We also see that ID2 has significantly
improved ∆ and Σ compared to the current infrastructure. However, ID3 has the most significant improvement in terms of both ∆ and Σ compared to the current
infrastructure. From Figure 3, we see that the respective Σ of ID2 and ID3 are not statistically different from
each other, so either ID can provide the best Σ. Based
on the two metrics, however, we say that ID3 can provide the improvement that we desire for the vehicular
flow.

Figure 2: The respective mean ∆ of the various IDs with DMRT grouping. Means with the
same letter are not statistically different from
each other by DMRT.

4.4 Effect of Increased Traffic Volume
Figures 4 and 5 show the respective effect of increasing V
from 10% to 100% on ∆ and Σ while under ID3. As V
is increased one percent at a time, ∆ increases linearly
with R2 = 0.98 at a rate of 1.03 s while Σ decreases
linearly with R2 > 0.99 at a rate of 0.15 km/h. At this
linear rate of decrease, all vehicles will come to a halt
(i.e., Σ = 0) at V ≈ 180%.

5. CONCLUSION
In this paper, we presented a multi-agent-based
microsimulation of a real-world vehicular traffic based
on the existing psycho-physical car following model

Figure 3: The respective mean Σ of the various IDs with DMRT grouping. Means with the
same letter are not statistically different from
each other by DMRT.
for simulating the movement of vehicles along the
lane [6, 8] and a rule-based algorithm for simulating
the lane-changing behaviors of the vehicles [1, 5].
We simulated a real-world infrastructure involving six
routes merging into a three-road fork via an unsignalized roundabout. This roundabout is the BR located in
Taguig City, Metro Manila. From the observed distribution of vehicles per route and per vehicular type, we
simulated the vehicular traffic of the current situation
and collected the τ metric. Based on the ANOVA, we
found out that τo and τs are not statistically different
from each other based on the F statistics. This means
that our simulation model can represent the real-world.
We then simulated three proposed infrastructure developments (ID1, ID2, and ID3) and found out that ID3
offers the least ∆ and at the same time the fastest Σ.
We increased the V from 10% to 100% under ID3 and
found out that ∆ increases linearly with V, while Σ
decreases linearly. Based on these results, we provide
the following conclusions:
1. Our microsimulation model can stastically represent the real-world;
2. The best ID for BR is ID3; and
3. Under ID3, V has a positive linear effect on ∆ and
a negative linear effect on Σ.

6. RECOMMENDATIONS
We recommend to the stakeholders within BR, such
as the City Government of Taguig through its TMO,
as well as the Engineering Section of the MMDA, the
residents of Lower Bicutan, the public and private

Table 6: The ANOVA table comparing the different IDs in terms of ∆ using the F statistics. SOV
means Source of Variation and DF means Degree of Freedom
SOV
DF Sum of Squares Mean Square
F
α∆
F
Replication
9
2,328.63
258.74
5.12
0.0004
ID
3
3,556.14
1,185.38 23.47 < 0.0001
Error
27
1,363.39
50.50
Total
39
7,248.15

Table 7: The ANOVA table comparing the different IDs in terms of Σ using the F statistics. SOV
means Source of Variation and DF means Degree of Freedom
SOV
DF Sum of Squares Mean Square
F
αΣ
F
Replication
9
14.73
1.64
5.63
0.0002
ID
3
11.83
3.94 13.56 < 0.0001
Error
27
7.85
0.29
Total
39
34.41

Figure 4: The respective mean ∆ of ID3 when V
was increased from 10% to 100%. The line represents the regression with ∆ = 1.03 × V + 72.64
(R2 = 0.98).

Figure 5: The respective mean Σ of ID3 when V
was increased from 10% to 100%. The line represents the regression with Σ = −0.15 × V + 26.93
(R2 > 0.99).

employees within BR such as the DOST, the TUP-T,
and the PUP-T, and the private sector, particularly
SM Corporation, to help call for the implementation
of infrastructure developments described in ID3. We
caution, however, that we may be able to find a better
vehicular flow if hybrid improvements of ID1 through
ID3 are considered. We further caution that it may
be cheaper to implement signalized traffic schemes
together with the ID. As we have mentioned before,
we have already reported our findings with various
traffic schemes elsewhere [12], while we are currently
working [2] on the interactive effects of ID and traffic
schemes to both ∆ and Σ.

3. Evaluation of evacuation plans for Tsunamivulnerable shoreline communities;
4. Designing re-routable and adaptive traffic plans
for time-responsive traffic schemes; and
5. Incorporating the “greedy” behavior of Jeepney
and tricycle drivers in the car-following model.

8. ACKNOWLEDGMENT

1. Studying the effects of installing railings in busy
sidewalks;

This research effort is funded by the DOST-SEI Graduate Scholarship Program with M.G. Arada and
M.F. Tataro as scholar-grantees taking the Master
of Information Technology graduate program in the
University of the Philippines Los Baños under the
research supervision of J.P. Pabico. We thank the PTV
Group [11] for allowing us to use their visual simulation
framework free of charge for academic and research
purposes.

2. Pedestrian dynamics for a very, very large academic building beyond its designed full capacity;

The following are the respective contributions of the
authors:

7. FUTURE WORK
Several scientific research efforts have already forked out
from this initial investigation. These are:

1. M.G. Arada implemented the computational solution, conducted the real-world observation, the
computational experiments and the statistical
analyses, and prepared the final manuscript;
2. M.F. Tataro implemented the computational solution and conducted the real-world observation and
computational experiments; and
3. J.P. Pabico formulated the computational solution
to the problem, conducted and interpreted the statistical analyses, and edited the final manuscript.

[13] The City Government of Taguig. The Official Website of Taguig City, Metro Manila, Philippines,
2013.
[14] Qi Yang and Haris N. Koutsopoulos. A microscopic
traffic simulator for evaluation of dynamic traffic
management systems. Transportation Research
Part C: Emerging Technologies, 4(3):113–129,
1996.

All authors declare no conflict of interest.

References
[1] Talal Al-Shihabi and Ron R. Mourant. Toward
more realistic driving behavior models for
autonomous vehicles in driving simulators. Transportation Research Record, 1843:41–49, 2003.
[2] Marian G. Arada, Merly F. Tataro, and Jaderick P.
Pabico. Simulating the interactive effects of infrastructure improvements and traffic schemes to vehicular travel time in the bicutan roundabout. Transactions of the National Academy of Science and
Technology, Philippines, 36(1), 2014. (accepted).
[3] Joe G. Bared and Evangelos I. Kaisar. Median
u-turn designs as an alternative treatment for left
turns at signalized intersections. Institute of Transportation Engineers Journal, 72, 2002.
[4] Robert E. Chandler, Robert Herman, and
Elliott W. Montroll. Traffic dynamics: Studies in
car following. Operations Research, 6(2):165–184,
1958.
[5] Martin Fellendorf and Peter Vortisch. Microscopic
traffic flow simulator VISSIM. In Jaume Barcelo,
editor, Fundamentals of Traffic Simulation, pages
63–94. Springer, New York, 2010.
[6] Hans-Thomas Fritzsche. A model for traffic simulation. Transportation Engineering Contribution,
5:317–321, 1994.
[7] Denos C. Gazis, Robert Herman, and Richard W.
Rothery. Nonlinear follow the leader models of
traffic flow. Operations Research, 9(545–567), 1961.
[8] Peter G. Gipps. A behavioural car-following model
for computer simulation. Transportation Research
Part B: Methodological, 15(2):105–111, 1981.
[9] Google Corporation.
Google
http://www.googlemap.com.

[12] Merly F. Tataro, Marian G. Arada, and Jaderick P.
Pabico. Selecting the best traffic scheme for the
bicutan roundabout using microsimulations of multiple driver-agents. In Proceedings of the 6th UPLB
CAS Student-Faculty Research Conference, NCAS
Auditorium, UPLB, 2013. CAS Standing Committe on Research.

Map,

2013.

[10] Robert Herman and Richard W. Rothery. Carfollowing and steady state flow. In J. Almond,
editor, Proceedings of the 2nd International Symposium on the Theory of Road Traffic Flow, pages
1–13, 1963.
[11] PTV Group. The Official Website of PTV Group,
2013. http://company.ptvgroup.com/en.

