Next Article in Journal
Spatial Analysis of Industrial Benzene Emissions and Cancer Incidence Rates in Texas
Next Article in Special Issue
Risk Assessment in Urban Large-Scale Public Spaces Using Dempster-Shafer Theory: An Empirical Study in Ningbo, China
Previous Article in Journal
How to Stop Victims’ Suffering? Indirect Effects of an Anti-Bullying Program on Internalizing Symptoms
Previous Article in Special Issue
Risk Riding Behaviors of Urban E-Bikes: A Literature Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Investigation on the Injury Severity of Drivers in Rear-End Collisions Between Cars Using a Random Parameters Bivariate Ordered Probit Model

The Key Laboratory of Road and Traffic Engineering, Ministry of Education Tongji University, Shanghai 201804, China
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(14), 2632; https://doi.org/10.3390/ijerph16142632
Submission received: 25 June 2019 / Revised: 22 July 2019 / Accepted: 22 July 2019 / Published: 23 July 2019

Abstract

:
The existing studies on drivers’ injury severity include numerous statistical models that assess potential factors affecting the level of injury. These models should address specific concerns tailored to different crash characteristics. For rear-end crashes, potential correlation in injury severity may present between the two drivers involved in the same crash. Moreover, there may exist unobserved heterogeneity considering parameter effects, which may vary across both crashes and individuals. To address these concerns, a random parameters bivariate ordered probit model has been developed to examine factors affecting injury sustained by two drivers involved in the same rear-end crash between passenger cars. Taking both the within-crash correlation and unobserved heterogeneity into consideration, the proposed model outperforms the two separate ordered probit models with fixed parameters. The value of the correlation parameter demonstrates that there indeed exists significant correlation between two drivers’ injuries. Driver age, gender, vehicle, airbag or seat belt use, traffic flow, etc., are found to affect injury severity for both the two drivers. Some differences can also be found between the two drivers, such as the effect of light condition, crash season, crash position, etc. The approach utilized provides a possible use for dealing with similar injury severity analysis in future work.

1. Introduction

With the increase of vehicle miles/kilometers travelled, traffic crash has become one of the main factors that cause human injury and death along with huge property damage. According to the World Health Organization, nearly 1.25 million people died in traffic crashes worldwide annually and each fatality, on average, causes a loss of about 1.42 million dollars [1]. Injury crashes become a major concern of researchers, policy-makers, and the public [2]. There are various types of traffic crash such as rear-end, head-on, forward impact, side wipe, etc., among which rear-end crash results in a much great proportion of severe injuries or fatalities [3]. For instance, rear-end crashes account for 30% of all injuries and 29.7% of all property damage in the USA [4]. As one of the main concerns of the rear-end collision, drivers’ injury severity is influenced by many contributing factors, such as roadway alignments, environment characteristics, driver characteristics, traffic flow, etc. Accordingly, it is necessary to identify the factors related to rear-end crashes and having a comprehensive knowledge of these factors’ potential positive or negative effects is essential to prevent crashes and reduce injury severities [5].
To investigate the contributing factors of rear-end crash, vehicle laboratory crash test [6], numeric simulation [7,8], and field crash data analysis are three effective tools. As a direct and well-targeted method, vehicle laboratory crash test usually uses the scrapped vehicles and human body models to conduct vehicle collision tests. Such kinds of tests are non-repeatable and expensive; thus, fewer samples can be obtained. Conversely, the numerical simulation can be used to explore the injury mechanism by setting different parameter values and combinations of significant factors such as collision angle, impact speed, vehicle types, etc. The numerical simulation is a powerful tool; however, it needs many refined modeling and complex mechanics analysis, and some natural conditions are usually neglected.
As the crash data recorded in the real world become more and more detailed and specific, the field crash data analysis has been the most widely used method to analyze human injury severity in rear-end crashes. As pointed out by Chen et al. [9], these naturally collected data usually include almost all crash-related perspectives like human factors (age, seat position, sex, fatigue, alcohol usage, and so on), vehicle (speed, vehicle type, weight, and so on), and environment (light condition, road surface, weather, and so on). These data reflect the real state of crashes pretty well. Moreover, Wang et al. developed a quasi-vehicle-trajectory-based method to find crash contributing factors, presenting a possible use of real-time data [10]. Based on various methodological and analytical techniques in big data analysis, many prediction models are also obtained and applied, including multinomial logit models, latent class logit models, random parameter (mixed) logit models, ordered logit/probit models, support vector machine, and so on [11,12,13].
In terms of specific models for investigating injury severity in previous studies, there are also abundant researches. Yuan et al. developed a binary logistic regression model to predict occupant injury severity. This model identified corresponding affecting factors in rear-end crashes involving trucks as the front vehicle [14]. A multinomial logit model was utilized by Konstantina et al. to investigate severity outcomes of farm vehicle crashes and it was found that antiquated farm vehicles were more likely to suffer serious injuries or fatalities in crashes under insufficient light conditions [15]. Via the same approach, the corresponding factors were identified to classify four severity outcomes based on the most severe injuries in crashes [16]. Substantial difference between the impacts of variables on the driver-injury severity in single- and multi-vehicle accidents had been found by a mixed logit model developed by Chen and Chen [17].
Recently, some innovative methods have also been used in analyzing injury severity. Kunt et al. compared a genetic algorithm and an artificial neural network. The results demonstrated that the artificial neural network provides the best prediction [18]. In [3], Chen et al. studied the driver injury severity in rear-end crashes with the use of a multinomial logit model-Bayesian network hybrid approach, which performs reasonably well. The results indicated that the factors including inferior lighting conditions, windy weather, etc. could significantly increase driver injury severities in rear-end crashes. Other related researches can also be found and they present detailed insights into injury severity analysis [19,20,21,22,23,24,25,26,27].
Abundant as the models in assessing the injury severity seem, each model often has its own limitations. For instance, the genetic algorithm, artificial neural network, and Bayesian network have some deficiencies in interpreting the relationship between a certain factor and the injury severity. Moreover, injury severity has its natural character as severity outcomes are ordered from a low level to a high level naturally, and this should be noted in the model. After decades of exploration, the ordered probit model has been one of the most common approaches used in accident severity studies. In [28], Chandler et al. demonstrated an ordered probit model to examine the impact of various factors on injuries severity to passenger car occupants involved in truck-car collisions. Abdel-Aty developed several ordered probit models to examine driver injury severity for crashes at roadway sections, signalized intersections, and toll plazas in Central Florida, respectively [29]. Alcohol, lighting conditions, and the existence of a horizontal curve were found to be of significance in the roadway sections’ model. An ordered probit model was developed to analyze the injury severity of wrong-way driving crashes, which also found the lighting condition to be significant to injury severity [30].
In this paper, instead of the most severe injury, we focused mainly on comparing the injury severity of the two drivers in passenger cars involved in the same rear-end crash. There exists a problem if we simply use the ordered probit model, since the two drivers share some common conditions such as road surface condition, but the model allows only one dependent variable. Moreover, there may be some unobserved factors causing correlation. To address this correlation, the bivariate ordered probit model was proposed and two dependent variables were designed for each crash. Furthermore, compared with the fixed-parameter model, the random-parameter regression model is often adopted, because the unobserved heterogeneity issues can be addressed by allowing parameters to vary across observations [2,5,31].
The remainder of this paper is organized as follows: in Section 2, the random parameter bivariate ordered probit model adopted in this paper is introduced in detail; in Section 3, a brief description of the rear-end crashes is provided; Section 4 presents the results and discussion of the model; Finally, Section 5 provides the summaries and conclusions of the work.

2. Methodology

As previously stated, the injury severity levels of two drivers involved in the same rear-end crash are typically correlated considering they usually share the same lighting condition, road alignment, road surface condition, and, especialy, other unobserved factors. To address such kinds of possible correlation problems, the bivariate ordered probit model, which is a hierarchical system of two equations, can be employed to model a simultaneous relationship of two response variables [22].
Suppose y i j is the observed injury severity; (i, j) is the index indicating the two drivers involved in the same rear-end crash, where i (i = 1, 2, 3, …, n) refers to the crash number and j refers to the driver number in a certain crash i (j = 1 for the rear vehicle, 2 for the front vehicle). Thereafter, the latent (unobserved) injury severity propensities of the two drivers match their actual injury severity, as presented in the following equations [7]:
y i , j   =   1 = k ,       i f   μ j   =   1 ,   k 1 < y i ,     j = 1 * < μ j   =   1 ,   k
y i , j   =   2 = l ,         i f   μ j   =   2 , l 1 < y i ,     j = 2 * < μ j   =   2 ,   l
where μ j , k 1 , μ j , k , μ j , l 1 , μ j , l are thresholds or cut-off values used to determine observed injury severity levels of both two drivers, their values are relative to their corresponding injury factors in crash i . Additionally, k (k = 0, 1, 2, …, K) and l (l = 0, 1, 2, …, L) represent ordinal categories of injury severity sustained by each driver. y i ,     j   =   1   * and y i ,     j = 2   * can be calculated using real data as follows:
y i ,     j = 1 * = β 0 , j = 1 + β 1 , j = 1 * x 1 i , j = 1 + β 2 , j = 1 * x 2 i , j = 1 + + ε i , j = 1 = β 1 X i , j = 1 + ε i , j = 1
y i ,     j = 2 * = β 0 , j = 2 + β 1 , j = 2 * x 1 i , j = 2 + β 2 , j = 2 * x 2 i , j = 2 + + ε i , j = 2 = β 2 X i , j = 2 + ε i , j = 2
where X i j is the variable vector; β j is a parameter vector remaining to be estimated; ε i j represents the random components that capture all unobserved factors associated with two involved parties, which is assumed to follow a bivariate normal distribution as follows:
( ε i , j = 1 ε i , j = 2 ) ~ N [ ( 0 0 ) ,   ( 1     ρ ρ     1 )   ]
where ρ is the estimated correlation parameter between ε i , j   = 1 and ε i , j = 2 . Therefore, the joint probability when the injury severities of two drivers in a rear-end crash are k and l , which, respectively, can be expressed as follows:
Pr ( y i , j = 1 = k ;   y i , j = 2 = l )                 = Pr ( μ j = 1 ,   k 1 < y i ,     j = 1 * < μ j = 1 ,   k ;   μ j = 2 , l 1 < y i ,     j = 2 * < μ j = 2 ,   l )                 = Pr ( μ j = 1 ,   k 1 < β 1 X i 1 + ε i 1 < μ j = 1 ,   k ;   μ j = 2 , l 1 < β 2 X i 2 + ε i 2 < μ j = 2 ,   l )                 = Pr ( μ j = 1 ,   k 1 β 1 X i 1 < ε i 1 < μ j = 1 ,   k β 1 X i 1 ;   μ j = 2 , l 1 β 2 X i 2 < ε i 2 < μ j = 2 ,   l β 2 X i 2 )       = Φ 2 ( μ j = 1 ,   k β 1 X i 1 , μ j = 2 ,   l β 2 X i 2 ; ρ ) Φ 2 ( μ j = 1 ,   k 1 β 1 X i 1 , μ j = 2 ,   l β 2 X i 2 ; ρ )             Φ 2 ( μ j = 1 ,   k                       β 1 X i 1 , μ j = 2 ,   l 1 β 2 X i 2 ; ρ )             + Φ 2 ( μ j = 1 ,   k 1 β 1 X i 1 , μ j = 2 ,   l 1 β 2 X i 2 ; ρ )
where Φ2( ) is the standard bivariate normal cumulative distribution function.
Due to the characteristic of the bivariate ordered probit model, the signs of the parameters are of great interest. A positive sign of the parameter indicates a positive effect on the probability of the injury severity.
While bivariate ordered probit can address the problem of factors correlation between two drivers involved in the rear-end crash, this method assumes the parameters β 1 ,   β 2 to have a certain value neglecting the effect of unobserved heterogeneity of observations. This constraint on model parameters may lead to inconsistent and biased parameter estimation [2]. As previously stated, the random-parameter method can address the unobserved heterogeneity by allowing the parameters to vary across observations. Therefore, the random parameters bivariate ordered probit model can be derived by setting:
β i = β + γ i
where β i is the vector of specific parameters; γ i is the randomly distributed term which is normally distributed with a zero mean value and variance σ 2 .

3. Data Description

In this paper, the crash data from the National Automotive Sampling System (NASS) General Estimates System (GES) were obtained. Data for all car-car rear-end collisions occurring from the calendar year 2011 to 2015 were sampled as the basic dataset. These crash data recorded details of the crash time, driver information, environmental condition, roadway condition, vehicle characters, etc. Some data were excluded for those crashes in which the driver’s injury severity was unknown or serious errors and missing existed.
The final dataset included a total of 15,159 rear-end crashes involving 30,318 motor vehicle drivers. The original injury severity was coded on the KABCO scale and was regrouped due to the low proportion of fatal injuries and incapacitating injury levels. The new groups were presented as follows:
  • L1: No injury (O) and no apparent injury (O);
  • L2: Possible injury (C);
  • L3: Non-incapacitating evident injury (B), suspected minor injury (B);
  • L4: Incapacitating injury (A), suspected serious injury (A), fatal injury (K).
As previously stated, the differences and similarities of the two drivers in the same rear-end crash should be addressed. When analyzing the crash data, drivers were separated into two groups: the front vehicle and the rear vehicle. Table 1 presents a summary of the joint distribution of the injury severity for drivers in front vehicles and drivers in rear vehicles. It can be found that the distributions of injury severity for drivers in front vehicles and rear vehicles vary in some degrees. For instance, 33.4% of the drivers in front vehicles tend to suffer an injury while the value is 19.4% for rear vehicles. Moreover, the proportions of L2, L3, L4 for rear vehicles are all smaller than those of front vehicles. This indicates that the drivers in front vehicles seem to suffer more serious injury.
A large number of explanatory variables were examined and possible related variables are presented in Table 2 and Figure 1 (some variables with complex classifications are depicted in Figure 1, the percentage curves are also presented). It is interesting to note some substantial differences between the two drivers. For instance, the drivers in rear vehicles tend to have a higher proportion of alcohol involved and drug use. Since prior work has identified many injury-related variables, investigation on the differences and similarities between the two drivers involved in the same rear-end crash needs further efforts and the quantification of variables’ effect on injury severity needs to be addressed by the statistical model.
To improve the process of modeling, the explanatory variables were coded into dummy variables (0/1) when conducting model estimations. For instance, the crash season can be described using three binary indicates (1,0,0 for Spring; 0,1,0 for Summer; 0,0,1 for Autumn and 0,0,0 for Winter).

4. Model Estimation and Discussion of Results

To start the estimation process of the random parameter bivariate ordered probit model, the initial value of variable coefficients had to be set up. Therefore, two independent ordered probit models were developed respectively after separating the dataset into two parts (one for the driver in the front vehicle and another for the driver in the rear vehicle). All variables that were significant at a 90% confidence level in two separate models were remained for the subsequent analysis of random parameters bivariate ordered probit model.
Table 3 presents the results of the random parameters bivariate ordered probit models. Parameter estimates, p-values, standard errors, and z value of the estimates were included. Compared with the results of two separately ordered probit models, substantive improvements can be found in the random parameters bivariate ordered probit model, such as the significance level of variables and the z value. This indicates that the random parameters bivariate ordered probit model indeed addressed the issue of correlation between the outcomes in the same crash and the unobserved heterogeneity across the observations.
As presented in Table 3, the final model includes all significant variables and some of them are set as fixed parameters or random parameters according to constant estimation (the parameter without a significant standard deviation is set as a fixed parameter). The value of ‘Percent observations’ indicates the distribution of the value of random parameter (above zero and below zero, respectively), which is determined by the value of mean and standard deviation.
From Table 3, the coefficient value of younger drivers (age 24 and below) was negative for both the two parties indicating that younger drivers generally sustained lower level of injury severity. This result was consistent with previous studies of Abay et al. and Chiou et al. and this may be because of the physiological or driving behavior differences such as reaction time [32,33]. Moreover, drivers in the rear vehicle at an age between 25 and 63 were found to suffer a lower level of injury severity, while the result of drivers in the front vehicle went in the opposite way.
As for the driver’s gender, male drivers in both two vehicles seemed to suffer less severe injuries. However, the effect of gender seemed to be mixed as the male driver had been found to be less injured in a previous work by Russo et al. [11] and more at risk for injuries in another work by Chen et al. [34]. In fact, the ‘Percent observations’ of driver gender may be used to account for this mixed effect as coefficient’s value has a distribution of 20.05% versus 79.95% (above zero versus below zero) for the front vehicle, 14.01% versus 85.99% for the rear vehicle. It was also found that drivers tended to sustain more severe injuries when the airbag was not deployed or the seat belt was not in use. This is consistent with common sense and highlights the importance of seat belt and airbag checking (occupant protection systems).
With positive coefficient values, injury severities tended to increase with alcohol or drug use by the driver in the rear vehicle. However, this effect had not been found for the driver in front vehicle. This difference may be explained as follows: drivers who were involved in alcohol or drug were usually found to drive at a relatively higher speed, and speeding was one of the main causes of rear-end crashes with severe injuries.
From the aspect of vehicle information, the front vehicle whose manufacturer’s model year was before 2001 (1995 and before, from 1996 to 2000) and the rear vehicle whose manufacturer’s model year was before 1996 (1995 and before) were more likely to bring the drivers a serious injury. This effect had been found in a previous study, which showed that old vehicles had a significantly greater chance of being involved in an injury crash [35]. This may be explained by the lack of fine maintenance and advanced safety design on old vehicles. This factor has a mixed effect on injury severity for the rear vehicle according to the value of ‘Percent observations’.
Interesting results can be found regarding the effect of crash time. For drivers in the front vehicles, the injury severity of crashes that happened on Mondays tended to be higher. For drivers in the rear vehicles, more severe injuries might be found in crashes happened on a summer day. For both front and rear vehicles, drivers tended to suffer a severe injury during nighttime (beyond 7:00–18:00). Related conclusions about time factors’ effect on injury severity had also been obtained in previous studies. For instance, lower level of injury had been found in the crashes during the winter month and relative researches attributed this to lower speed and more cautious drivers on a winter day [11]; this may explain why lower injury severity would be sustained by drivers in rear vehicles on a snowy road surface condition as shown in Table 3. Summer season and nighttime had been found to be significantly associated with the injury severity, causing a positive effect on the possibility of the injuries [5]. Moreover, both drivers seemed to suffer a more severe injury when light level that existed at the time of the crash is daylight.
The traffic way related factor (VTRAFWAY), which focused on the value of coefficients only, it can be concluded that the level of injury severity for both the two drivers in rear-end crashes may increase under a two-way divided traffic flow. Moreover, the injury severity for drivers in front vehicles may also become more serious at the entrance or exit ramp where the sight distance is usually limited. As for the influence of traffic flow, drivers in rear vehicles may suffer a more severe injury at the one-way traffic flow. Interestingly, this factor’s effect is not absolutely positive or negative for the front vehicle according to the result of ‘Percent observations’ for this parameter. It is much better to say the effect is more possible to be positive (over a 50% probability).
For lane-related factors, the injury outcome of drivers in rear vehicles tends to be higher within a straight lane. The crash position reflected the regional characteristics of the crashes and the crash position had also been found to be related to the injury severity as presented in Table 3. A converse effect of crash position was shown on two drivers as the drivers in front vehicles may sustain a more severe injury in the Midwest of the country and the drivers in rear vehicles may suffer a lower level of injury in the Midwest and South of the country (the classification of region can be found in GES Analytical User’s Manual, for instance, Midwest includes OH, IN, IL, MI, WI, MN, ND, SD, NE, IA, MO, and KS of the United States). This phenomenon indicates that there may be a region difference on injury severity of the rear-end crash.
Apart from addressing the issue as to unobserved heterogeneity across observations, the random parameter bivariate ordered probit models also quantifies the correlation in the error terms between the models for each pair of crash-involved drivers by estimating the correlation parameter (𝜌). This index actually reflects the correlation in unobserved factors affecting injury severity in each rear-end crash, which may not be included in the recorded data. As presented in Table 3, the value of 𝜌 is positive suggesting that unobserved factors tend to jointly increase (or decrease) the level of injury sustained by two drivers involved in the same crash.

5. Conclusions

Rear-end crashes have become one of the main factors of human injury along with huge property damage; thus, it is necessary to figure out the possible cause of such kind of traffic accidents. This study examined the degree of injury sustained by drivers involved in the car-car rear-end crashes and the significant factors affecting injury severity. To address the within-crash correlation between the two drivers involved in the same crash and varying effects across observations, a random parameters bivariate ordered probit model was developed concerning the ordered probit model. Specifically, this model used two equations to describe injury severity for two drivers and proposed a correlation parameter to demonstrate correlation in unobserved factors between two drivers. The value of the correlation parameter suggests that these unobserved factors throw a joint increase (or joint decrease) effect on the degree of injury severity of the two drivers in the same crash. This indicates it is necessary to account for this within-crash correlation in model developing instead of heterogeneity in parameter effects only.
Specific effects of all significant variables on the injury for two drivers were also investigated and presented in this paper. For both vehicles, older age and female drivers tended to sustain a greater injury in rear-end crashes. The injury severity level with airbag or seat belt not in use becomes relative higher indicating the importance of a guaranteed occupant protection system. Furthermore, an older vehicle, nighttime, two-way traffic, and a light condition of daylight also caused positive effects on the injury severity for both two drivers. Although many common variables and consistent effects can be found, there are also some differences between the two drivers. For the drivers in rear vehicles, alcohol or drug use show positive effects on injury severity suggesting a necessary restrict on drivers’ alcohol or drug use. However, this effect has not been found on drivers in front vehicles. The outcome of injury for drivers in rear vehicles in a summer day or on a snowy road surface also increased, indicating a time-relation with the crash injury. Similar results have been presented about these effects on injury severity in previous studies [5,11,32,33,34]. One interesting result is that converse effects by crash position on two drivers can be found: when crashes occur in Midwest of the USA, drivers in front vehicles are prone to suffer a more severe injury, while drivers in rear vehicles sustain less injury. It has to be addressed that the paper presents the ‘Percent observations’ for random parameters, which may be used to explain the mixed effect for a certain factor.
Since rear-end crashes are related to many factors, this study gives an insight into several significant variables. Abundant as the rear-end crash data seem, the outcome of the current study is still limited with a certain dataset, which may cause a certain over-estimation. Exhaustive as the consideration of possible variables is, some potentially crucial information such as traffic volume has been neglected due to the lack of available data. However, the outcome of this study can still be used in some possible aspects like developing effective driver training, setting vehicle’s using year, offering risk information on road management, etc. There are also several promising extending of this work, e.g., examining different collision type, addressing temporal correlation in injury severity, or the injury of occupants in the same vehicle.

Author Contributions

Conceptualization, M.S. and X.M.; Data curation, M.S.; Formal analysis, M.S.; Funding acquisition, F.C.; Investigation, M.S.; Methodology, M.S.; Project administration, F.C.; Resources, X.M.; Supervision, F.C.; Validation, M.S. and X.M.; Visualization, F.C.; Writing—original draft, M.S.; Writing—review & editing, F.C. and X.M.

Funding

This research was jointly supported by Project 51578417 and U1733113 sponsored by the National Natural Science Foundation of China, and Project 18YF1424600 sponsored by Shanghai Sailing Program.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. World Health Organization. Global status report on road safety 2015. Inj. Prev. 2015, 15, 286. [Google Scholar]
  2. Dong, C.; Clarke, D.B.; Nambisan, S.S.; Huang, B. Analyzing injury crashes using random-parameter bivariate regression models. Transp. A Transp. Sci. 2016, 12, 794–810. [Google Scholar] [CrossRef]
  3. Chen, C.; Zhang, G.; Tarefder, R.; Ma, J.; Wei, H.; Guan, H. A multinomial logit model-Bayesian network hybrid approach for driver injury severity analyses in rear-end crashes. Accid. Anal. Prev. 2015, 80, 76–88. [Google Scholar] [CrossRef] [PubMed]
  4. Singh, S. Driver attributes and rear-end crash involvement propensity; DOT HS 809 540; National Technical Information Service: Springfield, Virginia, VA, USA, March 2003.
  5. Zhang, K.; Hassan, M. Identifying the Factors Contributing to Injury Severity in Work Zone Rear-End Crashes. J. Adv. Transport. 2019, 2019, 9. [Google Scholar] [CrossRef]
  6. Faller, R.K.; Magdaleno, J.A.; Post, E.R. Full-Scale Vehicle Crash Tests on the Iowa Retrofit Concrete Barrier Rail; Transportation Research Report TRP-03-015-88; University of Nebraska-Lincoln: Nebraska, NE, USA, January 1989. [Google Scholar]
  7. Xu, J.; Li, Y.; Lu, G.; Zhou, W. Reconstruction model of vehicle impact speed in pedestrian–vehicle accident. Int. J. Impact Eng. 2009, 36, 783–788. [Google Scholar] [CrossRef]
  8. Xu, J.; Shang, S.; Qi, H.; Yu, G.; Wang, Y.; Chen, P. Simulative investigation on head injuries of electric self-balancing scooter riders subject to ground impact. Accid. Anal. Prev. 2016, 89, 128–141. [Google Scholar] [CrossRef]
  9. Chen, F.; Chen, S.; Ma, X. Analysis of hourly crash likelihood using unbalanced panel data mixed logit model and real-time driving environmental big data. J. Saf. Res. 2018, 65, 153. [Google Scholar] [CrossRef]
  10. Wang, L.; Abdel-Atym, M.; Ma, W.; Hu, J.; Zhong, H. Quasi-vehicle-trajectory-based real-time safety analysis for expressways. Transport. Res. C-Emer. 2019, 103, 30–38. [Google Scholar] [CrossRef]
  11. Russo, B.J.; Savolainen, P.T.; Iv, W.H.S.; Anastasopoulos, P.C. Comparison of factors affecting injury severity in angle collisions by fault status using a random parameters bivariate ordered probit model. Anal. Methods Accid. Res. 2014, 2, 21–29. [Google Scholar] [CrossRef]
  12. Li, Z.; Liu, P.; Wang, W.; Xu, C. Using support vector machine models for crash injury severity analysis. Accid. Anal. Prev. 2012, 45, 478–486. [Google Scholar] [CrossRef]
  13. Ma, C.; Hao, W.; Xiang, W.; Yan, W. The Impact of Aggressive Driving Behavior on Driver-Injury Severity at Highway-Rail Grade Crossings Accidents. J. Adv. Transport. 2018, 58, 1–10. [Google Scholar] [CrossRef]
  14. Yuan, Q.; Lu, M.; Theofilatos, A.; Li, Y.-B. Investigation on occupant injury severity in rear-end crashes involving trucks as the front vehicle in Beijing area, China. Chin. J. Traumatol. 2017, 20, 20–26. [Google Scholar] [CrossRef] [PubMed]
  15. Konstantina, G.; Kinzenbaw, C.R.; Shauna, H.; Neal, H. An empirical analysis of farm vehicle crash injury severities on Iowa’s public road system. Accid. Anal. Prev. 2010, 42, 1392–1397. [Google Scholar]
  16. Manner, H.; Wünsch-Ziegler, L. Analyzing the severity of accidents on the German Autobahn. Accid. Anal. Prev. 2013, 57, 40–48. [Google Scholar] [CrossRef] [PubMed]
  17. Chen, F.; Chen, S. Injury severities of truck drivers in single- and multi-vehicle accidents on rural highways. Accid. Anal. Prev. 2011, 43, 1677–1688. [Google Scholar] [CrossRef] [PubMed]
  18. Kunt, M.M.; Aghayan, I.; Noii, N. Prediction for traffic accident severity: Comparing the artificial neural network, genetic algorithm, combined genetic algorithm and pattern search methods. Transport 2011, 26, 353–366. [Google Scholar] [CrossRef]
  19. Yan, X.; Radwan, E. Analyses of Rear-End Crashes Based on Classification Tree Models. Traffic Inj. Prev. 2006, 7, 276–282. [Google Scholar] [CrossRef]
  20. Das, A.; Abdel-Aty, M.A. A combined frequency–severity approach for the analysis of rear-end crashes on urban arterials. Saf. Sci. 2011, 49, 1156–1163. [Google Scholar] [CrossRef]
  21. Wang, W.; Liu, C.; Chen, D. Predicting driver injury severity in freeway rear-end crashes by support vector machine. In International Conference on Transportation, Changchun, China 16–18 December 2011; IEEE: Piscataway, NJ, USA, 2011. [Google Scholar]
  22. Savolainen, P.T.; Manneringm, F.L.; Lord, D.; Quddus, M.A. The statistical analysis of highway crash-injury severities: A review and assessment of methodological alternatives. Accid. Anal. Prev. 2011, 43, 1666–1676. [Google Scholar] [CrossRef]
  23. Wu, Y.; Abdel-Aty, M.; Cai, Q.; Lee, J.; Park, J. Developing an algorithm to assess the rear-end collision risk under fog conditions using real-time data. Transport. Res. C-Emer. 2018, 87, 11–25. [Google Scholar] [CrossRef]
  24. Zhao, P.; Lee, C. Assessing rear-end collision risk of cars and heavy vehicles on freeways using a surrogate safety measure. Accid. Anal. Prev. 2018, 113, 149. [Google Scholar] [CrossRef] [PubMed]
  25. Xie, Y.; Zhao, K.; Huynh, N. Analysis of driver injury severity in rural single-vehicle crashes. Accid. Anal. Prev. 2012, 47, 36–44. [Google Scholar] [CrossRef] [PubMed]
  26. Ma, X.; Chen, S.; Chen, F. Multivariate space-time modeling of crash frequencies by injury severity levels. Anal. Methods Accid. Res. 2017, 15, 29–40. [Google Scholar] [CrossRef]
  27. Zeng, Q.; Gu, W.; Zhang, X.; Wen, H.; Lee, J.; Hao, W. Analyzing freeway crash severity using a Bayesian spatial generalized ordered logit model with conditional autoregressive priors. Accid. Anal. Prev. 2019, 127, 87–95. [Google Scholar] [CrossRef] [PubMed]
  28. Chandler, S.; Duncan, A.J.K.; Forrest, M.; Khattak, A.J.; Council, F.M. Applying the ordered probit model to injury severity in truck-passenger car rear-end collisions. Transport. Res. Rec. 1998, 1635, 63–71. [Google Scholar]
  29. Abdel-Aty, M. Analysis of driver injury severity levels at multiple locations using ordered probit models. J. Saf. Res. 2003, 34, 597–603. [Google Scholar] [CrossRef]
  30. Jalayer, M.; Shabanpour, R.; Pour-Rouholamin, M.; Golshani, N.; Zhou, H. Wrong-way driving crashes: A random-parameters ordered probit analysis of injury severity. Accid. Anal. Prev. 2018, 117, 128–135. [Google Scholar] [CrossRef]
  31. Mannering, F.L.; Bhat, C.R. Analytic methods in accident research: Methodological frontier and future directions. Anal. Methods Accid. Res. 2014, 1, 1–22. [Google Scholar] [CrossRef]
  32. Abay, K.A.; Paleti, R.; Bhat, C.R. The joint analysis of injury severity of drivers in two-vehicle crashes accommodating seat belt use endogeneity. Transport. Res. B-Meth. 2013, 50, 74–89. [Google Scholar] [CrossRef] [Green Version]
  33. Chiou, Y.-C.; Hwang, C.-C.; Chang, C.-C.; Fu, C. Reprint of “Modeling two-vehicle crash severity by a bivariate generalized ordered probit approach”. Accid. Anal. Prev. 2013, 61, 97–106. [Google Scholar] [CrossRef]
  34. Chen, H.; Cao, L.; Logan, D.B. Analysis of risk factors affecting the severity of intersection crashes by logistic regression. J. Crash Prev. Inj. Control 2012, 13, 300–307. [Google Scholar] [CrossRef] [PubMed]
  35. Blows, S.; Ivers, R.Q.; Woodward, M.; Connor, J.; Ameratunga, S.; Norton, R. Vehicle year and the risk of car crash injury. Inj. Prev. 2003, 9, 353. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Summary statistics for driver-related, vehicle-related and other independent variables.
Figure 1. Summary statistics for driver-related, vehicle-related and other independent variables.
Ijerph 16 02632 g001
Table 1. Summary of injury severity for rear-end crash-involved drivers.
Table 1. Summary of injury severity for rear-end crash-involved drivers.
Injury Severity of Driver in the Front VehicleInjury Severity of Driver in the Rear VehicleTotal
L1L2L3L4
No injury (L1)872357658520910,093 (66.6%)
Possible injury (L2)1683548220802531 (16.7%)
Non-incapacitating injury (L3)1374188264471873 (12.3%)
Incapacitating or fatal injury(L4)4376059106662 (4.4%)
Total12,217 (80.6%)1372 (9.0%)1128 (7.4%)442 (2.9%)15,159 (100%)
Table 2. Summary statistics for driver-related, vehicle-related and other independent variables.
Table 2. Summary statistics for driver-related, vehicle-related and other independent variables.
VariablesDescription and ClassificationDriver in the Front VehicleDriver in the Rear VehicleTotal
L1L2L3L4L1L2L3L4
HOURDay: 7:00~18:0084372098153347210,193113089032725,080
Night: else165643334019020242422381155238
AGEAge 24 and below23034333108746324383661318700
Age between 25 and 63677018511395499667379063125818,867
Age 64 and above102024716876912144131532751
SEXMale4827952637267631257445418814,211
Female526615791236395590579867425416,107
ALCOHOLAlcohol involved22848275504837452
Alcohol not involved10,0712523186965411,9421322108040529,866
DRUGSDrugs use712470202623153
No drugs use10,0862530187165812,1471352110241930,165
AIR_BAGAirbag not deployed26015716512418965885271743891
Airbag deployed98332374170853810,32178460126826,427
REST_USESeat belt equipment not in use62352417104305251375
Seat belt equipment in use10,0312496184964512,1131342107639129,943
VALIGNRoadway alignment: Straight95392440178762011,5471326109242128,772
Curve5549186426704636211546
REGIONNortheast21646429014224733571041046076
Midwest17973135121392306142225885522
South497111551122339607468262520615,174
West1161421149421364191174443546
SPEEDRELDriver’s speed related10917111024923342811153369
Driver’s speed not related9984251418626529725103884732726,949
LGTCONDaylight7785194414444439373106786331323,232
Dark195451437719124312562311186072
Dawn or Dusk3547352284134934111014
WEATHERClear716418361271482862797082233421,506
Rain115425419368139614695323338
Snow, Freezing Rain or Drizzle1353324917411106402
Fog, Smog, Smoke and so on258303033072
Cloudy16154003821031990242198705000
INT_HWYCrash occurred on an interstate way144236029916917651981921154540
Crash not occurred on an interstate way86512171157449310,452117493632725,778
MDLYRManufacturer’s model year: before 199536512510339873100119301754
1996–20001113361322872303299208854778
2001–2005240066949420536624063391238298
2006–2010372681859419335793802961379723
2011–201524895583601381800187166675765
BDYTYPThe car has four doors7152185713734918834101378712921,636
The number of car’s doors is not four294167450017133833593413138682
WRKZONEWithin the boundaries of a work zone191404612228192814578
Not within the boundaries of a work zone99022491182765011,9891353110042829,740
Table 3. Results of random parameter bivariate ordered probit models.
Table 3. Results of random parameter bivariate ordered probit models.
VariablesCoefficientStandard ErrorZ ValuePercent Observations
Above 0Below 0
Driver in the Front Vehicle
Constant−0.2028 ***0.0677−3.0
▲ AGE: ≤24−0.3182 ***0.0269−11.82.74%97.26%
  Standard deviation0.1657 ***0.02387.0
▲ SEX: Male−0.4426 ***0.0223−19.920.05%79.95%
  Standard deviation0.5292 ***0.017729.8
▲ AIRBAG: Not deployed0.8767 ***0.044919.599.95%0.05%
  Standard deviation0.2692 ***0.04326.2
▲ SEAT BELT: Not in use0.5594 ***0.10255.5
  SPEEDREL: Speed related−0.3872 ***0.1163−3.3
▲ MDLYR: 2001–2005−0.1392 ***0.0320−4.37.08%92.92%
  Standard deviation0.0945 ***0.02034.7
▲ MDLYR: 2006–2010−0.3648 ***0.0307−11.90.96%99.04%
  Standard deviation0.1559 ***0.01788.8
▲ MDLYR: 2011–2015−0.3569 ***0.0331−10.80.23%99.77%
  Standard deviation0.1255 ***0.02165.8
BDYTYP: Four doors0.0839 ***0.02323.676.42%23.58%
  Standard deviation0.1164 ***0.01219.6
WEEK: Monday0.0592 **0.02912.0
▲ HOUR: 7:00~18:00−0.1753 ***0.0370−4.70.69%99.31%
  Standard deviation0.0713 ***0.01146.3
▲ VTRAFWAY: Two-way, divided no median0.1414 ***0.02924.968.44%31.56%
  Standard deviation0.2967 ***0.024212.3
VTRAFWAY: Entrance/Exit Ramp0.0971 *0.05881.757.53%42.47%
  Standard deviation0.4995 ***0.05988.4
REGION: Midwest0.0876 ***0.02733.256.75%43.25%
  Standard deviation0.5032 ***0.025719.6
▲ Light Condition: Daylight0.1237 **0.05772.2100%0
  Standard deviation0.0306*** 0.01182.6
Light Condition: Dark0.1215 **0.06082.0
μ 1 0.6191 ***0.011354.8
μ 2 1.4814 ***0.020273.3
Driver in the rear vehicle
Constant|−0.7899 ***0.0921−8.580100%
  Standard deviation0.2067 ***0.012316.78
▲ AGE: ≤24−0.4894 ***0.0454−10.798.85%91.15%
  Standard deviation0.3615 ***0.022116.36
AGE: 25–63−0.2780 ***0.0426−6.530100%
  Standard deviation0.0754 ***0.01664.54
▲ SEX: Male−0.3822 ***0.0261−14.6714.01%85.99%
  Standard deviation0.3554 ***0.019218.49
ALCOHOL: Involved0.2452 ***0.06883.56
DRUGS: Involved0.5898 ***0.11105.3181.06%18.94%
  Standard deviation0.6645 ***0.11285.89
▲ AIRBAG: Not deployed1.0486 ***0.026739.3
▲ SEAT BELT: Not in use1.1682 ***0.081814.2891.77%8.23%
  Standard deviation0.8390 ***0.09069.27
MDLYR: 1995–2000−0.1192 **0.0531−2.256.18%93.82%
  Standard deviation0.0774 ***0.02812.75
▲ MDLYR: 2001–2005−0.1669 ***0.0506−3.317.36%82.64%
  Standard deviation0.1782 ***0.02267.88
▲ MDLYR: 2006–2010−0.2833 ***0.0517−5.4825.78%74.22%
  Standard deviation0.4353 ***0.024417.81
▲ MDLYR: 2011–2015−0.2049 ***0.0554−3.70
SEASON: Jun, Jul, Aug0.0933 ***0.02913.2094.74%5.26%
  Standard deviation0.0575 **0.02462.34
▲ HOUR: 7:00~18:00−0.2304 ***0.0451−5.11
VSURCOND: Snow−0.6991 ***0.2430−2.88
  Standard deviation0.5531 **0.23442.36
▲ VTRAFWAY: Two-way, divided no median0.1426 ***0.036-3.9798.98%1.02%
  Standard deviation0.0613 **0.02962.07
VTRAFWAY: Two-way, divided, median0.2466 ***0.03038.14
VTRAFWAY: One-way0.2140 ***0.08152.62
VALIGN: Straight0.1703 ***0.06092.8
REGION: Midwest−0.1221 ***0.0375−3.25
REGION: South−0.0916 ***0.0302−3.04
▲ Light Condition: Daylight0.1929 ***0.04154.65
μ 1 0.5090 ***0.012939.56
μ 2 1.3173 ***0.024354.17
ρ (correlation parameter)0.2440 ***0.014317.06
Final log-likelihood−23,365
Note: ***, **, * mean Significance at 1%, 5%, 10% level. ▲ means this variable is significant for both drivers.

Share and Cite

MDPI and ACS Style

Chen, F.; Song, M.; Ma, X. Investigation on the Injury Severity of Drivers in Rear-End Collisions Between Cars Using a Random Parameters Bivariate Ordered Probit Model. Int. J. Environ. Res. Public Health 2019, 16, 2632. https://doi.org/10.3390/ijerph16142632

AMA Style

Chen F, Song M, Ma X. Investigation on the Injury Severity of Drivers in Rear-End Collisions Between Cars Using a Random Parameters Bivariate Ordered Probit Model. International Journal of Environmental Research and Public Health. 2019; 16(14):2632. https://doi.org/10.3390/ijerph16142632

Chicago/Turabian Style

Chen, Feng, Mingtao Song, and Xiaoxiang Ma. 2019. "Investigation on the Injury Severity of Drivers in Rear-End Collisions Between Cars Using a Random Parameters Bivariate Ordered Probit Model" International Journal of Environmental Research and Public Health 16, no. 14: 2632. https://doi.org/10.3390/ijerph16142632

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop