Assessment and prediction of road accident injuries trend using time-series models in Kurdistan
© The Author(s) 2018
Received: 17 July 2017
Accepted: 12 February 2018
Published: 9 March 2018
Road traffic accidents are commonly encountered incidents that can cause high-intensity injuries to the victims and have direct impacts on the members of the society. Iran has one of the highest incident rates of road traffic accidents. The objective of this study was to model the patterns of road traffic accidents leading to injury in Kurdistan province, Iran.
A time-series analysis was conducted to characterize and predict the frequency of road traffic accidents that lead to injury in Kurdistan province. The injuries were categorized into three separate groups which were related to the car occupants, motorcyclists and pedestrian road traffic accident injuries. The Box-Jenkins time-series analysis was used to model the injury observations applying autoregressive integrated moving average (ARIMA) and seasonal autoregressive integrated moving average (SARIMA) from March 2009 to February 2015 and to predict the accidents up to 24 months later (February 2017). The analysis was carried out using R-3.4.2 statistical software package.
A total of 5199 pedestrians, 9015 motorcyclists, and 28,906 car occupants’ accidents were observed. The mean (SD) number of car occupant, motorcyclist and pedestrian accident injuries observed were 401.01 (SD 32.78), 123.70 (SD 30.18) and 71.19 (SD 17.92) per year, respectively. The best models for the pattern of car occupant, motorcyclist, and pedestrian injuries were the ARIMA (1, 0, 0), SARIMA (1, 0, 2) (1, 0, 0)12, and SARIMA (1, 1, 1) (0, 0, 1)12, respectively. The motorcyclist and pedestrian injuries showed a seasonal pattern and the peak was during summer (August). The minimum frequency for the motorcyclist and pedestrian injuries were observed during the late autumn and early winter (December and January).
Our findings revealed that the observed motorcyclist and pedestrian injuries had a seasonal pattern that was explained by air temperature changes overtime. These findings call the need for close monitoring of the accidents during the high-risk periods in order to control and decrease the rate of the injuries.
Road traffic accidents are associated with high mortality, severe injuries, and considerably high economic losses globally . The road traffic accidents account for 12% of the burden of disease and to the high rates of unintentional injuries . About 1.2 million deaths and more than 50 million injuries occur annually. Therefore, road traffic accident is one of the most common causes of mortality and disability in the world [3, 4]. Evidences from Iran on the disability-adjusted life year (DALY) index indicated that road traffic accident was ranked ninth in 1990, and is expected to be the third by 2020 .
The pedestrians, cyclists, and motorcyclists were particularly at high-risk of road traffic accidents and more than 90% of the deaths from the accidents occurred in low-income countries . The death rate in the developing countries has been on the rise while it was continuously decreasing in the developed countries. Currently, the car ownership in most Asian countries has increased twice faster than the per capita income , and this might have contributed to the increasing rates of accidents in the region.
The death rate of road traffic accidents in Iran in 2002 which was at 44 per 100,000 populations was more than twice the global average (19 per 100,000) . The annual DALY loss from road traffic accidents was more than 1,200,000 per year . This event accounted for 57% of the DALY lost and more than 80% of the victims were among males of 15–29 years .
Since the past seven decades, the deaths from road traffic accidents in Iran have been increasing. The mortality trend during the years 2001–2005 reached as high as 28,000 deaths per year . Despite the reduction in mortality from 38 per 100,000 in 2006 to 31 per 100,000 in 2008, the number of injured people from road traffic accident has been increasing overtime [11, 12]. Thus, the mortality of road traffic accidents in Iran remained among the highest rates in the world . Despite all interventions to reduce traffic accident injuries, even fatal injuries have still a high incidence rate in Iran .
The geographic, environmental, demographic, and individual factors can contribute to the occurrence of road traffic accidents. Hence, access to the information about road traffic accidents in a given context is significant to generate evidence to contribute to the prevention and control of context-specific accidents . The objective of this study was modeling the frequency of injuries due to road accidents in Kurdistan Province, located in the west of Iran, from March 2009 to February 2015 and to predict injuries until February 2017 using time series models.
Study design and population
This cross-sectional study was conducted to analyze time-series observations of traffic accident injuries from March 2009 to February 2015 in Kurdistan Province. The data was obtained from the Accidents’ Unit of the Health Deputy in the Kurdistan University of Medical Sciences. The unit regularly gathers data on road traffic accident injuries from hospitals, police stations, forensic medicine, and road organization. In order to ensure the quality of the collected data, any duplicate or redundant information concerning the injuries was cleaned. Then, the data on the injured persons was divided into three separate categories: pedestrian injuries, motorcyclist injuries, and car occupant injuries and each category were analyzed separately.
Where φ p (β) = (1 − φ1β − … − φ p β p ), θ q (β) = (1 − θ1β − … − θ q β q ), Φ P (β S ) = (1 − Φ1β S − … − Φ P β Ps ) and Θ Q (β s ) = (1 − Θ1β S − … − Θ Q β Qs ). In these extensions, φ and θ are the ordinary autoregressive and moving average parameters; and, Φ and Θ are the seasonal autoregressive and moving average parameters. Besides, β represents the time-delay parameter so that β r Z t is equal to Zt − r.
The Box-Jenkins modeling involved several steps, including model identification, parameter estimation, checking, and prediction . The first step included testing for “white noise” in the data using the Portmanteau and Bartlett test. The stationarity in the seasonal and non-seasonal series was checked by Dickey-Fuller test. For the SARIMA model construction, the order of (p, d, q) (P, D, Q) was distinguished based on the ACF and PACF plots. The cutoff point for the ACF plot and the slow decay in the PACF plot were used to identify the correct parameters. The model parameters were estimated using ordinary least-squares approach and t-test was considered to decide on the significant difference . Based on the Box-Jenkins approach, the most parsimonious SARIMA model that requires the least moving average and autoregressive parameters should be selected.
The Akaike Information Criteria (AIC) and Schwartz Bayesian Criteria (SBC) were calculated to evaluate the goodness of fit for each model . These indicators evaluated the model fitness based on the likelihood model and a number of parameters. The smaller the size, the better was the model. The absence of process autocorrelation and partial autocorrelation in the residuals are also indicators of the model’s goodness of fit. Finally, the fitted model was used to predict the trend of road accident injuries from February 2015 to February 2017 in the context. All the analyses and the forecast were computed using the R statistical software package. The statistical significance was decided at p < 0.05.
Frequencies of pedestrians’, motorcyclists’, and car occupants’ injuries from March 2009 to February 2015
Frequency of injuries per year (n)
Type of accident
Seasonality of pedestrian accident injuries
Parameter estimations and goodness of fit measures for pedestrians’, motorcyclists’, and car occupants’ injuries from Mar 2009 to Feb 2015
Car occupants’ injuries
Predicted frequencies of pedestrians’, motorcyclists’, and car occupants’ injuries from March 2015 to February 2017
Car occupants (n)
Seasonality of motorcyclist’s accident injuries
The residuals’ ACF, PACF, and trend plots exhibited a white noise pattern and the goodness of the fitted model (Fig. 2e–g). Finally, the frequency of motorcyclists’ accident injuries was predicted until February 2017 using the above mentioned Eq. (2) (Fig. 2d, Table 3).
Car occupants’ injuries in car accident
The model’s goodness of fit was confirmed based on the residuals’ behavior as shown in Fig. 3e–g, showed the predicted frequency of the car occupants’ injuries for the 24 months period using the model which was presented in Eq. 3 is shown in Table 3 and Fig. 3d.
This study applied the time-series modeling to determine and predict road traffic accident injuries patterns in Kurdistan Province, Iran. Time-series analysis was a widely used method for modeling phenomena . Several studies also applied this model to predict the road traffic accident fatal injuries [12, 18, 19] and were proven effective for determining accident prevention strategies in a short time . The remarkably high frequency of road traffic injuries (n = 43, 120 cases) observed in our study in about 6 years indicated the public health importance of the issue. Despite the report of a descending trend of deaths from road traffic accident , our findings indicated an increasing trend in the frequency of non-fatal injuries. Another study in Kermanshah Province also indicated a seasonally descending trend of deaths between the years 2006 and 2013 . More than 67% occupant injuries and about 12% pedestrian injuries were the maximum and minimum frequency of road traffic injuries. The predicted frequency of car occupant injuries during the period 2015 to 2016 revealed a high rate of injuries. These findings were not consistent with the findings of another study during the years 2008 to 2014 which reported a decrease trend of the car occupants’ injuries . The reductions in the deaths from accidents can be the consequence of extensive and coordinated efforts of several organizations that bold the importance of some interventions and strict monitoring the implementation of the interventions. The interventions included the mandatory use of seat belt by drivers and passenger vehicles, helmet use by motorcyclists, improving the quality of vehicles and roads, and public education programs aimed at prevention of accidents [23, 24]. These interventions can reduce the occurrence of death from the road traffic accidents, but may not completely prevent injuries.
Despite the reducing trend in the death rate from road traffic accidents, evidence showed that Iran was one of the countries with high incidence of road traffic accidents . There were 23,249 fatal road traffic injuries in 2013 in Iran  which was more than seven times the report from Germany (3648) and about six times more than that in Turkey (4045) . The extremely high incidence of road accidents in Iran may relate to geographic factors, low cost of gasoline, poor public transportation system, rapidly growing automotive industry, poor quality of cars and motorcycle equipment, and unsafe roadway networks .
The trend of motorcycle injuries in this study was extremely increasing overtime. The high price of cars and the heavy traffic jams, might have contributed to the increased number of individuals who used motorcycle for transportation. Thus, the increasing trend of injuries might be associated with an increased number of motorcycle riders. Epidemiologic studies indicated that most accidents due to motorcyclists were related to younger men and the low-income groups [26, 27]. However, there is concrete evidence concerning the factors related to motorcyclist accidents in Iran . Hospital-based evidence from USA indicated an increasing number of motorcyclist accident during the years 2001 to 2008 regardless of their age . Others reported forward inattention, yield non-conformity, and moving in the opposite direction as the most common causes of motorcyclist accidents . The helmet use was one of the effective ways of decreasing the likelihood of death and head injuries among motorcyclist . According to the 2010 WHO report, only 30% motorcycle riders and 10% of motorcycle passengers use a helmet . In a study conducted in Kerman Province, it was also found that from total of 2880 motorcyclists, only 324 (11.2%) used a helmet and 5% of them were passengers .
As you can see in Figs. 1, 2, and 3, our predicted trends in the number of injuries after February 2015 is so smooth, this case may be due to small autocorrelation and partial autocorrelation coefficients of the fitted models (Table 2). Because small-value coefficients make the process of injuries similar to a random walk model which it returns a flat forecast function.
Our study followed 12-month seasonal pattern of injuries and most of the car occupants’ injuries were predicted to occur in August and September. The evidences from other studies also showed similar findings [32, 33]. A study in Esfahan Province found that most of the motorcycle accidents occurred during the August and the least occurred during the December , while a report from a meta-analysis indicated that the majority of the motorcycle accidents in Iran occurred during the first half of the year in the evening . The large number of motorcycle injuries during August in our study might be related to the low helmet use during the hot season. Furthermore, it might be related to the high mobility of people during summer season usually for holidays .
Kurdistan Province is geographically mountainous, and the road networks connect it with other provinces in the northern, western, and southern parts of Iran. A highly busy roadway can also increase the risk of human factors that lead to road traffic accidents. The road traffic accidents in a mountainous province can be highly fatal . Thus, the road traffic accident prevention intervention programs and periodic evaluation of those actions can help reduce the incidence of accidents.
Our findings indicated increasing trends of the car occupants’ and motorcyclists’ injuries in Kurdistan Province during the period of 2009 to 2017. The majority of the road traffic injuries occurred during the summer season. The decision-makers and planners in Iran can draw advantage of the findings for planning the prevention and control measures of road traffic injuries in the context.
This work was supported by a research grant from the Social Determinants of Health Research Center, Kurdistan University of Medical Sciences, Sanandaj, Iran (grant number IR.MUK. REC.1395.88).
This study was supported by Kurdistan University of Medical Sciences, Deputy of Research (Project 1395/88).
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
MP and BN carried out the study conception, design, and implementation. AK contributed to the analyzing and interpretation of manuscripts. AW and SR provided assistance in the interpretation and writing of the manuscript. SN contributed to the design and preparation of the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
This study was approved by the ethics committee of Kurdistan University of Medical Sciences (IR.MUK. REC.1395.88).
Consent for publication
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hesari A, Esmaeli A. Estimates of deaths from traffic accidents on life expectancy at birth and the financial burden it (2002). Health Inf Manag. 2004;1(2):27–35.Google Scholar
- Marasy MR, Tabar IM. The burden of road traffic injuries in Isfahan, Iran in 2010. J Kerman Univ Med Sci. 2013;20(5):505–19.Google Scholar
- Yazdani CJ, Ahmadi BE, Ghadami M. Mapping of mortality rate in suburban accidents, Mazandaran Province, 2007-2010. J Mazandaran Univ Med Sci. 2012;22(97):50–8.Google Scholar
- Ghorbani A, Nabavi fard H, Khoshhal M, Hosseini H. Costs imposed on the effects of mortality due to traffic accidents (Sabzevar). Traffic Manag Stud. 2011;20:49–58.Google Scholar
- Ayatollahi SH, Hassanzadeh J, Ramezani A. The burden of traffic accidents in South Khorasan Province, Iran in 2005. Iran J Epidemiol. 2009;4(3):51–7. [In Persian]Google Scholar
- Herman J, Ameratunga S, Jackson R. Burden of road traffic injuries and related risk factors in low and middle-income Pacific Island countries and territories: a systematic review of the scientific literature (TRIP 5). BMC Public Health. 2012;12:479.View ArticlePubMedPubMed CentralGoogle Scholar
- Mohan VR, Sarkar R, Abraham VJ, Balraj V, Naumova EN. Differential patterns, trends and hotspots of road traffic injuries on different road networks in Vellore district, southern India. Tropical Med Int Health. 2015;20(3):293–303.View ArticleGoogle Scholar
- Naghavi M, Shahraz S, Bartels D, Puthenpurakal JA, Motlagh ME. Adverse health outcomes of road traffic injuries in Iran after rapid motorization. Arch Iran med. 2009;12(3):284–94.PubMedGoogle Scholar
- Naghavi MA, Abolhassani F, Pourmalek F, Jafari N, Moradi LM, Eshrati B. The burden of disease and injury in Iran 2003. Iran J Epidemiol. 2008;4(1):1–19. (Persian)Google Scholar
- Moradi A, Kh R. Trend of traffic accidents and fatalities in Iran over 20 years (1993-2013). J Mazandaran Univ Med Sci. 2014;24(118):186–97.Google Scholar
- Ghadirzadeh MR, Shojaei A, Khademi A, Khodadoost M, Kandi M, Alaeddini F, et al. Status and trend of deaths due to traffic accidents from 2001 to 2010 in Iran. Iran J Epidemiol. 2015;11(2):13–22.Google Scholar
- Bakhtiyari M, Mehmandar MR, Riahi SM, Mansournia MA, Sartipi M, Bahadorimonfared A. Epidemiologic pattern of fatal traffic injuries among Iranian drivers; 2004–2010. Iran J Public Health. 2016;45(4):503–14.PubMedPubMed CentralGoogle Scholar
- Moradi A, Rahmani K, Hoshmandi-Shoja M, Rahimi-Sepehr H, Khorshidi A. An overview of the situation of traffic accidents in Iran in comparison with other countries. Iran J Forensic Med. 2016;22(1):45–53.Google Scholar
- Duenas C, Fernandez MC, Canete S, Carretero J, Liger E. Stocastic model to forecast ground level ozone concentration at urban and rural areas. Chemosphere. 2005;61(10):1379–89.View ArticlePubMedGoogle Scholar
- Lin Y, Chen M, Chen G, Wu X, Lin T. Application of an autoregressive integrated moving average model for predicting injury mortality in Xiamen. China BMJ Open. 2015;5(12):1–8.Google Scholar
- Box GEP, Jenkins GM, Reinsel GC, Ljung JM. In: David JB, Noel AC, Garrett MF, editors. Time series analysis: forecasting and control. New Jersey: Wiley; 2016. p. 179–209. Academic.Google Scholar
- Mansouri F, Khanjani N, Rananadeh KL, Pourmousa R. Forecasting air pollutant situation using the time series models in Kerman, Iran. Sci J School Pub Health Ins Pub Health Res. 2013;11(2):75–86.Google Scholar
- Bahadorimonfared A, Soori H, Mehrabi Y, Rahmati RM, Esmaili AR, Salehi M, et al. Trends of fatal road traffic injuries in Iran (2004–2011). PLoS One. 2013;8(5):1–5.View ArticleGoogle Scholar
- Bahadori MA, Soori H, Mehrabi Y, Rahmati RM, Esmaili AR, Salehi M, et al. A model for prediction of on the rate of mortality due to road traffic accidents in Iran. Res Med. 2013;36(5):7–11.Google Scholar
- Mehmandar MS, Soori H, Mehrabi Y. Predicting and analyzing the trend of traffic accidents deaths in Iran in 2014 and 2015. Int J Crit Illn Inj Sci. 2016;6(2):74–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Zolala F, Haghdoost AA, Ahmadijouybari T, Salari A, Bahrampour A, Baneshi MR, et al. Forecasting the trend of traffic accident mortality in West Iran. Health Scope. 2016;5(3):1–5.View ArticleGoogle Scholar
- Yousefzadeh CS, Ranjbar TF, Malekpouri R, Razzaghi A. A time series model for assessing the trend and forecasting the road traffic accident mortality. Arch Trauma Res. 2016;5(3):1–6.Google Scholar
- Rasouli MN, Nouri M, Zarei MR, Saadat S, Rahimi MV. Comparison of road traffic fatalities and injuries in Iran with other countries. Chin J Traumatol. 2008;11(3):131–4.View ArticlePubMedGoogle Scholar
- Soori H, Royanian M, Zali AR, Movahedinejad A. Study of changes on road traffic injury rates, before and after of four interventions by Iran traffic police. Pajoohandeh J. 2009;14(1):15–20.Google Scholar
- Soori H, Iranfar M. Road traffic status in the world and Iran: review of results from the World Health Organization. J Saf Prom Inj Prev. 2013;1(2):53–62.Google Scholar
- Mirzaei M, Mirzadeh M, Shogaei-Far H, Mirzaei M. Trends in road traffic deaths in Yazd, Iran, 2004 - 2010. Arch Trauma Res. 2016;5(2):1–6.Google Scholar
- Sadeghi 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):1–10.Google Scholar
- Jackson TL, Mello MJ. Injury patterns and severity among motorcyclists treated in US emergency departments, 2001–2008: a comparison of younger and older riders. Inj Prev. 2013;19(5):297–302.View ArticlePubMedGoogle Scholar
- Khorshidi A, Ainy E, Soori H. Epidemiological pattern of road traffic injuries among Iranian motorcyclist in 2012. J Saf Promot Inj Prev. 2016;4(1):47–54.Google Scholar
- Liu BC, Ivers R, Norton R, Boufous S, Blows S, Lo SK. Helmets for preventing injury in motorcycle riders. Cochrane Database Syst Rev. 2008;23:1.Google Scholar
- Mokhtari AM, Samadi S, Hatami SE, Jalilian H, Khanjani N. Investigating the rate of helmet use and the related factors among motorcyclist in Kerman between 1391–92. J Saf Prom Inj Prev. 2014;2(3):209–14.Google Scholar
- Heydari SH, Hoseinzadeh A, Sarikhani Y, Hedjazi A, Zarenezhad M, Moafian G, et al. Time analysis of fatal traffic accidents in Fars Province of Iran. Chin J Traumatol. 2013;16(2):84–8.PubMedGoogle Scholar
- Hasanzadeh J, Moradinazar M, Najafi F, Ahmadijouybary T. Trends of mortality of road traffic accidents in Fars Province, southern Iran, 2004–2010 Iranian. Aust J Public Health. 2014;43(9):1259–65.Google Scholar
- Mohammadian M, Hajare A, Mohammadian HA. Incidence trends of injury and mortality from traffic accidents in urban and suburban areas of Isfahan Province during 2002-2010. J Police Med. 2014;3(1):1–10.Google Scholar