Abstract

In the corresponding research available, the safety impact remains controversial in implementing signal coordination on arterials, which calls for an in-depth exploration with the appropriate statistical methods. Based on the traffic data from Ann Arbor City (Michigan, USA), the paper proposes a safety evaluation model considering the multiple heterogeneities. In terms of arterials with the coordinated signalization, modeling results imply that (1) the multivariate heterogeneity shows the strongest interaction on crash frequency, followed by the spatiotemporal and structural heterogeneities, and (2) the spatial variation is unrelated to the temporal change among crashes in the denoted traffic analysis zones (TAZs). In an attempt to alleviate the coupled crash risks along the coordinated arterials, the study emphasizes the necessity of dividing the subcontrol traffic areas in real time according to the correlative degree of crash distribution. Meanwhile, the modeling framework with multiple heterogeneities can be applied for the safety analysis of other urban roads.

1. Introduction

Arterials obtain growing concerns with the increase of the major traffic flow among urban road networks. For ensuring operational efficiency, the control technique of signal coordination was widely applied along the arterials in recent decades [13]. By coupling the intersections simultaneously along arterials, signal coordination generates the “green wave” by optimizing the signal attributes of phase, timing, and cycle length [4]. The advantages of coordinated signalization are verified from several aspects, such as increasing the traffic output, alleviating the vehicular delay, and decreasing the phenomenon of “stop and go” [3].

The previous study focuses mainly on the efficiency improvement of arterials with coordinated signalization. It should be noted that the efficiency and safety impacts are not consistent. Signal coordination may induce safety hazards along arterials with the efficiency optimization. Since safety is the precondition of sustainable operation, it is imperative to provide an in-depth safety exploration of coordinated arterials with appropriate models.

The safety evaluation model is a useful tool to identify the crash risks [5, 6]. During the modeling procedure, the heterogeneity should be nonnegligible on account of the connected traffic flow along signal-coordinated arterials [7]. Heterogeneous correlations, posing an important effect on safety analysis results, stand for the interference involved in the data source of coordinated arterials (Figure 1). From the micro-aspect (i.e., on the roadway segments/or at the signalized intersections), there is the multivariate heterogeneity between injury and property damage only (PDO) crashes, since crashes from different severities may show the diverse contributors [8]. Meanwhile, signal coordination results in a strong association among vehicles with the “green wave,” thus forming the structural heterogeneity between crashes in each traffic analysis zone (TAZ) at the meso-viewpoint [9]. Regarding the macrolevel (i.e., along the whole coordinated arterials), it calls for the need to distinguish the specific form of spatiotemporal heterogeneity among the TAZs [10].

The innovation of the paper is to provide the comprehensive safety analysis of signal-coordinated arterials with multiple heterogeneities. The paper begins with the literature review in Section 2. Section 3 divides the signal-coordinated arterials into roadway segments, signalized intersections, and TAZs. Each corresponding area retrieves the data of road design features, traffic characteristics, and crash counts. Section 4 describes seven statistical models, including Poisson log-normal (PLN) model, multivariate Poisson log-normal (MPLN) model, joint multivariate Poisson log-normal (JMPLN) model, joint multivariate Poisson log-normal conditional autoregressive (JMPLN-CAR) model, JMPLN-CAR model with fixed spatial and linear temporal variation, JMPLN-CAR model with fixed spatial and linear-quadratic temporal variation, and JMPLN-CAR model with time-varying spatial interaction. To compare the predictable performance of models, the paper adopts three kinds of criteria, i.e., classical criteria (mean absolute deviance (MAD), the weighted mean absolute deviance (WMAD), mean squared prediction error (MSPE), and weighted mean squared prediction error (WMSPE)), prevalent criteria (deviance information criterion (DIC) and total deviance information criterion (TDIC)), and emerging criteria (log pseudo marginal likelihoods (LPML) and total log pseudo marginal likelihoods (TLPML)). In addition, variation proportion is feasible to reveal the modeling fitness explained by the heterogeneous correlations. With the evaluated results (Section 5), safety improvement suggestions are outlined in Section 6.

2. Literature Review

2.1. Efficiency Optimization of Signal-Coordinated Arterials

The research of signal coordination draws surging attention from scholars. In contrast to previous control modes (e.g., the fixed period and isolated signal), the coordinated signalization systematically integrates arterials as a whole with advanced approaches (e.g., the coevolutionary analysis, sensor technology, automatic control, and information theory) [11], thus significantly improving the movement efficiency of vehicles.

Specifically, the efficiency optimization of signal coordination had two aspects: (1) setting the green bandwidth to increase the traffic output [12] and (2) quantifying the performance index, such as reducing the vehicular delay, decreasing the queue length, considering the product of traffic flow, and improving the traveling speed [13]. To achieve the mentioned optimization objects, several algorithms were provided of the Max-Band control [14], forward-backward [15, 16], fuzzy control, game theory [17], ant colony [18], hierarchical and distributed control [19], and linear programming of mixed integers [14].

In summary, the previous literature focuses mainly on the efficiency improvement of signal coordination. It needs to provide the corresponding safety evaluation since the efficiency and safety impacts are not consistent.

2.2. Safety Evaluation of Signal-Coordinated Arterials

Signal coordination is able to optimize the moving efficiency of vehicles, while its safety impact is controversial in the current studies. One viewpoint is that the safety status improves significantly after applying the signal coordination on arterials. Carter and St-Onge [20] verified the 4.9% crash reduction of Phoenix city (Arizona, USA) under the operation of signal-coordinated arterials. Meanwhile, the unsafe behaviors were decreased in the study of Khattak et al. [21] due to the existence of consecutively vehicular platoons with signal coordination. Moreover, scholars identified the decrease of pedestrian, right-turn, and right-angle crashes with the control of signal coordination in the before-after study [22, 23]. An explanation was that signal coordination mitigates the phenomenon of “stop and go,” therefore decreasing the possibility of crash occurrence [24, 25].

In contrary to the above findings, there is a conclusion that signal coordination incurs safety hazards along arterials. For instance, crashes occurring at the signalized intersections were migrated to the roadway segments due to the performance of signal coordination [9]. To alleviate the vehicular delay, signal coordination was also allocated the high-speed limit on arterials. However, this countermeasure was potentially correlated with severe injury outcomes [26]. When vehicles entered the signalized intersections at the end of “green wave,” they might execute the unsafe behaviors of speeding and frequent lane changing [27]. In addition, there was a longer vehicular delay on the adjacent uncoordinated arterials, where the drivers became inpatient and might engage in “road rage” [8].

In conclusion, the safety study of signal coordination deserves an in-depth exploration with effective evaluation models.

2.3. Statistical Models of Crash Frequency considering the Heterogeneous Correlations

There were a number of classical safety evaluation models of crash frequency, such as the Poisson regression model (the primary Bayesian counting model) [28], negative binominal/Poisson-Gamma regression model (solving the issue of overdispersion) [29], and zero-inflated Poisson regression model (considering the existence of zero crashes in units) [6]. Based on the traditional approaches, several models revealed the heterogeneity to ensure the stability of regression results, including the random effects model, finite mixture/latent class model, Markov switching model, and generalized linear model [30]. Alternatively, the random effects model is more flexible to reflect the heterogeneous differentiation as the random term [30, 31]. With the method of maximum likelihood estimation, the model could acquire a better goodness-of-fit for variables.

2.3.1. Multivariate Heterogeneity

Crashes of different types may show the diverse contributors [32, 33]. For achieving the safety analysis at the microscopic aspect, crashes are usually classified into several levels as injury, PDO, and fatality. The multivariate random effects model is able to reflect heterogeneous correlations among crashes with different severities [30, 34].

2.3.2. Structural Heterogeneity

The structural heterogeneity exists between crashes at units with different road design and traffic flow features (e.g., roadway segments and signalized intersections along the coordinated arterials) [9]. Scholars concluded that the random effects model suffices to address the structural heterogeneity [30, 35]. However, the structural relationship between the units was not truly formulated. With the joint random effects model, Fan et al. [9] recently rebuilt the structural association between the roadway segments and signalized intersections.

2.3.3. Spatiotemporal Heterogeneity

In contrast to the structural correlation, spatial heterogeneity derives from crashes in units with similar characteristics [36, 37]. The analysis of spatial heterogeneity has attracted a wide popularity in the aspects of intersections [38], segments [36, 39], TAZs [33, 40], corridors [41], counties [26], and states/provinces [42]. Generally, the spatial heterogeneity in the random effects model was reflected as the distance function of homogeneous units [10].

There are the analyses of temporal variation among multiple years of crash data. For instance, Cheng et al. [43] explored the spatial and temporal heterogeneities (i.e., spatiotemporal heterogeneity) in the same model, where the specific forms included the fixed spatial and linear temporal association, fixed spatial and linear-quadratic temporal association, and time-varying spatial association.

3. Data

3.1. Description of Signal-Coordinated Arterials

The data records six arterials with the signal coordination from Ann Arbor City (Michigan State, USA) [44], located in Eisenhower, Packard, Plymouth, Stadium, State, and Washtenaw Roads. Since 2009, the signal-coordinated arterials apply the Adaptive Control System (ATCS), which poses the advantages of presenting the flexible signal plan, reflecting the timely traffic status, and providing the reasonable interference of signal coordination [45]. More specifically, ATCS detects the real-time traffic status (i.e., occupancy, flow, and congestion) with the use of monitors installed upstream of signalized intersections. The detected data serves to predict the vehicular queuing length, thus optimizing the corresponding signal attributes. It is mentioned that the relative signal attributes (cycle, timing, and phase) could not be directly acquired from the arterials due to their changeable features with the time variation.

Along the coordinated arterials, we investigated 38 roadway segments and 38 signalized intersections. For the purpose of differentiating the specific location, each signalized intersection locates in the area consisting of the central district and safety affecting region (i.e., 250 feet away from the stop bar) (Figure 1) [36, 39]. The roadway segment lies in the section between two adjacent intersections. Every TAZ is formulated by the signalized intersection and its adjacent roadway segment. With the distinguished units, the paper extracts the dependent variables, road design features, and traffic characteristics, respectively (Table 1).

Google Earth 2014 is capable of pinpointing the latitude and longitude intervals for the 38 roadway segments (or signalized intersections). With the geographical information, nine-year (2009–2014) crash data were retrieved from the Michigan Department of Transportation (MDOT). The total crash count, excluding the crashes related to animals, nonmotor vehicles (e.g., bicycles and motorcycles), and pedestrians, is 5613 along the coordinated arterials composed by 1996 crashes of roadway segments and 3647 crashes of signalized intersections. There are 1576/2989 PDO crashes and 387/657 injury crashes occurring on the segments/intersections. The fatal crashes are filtered due to the small sample size.

Road design variables, revealing operational requirements of the signal-coordinated arterials, are obtained by the road survey with Google Earth 2014. Meanwhile, the paper acquires the related traffic characteristics with the software of ArcGIS and the traffic volume map of the Southeast Michigan Council of Governments (SEMCOG) [46]. The related traffic variables suffice to describe the development status of arterials under the coordinated signalization.

3.2. Data Processing

According to the conclusion of Alarifi et al. [39], multicollinearity exists among the variables. The potential disturbance involved in the data may lead to unstable results during the modeling process. For alleviating the multicollinear effect, we employ the techniques of product-moment and Spearman rank-order [47]. When the correlation coefficient is larger than 0.6 between two estimated variables, only the more important one (i.e., the variable reflecting the characteristic of arterials more in detail) will remain. The filtered variables are in italics in Table 1.

4. Methodology

4.1. Model Configuration

There are seven models in the paper, including four existing models (i.e., PLN model, JMPLN model, JMPLN-CAR model, and MPLN model) and three newly formulated models (i.e., JMPLN-CAR model with fixed spatial and linear temporal variation, JMPLN-CAR model with fixed spatial and linear-quadratic temporal variation, and JMPLN-CAR model with time-varying spatial interaction). Detailed descriptions are summarized in Table 2.

According to Table 2, the comparison of models 1-2 verifies the existence of multivariate heterogeneity between PDO and injury crashes. Similarly, the paper could compare models 2-3, 3-4, and 5–7 to explore the structural, spatial, and spatiotemporal heterogeneities, respectively.

4.1.1. Poisson Log-Normal Model (Model 1)

The framework of model 1 is specified as [34]where stands for the crash count of unit ( = 1, 2, …, ;  = 38 is the total number of units) for severity level ( = 1, 2, …, ;  = 1 for injury and  = 2 for PDO) that belongs to structure ( = 1 for roadway segment and  = 2 for signalized intersection). represents the expectation of .

The probability of crash counts satisfies the formula as follows:

further follows the distribution aswhere is the intercept of severity level and structure ; reflects the regression vector related to the explanatory vector ; stands for the dichotomous variable, equaling to the length of segment if  = 1 and equaling to 0 if  = 2; ; and [8].

4.1.2. Multivariate Poisson Log-Normal Model (Model 2)

We havewhere indicates the multivariate-heterogeneous effect between injury and PDO crashes [8]. Also, is the corresponding matrix of variance-covariance denoted bywhere is the inverse matrix of ; represents the scale matrix of Jonathan et al. [48]; and is the freedom degree, assigned as 2 according to the classification of structure . Model 2 further discusses the multivariate heterogeneity on roadway segments (or at signalized intersections) in the microscopic aspect compared to model 1.

4.1.3. Joint Multivariate Poisson Log-Normal Model (Model 3)

Compared to model 2, model 3 incorporates the structural-heterogeneous correlation in each TAZ from the mesoscopic level. Inspired by the study of Fan et al. [9], model 3 follows the formula aswhere is the vector of regression parameters related to the explanatory vector of TAZ and represents the term of structural heterogeneity between each segment and its adjacent intersection [9]. More specifically, is the structural heterogeneity from the intersections acting on the segments, while  = () stands for the structural heterogeneity from the segments acting on the intersections; ( = 1, 2, 3) reflects the estimated coefficient. satisfies and is the precision parameter [33].

4.1.4. Joint Multivariate Poisson Log-Normal Conditional Autoregressive Model (Model 4)

Model 4 further considers the spatial-heterogeneous correlation among the TAZs. The corresponding equation is shown aswhere stands for the CAR term of spatial relationship. The specific distribution of is assigned aswhere is the set derived from [39]; reflects the CAR precision parameter assumed to be the gamma hyper-prior (0.01, 0.01); and is represented by the matrix , where each entry is the spatial heterogeneity between units and . The formula of is given aswhere is shown as the geometry centroid length and (>0) indicates the decline rate of distance [49].

Derived from model 4, models 5–7 identify the different types of spatiotemporal-heterogeneous correlations in Sections 4.1.54.1.7, respectively.

4.1.5. Joint Multivariate Poisson Log-Normal Conditional Autoregressive Model with Fixed Spatial and Linear Temporal Variation (Model 5)

where is the crash count of divided by time ( = 1, 2, …, ); reveals the spatiotemporal heterogeneity, consisting of the fixed spatial correlation and the linear temporal variation ; is the linear trend of time; is the average variation among all units; and follows the similar distribution as .

In model 5, it is emphasized that the spatial variation ( and ) is independent of the temporal changing () [43].

4.1.6. Joint Multivariate Poisson Log-Normal Conditional Autoregressive Model with Fixed Spatial and Linear-Quadratic Temporal Variation (Model 6)

We have

Model 6 addresses the fixed spatial and linear-quadratic temporal variation () based on model 5. In equation (12), reflects the quadratic temporal variation; and is the coefficient corresponding to the quadratic time .

and are unrelated to and , which indicates that the temporal trend does not affect the spatial connectivity and vice versa [10].

4.1.7. Joint Multivariate Poisson Log-Normal Conditional Autoregressive Model with Time-Varying Spatial Interaction (Model 7)

Model 7 presents the spatiotemporal interaction between the spatial and temporal heterogeneity in reality [10]:where the time-varying spatial interaction satisfies the distribution as . is the variance-covariance matrix expressed as is the inverse matrix of ; is the scale matrix defined by Aguero-Valverde and Jovanis [50]; and stands for the degree of freedom [10].

4.2. Bayesian Inference and Model Comparison

Three types of methods are feasible to identify the differentiation among the models, including the classical criteria (i.e., WMAD/MAD and WMSPE/MSPE), pervasive criteria (i.e., TDIC/DIC), and emerging criteria (i.e., TLPML/LPML). The summary of the criteria is presented in Table 3.

4.2.1. Mean Absolute Deviance and Mean Squared Prediction Error

The classical approaches of and are employed to compare the model fitness [41, 52, 53]. Regarding the structure , they are separately presented in the following equations:

For each TAZ, the weighted (i.e., ) is shown as [52]where and are related to roadway segments and signalized intersections, respectively.

In terms of and , the calculation of is similar to that of .

4.2.2. Deviance Information Criterion

The formula of for structure is described as [36, 39, 40, 53, 54]where represents the posterior mean deviance and shows the complexity measurement of the model.

Total deviance information criterion () of total TAZs is formulated as

4.2.3. Log Pseudo Marginal Likelihoods

of each structure satisfies the formula as follows [10, 43, 51]:where is the conditional predictive ordinate given by

In equation (21), is the conditional predictive distribution.

In terms of TAZs, the total is denoted as follows:

4.3. Variation Proportion

The paper utilizes , , and to reveal the variation proportion explained by the multivariate, structural, and spatiotemporal correlations (i.e., , H, and ST), respectively. The specific formula is shown aswhere is the standard deviation. The lower value of ( = 1, 2, 3) indicates the relatively weaker heterogeneous correlation corresponding to the total sum of heterogeneities [34].

5. Results

The paper conducts the modeling estimation with the software of WinBUGs 1.4.3 [55], where the total iterations perform as 30,000 containing discarded 10,000 iterations at the burn-in-step.

5.1. Goodness-of-Fit for Models

Tables 4 and 5 provide the predictive performance of models 1–4 and 5–7, respectively. According to Table 4, model 4 with the lowest value (888.565) prefers models 1–3 (separately with the value of 903.119, 926.193, and 973.746). However, the evaluation result of and is inconsistent with the counterpart of . It is because and fail to account for the complexity of models in contrast to (Table 3).

Table 5 adds the approach of to further quantify the division of temporal dimension in models 5–7. Model 5 notably rules out models 6-7 as associating with the highest value (−65.750 > −66.890 > −75.470) and the lowest value (4780.790 < 4785.050 < 4940.530). However, the estimate of and generates the biases compared to and according to Table 3.

With the above criteria, it is difficult to directly compare the fitness between models 1–4 and 5–7 due to the increase of temporal dimension in latter models. Therefore, the paper utilizes the method of variation proportion to complete the model comparison (Table 6). It is shown that the values of , , and are greater than 0, which verifies the existence of multivariate, structural, and spatiotemporal heterogeneities. By incorporating the multiple-heterogeneous correlations simultaneously, models 5–7 implicitly obtain a better predictive performance compared to models 1–4. Combined with the modeling results in Table 5, it can be summarized that model 5 poses a better fitness than other models. The parameter estimation and discussions below are based on model 5.

Meanwhile, the value of (0.764) is larger than (0.166) and (0.070) (Table 6). It implies that multivariate heterogeneity has the most powerful impact on the assessment results, followed by the spatiotemporal and structural heterogeneities [34, 40, 53].

5.2. Parameter Estimation

Presented in Table 7 is the regression result of models 5–7. The corresponding results of models 1-2 and 3-4 are shown in the Supplementary Materials section in Tables S1 and S2, respectively. Meanwhile, the point range plots of Table 7 are provided in the Supplementary Materials section in Figures S1S5. Based on the results of roadway segments in model 5, several contributors are positively associated with the occurrence of injury and PDO crashes, including the median opening width, lane density, posted speed limit, and log AADT. Meanwhile, the principal arterial is assumed to cause more crashes than the minor arterial according to the modeling result of functional roadway classification. In terms of the TAZs, the estimation results are summarized as (1) the increase of average AADT causing more crashes; (2) the higher average posted speed limit decreasing the possibility of crash occurrence; and (3) the length of TAZ not being a significant indicator of crash frequency. Regarding the signalized intersections, four variables (i.e., NETMA, NETMI, number of legs, and log AADT on major road) are the positively statistical evidence impacting on the number of injury and PDO crashes at the significance level of 95%. PLMA, PLMI, PMMA, and PMMI increase the crash count. The variation of ASD is unrelated to both crash occurrences.

6. Conclusion and Discussion

Previous research lacks the comprehensive heterogeneity consideration for the signal-coordinated arterials. This paper aims to build the safety evaluation model by incorporating multiple heterogeneities simultaneously, including the multivariate heterogeneity (at the microscopic level), structural heterogeneity (from the mesoscopic aspect), and spatiotemporal heterogeneity (in the macroscopic viewpoint). The model serves to identify the significant contributors to crash frequency, thus achieving the accurate safety analysis of arterials with the coordinated signalization.

Tables 5 and 6 show that models 5–7 perform better than models 1–4 due to the existence of multivariate, structural, and spatiotemporal heterogeneities. Signal coordination forms the “green wave” to impose the high-speed limit for the consecutively vehicular platoon, where the vehicles are densely correlated. Thus, the coupled vehicles connect as a whole generating the multivariate, structural, and spatiotemporal heterogeneities. The conclusion verifies the necessity of proposing the safety evaluation model with the multiple heterogeneities. In addition, it is imperative to provide the safety improvement countermeasures to alleviate the correlated crashes caused by the multiple heterogeneities, such as (1) adding the content of signal-coordinated arterials during the driving license test to enhance the drivers’ theoretical and practical cognition; (2) setting the warning signs at the entrance of coordinated arterials [9]; (3) equipping the variable message sign system (VMSS) to present the traffic information for vehicles in real time [56]; (4) providing vehicles with the advanced technologies of autonomous driving, connected vehicle, and artificial intelligence [5759], thus achieving the accurate crash control and eliminating the interference of humans.

Model 5 is superior to models 6-7 (Table 5), implying that there exists the fixed spatial and linear temporal variation (Section 4.1.5). Specifically, it indicates that the variation of spatial heterogeneity is unrelated to the temporal changing among the TAZs. It is possible that, with the aid of signal coordination, the stable spatial correlations on signalized intersections evenly affect the whole TAZs [8]. The signal coordination allows the vehicular platoon to pass through the intersections with the invariant headways, while the invariant headways are associated with the stable distribution of platoon, resulting in the unchanged spatial heterogeneity. Given the consistent spatial association, it is necessary to divide the subcontrol traffic zones based on the correlative degree of vehicles in real time. With the response to the dynamic traffic demand, the approach can avoid the “green wave” being too long or too short, thus ensuring the stability of vehicular movement and alleviating the unsafe behaviors (e.g., acceleration and “stop and go”) [60, 61].

Meanwhile, it is found that multivariate heterogeneity obtains a stronger impact on the evaluation results than the structural and spatiotemporal heterogeneities (Table 6). In view of the differentiation among the multivariate, structural, and spatiotemporal heterogeneities, appropriate suggestions are recommended for divers, such as (1) paying more attention to the vehicular distribution in real time, thus dynamically responding to the correlative criterion; (2) maintaining the driving speed and headway to ensure the awareness of safety risks; (3) adopting the actuated advance warning system (ADWS) [62]; and (4) equipping the on-board sensors in vehicles to retrieve a large amount of data (e.g., the acceleration, application status of vehicular components, and the distance between other vehicles) with the video recording technology [63].

For the roadway segments (Table 7), median opening width is positively correlated with injury and PDO crash counts. A feasible explanation is that the high-speed vehicular platoon migrates from the signalized intersections to the adjacent segments due to the implementation of signal coordination. Without the signalized control, the aggregated vehicles on the segments pose relatively short headways [8, 9]. Under the scenario, the wider median opening may incur more risky driving behaviors (e.g., the simultaneous U-turn of vehicles), thus resulting in the phenomenon of traffic jams more easily. Similarly, the increase of lane density, posted speed limit, log AADT, and functional roadway classification causes safety hazards along the roadway segments [39, 6466]. In order to alleviate the crash hazards of roadway segments along the coordinated arterials, the countermeasures are provided: (1) reducing the width of median opening, evenly decreasing the setting of median opening; (2) implementing the low-speed limit on the segments; (3) applying the vibrating deceleration pavement markings to control the vehicular speed; and (4) maintaining a sufficient sight distance for the drivers.

At the signalized intersections (Table 7), the higher possibility of injury and PDO crash occurrences is due to the increase of NETMA, NETMI, the number of legs, and log AADT on the major road. Meanwhile, PLMA, PLMI, PMMA, and PMMI result in more crash counts. The above findings are in accordance with the previous literature [8, 39, 6769]. The variable of ASD is independent of crash count. As opposed to the minor road, signal coordination allocates the adequate “green wave” on the major road with a high-speed limit [9, 41]. In the presentation of decreasing the ASD (or appropriately increasing the posted speed limit on the minor road), the crash variation of whole intersections is insignificant due to the limited traffic flow and “green bandwidth” on the minor road.

Although the study considers the potential disturbance involved in the data, the issue of multicollinearity may not be completely resolved. The safety evaluation of the signal-coordinated system needs to be further explored by employing more effective methods to address the multicollinearity among variables, e.g., testing each variable in the safety evaluation model and then excluding ones to be highly correlated.

Data Availability

The original accident data used to support the findings of this study are currently under embargo. Requests for data, 12 months after publication of this article, will be considered by the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The work was supported by the National Natural Science Foundation of China (NSFC-71771191), the Doctoral Research Foundation of Taiyuan University of Science and Technology (20202048), the Shanxi Scholarship Council of China (2021-136), and the Natural Science Foundation of Zhejiang Province (LQ21G010002). The authors would like to thank He Huang (School of Transportation and Logistics Engineering, Taiyuan University of Science and Technology) for the generous help.

Supplementary Materials

The following is the concise description for supplementary files. Tables S1 and S2 present the regression results of variables related to models 1-2 and 3-4, respectively. Meanwhile, Figures S1–S5 provide the point range plots of Table 7 from Section 5.2. (Supplementary Materials)