#### Abstract

Identification of key indicators that cause safety challenges and vulnerable roads is crucial for improving traffic safety. This paper, therefore, entails to the development of a novel multiphase multicriteria decision-making (MCDM) model to evaluate the vulnerability of urban roads for traffic safety. This was conducted by using data from 17 important roads of a South African city and combining several methods such as CRiteria Importance through Intercriteria Correlation (CRITIC), data envelopment analysis (DEA), and measurement of alternatives and ranking according to compromise solution (MARCOS). Taking the elements of the DEA method, two new approaches for calculating the weights of criteria, the DEA-1 and DEA-2 models, were formed and integrated with the CRITIC method to obtain the final values of criteria weights. The MARCOS method was applied to evaluate 17 alternatives, for each direction separately. The aim of developing such a model is to use the advantages of obtaining objectivity of criteria weights through linear programming and correlation of values of the collected data. Also, the MARCOS method, as one of the newer and applicable methods, provides additional significance. Extensive sensitivity analyses were conducted to validate the model. The findings suggest that there are a certain number of roads that have a high level of safety for both directions, as well as a group of risky roads, which need traffic improvement measures. Thus, the results indicate that the model is sensitive to various approaches and can prioritize vulnerable roads comprehensively based on which safety measures can be taken.

#### 1. Introduction

The occurrence of road accidents is a recurring challenge across the world. The problem is more pervasive in the cities of developing countries such as South Africa. The majority of the people in the South African cities live in the sprawling suburbs and usually travel daily to various parts of the cities and towards the city centres to work and for other civic activities. The cities are characterized by the location of important public and commercial activities on the roadsides or close to road intersections and squares and having predominantly mixed land-use systems in and around the central areas. Due to the lack of a significant and efficient public transportation system, the majority of people travel by their individual cars or shared public taxis. Moreover, access roads linking suburban residential areas and arterial roads to provide higher accessibility also create more conflict points on these roads. The combined effect of various land-use, activity, road, and traffic-related challenges makes the urban roads susceptible to the incidence of traffic accidents. The occurrence of higher vehicular traffic accidents has been experienced not only on the principal thoroughfares of the cities but also on the arterial roads and intersections located in the suburban areas. For instance, according to the reports from the police stations of a city, more than 45% of the accidents that happen in the city are observed on the urban arterial roads that connect different suburbs of the city. However, it is also observed that some of the roads are more prone to accidents than others and need adequate interventions.

Improvement of road safety and the reduction of vehicular accidents have been considered as a prime focus in the cities. The road and traffic departments and municipalities are undertaking remedial measures such as improving road conditions, providing road pavement marking and signage, improving road traffic management systems such as signals, implementing traffic regulations, creating traffic safety awareness programs, and traffic safety audits. However, with limited resources, it becomes difficult to take up remedial measures on each road across the city, and therefore, prioritization of vulnerable roads and areas of accidents is necessary.

In the context of South African cities, most of the studies are aimed at either establishing the linkage between different road, traffic, environmental or human factors with the occurrence of the accidents on the roads or developing remedial measures to improve road safety. However, identification of vulnerable roads and prioritization of these roads based on road geometry, traffic-related parameters, and their linkage with the occurrence of accidents have been undermined. Moreover, most of the studies that have been conducted were based on the use of conventional statistical methods, and studies by using advanced methods are scarce. Since remedial measures to improve road safety demand prioritization or ranking of vulnerable roads in addition to exploring the causes of the occurrence of accidents, there is a need for robust, advanced, and flexible methods that could accommodate various parameters effectively and offer reliable predictions.

Therefore, the objective of the study is to develop an appropriate and robust method that can evaluate and rank vulnerable roads for the occurrence of vehicular traffic accidents effectively. For this purpose, a novel multiphase model consisting of several methods, CRiteria Importance through Intercriteria Correlation (CRITIC), data envelopment analysis (DEA), and measurement of alternatives and ranking according to compromise solution (MARCOS), was developed by using the data from 17 important arterial roads (U3-minor arterial roads) (COTO, TRH26, 2012) of Bloemfontein city in South Africa. The paper presents a study in the city of Bloemfontein related to traffic safety in 17 different streets. Data for both directions of movement (N and S) were collected, and a new integrated multiphase model was developed in order to obtain the most relevant indicators of traffic safety in the city. Six criteria were defined. Five of them represented input parameters in the DEA model, while the sixth criterion represented an output parameter, and the values of criteria weights were determined by developing two new approaches DEA-1 and DEA-2. After that, the values obtained with the CRITIC method were integrated, and the final values of criteria weights were calculated. They show that the number of traffic accidents per km is the most important factor, which is understandable. After that, implementation of the stated weights of criteria in the MARCOS method for ranking streets has been performed. Fulfilling the seated objective and the creation of an integrated model where the DEA method was used in one innovative way can be manifested as the main contribution of the study.

The rest of the paper is structured in the following manner. A succinct account of the literature concerning the traffic crash prediction models (CPMs), specifically the use of CRITIC, DEA, and MARCOS, and demands for hybrid models are given in the next section. This is followed by the methodology, which offers an account of the research flow for the development of the model and description of the CRITIC, DEA, and MARCOS methods. Case study, data collection, and determination of weights by using CRITIC, DEA-1, and DEA-2 models and final weights by the combination of these models as well as evaluation of alternatives by the use of the MARCOS model are presented in the following section. The next section includes the multiphase sensitivity analysis which was conducted by using a reverse rank matrix, changing criteria weights, and comparing with other MCDM models to validate the developed hybrid model. Finally, the implication of the model in terms of ranking of the streets concerning traffic safety is presented in the conclusion section.

#### 2. Literature Review

Many traffic accident prediction models are available and used to estimate the occurrence of the accidents on a road section because of road and traffic-related parameters. These models enable establishing quantitative relations between road and traffic characteristics and traffic accidents, as well as the influence of the application of safety interventions on a road section [1]. These models include both traditional statistical or mathematical methods and advanced models based on artificial intelligence models such as multicriteria decision-making models. The conventional assessment methods include relative accident rate analysis, time series analysis, and regression analysis [2, 3]. These methods face a common challenge—the assessment results disagree due to the often different parameter choice [4]. Similarly, weight evaluation methods such as the analytic hierarchy process (AHP), expert scoring method (Delphi), principal component analysis method, eigenvalue method, and gray correlation method are used. The general problem of these methods is that they do not consider the internal correlation and inconsistency of indicators, and the outcomes obtained are often fairly different from the real ones [5]. Consequently, advanced MCDM models such as DEA, MARCOS, TOPSIS, and CRITIC methods that consider the uncertainties and relationship between qualitative concepts and quantitative concepts, which could assist in justifiable, explainable, and transparent decision-making, are increasingly used [3].

Furthermore, integrated model based on the q-ROULSs and EDAS method [6] and OHS risk assessment model by integrating picture fuzzy sets (PFSs) and alternative queuing method (AQM) [7] were used to evaluate and rank the risk of occupational hazards. In another recent study, association rule mining technique, as well as correlating various attributes to the severity of the accident, has been used to identify accident spots [8]. Apart from these models, although not directly linked to traffic safety or OHS, a new method based on double hierarchy hesitant linguistic term sets (DHHLTSs) and alternative queuing method (AQM) to evaluate the satisfaction of the rail transit network under a large group environment was developed recently [9].

Data envelopment analysis (DEA) was used by various scholars for analyzing traffic safety aspects. It has been used to prioritize traffic crash sites [10], examine operational efficiency related to the traffic safety [11], assess the urban safety on urban roads, Runde et al., prioritize road safety needs [12], develop road safety policies, and so on. DEA is able to consider multiple inputs and outputs and does not require a functional form, which relates to inputs and outputs. It can optimize each observation and can compare it against best practice observations. This method is also different from other ranking methods as it can add standard errors of crash modiﬁcation factor and crash costs in the selecting process, as well as the average values (Sadeghi & Moghaddam, 2016). However, DEA has certain limitations that include the following: it can only calculate relative efficiency measures; and since it is a nonparametric technique, statistical hypothesis tests are difficult [13]. Also, when the number of inputs and outputs is more than the number of decision-making units, DEA models will not be able to separate DMUs [14]. Consequently, DEA, in recent years, is integrated with other methods such as artificial neural networks [15], analytical hierarchy process [16], analytic network process [17], TOPSIS [13, 18], and fuzzy logic [15, 16] to convert the qualitative variables into quantitative equivalences.

The MARCOS method has also been used to evaluate the road infrastructure on a road section to rank the accident risk on road sections. The ability to consider fuzzy reference points through the fuzzy ideal and fuzzy anti-ideal solution at the very beginning of model formation is a major advantage of this method. Also, the ability to determine the degree of utility concerning both set solutions more precisely, proposing new ways to determine utility functions and their aggregation, and the possibility to handle a large set of criteria and alternatives are the other advantages [19].

The CRITIC model determines the objective weights that are premised on the quantification of two fundamental notions of multicriteria decision-making (MCDM) such as the contrast intensity and the conflicting character of the evaluation criteria [20]. The extraction and exploitation of these two fundamental features which are stored as intrinsic information in the data defining the multicriteria problem benefit the decision-making process. In this method, the objective weights derived incorporate both the contrast intensity of each criterion and the conflict between the criteria. Also, the contrast intensity of the criteria is considered by the standard deviation. The conflict between them is assessed by the correlation coefficient [21]. This method comprehensively measures the objective weight by the contrast intensity and conflict between different indexes, which is argued to be more scientific and reasonable. However, studies on the application of the CRITIC method in traffic safety are scarce.

Each of the models discussed has limitations and advantages, and the performance of each model depends on the information available and parameter selection [3, 4]. As a result, often, the models are used in combination with other relevant models, or in other words, hybrid models are popular. Therefore, in this study, the integration of DEA, MARCOS, and CRITIC models has been made, which is a major contribution to the body of the knowledge in the analysis of road traffic safety.

#### 3. Methodology

##### 3.1. A Developed Model for Traffic Safety Evaluation

Figure 1 presents a whole diagram of the research flow which is divided into four main phases that are interdependent. The first and second phase have four steps each. The third phase has seven activities, while the last, fourth phase includes a sensitivity analysis that needs to be performed through seven steps, too.

The model presented in Figure 1 encompasses a complete research flow consisting of a total of 22 activities that are causally related. The first phase is the recognition of the need for research, and the area of traffic safety is always present and important for the study of key indicators that affect output parameters, which relate to the number of traffic accidents. After the need for research in the city of Bloemfontein has been determined, the next step is field exploration and then the formation of a database, sorting the data according to the needs of research. The second phase involves determining adequate criteria based on which the evaluation of alternatives will be performed. Six criteria are defined. Five of them refer to input parameters, and the sixth is the output parameter in DEA. After that, data are sorted for each of 17 streets in which the research was conducted. It is important to note that a double MCDM model has been formed, which includes separate data collection and evaluation of streets for both directions (north and south). The third phase is the creation of a new integrated multicriteria model consisting of several steps. CRITIC [22, 23] and DEA methods were used to determine the weight values of criteria. The DEA model was used to determine the values of the criteria in two ways, which is presented in more detail below. Based on their individual values, averaging was performed, the final values were obtained, and they were included in the MARCOS method [24–26] when evaluating alternatives. The last, fourth phase is the verification of previously obtained results through a multiphase sensitivity analysis. Validity tests refer to the (a) formation of reverse rank matrices in which the worst alternative is eliminated from further calculation, (b) change of criterion weight in a dynamic environment, (c) evaluation of alternatives based on individually obtained values of criteria weights using CRITIC, DEA-1, and DEA- 2 methods; in addition, a scenario in which all criteria have equal significance was created, (d) comparative analysis of the obtained results with five other MCDM methods, (e) calculation of SCC for all ranks obtained by applying different approaches, and (f) calculation of the standard deviation of ranks obtained by applying different approaches.

##### 3.2. CRITIC Method

Diakoulaki et al. [20] introduced the CRITIC method as a tool for determining the objective weights of criteria in MCDM problems. Steps of this method are presented as follows: Step 1: formation of the decision matrix (*X*): Step 2: normalization of the initial decision matrix depending on the criterion type: Step 3: calculation of symmetric linear correlation matrix *r*_{ij}: Step 4: determination of objective weights:

##### 3.3. DEA Method

The DEA-based road safety (DEARS) model proposed by Shen et al. [27] is as follows:

The decision-making unit (DMU) consists of *m* input parameters, while *s* represents output parameters taking into account the weights of the parameters denoted by .

In this paper, the DEA method is applied to calculate the weights of the criteria in one part of the model. In order to be able to accurately determine the weights after obtaining them by applying equation (6), the following two equations (7) and (8) are formed for the DEA-1 model and (9) and (10) for the DEA-2 model:where *n* represents the total number of DMUs. In essence, in this equation, all the weights of the criteria are averaged using the arithmetic mean.

After that, the final values according to the DEA-1 model are obtained as follows:where *V*_{j} represents the mean values of weights from the previous equation.where all individual weights are divided by their sum.

Final values according to the DEA-2 model are obtained by applying equation (10). In essence, in this equation, all the weights of the criteria are averaged using the arithmetic mean. Also, *n* represents the total number of DMUs.

##### 3.4. MARCOS Method

The MARCOS method developed by Stević et al. [24] consists of the following steps [28]: Step 1: formation of an initial decision-making matrix. Step 2: formation of an extended initial matrix. In this step, the extension of the initial matrix is performed by defining the ideal (AI) and anti-ideal (AAI) solution. The anti-ideal solution (AAI) is the worst alternative, while the ideal solution (AI) is an alternative with the best characteristic defined by applying equations (12) and (13): where *B* represents a benefit group of criteria, while *C* represents a group of cost criteria. Step 3: normalization of the extended initial matrix (*X*). The elements of the normalized matrix are obtained by applying equations (14) and (15): where elements and represent the elements of matrix *X*. Step 4: determination of the weighted matrix : Step 5: calculation of the utility degree of alternatives *K*_{i} by applying equations (17) and (18): where (*i* = 1, 2, ..., *m*) represents the sum of the elements of weighted matrix *V*: Step 6: determination of the utility function of alternatives *f*(*K*_{i}) defined by the following equation: where represents the utility function in relation to the anti-ideal solution, while represents the utility function in relation to the ideal solution. Utility functions in relation to the ideal and anti-ideal solution are determined by applying equations (21) and (22): Step 7: ranking the alternatives. Ranking of the alternatives is premised on the final values of utility functions. It is desirable that an alternative has the highest possible value of the utility function.

#### 4. Case Study

The arterial roads of Bloemfontein city in South Africa (Figure 2) were chosen for collecting data for this research. Bloemfontein is a typical middle-sized sprawling city with a population of more than 500,000 and functions as the capital of the Free State province. The majority of the people live in suburban areas, thereby necessitating large-scale intracity vehicular travel.

The city has a hierarchical road system with principal arterials (U1 and U2), minor arterial (U3), collector roads (U4), local streets (U5), and cul-de-sacs (COTO, TRH26, 2012). The road network constitutes a general gridiron and loop pattern and the integration of the two. The U3 roads, which pass through the suburban areas, connect to the U1 and U2 roads of the city. Also, these U3 roads act as the major thoroughfares in the suburban areas offering linkage to the collector roads and the local streets to offer access to the people.

The road surfaces of most of the roads are paved and well maintained. The principal mode of travel in the city is personal/individual-driven motor cars. About 60% of people travel by their private vehicles for daily commute. However, according to police reports, a significant share (about 45%) of the total accidents, which occur in and around the city, is observed on these U3 arterial roads. Thus, the importance of the arterial roads (U3) and the high incidence of traffic accidents on these roads of the city warranted this investigation.

##### 4.1. Survey and Data Collection

A survey research method was used to collect data. Data were collected from both primary and secondary sources. Three different types of survey, household surveys, traffic surveys, and road geometrical parameter survey, were made. The household sample survey was conducted among 410 road users (from households located in six important suburbs of the city. The suburbs include Fichardt Park, Pellissier, Universitas, Gradenia Park, Wilgehof, and Langenhoven Park. The suburbs were chosen based on their functional importance, population, vehicle ownership, accessibility, the complexity of the road network, the location of arterial thoroughfares inside the residential area, and the occurrence of the number of accidents. The sample size is adequate (≥384) for a population 500,000 at a 95% confidence level and worst-case percentage of 50%). This household survey was conducted to understand the perceptions of road users on the factors that cause accidents on urban roads. A systematic stratified random sampling process was used to conduct the survey. Traffic and road surveys were conducted on important road sections. A traffic survey was conducted at important road sections of the selected 17 U3 roads passing through the selected suburban areas to observe the volume (ADT), speed, and traffic control and management scenarios. The roads selected for this investigation include Jan Spies Street, Totius Street, Wynand Mouton Drive, De Bruyn Street, Paul Kruger Drive North, Paul Kruger Drive South, Volkspele Drive, Pellissier Drive, Paul Kruger Drive (E), Gardenia Avenue, Van Schalkwyk Street, Haldon Drive, Edeling Street, Stals Road, Benade Drive (N), Benade Drive (S), and Eric Rosendorf Street. These U3 roads connect different suburbs and important activity centres of the city. Generally, these roads are flexible paved one-way roads having two or more lanes with or without medians. The junctions of the roads are controlled by automated signalling systems. The traffic survey was conducted uninterruptedly for sixteen hours a day (6.00 am to 22.00 hours) for a week that includes both weekdays and weekends. Similarly, physical road surveys were conducted at important road sections of these roads to assess the current status of the road geometry parameters (road width, shoulder width, lane width, number of lanes, curbs, curvature, median width, gradient, sight distance, and road surface condition) and their influence on the incidence of the accidents. The households’ surveys and road and traffic surveys were conducted during the year 2017. The average accident data on different roads over the period 2010 to 2017 were collected from published and unpublished documents and the police station records of the city. The surveys and assessment of the incidence of accidents were conducted on both directions represented by N and S, indicating the two different ways of these roads.

##### 4.2. Determining Criteria Weights Using the CRITIC Method

Step 1: development of the decision matrix (*X*) that is shown in Table 1. The initial matrix is shown for both directions since, as already mentioned, the data are different by directions, which means that we practically have two initial matrices. Step 2: normalization of the initial matrix is shown in Table 2. Normalization for cost criteria is performed using equation (2), for example, Normalization for benefit criteria is performed using equation (3), for example, Step 3: the symmetric linear correlation matrix (

*m*

_{ij}) is shown in Table 3. Step 4: determination of the objective weights of criteria, shown in Table 4.

After applying all the steps of the CRITIC method and the weighted normalized matrix shown in Table 5, it can be seen that the third criterion, road width (*m*), is the most significant with a value of 0.220 for the N-direction and a value of 0.204 for the S-direction. It is interesting to note that the sixth criterion, the average number of accidents (per km), is the second most significant and has an identical value of 0.187 for both directions. Essentially, the differences in the values of the criteria are not dominant in any case.

##### 4.3. Determining Criteria Weights Using DEA-1 and DEA-2 Models

By applying equation (6), the individual values of the criteria for each DMU are obtained. An example of obtaining the values of criteria for DMU17-N is shown in Appendix. After solving this model shown in Appendix in Lingo 17 software, the following values are obtained: = 0.005, = 0.000, = 0.095, = 0.000, = 0.000, and = 0.111. It is necessary to form such a model for all 17 DMUs in order to obtain individual values. By solving all 17 models, the weights shown in Table 6 are obtained.

By applying equations (7) and (8), the values of the criteria according to the DEA-1 model are obtained as follows:

By applying equations (9) and (10), the values of the criteria according to the DEA-2 model are obtained as follows:

The previous calculation refers to the N-direction, while the values for the S-direction are obtained in an identical way, which together with the previous calculations and the final values are shown in Figures 3 and 4.

**(a)**

**(b)**

**(a)**

**(b)**

##### 4.4. Determining the Final Criteria Weights

The final values of the criteria are obtained by averaging the weight values of the criteria obtained by individual methods: CRITIC, DEA-1, and DEA-2. Figure 3 presents the individual and final values of the criteria for the N-direction (Figure 3(a)) and the S-direction (Figure 3(b)).

From Figure 3, it can be seen that there is a difference in applying different approaches and that different values are obtained. This is the reason for their integration when obtaining the final values. Using the CRITIC method, the following values are obtained: = 0.145, = 0.170, = 0.220, = 0.122, = 0.157, and = 0.187 for the N-direction and = 0.159, = 0.182, = 0.204, = 0.120, = 0.147, and = 0.187 for the S-direction. In the previous section of the paper, the values obtained using DEA models are explained. The DEA-1 model for the N-direction implies the values = 0.013, = 0.0003, = 0.139, = 0.154, = 0.221, and = 0.472 and for the S-direction, = 0.0003, = 0.0001, = 0.031, = 0.300, = 0.323, and = 0.416. Applying a different approach, i.e., the DEA-2 model, the values of the criteria for the N-direction are = 0.014, = 0.000, = 0.144, = 0.054, = 0.274, and = 0.515, while for the S-direction, they are = 0.0003, = 0.0001, = 0.044, = 0.086, = 0.414, and = 0.456. Observing the range of weight values of the criteria, including the final values, the largest deviation is for the sixth criterion with a value of 0.328 for the N-direction. There is a slightly smaller deviation for the sixth and fifth criterion for the S-direction. For other criteria, the range of values is smaller.

The final values of the criteria further included in the MARCOS method are = 0.057, = 0.057, = 0.168, = 0.110, = 0.217, and = 0.391 for the N-direction and = 0.053, = 0.061, = 0.093, = 0.145, = 0.295, and = 0.353 for the S-direction. Figure 4 also shows the ranks, i.e., the significance of all criteria according to each individual model and for the final values. The most significant criterion according to DEA-1 and DEA-2 models and the final values is the average number of accidents (per km), which is understandable, while the same criterion is the second according to the CRIITIC method. The second most significant criterion is the number of access (per km) for all models, except for the already mentioned CRITIC, where it is in the fourth position. The third most significant criterion is road width (m) according to the final values and the DEA-2 model, while by applying DEA-1, it is in the fourth place, i.e., in the first place by applying the CRITIC method. The fourth criterion, the number of lanes, is positioned in the fourth place, while the last two places occupy the first and second criterion, respectively. The listed ranks and significance refer to the left side of the figure which implies the N-direction. When it comes to the significance of the criteria for the S-direction, it can be concluded that the ranks are more correlated than for the N-direction. The difference in the final values in terms of ranks is only for the third and fourth criterion, which change positions with each other.

##### 4.5. Evaluation of Alternatives Using the MARCOS Method

The following section of the paper presents the evaluation of alternatives using the MARCOS method. For each step of the method, an example of the calculation is given in order to facilitate the monitoring of the obtained results. Table 7 shows an extended initial decision matrix obtained when Table 1 is expanded with an ideal (AI) and an anti-ideal (AAI) solution depending on the type of criteria. It is important to emphasize that the third and fourth criterion are of the benefit type (orientation max), while the others are of the cost type (orientation min).

Applying equation (14), the normalized values for cost criteria are obtained, e.g., , and for benefit criteria using equation (15), , and a complete normalized matrix, shown in Table 8, is obtained.

The next step is weighting the normalized matrix using equation (16) by multiplying all the values of the normalized matrix with the values of the criteria. The weighted normalized matrix is shown in Table 5.

Applying equations (17)–(22), the final results of Table 9 are obtained using the MARCOS method. The results are obtained as follows.

By applying equation (19), all the values (by rows) for alternatives are summarized as follows:

Similarly, the values for the remaining alternatives are obtained for both directions.

By applying equation (17), the utility degrees in relation to the anti-ideal solution are calculated. The example of the calculation iswhile using equation (18), the utility degrees in relation to the ideal solution are obtained, e.g.,

The utility function in terms of the anti-ideal solution is obtained using equation (21) as follows:while the utility function in terms of the ideal solution is obtained using equation (22) as follows:

Finally, the utility function of alternative A1 is obtained by applying equation (20):

After applying the whole methodology, the results for both directions are summarized in Table 9. When it comes to the rank, i.e., determination of the traffic safety level in 17 streets of the city of Bloemfontein, it can be noticed that alternative A11, i.e., Van Schalkwyk Street, is the safest with a value of 0.794 for the N-direction. Volkspele Drive is the second in terms of safety, with Eric Rosendorf Street and Gardenia Avenue in the third and fourth place, respectively. Another alternative A13, i.e., Edeling Street, can be put into a group of safe streets comparing all alternatives. The group of most risky streets includes Wynand Mouton Drive, Paul Kruger Drive North, Benade Drive (N), Totius Street, and Paul Kruger Drive South with a range of values from 0.341 to 0.421.

As already mentioned in the previous section of the paper, data have been collected for both directions, and the evaluation is performed accordingly, so it is evident that there are differences in terms of safety of individual streets. When it comes to the S-direction, the group of the five safest streets includes the following, respectively: Eric Rosendorf Street, Volkspele Drive, Stals Road, Gardenia Avenue, and Pellissier Drive. The group of least safe streets for the S-direction includes Wynand Mouton Drive, Paul Kruger Drive North, Paul Kruger Drive (E), and Benade Drive (N).

Three streets that belong to the group of safe regardless of the direction are Eric Rosendorf Street, third and first positions, Volkspele Drive, second place for both directions, and Gardenia Avenue which is in the fourth place. In terms of the most risky streets for both directions, Wynand Mouton Drive is in the last place, Paul Kruger Drive North in the 16th place, Benade Drive (N) in the 15th and 14th place, respectively, and Totius Street in the 14th and 13th place, respectively.

#### 5. Sensitivity Analysis

Verification of previously obtained results was performed through a multiphase sensitivity analysis. Validity tests refer to the (a) formation of reverse rank matrices in which the worst alternative is eliminated from further calculation, (b) change of criterion weight in a dynamic environment, (c) evaluation of alternatives based on individually obtained values of criteria weights using CRITIC, DEA-1, and DEA-2 methods; in addition, a scenario is created in which all criteria have equal significance, (d) comparative analysis of the obtained results with five other MCDM methods, (e) calculation of SCC for all ranks obtained by applying different approaches, and (f) calculation of standard deviation of ranks obtained by applying different approaches.

##### 5.1. Reverse Rank Matrix

In this part of the sensitivity analysis, the worst alternative is eliminated through 16 scenarios. In each of the scenarios, the last-ranked alternative is eliminated, and a full calculation is performed again, until one alternative remains, as is the case in scenario 16.

Figure 5 presents the results of the ranks of alternatives through previously formed scenarios for both directions. As can be seen, there is no change in ranks in any case, which means that the model is stable, and the results are valid in such conditions.

##### 5.2. Changing Criteria Weights

In this phase of the validation test, the impact of changing the three most important criteria *C*6, *C*5, and *C*3 for the N-direction and *C*6, *C*5, and *C*4 for the S-direction on ranking results was analyzed. Applying equation (33), a total of 18 scenarios were formed.

The sixth criterion was changed in scenarios *S*1–*S*6, criterion *C*5 was changed in scenarios *S*7–*S*12, and criterion *C*3 was changed in scenarios *S*13–*S*18 for the N-direction. In equation (33), represents the new value of criteria *C*1–*C*5 for scenarios *S*1–*S*6 and then *C*1–*C*4 and *C*5 for scenarios *S*7–*S*12, i.e., *C*1-*C*2, *C*4, and *C*5-*C*6 for scenarios *S*13–*S*18. represents the corrected value of criteria *C*6, *C*5, and *C*3, respectively, by groups of scenarios, represents the original value of the criterion considered, and represents the original value of the criterion whose value is reduced, in this case, *C*6, *C*5, and *C*3.

For the S-direction, the difference is that, in scenarios *S*13–*S*18, the value of criterion *C*4 was changed instead of *C*3.

In all scenarios, the value of criteria was reduced by 15%, while the values of the remaining criteria were proportionally corrected by applying equation (33). After forming 18 new vectors of the weight coefficients of the criteria (Table 10), new results were obtained, as presented in Figure 6.

**(a)**

**(b)**

Figure 6(a) presents the ranks of alternatives for the N-direction due to the formation of 18 scenarios with new criteria weights. As can be seen, the rank changes due to the simulation of different values of the criteria, which means that the results are sensitive to the significance of the criteria. The biggest changes compared to the initial rank (SET 0) occur in the sixth set when alternative ST11 falls from the first to the third place. The consequence of such a change is in the value of the sixth criterion, which is only 0.039. If we take into account that, in real conditions, the sixth criterion should not have such a low value, it remains the best alternative option. In the sixth scenario, dominant criteria are the fifth and the third, which consequently influence the sixth-ranked alternative to become the first compared to the initial ranking. In other scenarios, there is also a change in ranks, but they are negligible compared to the above scenario.

When it comes to changes of ranks for the S-direction, Figure 6(b), smaller changes can be noticed, i.e., greater stability of the initial rank. The best three alternatives ST17, ST7, and ST14 do not change the ranks of the alternatives at all regardless of the criteria values. In addition, there are fewer changes through the other scenarios, with the exception of SET11 and SET12 in which changes occur by four positions, which refer to ST4 which moves from the tenth position to the sixth. Taking into account the large changes in the significance of the criteria, it can be concluded that the changes in the ranks are not huge, which will be proven later through the calculation of the correlation of ranks.

In addition to creating dynamic matrices in which the values of the criteria weights change, this part of the sensitivity analysis shows a comparison of the ranks of the proposed model with CRITIC-MARCOS (CM), DEA-1-MARCOS, and DEA-2-MARCOS models and an approach with equal criteria weights, which is shown in Figure 7.

**(a)**

**(b)**

Figure 7(a) presents the ranks for the N-direction under different approaches to determine the weight values of the criteria. The first two best alternatives retain their positions regardless of the approach applied, while ST17 falls from the third position to the fourth using the DEA-1-MARCOS and DEA-2-MARCOS models, respectively. The biggest changes are related to the change of rank by two places compared to the initial results. When it comes to the S-direction, Figure 7(b), the results are similar, with three best-positioned alternatives retaining their places for all applied approaches. In addition, it can be concluded that the model is sensitive to the application of different approaches, but the results are highly correlated.

##### 5.3. Comparison with Other MCDM Methods

In this section of the paper, a comparative analysis is performed with five other MCDM methods: ARAS—additive ratio assessment [29], MABAC—multiattributive border approximation area comparison [30, 31], SAW—simple additive weighting method [32], WASPAS—weighted aggregated sum product assessment [33, 34], and EDAS—evaluation based on distance from average solution [35]. Figure 8 shows the results of a comparative analysis for the N-direction, while Figure 9 shows a comparative analysis for the S-direction.

By applying the MABAC method, alternative ST11 takes the second position, while by applying all other methods, it represents the best solution. Essentially, the two best alternatives, ST11 and ST7, change their ranks with each other when calculating with the MABAC method. Alternative ST17 is in the third position using all MCDM methods. When it comes to the fourth and fifth position, the situation is the same as for the previously mentioned changes where ST10 and ST13 change their ranks with each other in the MABAC method. The biggest changes are in relation to the MABAC method, where one alternative can change its position by two ranks, while the application of other methods leads to a shift of rank by only one position in some cases.

A comparative analysis for the S-direction shows that the completely observed model has minor deviations in the ranks since, e.g., the first six positioned alternatives, ST17, ST7, ST14, ST10, ST8, and ST11, do not change their positions regardless of the applied method. When we observe individual rank deviations, there is a much different situation since there are deviations by three positions using the ARA method when the ST12 alternative moves from the seventh to the tenth position. An even more drastic change relates to the application of the MABAC method and the ST12 alternative when it changes its place by five positions and comes in the 12th place. Regarding other alternatives and methods, changes in rankings compared to the initial results are practically negligible.

##### 5.4. Calculation of SCC and ST_{DEV} for All Parts of the Sensitivity Analysis

This part of the sensitivity analysis refers to the calculation of the Spearman correlation coefficient for all previously applied approaches as well as the calculation of the standard deviation shown in Figure 10.

Observing the results shown in green, it can be noticed that the largest deviation of 2.251 is for alternative ST12 which changes its position by three to five places, and it refers to a comparative analysis for the S-direction. Then, the fourth alternative ST4 shows the largest subsequent deviation in ranks of 0.753 which results in a change by one to two positions. Alternative ST13 has a deviation of 0.548 because it changes its position by one place in the calculation of three methods (MABAC, WASPAS, and EDAS). Alternatives ST1, ST6, and ST16 show a deviation of 0.516 because they change by one position in two cases. Alternative ST5 has a deviation of 0.408 because only by applying the EDAS method, it changes its rank by one position. All remaining alternatives have a standard deviation of zero, which means that they retain their ranks by applying all methods. When it comes to deviations for the N-direction, they are generally individually smaller because the largest standard deviation is 0.753 for ST1, ST4, and ST12, which means that there is a change of rank by two positions. Alternatives ST3, ST5, and ST16 have a ST_{DEV} of 0.516, which means a change of rank by one position. Other alternatives have deviations of 0.408 or zero. Observing the results obtained by applying different approaches to determine the weights of the criteria, larger deviations are noticed for both directions, in a range of 0–1.643 for S1 and 0–1.517 for N1.

Table 11 and Figures 11 and 12 show the Spearman and WS [36] correlation coefficient of all ranks obtained by different models and approaches through a whole sensitivity analysis.

When it comes to the ranks for a comparative analysis of the obtained results, the calculated SCC shows the following correlations: MARCOS and SAW methods have a complete correlation, SCC = 1.000 for both directions. For the N-direction, the MARCOS and WASPAS methods have an almost full correlation of 0.998 resulting from the change of alternatives ST1 and ST12 by one position. ARAS with WASPAS and EDAS for the N-direction and ARAS with WASPAS and MABAC with EDAS for the S-direction have the same correlation value. SAW and ARAS for the N-direction have SCC = 0.995. MARCOS with EDAS for the N-direction and MABAC with WASPAS for the S-direction have the correlation value of 0.993. ARAS and MABAC and WASPAS and EDAS for the S-direction have a correlation of 0.990. MABAC and WASPAS for the N-direction and ARAS and EDAS for the S-direction have SCC of 0.988. MARCOS has the lowest correlation with MABAC for both directions, 0.980 and 0.963, respectively. Since the SAW method has a full correlation with MARCOS, it means that it has the same correlation with other models as the MARCOS method does. Taking into account the values of SCC according to Subotić et al. [37] and Muravev and Mijić [38], it can be concluded that the ranks of all alternatives according to a comparative analysis are very high, i.e., almost completely correlated. After additional calculation of WS coefficient (marked with grey color), it can be concluded that ranks have extremely high correlation.

Comparing the correlation values of the MARCOS method with approaches in which the weight values of the criteria have been calculated, it can be noticed that there is an extremely high correlation which ranges from 0.966 to 0.988. MARCOS has the highest correlation with DEA-1-MARCOS with a value of 0.988, while with DEA-2-MARCOS, it has a slightly lower correlation with a value of 0.980. When calculating with equal criteria weights and inserting them in the MARCOS method, a correlation of 0.973 compared to the initial results is obtained. The results have the lowest correlation with comparisons by the CM model (0.966). These correlations refer to the N-direction, while for the S-direction, the correlation is quite lower and ranges from 0.848 to 0.978.

Figure 11 shows the SCC and WS values for the ranks through 18 different scenarios in which the changes of criteria weights for the N-direction were made. The average value of SCC is 0.983, which means an extremely high correlation. The average value of WS is 0.998, which means that ranks have almost full correlation. Although the model is sensitive to changes in the weights of the criteria and a change in ranks has been established based on the calculated SCC and WS, it can be concluded that changes in ranks are not large.

Figure 12 shows the SCC and WS values for the ranks through 18 different scenarios in which the changes of the criteria weights for the S-direction were made. The results show similarities with the N-direction since the average correlation value, in this case, is 0.989 for SCC and 0.998 for WS coefficient, which is a very high correlation.

#### 6. Conclusion

Implementing the developed model, the following results are obtained: alternative A11, i.e., Van Schalkwyk Street, is the safest with a value of 0.794 for the N-direction. Volkspele Drive is the second in terms of safety, with Eric Rosendorf Street and Gardenia Avenue in the third and fourth place, respectively. Another alternative A13, i.e., Edeling Street, can be put into a group of safe streets comparing all alternatives. The group of most risky streets includes Wynand Mouton Drive, Paul Kruger Drive North, Benade Drive (N), Totius Street, and Paul Kruger Drive South whose values range from 0.341 to 0.421. As already mentioned in the paper, data were collected for both directions, and the evaluation was performed accordingly, so it is evident that there are differences in terms of safety for individual streets. When it comes to the S-direction, the group of the five safest streets includes the following, respectively: Eric Rosendorf Street, Volkspele Drive, Stals Road, Gardenia Avenue, and Pellissier Drive. The group of least safe streets for the S-direction includes Wynand Mouton Drive, Paul Kruger Drive North, Paul Kruger Drive (E), and Benade Drive (N). Three streets that belong to the group of safe ones regardless of the direction are Eric Rosendorf Street, third and first positions, Volkspele Drive,second position for both directions, and Gardenia Avenue, which is in the fourth place. In terms of the most risky streets for both directions, Wynand Mouton Drive is in the last place, Paul Kruger Drive North in the 16th place, Benade Drive (N) in the 15th and 14th place, respectively, and Totius Street in the 14th and 13th place, respectively.

In order to validate the proposed model and the results obtained by its application, an extensive sensitivity analysis was performed consisting of (a) a reverse rank matrix, (b) a change of criterion weight in a dynamic environment, (c) an evaluation of alternatives based on individually obtained values of criteria weights, (d) a comparative analysis of obtained results with five other MCDM methods, (e) calculation of Spearman correlation coefficient (SCC) for all ranks obtained by applying different approaches, and (f) calculation of standard deviation of ranks obtained by applying different approaches. The validation results show that taking into account the large changes in the significance of the criteria, it can be concluded that changes in the ranks are not large, which was proven through the SCC calculation. Using DEA-1-MARCOS, i.e., DEA-2-MARCOS, and CM models, it can be concluded that the model is sensitive to the application of different approaches, but the results are highly correlated. A comparative analysis also confirmed the validity of the results, as well as the calculated SCC for all approaches and the calculated standard deviation.

Based on the findings from the model in terms of the ranking of vulnerable streets, the streets can be prioritized to take up road geometry and traffic-related intervention measures which could include applying the control measures to restrict the speed, optimal assignment of traffic volume during the peak hours, limiting the number of access streets, and the number of lanes, as well as providing adequate median width.

The scope of the study was confined to the development of the model based on road and traffic-related parameters, and the human and driver behaviour or weather-related parameters were not considered. However, the integration of such parameters to make the model more inclusive and holistic is the future scope of this research. However, at the current state, the model can enable assessment of the vulnerable streets adequately and prioritize them based on which traffic safety measures can be taken.

#### Appendix

#### Determining Criteria Weights Using DEA-1 and DEA-2 Models

An example of obtaining the values of criteria for DMU17-N is as follows:

#### Data Availability

The data used to support the findings of this study are included within this article. However, the reader may contact the corresponding author for more details on the data.

#### Conflicts of Interest

The authors declare no conflicts of interest.