Exploring the modeling and site-ranking performance of Bayesian spatiotemporal crash frequency models with mixture components

A full Bayesian multilevel approach for modeling interaction effects in single-vehicle crashes

Interaction effects constitute crucial crash attributes that can be classified into two distinct categories: spatiotemporal interactions and factor interactions. These interactions are rarely addressed systematically in modeling the severity of single-vehicle (SV) crashes. This study focuses on uncovering these crash attributes by designing a full Bayesian spatiotemporal interaction multilevel logit (STIML-logit) approach with heterogeneity in means and variances (HMV). Meanwhile, a nested Gaussian conditional autoregressive (CAR) structure is proposed to fit the spatiotemporal interaction component and its effectiveness is verified by calibrating four different interaction patterns. A standard multilevel logit (with and without HMV), a multilevel logit with HMV, and a spatiotemporal multilevel logit with HMV are constructed for comparison. Risk factors are decomposed into traffic environment factors (group level) and individual crash factors (case level) to construct a multilevel structure and to capture possible interactions between risk factors from different levels (cross-level factor interactions). We perform regression modeling utilizing SV crash cases covering 96 major urban roads in Shandong, China. The modeling results underscore several significant findings: (1) the STIML-logit with HMV demonstrates the best regression performance, suggesting that systematically dealing with the interaction effects and the HMV is a trustworthy modeling perspective; (2) crash models with the nested CAR outperform those with the traditional CAR and the result is supported by all the spatiotemporal statistical functions, highlighting the potential advantages of the nested structure; (3) all the environment factors maintain significant interactions with the case factors, highlighting that the contribution of the environment factors to crash injuries is not constant but is rather influenced by the specific case-related crash factors. The study introduces a promising regression architecture for modeling crash injuries and revealing subtle crash attributes.

Road Accident Hotspots on Jordan’s Highway Based on Geometric Designs Using Structural Equation Modeling

One of the primary objectives of transportation engineering is to increase the safety of road infrastructure. This study seeks to determine the relationship between geometric design parameters in relation to road accident criteria based on accident hotspots on Jordan’s Desert Highway. The road accident data (from 2016 to 2019) were collected from the Jordan Traffic Department. The spatial pattern of hotspots was identified using a GIS tool named Getis-Ord Gi* based on the severity index of road accidents. A topographic survey was conducted to investigate the road alignment and intersections at hotspot locations. The study utilized the Structural Equation Modeling (SEM) technique via SmartPLS to highlight the correlation between geometric designs in relation to road accidents. The hotspot analysis (Gits-Ord Gi) discovered 80 road accident hotspots along the highway. The study found that horizontal alignment and road intersections significantly impact road accidents in hotspot locations. Furthermore, vertical alignment has no effect on road accidents in hotspot areas. The study enhanced the comprehension of the factors associated with road geometrics and intersections that affect the occurrence of road accidents.

Modeling spatiotemporal interactions in single-vehicle crash severity by road types

INTRODUCTION

Spatiotemporal correlations have been widely recognized in single-vehicle (SV) crash severity analysis. However, the interactions between them are rarely explored. The current research proposed a spatiotemporal interaction logit (STI-logit) model to regression SV crash severity using observations in Shandong, China.

METHOD

Two representative regression patterns-mixture component and Gaussian conditional autoregression (CAR)-were employed separately to characterize the spatiotemporal interactions. Two existing statistical techniques-spatiotemporal logit and random parameters logit-were also calibrated and compared with the proposed approach with the aim of highlighting the best one. In addition, three road types-arterial road, secondary road, and branch road-were modeled separately to clarify the variable influence of contributors on crash severity.

RESULTS

The calibration results indicate that the STI-logit model outperforms other crash models, highlighting that comprehensively accommodating spatiotemporal correlations and their interactions is a recommended crash modeling approach. Additionally, the STI-logit using mixture component fits crash observations better than that using Gaussian CAR and this finding remains stable across road types, suggesting that simultaneously accommodating stable and unstable spatiotemporal risk patterns can further strengthen model fit. According to the significance of risk factors, there is a significant positive correlation between distracted diving, drunk driving, motorcycle, dark (without street lighting), and collision with fixed object and serious SV crashes. Truck and collision with pedestrian significantly mitigate the likelihood of serious SV crashes. Interestingly, the coefficient of roadside hard barrier is significant and positive in branch road model, but it is not significant in arterial road model and secondary road model.

PRACTICAL APPLICATIONS

These findings provide a superior modeling framework and various significant contributors, which are beneficial for mitigating the risk of serious crashes.

Modelling injury severity in single-vehicle crashes using full Bayesian random parameters multinomial approach

Single-vehicle (SV) crash severity model considering spatiotemporal correlations has been extensively investigated, but spatiotemporal interactions have not received sufficient attention. This research is dedicated to propose a superior spatiotemporal interaction correlated random parameters logit approach with heterogeneity in means and variances (STICRP-logit-HMV) for systematically characterizing unobserved heterogeneity, spatiotemporal correlations, and spatiotemporal interactions. Four flexible interaction formulations are developed to uncover the spatiotemporal interactions, including linear structure, Kronecker product, mixture-2 model, and mixture-5 model. Four candidate approaches-random parameters logit (RP-logit), RP-logit with heterogeneity in means and variances (RP-logit-HMV), correlated RP-logit-HMV (CRP-logit-HMV), and spatiotemporal CRP-logit-HMV (STCRP-logit-HMV)-are also established and compared with the proposed model. SV crash observations in Shandong Province, China, are employed to calibrate regression parameters. The model comparison results show that (1) the performance of the RP-logit-HMV model outperforms the RP-logit model, implying that capturing heterogeneity in the means and variances can strengthen model fit; (2) the CRP-logit-HMV model and the RP-logit-HMV model are comparable; (3) the STCRP-logit-HMV model outperforms the CRP-logit-HMV model, implying that addressing the spatiotemporal crash mechanisms is beneficial to the overall fitting of the crash model; (4) the STICRP-logit-HMV model performs better than the STCRP-logit-HMV model and this finding remains stable across different interaction formulations, indicating that comprehensively reflecting the spatiotemporal correlations and their interactions is a promising approach to model SV crashes. Among the four interaction models, the STICRP-logit-HMV model with mixture-5 component maintains the best fit, which is a recommended approach to model crash severity. The regression coefficients for young driver, male driver, and non-dry road surface are random across observations, suggesting that the influence of these factors on SV crash severity maintains significant heterogeneity effects. The research results provide transportation professionals with a superior statistical framework for diagnosing crash severity, which is beneficial for improving traffic safety.

Incorporating driving volatility measures in safety performance functions: Improving safety at signalized intersections

About 40 percent of motor vehicle crashes in the US are related to intersections. To deal with such crashes, Safety Performance Functions (SPFs) are vital elements of the predictive methods used in the Highway Safety Manual. The predictions of crash frequencies and potential reductions due to countermeasures are based on exposure and geometric variables. However, the role of driving behavior factors, e.g., hard accelerations and declarations at intersections, which can lead to crashes, are not explicitly treated in SPFs. One way to capture driving behavior is to harness connected vehicle data and quantify performance at intersections in terms of driving volatility measures, i.e., rapid changes in speed and acceleration. According to recent studies, driving volatility is typically associated with higher risk and safety-critical events and can serve as a surrogate for driving behavior. This study incorporates driving volatility measures in the development of SPFs for four-leg signalized intersections. The Safety Pilot Model Deployment (SPMD) data containing over 125 million Basic Safety Messages generated by over 2,800 connected vehicles are harnessed and linked with the crash, traffic, and geometric data belonging to 102 signalized intersections in Ann Arbor, Michigan. The results show that including driving volatility measures in SPFs can reduce model bias and significantly enhances the models' goodness-of-fit and predictive performance. Technically, the best results were obtained by applying Bayesian hierarchical Negative Binomial Models, which account for spatial correlation between signalized intersections. The results of this study have implications for practitioners and transportation agencies about incorporating driving behavior factors in the development of SPFs for greater accuracy and measures that can potentially reduce volatile driving.

Application of Extremely Randomised Trees for exploring influential factors on variant crash severity data

Crash severity models play a crucial role in evaluating the influencing factors in the severity of traffic crashes. In this study, Extremely Randomised Tree (ERT) is used as a machine learning technique to analyse the severity of crashes. The crash data in the province of Khorasan Razavi, Iran, for a period of 5 years from 2013 to 2017, is used for crash severity model development. The dataset includes traffic-related variables, vehicle specifications, vehicle movement, land use characteristics, temporal characteristics, and environmental variables. In this paper, Feature Importance Analysis (FIA), Partial Dependence Plots (PDP), and Individual Conditional Expectation (ICE) plots are utilised to analyse and interpret the results. According to the results, the involvement of vulnerable road users such as motorcyclists and pedestrians alongside traffic-related variables are among the most significant variables in crash severity. Results show that the presence of motorcycles can increase the probability of injury crashes by around 30% and almost double the probability of fatal crashes. Analysing the interaction of PDPs shows that driving speeds above 60 km/h in residential areas raises the probability of injury crashes by about 10%. In addition, at speeds higher than 70 km/h, the presence of pedestrians approximately increases the probability of fatal crashes by 6%.

An accelerated hierarchical Bayesian crash frequency model with accommodation of spatiotemporal interactions

Although spatial and temporal correlations of crash observations have been well addressed in the literature, the interactions between them are rarely studied. This study proposes a Bayesian spatiotemporal interaction (BSTI) approach for crash frequency modeling with an integrated nested Laplace approximation (INLA) method to greatly expedite the Bayesian estimation process. Manhattan, which is the most densely populated urban area of New York City, is selected as the study area. Hexagons are used as the basic geographic units to capture crash, transportation, land use, and demo-economic data from 2013 to 2019. A series of Bayesian models with various spatiotemporal specifications are developed and compared. The BSTI model with Type II interaction, which assumes that the structured temporal random effect interacts with the unstructured spatial random effect is found to outperform the others in terms of goodness-of-fit and the ability to reduce the dependency of residuals. It is also found that the unobserved heterogeneity is mostly attributed to the spatial effects instead of temporal effects. In addition, the BSTI Type II model also yields the lowest predictive error when the last year's data are used as the test set. The proposed BSTI approach can potentially advance safety analytics by achieving high prediction accuracy and computational efficiency while maintaining its interpretability on the effects of contributing factors and the unobserved heterogeneity.

A new spatial count data model with Bayesian additive regression trees for accident hot spot identification

The identification of accident hot spots is a central task of road safety management. Bayesian count data models have emerged as the workhorse method for producing probabilistic rankings of hazardous sites in road networks. Typically, these methods assume simple linear link function specifications, which, however, limit the predictive power of a model. Furthermore, extensive specification searches are precluded by complex model structures arising from the need to account for unobserved heterogeneity and spatial correlations. Modern machine learning (ML) methods offer ways to automate the specification of the link function. However, these methods do not capture estimation uncertainty, and it is also difficult to incorporate spatial correlations. In light of these gaps in the literature, this paper proposes a new spatial negative binomial model which uses Bayesian additive regression trees to endogenously select the specification of the link function. Posterior inference in the proposed model is made feasible with the help of the Pólya-Gamma data augmentation technique. We test the performance of this new model on a crash count data set from a metropolitan highway network. The empirical results show that the proposed model performs at least as well as a baseline spatial count data model with random parameters in terms of goodness of fit and site ranking ability.