Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015

Selecting the Best Traffic Scheme for the Bicutan Roundabout:
A Microsimulation Approach with Multiple Driver-Agents
Merly F. Tataro,1 Marian G. Arada,2 and Jaderick P. Pabico3
1,2

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

3

1

mftataro@pup.edu.ph, 2mgarada@pup.edu.ph, 3jppabico@up.edu.ph

Abstract – Understanding the dynamics of vehicular traffic is important to planning and
improving shared public roadways, not only to effectively and efficiently facilitate the
controlled movement of vehicles, but also to guarantee the safety of individuals in the
vehicles. Because data from real-world traffic is difficult to measure, and conducting
replicated experiments on humans has ethical and economical questions, selecting the
best scheme out of the proposed vehicular traffic controls has been confined only through
simulation.
We present the result of our microsimulation study on the effects of six traffic schemes
T={t0, t1, ..., t5} on the mean total delay time (Δ) and mean speed (Σ) of vehicles at the
non-signalized Bicutan Roundabout (BR) in Upper Bicutan, Taguig City, Metro Manila,
with t0 as the current traffic scheme being enforced and ti>0 as the proposed ones. We
present first that our simulation approach, a hybridized multi-agent system (MAS) with
the car-following and lane-changing models (CLM), can mimic the current observed
traffic scenario C at a statistical significance of α=0.05. That is, the respective absolute
differences of the Δ and Σ between C and t0 are not statistically different from zero. Next,
using our MAS-CLM, we simulated all proposed ti>0 and compared their respective Δ and
Σ. We found out using DMRT that the best traffic scheme is t3 (i.e., when we converted the
bi-directional 4-lane PNR-PNCC road into a bi-directional 1-lane PNR-to-PNCC and 3lane PNCC-to-PNR routes during rush hours). Then, we experimented on converting BR
into a signalized junction and re-implemented all t3 with controlled stops of S={15s, 45s}.
We found out that t3 with a 15-s stop has the best performance. Finally, we simulated the
effect of increased in vehicular volume V due to traffic improvement and we found out
that t3 with 15-s stop still outperforms the others for all increased in V={10%, 50%,
100%}.
Keywords – microsimulation, MAS, CLM, vehicular traffic, Bicutan Roundabout
I. INTRODUCTION
Transportation is an important aspect of both the economy and the lifestyle of the population it is
serving because it links the various destinations of people, goods, and services. A smoother and faster
vehicular and pedestrian traffic would mean people reaching their workplaces on time resulting into
higher number of productive hours. Increased productivity results into higher wages and profit which in
return encourage people to buy more vehicles. The problem of traffic congestion in our country is the
result mostly of having more automotive growth than the available road can accommodate. According to
the Metro Manila Development Authority’s (MMDA) report recently, vehicular speed drops to 5kph
during rush hours along Epifanio Delos Santos Avenue (EDSA). The MMDA also reports that major
roads are no longer sufficient to accommodate the consistently rising traffic volume, which is made worse
by the continuing use of outmoded traffic signals, worsening road conditions, and decreasing availability
of efficient public transport [1]. As a result, mere accident at an important intersection, for example, is
likely to paralyze many streets for hours [2]. To address this seemingly worsening problem, especially in
the Metro Manila area, the government has implemented various infrastructure projects including
widening of roads, and construction of overpasses and underpasses along major thoroughfares. One
1

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
among the many bottleneck areas in Metro Manila that is being considered for infrastructure
improvement is the Bicutan Roundabout (BR) in Taguig City. Both the MMDA and the City Government
of Taguig (CGT) are separately considering building infrastructure to improve the traffic condition in the
area [1,3].
BR is located at the Upper Bicutan area and has an inscribed diameter of 34m. 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) that lies along the ENE-WSW line and SLEX’s East Service Road (ESR)
that lies along the NEN-SWS line of the area (Figure 1). The nominal center of BR is approximately
100m WSW from DOST’s main gate along GPSA, 150m SWS of the ESR and San Martin de Porres fork,
and 106m ENE of PNR crossing. Jeepneys, buses, trucks, taxis, AUVs, motorcycles, tricycles, and
bicycles are the usual vehicles that pass the area.

Figure 1.

The area map of Upper Bicutan, Taguig City showing BR.

Because BR is a bottleneck, it has become common to see and experience traffic congestion in the area
especially at peak hours both in the morning and in the afternoon. Even though BR is non-signalized, the
local government unit has to employ traffic enforcers to manually intervene with controlling the traffic
flow during these hours to ease congestion in the area [3]. The enforcers also provide some sort of relief
to commuters when they see that city official are at least doing something [4]. The traffic enforcers react
to the traffic situation by implementing ad hoc solutions that only solve the various current situations,
instead of a programmed solution that can solve the big problem in the long run [5,6]. This problem is
expected to continue to worsen in the future due to the rapid growth of population that is aggravated
further by the increasing number of vehicles passing in the area [1].
In an earlier works by Arada, et al. [5,6], they simulated the effect of various proposed infrastructure
developments in BR on the mean total vehicular delay time (Δ) and on the mean speed (Σ). They found
out that the best infrastructure development (ID best) that minimizes Δ and maximizes Σ, out of the three
that are being considered, is the one that widens the eastbound lane of the GPSA from the PNR towards
the DOST gate. They were even successful in finding out the effect to Δ and Σ of increasing the vehicular
volume up to 100% in the BR area under the IDbest scenario. However, before we support to recommend
the implementation of the findings of Arada, et al. [5,6], we caution the reader that there may exists a
traffic scheme that can approximate the benefits of ID best in terms of Δ and Σ. If we can find this, the
government is posed to save a lot in construction costs and time. Thus, in this research effort, we state the
following hypotheses:
H0: (Null Hypothesis) Among the schemes that are available, there does not exist a traffic scheme that
can at least approximate the optimal Δ and Σ that the IDbest can provide.
H1: (Alternate Hypothesis) Among the schemes that are available, there does exist a traffic scheme
that can at least approximate the optimal Δ and Σ that the IDbest can provide.
2

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
To test these hypotheses, we used the methodology by Arada, et al. [5,6] in employing a hybridized
microsimulation with a multi-agent system (MAS) and the car-following and lane-changing models
(CLM) on the effect of six traffic schemes T={t0, t1, ..., t5} on the mean Δ and mean Σ. The t0 is the current
ad hoc scheme being implemented while ti>0 are the ones with suggested long-term programming.
Because data from real-world traffic is difficult to measure, and conducting replicated experiments on
humans has ethical and economical questions, selecting the best scheme out of the proposed vehicular
traffic controls can only be realistically done through simulation. Thus, first, we show statistically that
our MAS-CLM can simulate the real-world scenario, just like how Arada, et al. [5,6] did it. Second, we
use MAS-CLM to simulate various traffic schemes and compile their respective mean Δ and mean Σ.
Lastly, we simulate the effect of increased vehicular volume on Δ and Σ under the best traffic scheme
condition.
II. THE MAS-CLM MODEL
Our MAS-CLM model, described in detail in the works of Arada, et al. [5,6] and of Tataro, et al. [4],
uses a car-following model developed by Herman and Rothery [7]. Relating this model to our MAS, a
driver-agent is said to be in the following state when it is constrained by the speed of the moving vehicles
that is directly in front of it. Speeding at the driver-agent's desired speed may lead to collision under this
state. However, when the driver-agent is not constrained by another vehicle, then it is said to be in the free
state and travels, on the average, at its desired speed. The action of the driver-agent while in the following
state is defined by the driver-agent's acceleration, although in some models like the one developed by
Gipps [8] the action is dependent on the driver-agent's speed. Our MAS-CLM model uses the following
states to describe the driver-agent's behavior (the pertinent details of this model, including the
mathematical foundations of the equation of motion has already been discussed elsewhere (see Arada, et
al. [5,6]):
1. Free Driving State – The driver-agent at this state is unconstrained and tries to achieve its goal
speed vgoal, subject of course to pertinent legal speed limits being imposed at the portion of the
road. If the driver-agent's current speed vcurrent > vgoal, the driver-agent uses the normal deceleration
rate –anorm to slow down to vgoal. If vcurrent < vgoal, the driver-agent uses the maximum acceleration
rate amax to speed up to vgoal at the shortest possible time tmin = vgoal/amax. The anorm and amax are
parameters of the driver-agents and dependent on the vehicular type being simulated.
2. Normal Following State – The driver-agent under this state adjusts its speed with respect to the
speed of the vehicle directly in front of it. Its acceleration is defined by the asymmetrical GazisHerman-Rothery (GHR) model described by Chandler, et al. [9], Gazis, et al. [10], and Yang and
Koutsopoulos [11].
3. Emergency Deceleration State – The driver-agent while at this state tries to decelerate to avoid
colliding with the slow-moving vehicle directly in front of it.
III. METHODOLOGY
We present in this section the different activities we performed to test the hypotheses that we stated
above. We also present here the metrics that we used and the statistical analyses that we employed to
analyze our results.
A. Observation of the current traffic condition
We assigned twelve (12) research assistants to different identified points within the BR to record the
time of entrance and exit of the vehicle, and the type of vehicle that passes through. They were tasked to
record the plate numbers of the vehicles as well as to take videos of the same vehicles that pass through
certain points of entrance and exit within the BR. We conducted a series of actual observations during
peak and off peak hours to get an accurate data and determine the highest and the lowest traffic volume.
From these observations, we obtained the mean travel time of sampled vehicles, as well as the respective
3

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
distributions of each vehicle type. The following are the points of entrances and exits that we designated
as observation points:
1. Point A (PNR) is the route from PNR passing through the BR to Point B or C;
2. Point B (PNCC) is the route from PNCC passing through the BR to Point A or C; and
3. Point C (DOST) is the route from DOST passing through the BR to Point A or B.
B. Simulation of the current traffic condition
We conducted a replicated microsimulation of the vehicular ﬂow under the current BR using the data
on respective distributions of vehicles by type. We replicated the study to n = 10 and computed Δ and Σ.
We wanted to know if the simulated data t0 can statistically reproduce the current condition C, under the
same traffic scheme. Statistically, the respective differences of Δ and Σ between C and t0 must not be
different from zero at α = 0.05. We used analysis of variance statistics (ANOVA) to evaluate two
hypotheses as follows:
H2: (Null Hypothesis) There is no signiﬁcant difference between C and t0 at α = 0.05, i.e., |C – t0| ≡ 0,
where the symbol ≡ means statistically equal at 5% confidence level.
H3: (Alternate Hypothesis) There is signiﬁcant difference between C and t0 at the same α, i.e., at 5%
confidence, |C – t0| > 0.
An example visualization of one our microsimulation runs is shown in Figure 2. The visualization uses
the VisSim framework for transport planning, traffic engineering and traffic simulation [12].

Figure 2.

A screenshot of an example microsimulation run on BR using VisSim as a simulation
visualization framework. The red line across the road means that the vehicles coming
from that direction are stopped, while the green line means that the vehicles are
moving.

C. Microsimulation of the proposed traffic schemes
We created the respective replicated n = 10 microsimulations of the vehicular ﬂow using the vehicular
distribution in C under the five proposed traffic schemes t1, t2, ..., t5 as follows:
t1: Controlled one-way scheme for short-timed intervals of 60s time cycle: 30s stop those coming
from Point B followed by 30s stop those coming from Point C;
t2: Same as t1 but the time interval is 120s;
t3: Designation of one lane from Point A to Point B and three lanes in the opposite lane (from Point
B to Point A) during peak hours instead of the normal two-lane, two-way traffic scheme;
t4: Combination of t1 and t3; and
t5: Combination of t2 and t3.
We computed the respective means of Δ and Σ for all proposed t1, t2, ..., t5, including that of the t0 and
conducted ANOVA to ﬁnd out whether the respective means are signiﬁcantly different from each other at
α = 0.05. We, thus, state the following hypotheses (where Δ(ti) means the Δ under ti):
4

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
H4: (Null Hypothesis) The absolute pairwise differences |Δ(t0)−Δ(t1)|≡···≡|Δ(t1)−Δ(t2)|≡···≡|
Δ(t2)−Δ(t3)|≡···≡|Δ(t4)−Δ(t5)|≡0.
H5: (Alternate Hypothesis) At least one of the following respective pairwise absolute differences is
true: |Δ(t0)−Δ(t1)|>0, ..., |Δ(t1)−Δ(t2)|>0, ..., |Δ(t2)−Δ(t3)|>0, ..., or |Δ(t4)−Δ(t5)|>0.
We have a similar null H6 and alternate H7 hypotheses for Σ, where Σ(ti) means the Σ under ti:
H6: (Null Hypothesis) The absolute pairwise differences |Σ(t0)−Σ(t1)|≡···≡|Σ(t1)−Σ(t2)|≡···≡|Σ(t2)−Σ(t3)|
≡···≡|Σ(t4)−Σ(t5)|≡0.
H7: (Alternate Hypothesis) At least one of the following respective pairwise absolute differences is
true: |Σ(t0)−Σ(t1)|>0, ..., |Σ(t1)−Σ(t2)|>0, ..., |Σ(t2)−Σ(t3)|>0, ..., or |Σ(t4)−Σ(t5)|>0.
D. Signalized BR with controlled stops and increased volume
We experimented on converting BR into a signalized junction and re-implemented the best scheme
found above with controlled stops of 15s and 45s, making sure that the distribution by vehicular type is
preserved. Here, we wanted to know if the beneﬁts of the best proposed scheme will be carried with
controlled signal stops, and if so, up to how much increase? We then reran the respective replicated
microsimulations with respective increases of 10%, 50%, and 100% in vehicular volume V, again making
sure that the distribution by vehicular type is preserved. We wanted to know if the beneﬁts of the
proposed schemes will still be carried with the increase in V and if so, up to how much increase within the
0% to 100% range? 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 ﬂow due to
improved traffic schemes will attract more driver-agents, increasing V in the area, and thereby might
degrade the expected designed beneﬁts of the schemes in the long run.
IV. RESULTS AND DISCUSSION
A. Observed vs. Simulated Current Traffic Condition
Table 1 summarizes the statistics we computed after conducting ANOVA on C and t0. We see that the F
statistics αF = 0.6873 > α = 0.05, which prompts us to accept H2 and say that the difference between the
mean sampled observed data C and the mean simulated data t0 is not signiﬁcantly different from zero at α
≈ 0.05. We can, thus, say that 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.
Table 1.

ANOVA table comparing C and t0 using the F statistics. SOV means sources
of variations, and DF means degrees of freedom.
DF

Sum of
Squares

Replication

9

1,633.05

C vs. t0

1

Error
Total

SOV

Mean Square

F

αF

181.45

39.21

< 0.0001

0.80

0.80

0.17

0.6873

9

41.65

4.63

19

1,675.50

B. Comparison of the Traffic Schemes
Tables 2 and 3 show the respective ANOVA of Δ and Σ. We find that αF(Δ) < 0.0001 and αF(Σ) <
0.0001, which prompts us to accept H5 and H7 for Δ and Σ, respectively. This means that at least one of
the absolute pairwise differences between Δ(ti) and Δ(tj), and those between Σ(ti) and Σ(tj), ∀i ≠ j, are
significantly greater than zero for Δ and Σ, respectively.

5

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
Table 2.

ANOVA table comparing Δ at T (i.e., Δ(t0) vs. Δ(t1) vs. ... vs. Δ(t5)).
DF

Sum of
Squares

Mean Square

F

αF

5

37,354.14

7,470.83

319.01

< 0.0001

Error

54

1,264.61

23.42

Total

59

38,618.75

SOV
T

Figure 3 shows the pairwise comparison of the mean Δ and the mean Σ under various T={t0, t1, ..., t5}
using the Duncan's Multiple Range Test (DMRT). We can see from the table that the best traffic scheme is
t3, both in terms of Δ and Σ. From the figure, we can also see the relative mean Δ and mean Σ of IDbest
found by Arada, et al. [5,6]. Even though the data for ID best are not included in our DMRT, we can see
relatively that even both Δ(ti>0) and Σ(ti>0) outperform IDbest. We, thus, accept H1 and say in fact that all
traffic schemes that we considered here outperform ID best at both Δ and Σ.
Table 3.

ANOVA table comparing Σ at T (i.e., Σ(t0) vs. Σ(t1) vs. ... vs. Σ(t5)).
DF

Sum of
Squares

5

917.46

183.49

Error

54

20.49

0.38

Total

59

937.95

SOV
T

Figure 3.

Mean Square

F

αF

483.58

< 0.0001

Pairwise comparison of (a) mean Δ and (b) mean Σ over various traffic schemes T={t0,
t1, ..., t5}. Means with the same small letters (a) or the same capital letters (b) are not
significantly different by DMRT at α=0.05. ID best is the mean Δ and mean Σ of the best
infrastructure development found by Arada, et al. [5,6].

C. Effect of signalized t3
When we converted BR into a signalized junction and reimplemented t3 with controlled stops of 15s
and 45s (i.e., t3+15s and t3+45s), respectively, we found out that the metrics Δ and Σ improved statistically
over the non-signalized scheme (Figure 4). Specifically, Δ(t3+15s) is statistically better than both Δ(t3) and
Δ(t3+45s). Consequently, Δ(t3) is statistically better than Δ(t3+45s). This means that among the three
6

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
schemes, Δ is at minimum under t3+15s. Additionally, Σ(t3) and Σ(t3+15s) are statistically the same, as
well as Σ(t3) and Σ(t3+45s). However, Σ(t3+15s) is statistically better than Σ(t3+45s). This means that Σ is
at maximum under either t3 or t3+15s. Since the best Δ is under t3+15s and the best Σ is under either t3 or
t3+15s, we recommend that t3+15s be implemented in BR. Thus, we say that installing signal lights with
15s controlled stops within the scheme described by t3 will provide statistically significant improvement
to the traffic flow.

Figure 4.

Pairwise comparison of (a) mean Δ and (b) mean Σ over the traffic schemes t3, t3+15s,
and t3+45s. Means with the same small letters (a) or the same capital letters (b) are not
significantly different by DMRT at α=0.05.

D. Effect of increased V to a signalized t3
Figure 5 shows the effect to the metrics Δ and Σ of increasing V when the traffic scheme t3+15s is
implemented. The respective increases are 10%, 50%, and 100% of the observed volume. The Δ increased
when V was increased (Figure 5a), while Σ decreased when V was increased (Figure 5b). Increasing V to
10% increased Δ significantly for t3+15s, but the increase is still not enough to be greater than either t3 or
t3+45s (Figure 5a inset). At additional V of 10%, t3+15s is still a better scheme than t3, while t3 is better
than t3+45s. Equations 1 through 6 show the respective linear estimates of Δ(t3), Δ(t3+15s), Δ(t3+45s),
Σ(t3), Σ(t3+15s), and Σ(t3+45s), as V is increased up to 100%. In these equations, V(+) is the percent
increase in V, r2 is the coefficient of determination, and an asterisk and an “ns” over a coefficient
respectively means that the coefficient is significantly and not significantly different from zero via a twotailed t-test.
Δ(t3) = 4.58* V(+) – 2.69ns,
r2 = 0.97
(1)
*
(+)
*
2
Δ(t3+15s) = 0.47 V + 21.92 ,
r = 0.94
(2)
*
(+)
ns
2
Δ(t3+45s) = 0.91 V + 34.82 ,
r = 0.90
(3)
Σ(t3) = –0.15* V(+) + 18.43*,
r2 = 0.95
(4)
*
(+)
*
2
Σ(t3+15s) = –0.08 V + 21.02 ,
r = 0.96
(5)
*
(+)
ns
2
Σ(t3+45s) = –0.11 V + 18.47 ,
r = 0.85
(6)
Since all slope coefficients are significantly different from zero, we can say that further increasing V(+)
to 50% and to 100% increased Δ(t3), Δ(t3+15s), and Δ(t3+45s) linearly by 4.58s/% (Equation 1), 0.47s/%
(Equation 2), and 0.91s/% (Equation 3), respectively. It is clear to see that up to an increased V of 100%,
t3+15s is the best scheme with respect to minimizing the rate of increase for Δ.
In Figure 5b, increasing the additional V up to V(+)=100% linearly decreases Σ. The linear decrease for
Σ(t3), Σ(t3+15s), and Σ(t3+45s) are, respectively, 0.15kph/% (Equation 4), 0.08kph/% (Equation 5), and
0.11kph/% (Equation 6). With respect to minimizing the rate of decrease for Σ, the best traffic scheme is
t3+15s. Thus, it is recommended that a signalized version of t3 with stops of 15s be implemented in BR.

7

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015

Figure 5.

Effect on (a) mean Δ and (b) mean Σ on increasing V to 10%, 50%, and 100%. The
markers are mean data while the lines are the regression lines defined in Equations 1
through 6.

V. SUMMARY, CONCLUSION, RECOMMENDATION AND EXTENSION
In this paper, we presented our MAS-CLM which we used to simulate six traffic schemes in the
nonsignalized BR. We found out that our MAS-CLM can statistically simulate the real-world at α=0.05.
Using the same MAS-CLM, we simulated T={t0, t1, ..., t5} and collected Δ and Σ as metrics. DMRT at
α=0.05 shows that t3 outperforms the other schemes ti≠3. We were also able to show that the traffic
schemes ti>0 have better mean Δ and mean Σ than the one found in an earlier work involving a costly
infrastructure development [5,6]. We converted BR into a signalized roundabout and reimplemented all
t3, while we also considered increasing V from 10% up to 100%. We found out that t3 with 15-s stop-cycle
outperforms the other two schemes. Because ID best performed poorly compared to ti>0, we highly
recommend to either MMDA or to CGT to consider implementing t3 with 15-s stop-cycle instead of
developing the infrastructure described by ID best. We caution the reader, however, that we may be able to
find better mean Δ and mean Σ than the one we found in this study if we consider the combination effects
of infrastructure development and traffic scheme. Our study on the combination effects are already
underway as an extension of this work. We will report the result of this extension in our future paper.
V. ACKNOWLEDGMENTS, DISCLOSURE OF EARLIER PRESENTATION, AND AUTHOR CONTRIBUTIONS
This research endeavor is funded by DOST ASTHRDP Graduate Scholarship Program with MFT and
MGA as scholar-grantees taking MIT at UPLB under the research advisorship of JPP. VisSim is courtesy
of PTV Planung Transport Verkehr AG of Germany.
An earlier version of this paper was presented as an oral paper [13] during the 6th Annual University
of the Philippines Los Baños (UPLB) College of Arts and Sciences (CAS) Student-Faculty Research
Conference (SFRC 2013) held at the NCAS Auditorium, UPLB on 13 December 2013.
8

Asia Pacific Journal of Multidisciplinary Research Vol. xx No. yy, Month 2015
The following are the respective contributions of the authors: (1) MFT implemented the computational
solution, conducted the real-world observation, implemented the computational experiments and
statistical analyses, and prepared the final paper; (2) MGA implemented the computational solution,
conducted the real-world observation and the computational experiments, and edited the final paper; and
(3) JPP formulated the computational solution to the problem, conducted and interpreted the statistical
analyses, and edited the ﬁnal manuscript. All authors declare no conﬂict of interest.

[1]
[2]
[3]
[4]

[5]

[6]

[7]

[8]
[9]
[10]
[11]
[12]

[13]

REFERENCES
MMDA. 2012. The Official Website of the Metro Manila Development Authority
(www.mmda.gov.ph)
A.C. Robles. 2012. Commuters suffer. D+C Development and Cooperation/E+Z Entwicklung
und Zusammenarbeit (12):161 (www.dandc.eu).
CGT. 2012. The Official Website of the City of Government of Taguig (www.taguig.gov.ph)
M.F. Tataro, M.G. Arada and J.P. Pabico. 2013. A microsimulation study on the effects of traffic
schemes to vehicular travel time in a real-world busy roundabout. International Research
Conference in Higher Education (IRCHE 2013), The Manila Hotel, Manila, 03-04 October
2013.
M.G. Arada, M.F. Tataro and J.P. Pabico. 2013. Simulating the effects of various road
infrastructure improvements to vehicular traffic in a busy three-road fork, In A.A. Sioson (ed.)
Proceedings of the 11th National Conference on Information Technology Education (NCITE
2013), Villa Paraiso, Mambajao, Camiguin, 24-26 October 2013.
M.G. Arada, M.F. Tataro and J.P. Pabico. 2014. Simulating the effects of various road
infrastructure improvements to vehicular traffic in a busy three-road fork. Philippine Information
Technology Journal (submitted).
R. Herman and R.W. Rothery. 1963. Car-following and steady state flow, In J. Almond (ed.)
Proceedings of the 2nd International Symposium on the Theory of Road Traffic Flows, pp 113.
P.G. Gipps. 1981. A behavioural car-following model for computer simulation. Transportation
Research Part B: Methodological 15(2):105-111.
R.E. Chandler, R Herman and E.W. Montroll. Traffic dynamics: Studies in car following.
Operations Research 6(2):165-184.
D.C. Gazis, R Herman and R.W. Rothery. 1961. Nonlinear follow the leader models of traffic flow.
Operations Research 9:545-567.
Q. Yang and H.N. Koutsopoulos. 1996. A microscopic simulator for evaluation of dynamic traffic
management systems. Transportation Research Part C: Emerging Technologies 4(3):113-129.
N.E. Lownes and R.B. Machemehl. 2006. VISSIM: a multi-parameter sensitivity analysis, In L.F.
Perrone, B.G. Lawson, J. Liu and F.P. Wieland (eds.) Proceedings of the 38th Conference on
Winter Simulation (WSC'06), Portola Plaza Hotel, Monterey, CA, 03-06 December 2006, pp
1,406-1,413.
M.F. Tataro, M.G. Arada and J.P. Pabico. 2013. 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 (SFRC 2013), NCAS Auditorium, University of the
Philippines Los Baños, 13 December 2013, pp. 37-43.

9

