A Novel Approach in Evaluating The Impact of Vehic
A Novel Approach in Evaluating The Impact of Vehic
A Novel Approach in Evaluating The Impact of Vehic
models have been limited to examining the effect of which are categorized in this research as follows:
the roadway and environmental feature on a specif- Type 1 (skidding, skidding and overturning, skid-
ic accident type. These models usually use a wide ding & hitting objects), Type 2 (hitting objects in
range of logistic regression models [9, 10], assum- the carriageway), and Type 3 (run-off-road).
ing a single relationship between all of the explan- In order to find the indirect effect of the envi-
atory variables and the target variable. In contrast, ronmental conditions on AVC the roads have been
Haghighi and others have explored the nested re- divided into homogeneous segments. The segmen-
lationship between individual crash characteristics tation process included factors related to traffic
and environmental and roadway features [11]. He volumes, speed limits, roadside hazards, and the
has applied a multilevel ordinal logistic regression presence of horizontal curves, as recommended by
to address the hierarchical structure of accident data the Highway Safety Manual specifications [20] in
and its impact on accident severity outcome. order to obtain homogeneous road segments with
This research applies a multilevel multinomial distinguished environmental characteristics. The
logistic regression (MML) model developing fur- segmentation process has resulted in 576 segments
ther the methodology of Haghighi [11], to investi- with an average length of 5.3 kilometres, and aver-
gate the multilevel relationship between accident age accident content that equals 12.3 accidents per
types (ex. skidding, overturning, etc.) and roadway segment per year.
geometric and environmental conditions (ex. road- The AVC has been divided into five age groups:
side hazard, horizontal curve, traffic volume, speed 0-2, 3-5, 6-8, 9-11, and >11 years that are repre-
limit, etc.) variables and its impact on AVC, assum- sented in the following sections by AVC1, AVC2,
ing that the age of the vehicle strongly affects the AVC3, AVC4, and AVC5, respectively. The analysis
level of in-built safety technology of the vehicle and of AVC attribute treats every vehicle individually as
the adapted improvements in its safety designs [12]. a single input object being involved at least by one
injury or fatal accident.
2. DATA
The used dataset has been obtained from the UK 3. METHODOLOGY
Department of Transport website and has included Vehicle age, as a target variable, is conditioned
two types of data. The first dataset ([13-15]) has on the fact that an accident has happened. The out-
included accident information (i.e. location, date, comes of the AVC have been divided equidistant-
time, number of vehicles involved, road character- ly into five age groups (i.e. AVC1, AVC2, AVC3,
istics, presence of junctions, weather conditions), AVC4, and AVC5). Accident frequencies between
crashed vehicle characteristics (i.e. type, model, these randomly divided groups may vary for the
age, engine capacity, type of accident), [16] and given time period. Multinomial logistic regression
description of casualties (i.e. severity, age, gender, is an appropriate analytical approach in such situa-
type). The second dataset ([17-19]) has included tions when there is an outcome with more than two
traffic volume characteristics and speed limits for events that are not ordered, or where the ordered na-
different road segments. Additional roadway char- ture is not so clear. The analysis of the data showed
acteristics (i.e. horizontal curves, roadside hazard, a hidden influence of a hierarchical relationship be-
lane characteristics) have been identified with the tween explanatory variables. A two-level approach
help of the Google Earth and ArcGIS software. has been used to examine the nested relationship.
Based on the available information, homogeneous The description of the two-level model is given in
road segments have been generated during the first the following.
step of the analysis.
This study has concentrated on examining the Level-one modelling
impact of the nested relationship between various Assuming the data have a multinomial distribu-
environmental characteristics (i.e. roadway geomet- tion, the total outcome categories can be denoted as
ric characteristics, traffic volume, and speed limits) C, and each individual category can be indexed by
and accident types and their impact on AVC as a c. The probability of being in c-th category P(Y=c)
target variable applying the data of the UK motor- is πc(c=1,2,..,C). The probability of the outcome of
way roads between the years 2014-2016. The data each of the other category πc is compared against the
includes three variables describing accident types, reference category πC using the expected cumulative
log link function ηc. The general single-level model Table 1 – Two-level hierarchical structure between AVC and
to link the expected values of the outcome to the the explanatory variables
predicted values of ηc can be written as follows: Analysed object Evaluated attributes
Q
- Traffic volume
hc = log a r Cc k b 0^ c h +
r
/ X q b q^ c h (1) Roadway geomet- - Speed limit
q=1
Level-two rics and environ- - Horizontal curve
where: ηc is the log odd variable of being in a partic- mental features - Roadside hazard
ular category (c) versus the reference category C; β0 - etc.
is the intercept; Xq is the applied prediction variable - Crash type (shape)
Characteristics of
of predictor q; βq is the coefficient of predictor q. Level-one vehicle involved in
- Engine capacity
- Type of vehicle
The odd ratio can result from exponentiating the a single accident
- etc.
log odd coefficients (exp(ηc)). When the odd ratio
Vehicle age at the
is less than 1, the function value of the outcome Target
time of accident
Vehicle age at the time
becomes less than 0, while in case of an odd ratio variable of accident (AVC).
(AVC)
above 1, the function value becomes positive.
group j. Assuming there are no level-two predictors,
Level-two modelling (Developing multinomial the general form of the level-two model will be as
modelling) follows.
In the developed model, the multilevel model
treats the individual vehicle crashes as a member of b 0j^ c h = c 0j^ c h + u 0j^ c h (2)
a certain group. Each group corresponds to a specif- At level-two, more level-one intercepts (β0j(c))
ic road segment, including different environmental
can be modelled as a function of level-two variance
and roadway geometric characteristics. Since the in-
(uqj) and intercept (γ0) for different group j. The gen-
dividual crashes of a given group are likely to share
eral MML model is formed by combining both lev-
similar environmental characteristics, their causing
factors are assumed to be similar, more likely to els (Equations 1 and 2) in a single formula (Equation 3).
each other than to the causing factors of other indi- h ij^ c h = log b r ^ h l + b 0j^ c h
r ij^ c h
(3)
viduals from another group. Accordingly, the data ij C
insignificant variables. Tables 2 and 3 present the re- average accidents (in case of “off-road environ-
sulting random and fixed effect coefficients of the ment”, the effect of accident type on age is 0.7). If
model. the type of accident is hitting objects, then it is ex-
pected to include a newer vehicle (in case of “hit-
Table 2 – The random effect variance components of the CMM
model ting objects”, the effect of accident type on age is
-0.23). It needs to be mentioned that this is not sig-
Parameter Estimate Wald Z Sig. nificant at 95% confidence level (as the Sig. is 0.193
Residual 20.581 41.333 0.000 in Table 3). This finding is consistent with previous
studies. Boodlal and colleagues have suggested that
Intercept Variance
(Cluster ID)
0.086 1.025 0.031 crash occurrence and crash type can be partly ex-
plained by road geometric features, together with
The random effect coefficients represent the lev- vehicle characteristics [23]. Haghighi, Liu, Zhang,
el-two variance components of the intercepts. As & Porter have proved that geometric features can be
mentioned before, the inter-class correlation (ICC) applied to predict typical crash types [11].
describes the proportion of variance between groups. However, this outcome has to be interpreted very
In other words, it describes the characteristics of the carefully, since the risk of a vehicle characterized
nested relationship between the level-two explan- by a certain age group being involved in an acci-
atory variables (i.e. environmental conditions) and dent is also strongly influenced by the penetration of
the target variable (i.e. AVC). According to Table 2, the given age group compared to the whole vehicle
the residual parameter describes the variance due to fleet, and the traffic performance characterizing the
differences among individuals (i.e. AVC for each specific age group [24]. According to the penetra-
vehicle) within their respective units. In accordance tion-based comparison the oldest and the youngest
with the table, there is a significant variance with- age groups are significantly under-represented in
in the groups (Wald Z=41.33, p<0.001). Similarly, case of the investigated accident types [25, 26]. The
the intercept parameter indicates that the intercepts relative under-representation of newer cars in case
vary significantly across the different clusters (Wald of the analysed crash types seems to be rationally
Z=1.03, p=0.031). The Wald Z test provides a Z acceptable, since the new safety technologies can
statistic summarizing the ratio of the estimate to its strongly support the improvement of road safety
standard error. In other words, the developed CMM [27-29].
model can explain the hierarchical relation between Another interesting finding is the relationship
the hierarchical predictors and AVC. between the driver’s age and the vehicle age at the
The CMM has been built based on four predic- time of accident. The inverse relationship in Table 3
tors (i.e. driver age, accident type, vehicle type, (-0.11) indicates that older drivers are more likely to
and car-makers), all of which are statistically sig- be involved in new vehicles’ accidents while young
nificant. The estimated values of the predictors are drivers are more likely to be involved in old vehi-
presented in Table 3. cles’ accidents. It looks like the age compatibility
Table 3 provides the estimates of the fixed-effect between the car and its driver is more favourable to
coefficients. First, it can be seen that the intercept avoid accidents. With respect to vehicle type, newer
(β0=10.37) represents the average AVC when all trucks and buses are more likely to be included in
the predictor values are zero or at their reference accidents in comparison with cars. Regarding bus-
value. In the case of Accident type, two of which es, it has to be mentioned that it is not statistically
were found to be significantly related to the AVC. significant at a 95% confidence level (as Sig.=0.220
In the case of the first accident type (skidding, and in Table 3). This can be strongly influenced by the
skidding and overturning), accident risk is propor- commercial nature of these categories since com-
tional to AVC with positive estimated values (0.46 mercially operated vehicles are usually younger
and 0.35, respectively). In other words, accidents than the individually owned cars.
of older vehicles are more likely to occur due to Table 3 also estimated (with overall significant
skidding, and skidding and overturning. Similarly, value) the AVC for 44 different car makers but with
in case of the second accident type, if an accident different models. Among those, only 16 vehicles of
occurred in an off-road environment, the age of the same makers are statistically significant for the
vehicle is expected to be older than in the case of 95% confidence level. Figure 1 gives an indication of
0
-1
-2
Estimated value
-3
-4
-5
-6
-7
VAUXHALL
SEAT
FIAT
KIA
SKODA
BMW
ISUZU
MERCEDES-BENZ
TRIUMPH
HYUNDAI
AUDI
VOLVO
NISSAN
MERCEDES
CITROEN
FORD
the distribution of the model estimated value with target variable, considering it outstandingly relevant
respect to car-makers. The first to be noted from from the viewpoint of the accidents, since vehicle
Figure 1 is that all of the estimated values are neg- age actually represents the difference in the integrat-
ative, indicating an inverse relationship between ed safety systems in case of different vehicle gener-
the AVC and car-maker. In other words, newer cars ations. Accordingly, multilevel analysis in our case
have a higher probability of an accident in case all targets the investigation of the relationship between
the analysed car-makers. the different groups of crashed vehicles by identify-
In more detail, among the car-makers of Table 3, ing a hierarchical data structure that can contribute
the new models of Ford and Citroen have the least to using the advantages of the clustered dataset. This
predicted accident probability in comparison with multilevel approach is represented by a hierarchical
their old models (-1.78, -1.99, respectively). On the structure aimed to reveal the influence of group-level
contrary, the accident probability for newer models characteristics of data. A multinomial logistic regres-
from Isuzu and Mercedes-Benz were higher than sion model has been applied to compare the different
the older models that gave them higher negative categories of AVC outcomes using the cumulative
predicted values (-5.95, -5.08, respectively) in the odd log link function.
model. These probability values are influenced on Before the summary of the final results, it would
the one hand by the changes in the penetration and be good to introduce the limitation of the study brief-
on the other hand by the safety improvement efforts ly. It is important to emphasize that in most countries,
of the given producer. For example, if a producer’s a comprehensive statistical dataset is not available,
vehicles are more attractive in the present than in the which would be capable to describe the traffic per-
past, it causes an increased number of vehicles on formance of the vehicle types, categories, and age
the road from the certain producer, and thus, it can groups. It means that we cannot differentiate the cat-
increase accident probability. Besides if a producer egories based on the driven distances, which would
efficiently allocates higher resources on the safety be crucial to evaluate the accident risks and proba-
field, it can lead to safer cars and thus smaller recent bilities related to a specific category (e.g. age groups
accident probability values compared to the past. or vehicle types). In this case, if we characterize the
Due to the introduced factors and the relatively low effect of vehicle age on traffic safety without reliable
number of data used in the analysis for the different performance data, a relevant part of the methodology
car-makers the outcome has to be handled carefully. has to be based on assumptions and estimations. On
Furthermore, the variation in the odd ratio of the the other hand, we also need to mention that a data-
crash types between different AVC groups has proved base on the level of automation (especially focusing
the nested relationship between crash types and the on ADAS) related to specific vehicle categories and
environmental conditions related to different road age groups would strongly support the analysis of the
segments, which also support the applicability of the effect of automation on traffic safety.
multilevel approach. According to the results of the However, the model parameters of the current
analysis, older vehicles are more likely to be involved investigation have been estimated using road acci-
in traditional types of accidents. On the other hand,
dent data of the UK motorway network in the period
considering the penetration-based comparison, the
of 2014-2016. The first level of the model includes
oldest and the youngest age groups are significantly
parameters related to accident characteristics, while
under-represented in the case of the investigated ac-
the second level describes the environment and geo-
cident types. For instance, skidding type of accident
metric features of road segments where accidents oc-
has a higher probability of including casualties for
curred. The result has proven a significant influence
older vehicles compared to newer vehicles.
of the AVC by applying a hierarchical structure of
the considered explanatory variables. The values of
5. CONCLUSION the variance component have supported the assump-
In this paper, a multilevel multinomial logistic tion that the intercepts vary across the different AVC
regression has been applied to analyse the nested groups describing the correlation between the envi-
relationship between crash types, roadway geomet- ronmental features and different AVC groups.
ric, environmental conditions, and their impact on Finally, I would like to introduce the applicability
the change in the age of vehicles at the time of the of the results to other countries briefly. In this case,
accident (AVC). AVC has been investigated as the it seems to be reasonable to handle the issues of the
[9] Al-Ghamdi AS. Using logistic regression to estimate traffic-counts [Accessed 1 Aug. 2019].
the influence of accident factors on accident severi- [20] American Association of State Highway and Transporta-
ty. Accid. Anal. Prev. 2002;34(6): 729-741. Available tion Officials. Highway Safety Manual. 1st Edition; 2010.
from: https://www.sciencedirect.com/science/article/pii/ [21] Azen R, Walker CM. Categorical data analysis for the
S0001457501000732 [Accessed 3 Nov. 2019]. behavioral and social sciences; 2011.
[10] Yau KKW. Risk factors affecting the severity of single [22] Barcikowski RS. Statistical power with the group mean
vehicle traffic accidents in Hong Kong. Accid. Anal. as the unit of analysis. J. Educ. Behav. 1981;6(3):
Prev. 2004;36(3): 333-340. Available from: https://www. 267-285. Available from: https://www.jstor.org/sta-
ncbi.nlm.nih.gov/pubmed/15003577 [Accessed 3 Nov. ble/1164877?seq=1#metadata_info_tab_contents [Ac-
2019]. cessed 23 Nov. 2019].
[11] Haghighi N, Liu XC, Zhang G, Porter RJ. Impact of [23] Boodlal L, Donnell ET, Porter RJ, Garimella D, Le T,
roadway geometric features on crash severity on ru- Croshaw K, Himes C, Kulis P, Wood J. Factors Influenc-
ral two-lane highways. Accid. Anal. Prev. 2018;111: ing Operating Speeds and Safety on Rural and Suburban
34-42. Available from: https://www.ncbi.nlm.nih.gov/ Roads. Turner-Fairbank Highway Research Center. No.
pubmed/29169103 [Accessed 3 Nov. 2019]. FHWA-HRT-15-030, 2015.
[12] Greene D, Hossain A, Hofmann J, Helfand G, Beach [24] EurosStat. Available from: https://ec.europa.eu/eurostat/
R. Consumer willingness to pay for vehicle attributes: web/products-datasets/product?code=road_eqs_carage
What do we Know?. Transportation Research Part 2019.04.05. [Accessed 23 Nov. 2019].
A: Policy and Practice. 2018;118: 258-279. Available [25] Bair ST, Huang RJ, Wang KC. Can vehicle maintenance
from: https://www.sciencedirect.com/science/article/pii/ records predict automobile accidents?. Journal of Risk
S0965856417308546 [Accessed 3 Nov. 2019]. and Insurance. 2012;79(2): 567-584. Available from:
[13] Department for Transport. Road Safety Data 2014. https://onlinelibrary.wiley.com/doi/full/10.1111/j.1539-
[Online]. Available from: https://data.gov.uk/dataset/ 6975.2011.01433.x [Accessed 23 Nov. 2019].
cb7ae6f0-4be6-4935-9277-47e5ce24a11f/road-safe- [26] Knobloch K, Wagner S, Haasper C, Probst C, Krettek C,
ty-data [Accessed 1 Aug. 2019]. Otte D, Richter M. Sternal fractures occur most often in
[14] Department for Transport. Road Safety Data 2015. old cars to seat-belted drivers without any airbag often
[Online]. Available from: https://data.gov.uk/dataset/ with concomitant spinal injuries: Clinical findings and
cb7ae6f0-4be6-4935-9277-47e5ce24a11f/road-safe- technical collision variables among 42,055 crash vic-
ty-data [Accessed 1 Aug. 2019]. tims. The Annals of Thoracic Surgery. 2006;82(2): 444-
[15] Department for Transport. Road Safety Data 2016. 450. https://www.sciencedirect.com/science/article/pii/
[Online]. Available from: https://data.gov.uk/dataset/ S0003497506005935 [Accessed 23 Nov. 2019].
cb7ae6f0-4be6-4935-9277-47e5ce24a11f/road-safe- [27] Szalay Z, Tettamanti T, Esztergár-Kiss D, Varga I, Bar-
ty-data [Accessed 1 Aug. 2019]. tolini C. Development of a test track for driverless cars:
[16] Zöldy M, Zsombók I. Modelling fuel consumption and Vehicle design, track configuration, and liability con-
refuelling of autonomous vehicles. MATEC Web of Con- siderations. Periodica Polytechnica Transportation En-
ferences. 2018;235: 37. EDP Sciences. Available from: gineering. 2018;46(1): 29-35. Available from: https://
https://www.matec-conferences.org/articles/matecconf/ pp.bme.hu/tr/article/view/10753 [Accessed 23 Nov.
abs/2018/94/matecconf_hort2018_00037/matecconf_ 2019].
hort2018_00037.html [Accessed 1 Aug. 2019]. [28] Hulse LM, Xie H, Galea ER. Perceptions of autono-
[17] Department for Transport-b. GB Road Traffic Counts mous vehicles: Relationships with road users, risk, gen-
2014. [Online]. Available from: https://data.gov.uk/data- der and age. Safety Science. 2018;102: 1-13. Available
set/208c0e7b-353f-4e2d-8b7a-1a7118467acc/gb-road- from: https://www.sciencedirect.com/science/article/pii/
traffic-counts [Accessed 1 Aug. 2019]. S0925753517306999 [Accessed 23 Nov. 2019].
[18] Department for Transport-b. GB Road Traffic Counts [29] Zöldy M. Legal Barriers of Utilization of Autonomous
2015. [Online]. Available from: https://data.gov.uk/data- Vehicles as Part of Green Mobility. In: Burnete N, Varga
set/208c0e7b-353f-4e2d-8b7a-1a7118467acc/gb-road- B. (eds) Proceedings of the 4th International Congress of
traffic-counts [Accessed 1 Aug. 2019]. Automotive and Transport Engineering (AMMA 2018).
[19] Department for Transport-b. GB Road Traffic Counts Proceedings in Automotive Engineering. Springer, Cham;
2016. [Online]. Available from: https://data.gov.uk/data- 2019. Available from: https://doi.org/10.1007/978-3-
set/208c0e7b-353f-4e2d-8b7a-1a7118467acc/gb-road- 319-94409-8_29 [Accessed 28 Oct. 2019].