Investigation of the Hourly and Spatial Patterns of Traffic Offenses During March-April 2019 in Iran Using Bivariate Generalized Additive Models and Integrated Nested Laplace Approximation


avatar Mohammad Fayaz ORCID 1 , * , avatar Alireza Abadi ORCID 2 , 3 , avatar Alireza Razzaghi ORCID 4 , avatar Soheila Khodakarim ORCID 5 , avatar Mostafa Hosseini ORCID 6 , **

Department of Biostatistics, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran
Department of Community Medicine, Faculty of Medicine, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran
Social Determinants of Health Research Center, Shahid Beheshti University of Medical Sciences, Tehran, IR Iran
Road Traffic Injury Research Center, Tabriz University of Medical Sciences, Tabriz, IR Iran
Department of Biostatistics, School of Medicine, Shiraz University of Medical Sciences, Shiraz, IR Iran
Department of Epidemiology and Biostatistics, School of Public Health, Tehran University of Medical Sciences, Tehran, IR Iran
Corresponding Authors:

how to cite: Fayaz M, Abadi A, Razzaghi A, Khodakarim S, Hosseini M. Investigation of the Hourly and Spatial Patterns of Traffic Offenses During March-April 2019 in Iran Using Bivariate Generalized Additive Models and Integrated Nested Laplace Approximation. Int J High Risk Behav Addict. 2022;11(3):e118376.



The control, management, and prevention of driving accidents and risky driving are regarded as concerns for numerous countries, according to the World Health Organization. In this regard, many technologies, such as count stations, are recommended. They count traffic offenses, such as speeding and unsafe distance, hourly and daily, and have different patterns according to the hour of the day and the location.


This study aimed to investigate the risky driving behaviors according to traffic offenses in Iran and estimate their hourly and spatial patterns using generalized additive models (GAMs) and stochastic partial differential equation methods.

Materials and Methods:

There were 2,316 count data stations for one month within March-April 2019. This study estimated the hourly average of each traffic offense, Pearson’s and Spearman’s correlations, and the energy statistics for testing the bivariate normal distribution. There are five distributions, such as univariate Poisson, quasi-likelihood Poisson, Gaussian, location-scale Gaussian, and bivariate Gaussian in GAMs, to study the hourly patterns which were compared to the mean squared error (MSE) and correlation.


The hourly average of total vehicles and number of speeding and unsafe distance offenses per count station had positive skew distributions with mean values equal to 347 ± 456, 22.5 ± 44.2, and 65.9 ± 150, respectively. The correlation between traffic offenses in most provinces was significant, not large, and different. The GAM with the bivariate Gaussian distribution had the best performance according to the MSE and correlation. It revealed three hourly patterns for count predictions; the first was that speeding is higher than unsafe distances; the second was that unsafe distances are higher than speeding; the third was that speeding and unsafe distances do not have a specific pattern in some hours. The percentage of speeding was higher in the central, northeast, and southeast regions than in other parts of Iran, and the percentage of unsafe distances was higher for the north, northwest, west, and some parts of the southwest than in other parts of Iran, respectively.


The hourly pattern of traffic offenses exists and has a complex structure. The spatial pattern of traffic offenses shows the riskiest points in Iran.

1. Background

The traffic and speed cameras and count stations near the roads are some technologies developed to manage, control, and predict traffic status in different countries. These devices produce massive datasets hourly and daily, making them one of the primary resources for discovering patterns and relationships. For example, there are diverse indices, including the count of driving offenses and total vehicles based on their type, which could be considered (1, 2). According to the global status report on road safety by the World Health Organization, driving accidents and risky driving are among the remarkable causes of death worldwide (2).

This study investigated the hourly and spatial patterns of driving offenses using two advanced statistical methods, respectively. Firstly, the generalized additive model (GAM) extended the generalized linear model (GLM) idea by adding the smoothing functions, such as cubic regression splines for estimating, hypothesis testing, and confidence interval for coefficients. The GAM was enhanced to capture complex and nonlinear relationships between response variables and covariates (3-5). Secondly, the stochastic partial differential equation (SPDE) method is one of the techniques to study spatial variability with the integrated nested Laplace approximation (R-INLA) (6).

2. Objectives

The traffic offenses dataset from count stations is a source for studying driving behaviors in Iran. The speeding and unsafe-distance indices are the indicators of risky driving for the control and prevention for which governments worldwide have different laws (2). Firstly, this study estimated the hourly average of each traffic offense per count station national-wide and province-wise. It shows how many traffic offenses occurred on average in each count station hourly. Secondly, the correlations and bivariate distribution between speeding and unsafe distance province-wise offenses were evaluated by Pearson’s and Spearman’s correlation coefficients and bivariate normal test based on the energy statistics, respectively. Thirdly, this study modeled the nonlinear relationship between the number of traffic offenses and the time of a day for each province with GAMs separately. Finally, this study introduced the percentage of traffic offenses among all transported vehicles as a key index to study the risky behavior of drivers on roads and investigated their spatial variability with SPDE. Moreover, temporal variability for four different ranges of hours was evaluated.

3. Materials and Methods

3.1. Population and Dataset

The population of this study was all count data stations near the interprovincial roads of Iran, which are available online from Iran Road Maintenance and Transportation Organization on the website of the Ministry of Roads and Urban Development ( All days of the first Iranian month, Farvardin, were considered in this analysis (31 days within 3/21/2019 to 4/20/2019). The number of individuals traveling between provinces increases in the first and second weeks of this month due to the Iranian New Year holidays. The condition of the roads gets back to regularly slow during the third and fourth weeks of this month. Therefore, the risky driving behavior of most Iranians could be estimated in this month.

The count data stations record hourly and daily different indices, such as total vehicles and the number of speeding and unsafe distance offenses. There were 31 separated provincial datasets. The active count stations record 60 minutes an hour, and only 2,316 count data stations remain. The count of unsafe distance offenses is the total number of vehicles with a distance shorter than 2 seconds between them. The count of speeding is the total number of vehicles on the road with a speed higher than the speed control limit (

3.2. Statistical Analysis

3.2.1. Testing the Multivariate Normal Distribution

In order to study the bivariate distribution of these two traffic offenses, this study chose the energy test statistics that have the best performance among multivariate Gaussian distribution tests (7, 8). The observation Z ∈ Rd has a multivariate Gaussian distribution Nd (0, Id) with the mean vector of 0 and the variance-covariance matrix Id (9) as follows:

Equation 1.EZ-Z'd=2EZd=2Γd+12Γd2

The y1,…, yn are the standardized elements of the sample. The energy test statistics for the d-distribution standard normal are as follows:

Equation 2.nεn,d=nn2j=1nEyi-Zd-2Γd+12Γd2-1n2j,k=1nyj-ykd


Equation 3.Ea-Zd=2Γd+12Γd2+2πk=0-1kk!2kad2k+22k+12k+2Γd+12Γk+32Γk+d2+1

For computation, the sample of observation is standardized based on the mean and correlation matrix. The obtained test statistic nε ^n,d has limited distributions, such as nεn,d, with the rejection region on the upper tail. The nε ^n,d statistics are compared to the repeated energy statistics by the Monte Carlo sampling and 200 iterations from standardized normal with the same dimensions.

3.2.2. Correlation Tests

The correlation between the two random variables Y1 and Y2 or ρ12 was assessed with the maximum likelihood estimation of the Pearson product-moment correlation coefficient if their distribution was bivariate normal. When the joint distribution of the two random variables was not multivariate normal, the correlation was evaluated utilizing Spearman’s Rank correlation coefficient. In this test, the ranks of Y1 and Y2 are R1 and R2 for calculation, respectively (10).

3.2.3. Generalized Additive Model

The GAM is a GLM with a set of smoothing functions of covariates. The general form is as follows:

Equation 4.gμi=Aiθ+f1x1i+f2x2i+

where μi ≡ E (Yi) and Yi ~ EF (μi, ϕ). The response term is Yi with the exponential family, a mean of μi, and the scale parameters of ϕ. The Ai is related to the parametric part of the model, and θ represents a vector related to the parameters. The fi is the smoothing function of the covariate xk (4). This study considered and compared the five distributions of univariate Poisson, univariate quasi-Poisson, univariate Gaussian, univariate location-scale Gaussian, and bivariate Gaussian for Yi. There are two models for each distribution, with the response variable in models 1 and 2 being the number of speeding and unsafe distance offenses, respectively. The only covariate in the model is the hour of the day, and the smoothing function f1(x) is a cubic regression spline with k = 20. The mean squared error (MSE) and correlation between observations and fitted values are the criteria for choosing the best model (4, 5).

3.2.4. Spatial Data Analysis

The response variable was the percentage of total speeding and unsafe distances in the total traffic during one month for each count station, respectively. The count station location has some information that accounts for this model. The Yi refers to the percentage of each traffic offense at count station si and has a Gaussian distribution as follows:

Equation 5.Yi~ Nμi,σ2,i=1,,2316
Equation 6.μi=β0+Zsi

where Z(.) indicates a spatially structured random effect with zero-mean Gaussian process and Matérn covariance. The fitting method was the SPDE approach with R-INLA. The mesh was constructed with different margins to approximate Yi to a discrete Gaussian Markov random field. The models predict the lower and upper limits with 95% credible intervals, and the result is a map plot (6). All the statistical analyses were performed using R 4.0.2 ( with mgcv, energy, geoR, maps, rgdal, maptools, sf, viridis, and ggplot2 packages (11-18).

4. Results

Table 1 shows the results of descriptive statistics and p-values for the bivariate normal test based on the energy statistics and correlation coefficients. Diverse factors were associated with driving offenses, among which only the time of day as a spline function was considered in this study. Table 2 shows the results of fitting 1) univariate GAM with quasi-Poisson distribution and 2) bivariate Gaussian GAM (other models, including univariate Poisson, univariate Gaussian, and univariate location-scale Gaussian, in Appendix 1). These two models had the best performance among their Poisson-related and Gaussian-related families according to the MSE and correlation indices.

Table 1.

Descriptive Statistics a

ProvincesNumber of Count StationsHourly Average per Station (SD)P-Value of Bivariate Normal Test bCorrelation
SpeedingUnsafe DistancePearsonSpearman
Ardebil6412.2 (12.3)33.3 (29.2)< 0.050.16 (< 0.05)0.29 (< 0.05)
Isfahan12026.5 (35.4)30.3 (50.2)> 0.050.26 (< 0.05)0.50 (< 0.05)
Alborz4099.4 (197.6)320.8 (555.8)< 0.050.32 (< 0.05)0.49 (< 0.05)
Ilam485.5 (5.7)14.5 (17.4)> 0.050.21 (< 0.05)0.46 (< 0.05)
West Azarbayjan9017.9 (19.2)40 (28.6)< 0.050.25 (< 0.05)0.46 (< 0.05)
East Azarbayjan8215.4 (23.5)18.2 (63.6)< 0.05-0.01 (0.11)0.08 (< 0.05)
Bushehr5216.9 (14.8)53.5 (68)< 0.050.26 (< 0.05)0.30 (< 0.05)
Tehran10349.5 (79.3)259.6 (334.7)< 0.050.27 (< 0.05)0.41 (< 0.05)
Chaharmahal and Bakhtiari6223.7 (24.8)26 (21.3)> 0.050.46 (< 0.05)0.49 (< 0.05)
South Khorasan7220.5 (16.4)13.9 (12.7)> 0.050.51 (< 0.05)0.61 (< 0.05)
Razavi Khorasan13119.6 (31.5)70 (133.9)< 0.050.30 (< 0.05)0.48 (< 0.05)
North Khorasan4712.6 (17.4)32.5 (27.5)> 0.050.33 (< 0.05)0.26 (< 0.05)
Khuzestan10414.5 (17.8)34.7 (31.2)> 0.050.16 (< 0.05)0.41 (< 0.05)
Zanjan747.3 (8.7)40.5 (52.5)> 0.050.38 (< 0.05)0.48 (< 0.05)
Semnan5443.9 (49.3)33.1 (36.2)> 0.050.58 (< 0.05)0.69 (< 0.05)
Sistan and Balouchestan7026.2 (24)1.2 (1.3)> 0.050.50 (< 0.05)0.44 (< 0.05)
Fars12633.7 (46.1)74 (110.9)> 0.050.21 (< 0.05)0.35 (< 0.05)
Qazvin5714.8 (25.1)171.5 (268.2)> 0.050.25 (< 0.05)0.50 (< 0.05)
Qom7649 (74)82.4 (124.3)> 0.050.16 (< 0.05)0.20 (< 0.05)
Kordestan527.6 (9.9)26.9 (28)< 0.050.30 (< 0.05)0.37 (< 0.05)
Kerman10713.4 (13.4)19.6 (17.2)< 0.050.27 (< 0.05)0.33 (< 0.05)
Kermanshah729.8 (12.1)42.2 (47.9)< 0.050.32 (< 0.05)0.52 (< 0.05)
Kohkilouye and Boyerahmad287.1 (8.5)24.6 (19.7)< 0.050.34 (< 0.05)0.50 (< 0.05)
Golestan5213.1 (24.6)78.5 (109.9)< 0.050.00 (0.65)0.29 (< 0.05)
Gilan849.8 (15)193 (179.2)< 0.050.09 (< 0.05)0.30 (< 0.05)
Lorestan6211.7 (12.5)28.7 (31.8)< 0.050.38 (< 0.05)0.49 (< 0.05)
Mazandaran10223.7 (34.7)111.6 (164)> 0.050.13 (< 0.05)0.32 (< 0.05)
Markazi6430.6 (29.7)56.7 (64.2)< 0.050.36 (< 0.05)0.47 (< 0.05)
Hormozgan7227.5 (31.7)37.4 (37.9)< 0.050.49 (< 0.05)0.48 (< 0.05)
Hamedan7115.3 (20.3)44.4 (74.9)> 0.050.47 (< 0.05)0.69 (< 0.05)
Yazd7819.6 (32.6)19.8 (40.8)< 0.050.58 (< 0.05)0.53 (< 0.05)
Table 2.

Mean Squared Error and Correlation Between Observed Values and Predicted Values by Provinces in Three Generalized Additive Models

ProvincesUnivariate Quasi-PoissonBivariate Gaussian
SpeedingUnsafe DistanceSpeedingUnsafe DistanceSpeedingUnsafe DistanceSpeedingUnsafe Distance
West Azarbayjan1083.112514.080.230.31544.51712434.6030.25640.2809
East Azarbayjan1943.8115195.810.250.26839.9921988.4030.24960.3423
Chaharmahal and Bakhtiari1664.662716.230.30.411105.0881737.0310.31120.4543
South Khorasan1319.271133.260.330.32834.305873.050.33860.3355
Razavi Khorasan2241.1250497.250.180.221682.41640515.6890.17970.2487
North Khorasan936.375419.720.160.38767.0993592.4780.16470.4066
Sistan and Balouchestan2478.8612.470.340.31363.40110.2180.35280.3227
Kohkilouye and Boyerahmad339.825697.060.250.4238.1773585.0470.26140.4133

Figure 1 shows the predicted responses for the bivariate Gaussian model in Sistan and Baluchestan, and Gilan provinces, Iran. The patterns of these two plots were in contradiction. All figures for both univariate quasi-Poisson and bivariate Gaussian GAM models are presented in the appendix (Appendix 2). Figure 2 illustrates the predicted percentage of speeding and unsafe distance offenses with their spatial patterns in Iran. The two plots had a different pattern that indicated the percentage of speeding is higher in the central, northeast, and southeast regions than in other parts of Iran. The predicted percentage of speeding was higher than unsafe distance, indicating that speeding is high in Iran. The predicted percentage of unsafe distance with a range of 20 - 40% was higher for north, northwest, west, and some parts of southwest than other parts of Iran. These plots were produced with the triangulated mesh method. The other mesh method is available for comparison in Appendix 3.

Estimated number of speeding (blue) and unsafe distance (red) offenses in A, Gilan; and B, Sistan and Balouchestan, Iran
Estimated number of speeding (blue) and unsafe distance (red) offenses in A, Gilan; and B, Sistan and Balouchestan, Iran
Mean prediction of A, speeding; and B, unsafe distance rates (within 0 (yellow) and 1 (red))
Mean prediction of A, speeding; and B, unsafe distance rates (within 0 (yellow) and 1 (red))

The temporal and hourly patterns of the predicted percentage for traffic offenses are available in Appendix 4.1 and Appendix 4.2, respectively. There were four intervals of 6 hours, including 0 - 5, 6 - 11, 12 - 17, and 18 - 23. According to Appendix 4.1, the predicted percentage of speeding was the highest at 0 - 5. On the other hand, Appendix 4.2 shows that the predicted percentage of unsafe distance was the lowest at 0 - 5 and had the highest value at 12 - 17 (Appendix 4.1 and Appendix 4.2).

5. Discussion

The speed limit and speeding fines for intercity highways are different in European countries and Iran. As illustrated in Figure 3, Norway has the highest speeding fines (€711), and the Czech Republic has the lowest speeding fines (€19) among other countries, respectively. In some countries, the speeding limits are different based on the two lanes or other roads (Norway), winter and other seasons (Sweden and Finland), rainy and other weather conditions (Luxembourg and Italy), and free speed (Germany). The traffic fines for speeding in Iran within 2021 - 2022 is 2,100,000 Iranian Rial (IR) (equivalently about €44.30 if €1 = 47,377 IRR (i.e., official governmental NIMA, the Central Bank of the Islamic Republic of Iran exchange rate) and about €6.72 if €1 = 312,480 IRR (i.e., unofficial open market exchange rate) in Q1 2022). Nevertheless, this comparison of speeding fines is simple and naive according to the different purchasing power parity and gross domestic product between countries; therefore, it is suggested to compare them to new indices, such as Big Mac Index, in future studies (19). The Figure 3 dataset is made from and datasets.

Bar chart (green) of highway speed limit (km/h) and line chart (red) of speeding traffic fines (euro) by countries
Bar chart (green) of highway speed limit (km/h) and line chart (red) of speeding traffic fines (euro) by countries

A significant relationship has been reported between road traffic accidents (RTA) and time (e.g., the time of the day) in Yazd, Iran, during the New Year holidays and summer (20). It has been the other study for these holidays in the six most populous provinces of Iran, namely Fars, Khorasan Razavi, Tehran, Isfahan, Kerman, and Khuzestan, within 2011 - 2015, indicating that Fars and Khorasan Razavi, with attractive tourist sites, have different high RTA among others (21). The mortality rate due to traffic accidents is higher in Iran at midnight and summer (22, 23) and in spring and summer in Shiraz (24) than in other times. The present study showed that the speeding rate was higher from midnight to early morning throughout Iran. This finding might be due to the existence of black spots with low and not uniform lighting on the roads (25). It also suggests that risky driving, not darkness, is the main reason for accidents (26). However, there are some exceptions; for example, in the south of Iran (27) and Yasuj (28), the rate of accidents during the day is higher than at night. The provinces in the south of Iran are located in the warmest region in the country, with an average daily high temperature. The high rate of accidents during the day in these provinces might be due to heat stress on drivers. A similar finding was reported in a study in Saudi Arabia (29) and high ambient temperatures in Spain (30). The percentage of RTA in Fars, Isfahan, Ilam, and the southeast region of Iran has a nonlinear trend in 24 hours with different peaks of speeding and unsafe distance (26, 31-34).

The rate of unsafe distance offense is almost high in Kermanshah, Iran, and it can be added as a new risk factor of RDA in this province (35). The highest rate of mortality for drivers, passengers, and pedestrians has been reported during 13:00 - 18:00 in the west of Iran (36), 18:00 - 20:00 in Mashhad (37), and 16:00 - 18:00 in the southeast of Iran (38). The findings of other studies showed that most of the collisions occurred in the early hours of the night. A part of these collisions is due to poor visibility. Inadequate visibility has a key role in crashes involving pedestrians, motorcyclists, cyclists, and drivers (39). Moreover, crashes at dark hours cause severe injuries (40, 41). According to the findings of a meta-analysis study, the odds ratio of mortality in dark-hour crashes is 53% higher than in day-hour crashes (42). Consequently, traffic offenses might have a relationship with RTA-related mortality. In this regard, this study suggests adding traffic offenses statistics to the Iranian Integrated Road Traffic Injury Registry (43, 44) and RTA studies. The other risk factors are the spatial variations of traffic offenses and accidents, seatbelt and helmet status, gender, age group (45), and climatic conditions, such as fog in the north of Iran (46).

The traffic fines and risky driving in Iran are studied in different ways, including the relationship between the number of traffic offenses and fuel costs within 2011 - 2019 (47), the relationship between increasing traffic fines policies and the road traffic law enforcement (48), the prevalence and determination of speeding in Iran (49), the comparison study of traffic fines in Iran and other countries (50), risky driving fined by police in 2006 and 2007 in Tehran (51), the effect of cameras on speeding behavior of taxi drivers in two highways (52), and aggressive violations (e.g., “sound horn to indicate your annoyance”, “get angry, give chase”, and “aversion, indicating hostility”). Moreover, Iran and Great Britain, the Netherlands, and Finland are among the countries with higher speeding violations than other countries, such as Greek and Turkey (53). In addition, visual, perceptual, and cognitive capabilities and physiological condition of drivers (e.g., Barkley’s Attention Deficit Hyperactivity Disorder Screening Test, Risk Perception Questionnaire, Risk Taking Questionnaire, Sensation Seeking Scale Survey, and Driver Behavior Questionnaire), among other factors in SHRP2 naturalistic driving study, are assessed in the USA (54).

Advanced and sophisticated statistical methods are in demand for traffic-related datasets. The GAMs are among the statistical models that can be used for complex relationships, such as risky driving in Iran (4) (e.g., driving offenses near public places, such as airports (1)). The bivariate structure of the response can estimate the correlation and compare the traffic offenses between provinces at distinct times of the day (5). It also calculates the peak hours and 95% confidence interval with their pattern for each province. The Getis-Ord General G* statistic in geographically weighted regression models revealed that the hotspot for fatal pedestrian accidents is in Mazandaran, Iran, and it is more common in Yazd, East Azerbaijan, and Ardebil (55). Future studies can investigate clustering methods, statistical learning methods (56), functional data analysis (57), and GAM for location, scale, and shape techniques to estimate the exact distribution with many parameters and their estimation for the underlying distribution (58).

5.1. Conclusions

The present study concluded that the risky driving behaviors due to traffic offenses can be estimated straightforwardly at different times and locations and add new information about the time of the days and roads that have not registered or occurred any traffic accidents. In this regard, they are predictive models. The geographical status of the roads, such as mountains or deserts, is shown to be related to the type of traffic offenses. For example, speeding violations on desert roads are higher than mountain roads, and unsafe distance violations on mountain roads are higher than desert roads. The day-night, rush hours, and holidays are the main time-related factors for occurring traffic offenses. The future direction of this study is to investigate the relationships between the percentage of traffic offenses and traffic accident occurrence, climate status (e.g., raining, foggy, and sunny), and holidays and restrictions (e.g., coronavirus disease 2019 restriction) on all roads in different times of the day.

In highway safety research, crash modification factors and safety performance are introduced based on the traffic volume and road characteristics, and different statistical methods are proposed to estimate them (59, 60). Therefore, defining new and easy-to-compute indices is needed for future studies to measure and model the percentage of risky driving. This study had some limitations. Firstly, the police statistics have crime classification errors (e.g., some errors in detecting speeding and unsafe distance) and systematic errors (e.g., the failure of count stations in some hours) (61). Secondly, the statistics on traffic accidents are not publicly available.


  • 1.

    Fayaz M, Abadi AR, Khodakarim S, Hoseini MR, Razzaghi AR. The Data-Driven Pattern for Healthy Behaviors of Car Drivers Based on Daily Records of Traffic Count Data from 2018 to 2019 near Airports: A Functional Data Analysis. JP J Biostat. 2020;17(2):539-57.

  • 2.

    World Health Organization. Global status report on road safety 2018: summary. Geneva: World Health Organization; 2018, [cited 2022]. Available from:

  • 3.

    Hastie T, Friedman J, Tibshirani R. The Elements of Statistical Learning: Data Mining, Inference, and Prediction. New York, NY: Springer Series in Statistics; 2001.

  • 4.

    Wood SN. Generalized Additive Models: An Introduction with R. New York: CRC Press; 2017.

  • 5.

    Wood SN, Pya N, Säfken B. Smoothing Parameter and Model Selection for General Smooth Models. J Am Stat Assoc. 2016;111(516):1548-63.

  • 6.

    Moraga P. Geospatial Health Data: Modeling and Visualization with R-INLA and Shiny. New York: Chapman and Hall/CRC; 2019.

  • 7.

    Chen W, Genton MG. Are You All Normal? It Depends. Ithaca, NY: arXiv; 2020, [cited 2022]. Available from:

  • 8.

    Joenssen DW, Vogel J. A power study of goodness-of-fit tests for multivariate normality implemented in R. J Stat Comput Simul. 2012;84(5):1055-78.

  • 9.

    Székely GJ, Rizzo ML. A new test for multivariate normality. J Multivar Anal. 2005;93(1):58-80.

  • 10.

    Kutner MH. Applied linear statistical models. New York: McGraw-Hill Irwin; 2005.

  • 11.

    Rizzo ML, Szekely GJ, Rizzo MM. Package ‘energy’. 2021. Available from:

  • 12.

    Wood SN. mgcv: GAMs and generalized ridge regression for R. R News. 2001;1/2:20-5.

  • 13.

    Garnier S, Ross N, Rudis R, Camargo AP, Sciaini M, Cédric S. viridis - Colorblind-Friendly Color Maps for R. Newark, NJ: Sjmgarnier; 2021, [cited 2022]. Available from:

  • 14.

    Becker RA, Wilks AR, Ray Brownrigg R, Minka TP, Deckmyn A. maps: Draw Geographical Maps. Wien, Austria: The Comprehensive R Archive Network; 2018, [cited 2022]. Available from:

  • 15.

    Bivand R, Lewin-Koh N, Pebesma E, Archer E, Baddeley A, Bearman N, et al. maptools: Tools for Handling Spatial Objects. Wien, Austria: The Comprehensive R Archive Network; 2022, [cited 2022]. Available from:

  • 16.

    Ribeiro Jr PJ, Diggle P, Christensen O, Schlather M, Bivand R, Ripley B, et al. geoR: Analysis of Geostatistical Data. Wien, Austria: The Comprehensive R Archive Network; 2022, [cited 2022]. Available from:

  • 17.

    Bivand R, Keitt T, Rowlingson B, Pebesma E, Sumner M, Hijmans R, et al. rgdal: Bindings for the 'Geospatial' Data Abstraction Library. Wien, Austria: The Comprehensive R Archive Network; 2021, [cited 2022]. Available from:

  • 18.

    Wickham H. ggplot2: Elegant Graphics for Data Analysis. New York: Springer; 2009.

  • 19.

    Clements KW, Si J. Simplifying The Big Mac Index. J Int Financial Manag Account. 2017;28(1):86-99.

  • 20.

    Lotfi M, Montazeri M, Lashkardoost H, Shamsi F, Askari M, Hamedi E, et al. Road traffic accidents in Yazd province, Iran: A longitudinal study (2012–2016). Arch Trauma Res. 2018;7(2):68-72.

  • 21.

    Besharati MM, Azizi Bondarabadi M, Memariyan M, Tavakoli Kashani A. Patterns of road traffic fatalities in the six most populous provinces of Iran, 2011–2015. Arch Trauma Res. 2019;8(3):177-81.

  • 22.

    Sadeghi-Bazargani H, Ayubi E, Azami-Aghdash S, Abedi L, Zemestani A, Amanati L, et al. Epidemiological Patterns of Road Traffic Crashes During the Last Two Decades in Iran: A Review of the Literature from 1996 to 2014. Arch Trauma Res. 2016;5(3). e32985. [PubMed ID: 27800461]. [PubMed Central ID: PMC5078874].

  • 23.

    Shahbazi F, Soori H, Khodakarim S, Ghadirzadeh MR, Hashemi Nazari SS. Analysis of mortality rate of road traffic accidents and its trend in 11 years in Iran. Arch Trauma Res. 2019;8(1):17-22.

  • 24.

    Ghaem H, Hajipour M, Tababataee HR, Yadollahi M, Izanloo F. Time Series Analysis of Mortalities Resulting from Car Accidents in the Injured Individuals Hospitalized in Shiraz Shahid Rajaee Hospital During 2010 - 2016. Trauma Mon. 2018;23(1). e13573.

  • 25.

    Mohan A, Landge VS. Identification of Accident Black Spots on National Highway. Int J Civ Eng Technol. 2017;8(4):588-96.

  • 26.

    Mohammadi G. The pattern of fatalities by age, seat belt usage and time of day on road accidents. Int J Inj Contr Saf Promot. 2009;16(1):27-33. [PubMed ID: 19058047].

  • 27.

    Rakhshani T, Rakhshani F, Asadi ZS, Hadiabasi M, Khorramdel K, Zarenezhad M. Study of the pattern of mortality caused by Traffic Accidents (TAs) in The South of Iran. J Pak Med Assoc. 2016;66(6):644-9. [PubMed ID: 27339561].

  • 28.

    Rakhshani T, Kashfi M, Amirian I, Ebrahimi M, Hashemi Nazari S. Epidemiology of Fatal Road Traffic Accidents in Iran, Yasouj, 2014-2015. J Health Sci Surveill Syst. 2018;6(1):29-35.

  • 29.

    Nofal FH, Saeed AA. Seasonal variation and weather effects on road traffic accidents in Riyadh city. Public Health. 1997;111(1):51-5. [PubMed ID: 9033225].

  • 30.

    Basagana X, Escalera-Antezana JP, Dadvand P, Llatje O, Barrera-Gomez J, Cunillera J, et al. High Ambient Temperatures and Risk of Motor Vehicle Crashes in Catalonia, Spain (2000-2011): A Time-Series Analysis. Environ Health Perspect. 2015;123(12):1309-16. [PubMed ID: 26046727]. [PubMed Central ID: PMC4671248].

  • 31.

    Hasanzadeh J, Moradinazar M, Najafi F, Ahmadi-Jouybary T. Trends of Mortality of Road Traffic Accidents in Fars Province, Southern Iran, 2004 - 2010. Iran J Public Health. 2014;43(9):1259-65. [PubMed ID: 26175980]. [PubMed Central ID: PMC4500428].

  • 32.

    Mansouri M, Javad Kargar M. Analysis and Monitoring of the Traffic Suburban Road Accidents Using Data Mining Techniques; A Case Study of Isfahan Province in Iran. Open Transp J. 2014;8(1):39-49.

  • 33.

    Mohammadfam I, Karami Naserkhani R, Soltanian AR. The analysis of deaths caused by driving accidents in Ilam province, western Iran and the related factors by using the method of time series. Int J Occup Hyg. 2016;8(4):200-7.

  • 34.

    Khorshidi A, Ainy E, Hashemi Nazari SS, Soori H. Temporal Patterns of Road Traffic Injuries in Iran. Arch Trauma Res. 2016;5(2). e27894. [PubMed ID: 27703958]. [PubMed Central ID: PMC5037289].

  • 35.

    Izadi N, Khoram Dad M, Jamshidi P, Zanganeh AR, Shafiei J, Firouzi A. Epidemiological Pattern and Mortality Rate Trend of Road Traffic Injuries in Kermanshah Province (2009-2014). J Community Health Res. 2016;5(3):158-68.

  • 36.

    Hamzeh B, Najafi F, Karamimatin B, Ahmadijouybari T, Salari A, Moradinazar M. Epidemiology of traffic crash mortality in west of Iran in a 9 year period. Chin J Traumatol. 2016;19(2):70-4. [PubMed ID: 27140212]. [PubMed Central ID: PMC4897842].

  • 37.

    Sarbaz M, Kimiafar K, Khadem Rezaiyan M, Banaye Yazdipour AR. Epidemiology of transport accidents based on international statistical classification of diseases (ICD-10) in Mashhad, Iran. Int Electron J Med. 2018;7(1):23-9.

  • 38.

    Rad M, Martiniuk AL, Ansari-Moghaddam A, Mohammadi M, Rashedi F, Ghasemi A. The Pattern of Road Traffic Crashes in South East Iran. Glob J Health Sci. 2016;8(9):149-58. [PubMed ID: 27157159]. [PubMed Central ID: PMC5064071].

  • 39.

    World Health Organization. World report on road traffic injury prevention. Geneva: World Health Organization; 2004, [cited 2022]. Available from:

  • 40.

    Ackaah W, Apuseyine BA, Afukaar FK. Road traffic crashes at night-time: characteristics and risk factors. Int J Inj Contr Saf Promot. 2020;27(3):392-9. [PubMed ID: 32588731].

  • 41.

    Ramadani N, Zhjeqi V, Berisha M, Hoxha R, Begolli I, Salihu D, et al. Public Health Profile of Road Traffic Accidents in Kosovo 2010-2015. Open Access Maced J Med Sci. 2017;5(7):1036-41. [PubMed ID: 29362641]. [PubMed Central ID: PMC5771275].

  • 42.

    Yousefifard M, Toloui A, Ahmadzadeh K, Gubari MIM, Madani Neishaboori A, Amraei F, et al. Risk Factors for Road Traffic Injury-Related Mortality in Iran; a Systematic Review and Meta-Analysis. Arch Acad Emerg Med. 2021;9(1). e61. [PubMed ID: 34580659]. [PubMed Central ID: PMC8464012].

  • 43.

    Sadeghi-Bazargani H, Sadeghpour A, Lowery Wilson M, Ala A, Rahmani F. Developing a National Integrated Road Traffic Injury Registry System: A Conceptual Model for a Multidisciplinary Setting. J Multidiscip Healthc. 2020;13:983-96. [PubMed ID: 33061404]. [PubMed Central ID: PMC7520136].

  • 44.

    Marin S, Pourasghar F, Moghisi AR, Samadirad B, Haddadi M, Khorasani-Zavareh D, et al. Development and psychometric evaluation of data collection tools for Iranian integrated road traffic injury registry: Registrar-station data collection tool. Arch Trauma Res. 2019;8(3):170-6.

  • 45.

    Fathollahi S, Saeedi Moghaddam S, Rezaei N, Jafari A, Peykari N, Haghshenas R, et al. Prevalence of behavioural risk factors for road-traffic injuries among the Iranian population: findings from STEPs 2016. Int J Epidemiol. 2019;48(4):1187-96. [PubMed ID: 30843066].

  • 46.

    Khodadadi-Hassankiadeh N, Rad EH, Koohestani HS, Kouchakinejad-Eramsadati L. The Pattern of Road Accidents in Fog and the Related Factors in North of Iran in 2014-2018 ​. Durham, North Carolina: Research Square; 2020, [cited 2022]. Available from:

  • 47.

    Delavary M, Ghayeninezhad Z, Lavallière M. Evaluating the Impact of Increased Fuel Cost and Iran’s Currency Devaluation on Road Traffic Volume and Offenses in Iran, 2011–2019. Safety. 2020;6(4):49.

  • 48.

    Delavary Foroutaghe M, Mohammadzadeh Moghaddam A, Fakoor V. Impact of law enforcement and increased traffic fines policy on road traffic fatality, injuries and offenses in Iran: Interrupted time series analysis. PLoS One. 2020;15(4). e0231182. [PubMed ID: 32302374]. [PubMed Central ID: PMC7164613].

  • 49.

    Rahimi H, Hashemi Nazari SS, Soori H, Motevalian SA, Momeni E, Azar A. Traffic Police Effectiveness and Efficiency Evaluations, an Overview of Methodological Considerations. Arch Trauma Res. 2017;6(1). e36927.

  • 50.

    Safarzadeh M, Bagheri R. [Comparative studies of traffic fines by traffic police in Iran and other countries]. Rahvar. 2012;9(17):59-74. Persian.

  • 51.

    Shams M, Rahimi-Movaghar V. Risky driving behaviors in Tehran, Iran. Traffic Inj Prev. 2009;10(1):91-4. [PubMed ID: 19214883].

  • 52.

    Tavolinejad H, Malekpour MR, Rezaei N, Jafari A, Ahmadi N, Nematollahi A, et al. Evaluation of the effect of fixed speed cameras on speeding behavior among Iranian taxi drivers through telematics monitoring. Traffic Inj Prev. 2021;22(7):559-63. [PubMed ID: 34424783].

  • 53.

    de Winter JCF, Dodou D. National correlates of self-reported traffic violations across 41 countries. Pers Individ Differ. 2016;98:145-52.

  • 54.

    Antin J. Design of the In-Vehicle Driving Behavior and Crash Risk Study: In Support of the SHRP 2 Naturalistic Driving Study. Washington, DC: Transportation Research Board; 2011.

  • 55.

    Kavousi A, Moradi A, Soori H, Rahmani K. Environmental factors influencing the distribution of pedestrian traffic accidents in Iran. Arch Trauma Res. 2020;9(1):8-15.

  • 56.

    James G, Witten D, Hastie T, Tibshirani R. An introduction to statistical learning. New York, NY: Springer; 2013.

  • 57.

    Ramsay JO; Silverman. Functional Data Analysis. New York, NY: Springer; 2005.

  • 58.

    Stasinopoulos D, Rigby RA. Generalized Additive Models for Location Scale and Shape (GAMLSS) inR. J Stat Softw. 2007;23(7):1-46.

  • 59.

    Banks D, Persaud B, Lyon C, Eccles K, Himes S. Enhancing Statistical Methodologies for Highway Safety Research – Impetus from FHWA. McLean, VA; 2014. Contract No.: FHWA-HRT-14-081.

  • 60.

    Donnell E, Hanks E, Porter RJ, Cook L, Srinivasan R, Li F, et al. The Development of Crash Modification Factors: Highway Safety Statistical Paper Synthesis. McLean, VA; 2020. Contract No.: FHWA-HRT-20-069.

  • 61.

    Lohr S. Measuring Crime: Behind the Statistics. New York: Chapman and Hall/CRC; 2019.