A Novel CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS Model for Safety Evaluation of Road Sections Based on Geometric Parameters of Road

: Trends of globalization very often cause the emergence of phenomena that asymmetrically a ﬀ ect the overall sustainability of the transport system. In order to predict certain situations and potentially be able to manage the transport system, it is necessary to manage risk situations and tra ﬃ c safety in a timely manner. This study has conducted an investigation which implies deﬁning the level of safety of a total of nine sections of two-lane roads. The main aim of the paper is to create a new multiphase model consisting of CRITIC (The CRiteria Importance Through Intercriteria Correlation), Fuzzy FUCOM (Full Consistency Method), DEA (Data Envelopment Analysis), and Fuzzy MARCOS (Measurement Alternatives and Ranking according to the COmpromise Solution) methods for determining the level of tra ﬃ c safety on road sections under the conditions of uncertainty. In order for the created model to be adequately applied, eight parameters were created, and they were classiﬁed through four inputs and four outputs. To calculate the signiﬁcance of the inputs, the CRITIC method based on the symmetric correlation matrix was used, and taking into account the nature of the outputs, the Fuzzy FUCOM method based on averaged values using the fuzzy Bonferroni Mean (BM) operator was applied to determine their weights. To determine the degree of safety, the DEA model was created. After that, the Fuzzy MARCOS method was used in order to determine the ﬁnal ranking of the remaining ﬁve sections of the road network. Finally, the veriﬁcation of results was performed through three phases of Sensitivity Analysis (SA).


Introduction
Traffic safety and potential risk management is certainly one of the important goals in order to potentially achieve a safe, balanced, and sustainable performance of road transport. Considering current issues and the constant search for an increase in overall traffic safety, exploitation characteristics are taken into account in a large number of studies, and noteworthy models that contribute to the increase of traffic safety have been developed. Geometric road design on which road safety factors are based especially refers to the positioning of symmetric physical road elements, cost (efficiency) analysis, reduction of harmful impact on the environment, traffic volume, road accessibility, etc. A road Symmetry 2020, 12,2006 2 of 28 with flexible physical elements can significantly obstruct traffic volume. In order to obtain accurate estimates of the effects of increased traffic volume, it is necessary to analyze all potential factors that contribute to the optimization of traffic flow volume. One of the important factors that affects traffic volume refers to curvature, i.e., the change of horizontal curve radius. Additionally, a significant impact on traffic is the impact of longitudinal gradient/fall, increased/decreased number of access points, increased/decreased road width, increased/decreased speed limit, etc. The increase in traffic volume entails a higher probability of traffic accidents on two-lane roads, but still the cause of accidents primarily depends on road geometry. Risk management on two-lane roads with the aim of protecting human lives represents a possible impact on reducing the number of traffic accidents. By analyzing potential causes of traffic accidents, the geometric characteristics of roads are especially emphasized. Influential geometric parameters refer to the density of access points, minimum radii of horizontal curvature [1], longitudinal gradients, as well as other primary parameters. In addition to geometric parameters, many studies relate to traffic indicators on two-lane roads. There are different conflicting opinions in the scientific community about the individual influence of geometric characteristics on traffic safety on two-lane roads.
Taking into account the relevance and great importance of the area which is the subject of the research, it is possible to define several goals of the study. The most important goal is to create a novel CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model as a contribution to the overall field of multi-criteria decision-making. The next goal of the study is to determine the state of road traffic safety on the basis of exploitation parameters using the created integrated model. The last goal of the study refers to the integration of CRITIC-Fuzzy FUCOM methods as a partial model, which represents the integration of subjective-objective methods for determining the weights of inputs and outputs in the proposed model. All three goals describe the originality of the developed integrated model and its significance. The advantage of the developed integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model is the fact that through the multiphase model, each item is treated as an important and eliminated possibility of inadequate ranking compared to similar approaches as we stated further in the next section.
This research is divided into a total of six sections. In Section 2, a review of the situation in the field of road traffic safety is provided implying the use of various models and their impact on the mentioned field. Section 3 includes an overview of the research flow with detailed explanations in terms of methodology. The steps of the applied methods are given by defining a new extended linguistic scale for the evaluation of variant solutions when the criteria are of a cost type. After that, Section 4 consists of empirical research. At the beginning of the section, the created inputs and outputs are explained, then a detailed calculation of the integrated model is presented. In Section 5, a sensitivity analysis is created and divided into three parts: Checking the change in influence of the four most important criteria, comparative analysis with three MCDM methods in a fuzzy form, and changing the size of an initial matrix. Section 6 summarizes the contributions, results, and guidelines for continuing the research.

Literature Review
A large number of studies have analyzed the minimum radius of horizontal curvature as a negative influencing factor on potential possibilities of traffic accidents [2][3][4]. Based on the research in Connecticut on 655 road segments, each with a length of 1 km, without intersections and with control on major road approaches, variables that significantly affect the frequency of head-on crashes have been identified: Speed limit, sum of absolute change rate of horizontal curvature, maximum degree of horizontal curvature, and sum of absolute change rate of vertical curvature [5]. Based on the research conducted in Western Sweden, the statistical analysis has shown that road geometric characteristics have a significant impact on the accident ratio per million vehicle kilometer [6]. By the analysis of the report conducted in Texas [7], the regression analysis has found that the following variables affects the prediction of traffic accidents: Average Annual Daily Traffic (AADT), lane width, Symmetry 2020, 12, 2006 3 of 28 shoulder width, and section length. In addition, the research conducted on 12 rural two-lane roads in Egypt at locations that consisted of a tangent element and the succeeding horizontal curve was based on measuring vehicle flows at each study site. A capacity estimation based on extrapolation between the traffic flow and density was used in the study. The relationship between capacity loss and geometric characteristics is based on the minimum radius of the horizontal curve, as an independent variable [8]. Based on the research [9] conducted on two-lane roads in Valencia, it has been shown that traffic risk is particularly affected by access density, average sight distance, average speed limit, and the proportion of no-passing zones.
In laboratory conditions, the study was conducted on 36 drivers who participated in a driving simulator that tried to show the effect of three road configurations. The statistical analysis confirmed that the driver behavior was only affected by road-sections and geometric elements, but not the road configurations. Drivers do not change their behavior when barriers are not present [10]. An investigation [11] in southwest of Montana at 15.7 miles yielded results showing a link between accidents (sliding off the road) and a minimum radius of horizontal curvature, even a straight roadway at speeds above the posted speed limit, per the road design equation. Additionally, a higher accident rate was proven as a function of the difference between the free flow speed and the speed depending on the radius of the curvature. Based on a study [12] on rural two-lane roads, it has been determined that the impact on traffic accidents on such roads is influenced by section length, traffic volume, road width, lane width, pavement condition, average minimum radius of horizontal curve, and grades of vertical slopes. Moreover, the study has found that the reduction in the number of accidents due to different safety measures varies by crash severity type, highway class, level of deficiency on roads, and the extent to which the deficiency is addressed.
According to the conducted studies [4,13,14], the independent variable that influences the increase of traffic risk in the procedures of predicting traffic accidents on two-lane roads refers to the total traffic volume. Short-term effects of changes in traffic flow parameters (sudden short-term changes in traffic volume) can cause an increased frequency of traffic accidents. Moreover, refs. [2,3,[14][15][16][17], the influence of the number of access points on traffic risk prediction is especially emphasized, where the probability of head-on crashes increases with the increase in density of access points. The analysis of the longitudinal gradient (rise/fall) on the possible prediction of traffic accidents on two-lane roads also showed the negative impact of this factor through one of 12 factors [18]. In addition, the study did not determine the direct dependence of visibility, wind speed, and temperature on the number of accidents.
DEA, MCDM models, and other techniques are often applied in the field of traffic safety, which is confirmed by the following studies. However, the question is to what extent such studies can provide precise output parameters, since they are most often based on individual approaches. Alper et al. [19] applied a model based on an individual DEA to determine the efficiency of 197 municipalities. The model consists of two inputs and 14 outputs that are unequally divided into two groups. The results showed that due to the weights of the input, it is possible for a more efficient municipality to be ranked lower. The advantage of the developed integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model is precisely the fact that through the multiphase model, each item is treated as an important and eliminated possibility of inadequate ranking. The importance of integrating and improving the DEA model is emphasized by the authors [20] that have performed research, in which the traffic safety efficiencies were estimated using the double frontier DEA to reflect both the optimistic and pessimistic points of view. Moreover, they noted that the results obtained by the conventional DEA are not entirely reliable considering only the efficient frontier. This research improves the disadvantages of such model due to involving a few different approaches in one functional model. Due to the shortcomings of the application of only the classical DEA, the authors [21] have performed a determination of the efficiency of 31 provinces using the two-phase DEA model.
The MCDM model based on the AHP method has been applied in [22] for a determination influence of traffic factor interaction on the rate of traffic accidents. In order to achieve an optimization of multi  [24] as an improvement degree of traffic safety in a way that researchers have identified prioritized black spots. The PROMETHEE (Preference Ranking Organization METHod for Enrichment of Evaluations) method was used in [25] for the determination of urban road safety in Spain, while in [26] the authors used a novel approach for the assessment of pavement conditions. The study represents the MIVES model (The Integrated Value Model for Structural Assessment) [27], which includes the VF concept (Value Function) and AHP method. Moreover, MCDM methods such as fuzzy DEMATEL [28][29][30] can be used in the field of traffic safety. For example, in [31] the influence of cultural variables affecting traffic rules violation and causing traffic accidents have been investigated.
In the developed CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model, four different approaches were applied, treating each of their advantages equally. In that way, a multi-phase functional model for assessing the condition of road sections from a traffic safety aspect was created. Figure 1 shows the research flow with the created methodology for calculating the road sections safety degree of the road network. The overall methodology consists of four phases, i.e., 17 steps. Each phase is causally linked with the next phase by presenting input parameters for the continuation of the model, i.e., the next phase.  [25] for the determination of urban road safety in Spain, while in [26] the authors used a novel approach for the assessment of pavement conditions. The study represents the MIVES model (The Integrated Value Model for Structural Assessment) [27], which includes the VF concept (Value Function) and AHP method. Moreover, MCDM methods such as fuzzy DEMATEL [28][29][30] can be used in the field of traffic safety. For example, in [31] the influence of cultural variables affecting traffic rules violation and causing traffic accidents have been investigated.

Methodology
In the developed CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model, four different approaches were applied, treating each of their advantages equally. In that way, a multi-phase functional model for assessing the condition of road sections from a traffic safety aspect was created. Figure 1 shows the research flow with the created methodology for calculating the road sections safety degree of the road network. The overall methodology consists of four phases, i.e., 17 steps. Each phase is causally linked with the next phase by presenting input parameters for the continuation of the model, i.e., the next phase.

The First Phase
The constant need to analyze and improve this area is the first phase of the research, which includes determining road sections where the research was conducted. The total number of road sections is nine: Doboj Novi-Doboj (Poljice), Doboj (Poljice)-Border of Republic of Srpska-Federation of Bosnia and Herzegovina (RS-FB&H), Doboj-Border of RS-FB&H (Karuše), Rudanka-Doboj, Johovac-Rudanka, Klupe-Teslić (Barići), Obodnik-Klupe, Šešlije-Johovac, and Teslić (Barići)-Border of RS-FB&H. Based on expert consultations and monitoring of similar studies, input and output parameters, a total of eight, were defined. As inputs, a list of four criteria was formed: Road section length expressed in kilometers, Average Number of Access Points per kilometer (AAP/km), AADT, and slope of the road. The outputs are: Average number of traffic accidents with fatalities, average number of traffic accidents with severely injured persons, average number of traffic accidents with slightly injured persons, and average number of traffic accidents with material damage. In addition, the database for traffic accidents in the last 5 years have been considered. All inputs belong to a benefit group since they need to be maximized, while the outputs have a cost character. Four input criteria were selected and used as available data of importance in the performed research. Moreover, it is possible to select a larger number of input criteria, which in this case, was not used due to technical-operational limitations and traffic characteristics on measuring road sections. The number of access points was determined by visual measurement on the measuring sections.

The Second Phase
During this phase of the research, weights are calculated and the significance of inputs and outputs is determined using two MCDM methods, subjective and objective: CRITIC [32,33] and fuzzy FUCOM [34]. Taking into account that the nature of inputs allows the application of an objective method for their determination, the CRITIC method is applied, while the Fuzzy FUCOM method based on averaged values using the fuzzy Bonferroni mean operator is applied to calculate output weights. The subjective method is applied to determine the weights of outputs, since their nature requires it. As the last step of this phase, the weights obtained by applying the two mentioned methods are averaged.

CRITIC Method
The CRITIC method for the purpose of this research can be described as follows.
Step 1. Development of the decision matrix including four inputs for traffic safety evaluation (X).
where (x ij ) represents the characteristics of ith alternative (road sections) on jth criterion.
Step 2. Normalization of the decision matrix: Step 3: Determination of the symmetric linear correlation matrix (m ij ) of all inputs for traffic safety evaluation. Step 4: Estimation of both Standard Deviations (SD) of the traffic safety inputs and its correlation with other traffic safety criteria. The weight of the jth criterion (w j ) is obtained as follows: where C j is the amount of information contained in input j and is determined as follows: where σ is the SD of jth criterion and is the correlation coefficient between the two criteria.

Fuzzy Bonferroni Mean (BM) Operator
Since five experts participated in the research, the values of the interval fuzzy vector of weight coefficients are aggregated using a fuzzy Bonferroni aggregator [35], Equation (6).
where e represents the number of experts, while p, q ≥ 0 are a set of non-negative numbers.

Fuzzy FUCOM Method (FUCOM-F)
The following section presents the FUCOM-F algorithm through four steps.
Step 1. Forming a set of criteria including four outputs for traffic safety evaluation.
Step 2. Determination the ranking of criteria based on experts' preferences according to the importance of four outputs for traffic safety evaluation.
where k represents the rank of the last-ranked output.
Step 3. Comparisons of the outputs for traffic safety evaluation using triangular fuzzy numbers (TFNs) using a fuzzy linguistic scale. Since a specific significance of the outputs is defined for each expert, the application of the fuzzy Bonferroni aggregator results in an averaged significance of the outputs for traffic safety evaluation. Based on the averaged significance of the criteria, fuzzy comparative significance ϕ k/(k+1) is calculated by applying Expression (8): Therefore, a fuzzy vector of comparative significance of evaluation outputs for traffic safety evaluation is obtained, Expression (9): where ϕ k/(k+1) represents the significance of the output of C j(k) rank as compared to the output of C j(k+1) rank.
Step 4. Calculation of the optimal fuzzy weights. In the fourth step, the final values of the fuzzy weight coefficients of the outputs for traffic safety evaluation ( w 1 , w 2 , . . . , w n ) T are calculated. The final values of the weight coefficients should satisfy the conditions defined by Equations (10) and (11): ϕ k/(k+1) represents the comparative significance of the outputs for traffic safety evaluation C j(k) and C j(k+1) . After that, it is necessary to calculate the values of the weight coefficients of the outputs for traffic safety evaluation ( w 1 , w 2 , . . . , w n ) T which satisfy that with the minimization of χ. Based on the previously mentioned steps, the final nonlinear model ( w 1 , w 2 , . . . , w n ) T can be set.

Calculation of the Final Criteria Weights
This step involves dividing the criterion values calculated by using the CRITIC and Fuzzy FUCOM methods by two to obtain the final values of inputs and outputs with the sum of one.

The Third Phase
In this phase, it is necessary to apply the DEA model. Two models for determining the degree of traffic safety on the observed sections of the road are formed: Input oriented and output oriented models. Then, based on the obtained results, the alternatives with values greater than one are defined and eliminated from further calculation. In order to manage the safety of traffic, it is necessary to know the existing situation, which can include ranking dangerous sections of the road [36]. In order to perform the final ranking of the remaining sections of the road, the fuzzy MARCOS method is applied [37]. To apply the fuzzy MARCOS, the input and output values obtained using the methods described above are included. After applying the complete methodology, a list of road sections is created, starting from the safest to the riskiest one.

DEA Model
Here, the two DEA CCR models [38,39] formed consist of the input oriented inputs for traffic safety evaluation (max) and the output oriented model (min) outputs for traffic safety evaluation. The input oriented model is: In this model, the Decision-Making Unit (DMU) consists of m inputs for each road section xij, while s represents outputs for each road section yij.
The output oriented model is: Finally, Equation (15) is applied to obtain the road sections that enter the final ranking model.

Fuzzy MARCOS Method
The fuzzy MARCOS method was developed by Stanković et al. [37] and consists of the steps shown in Figure 2. The fuzzy MARCOS method was developed by Stanković et al. [37] and consists of the steps shown in Figure 2.  Step 1. Creating an initial fuzzy decision matrix consists of eight criteria (inputs and outputs) and five road sections.
Step 2. Expanding the previous matrix with the Anti-Ideal Solution (AAI) and the Ideal Solution (AID), which denotes the worst and best performance of road sections: Step 3. Normalizing the initial fuzzy decision matrix: Step 4. Weighting the normalized decision matrix: Step 5. Calculation of the S i matrix: Step 6. Calculation of the degree of usefulness K i for all road sections: Step 7. Calculation of the fuzzy matrix T i : Determination of the fuzzy number D: Step 8. Defuzzification of fuzzy numbers: Step 9. Determination of the utility functions f K i for all road sections: Step 10. Calculation of the final utility function for all road sections: Step 11. Ranking road sections.

The Fourth Phase
In the fourth phase of the created methodology, SA is performed in three parts. The first refers to checking the change in the significance of inputs and outputs: where W nα represents the reduced value of the criterion whose weight changes, W β represents the original value of the criterion under consideration, while W n denotes the original value of the criterion whose value increases. Then, a comparative analysis is done with three other MCDM methods: Fuzzy WASPAS [40], Fuzzy SAW [41], and Fuzzy TOPSIS [42]. At the end of checking the results, a calculation with a change in the size of the initial matrix is performed. Scenarios in which the worst ranked alternative is eliminated from the initial matrix are formed.

Defining Input-Output Parameters
Based on expert consultations and monitoring of similar studies, input and output parameters, a total of eight, were defined as explained in Section 3.1. Table 1 shows the values of all DMUs according to the input parameters.  Table 2 shows the values of all DMUs in relation to the output parameters.

Analysis of the State of Road Sections
The   Figure 4 shows AADT for the period from 2010 to 2015, and it can be seen that there are no significant deviations in the values of AADT in that period. A slight increase in AADT, which is higher than 4.0% per year, can be noticed on the Šešlije-Johovac section in some years. Figure 5 shows the overview of arithmetic means of the rise/fall (gradient) of the sections. Since the Obodnik-Klupe section stands out, by its length, it can be classified in the category of mountain sections of two-lane roads with an average gradient of over 1.0%. Some sections, despite the markedly large gradient that exceeds the value of 0.5%, do not have a significant length and cannot be classified as mountain sections. Typical plain sections are Šešlije-Johovac and Johovac-Rudanka. Figure 6 shows the total number of traffic accidents by type of accidents, with the DMU 4 -Rudanka-Doboj section standing out negatively, since there were also accidents with fatalities in the period from 2015 to 2019. On average, one person is killed per year on that section. When observing the stated time period, the number of traffic accidents with fatalities is 1, 2, 1, 1, and 0, respectively. This location, i.e., the road section also stands out for the total number of traffic accidents, which is 258, where five are with fatalities, 17 with severe injuries, 51 with minor injuries, and 185 with material damage. When it comes to DMU 1 , DMU 2 , DMU 8 , and DMU 9 , it is important to point out that there were no accidents with fatalities. Figure 3 shows the length and arithmetic mean of the number of access points of the previously stated nine measuring sections of two-lane roads. The two-lane road, Obodnik-Klupe, has the largest section length (20.13 km), and the first section, Doboj Novi-Doboj (Poljice), has the smallest with a length of 1.47 km. The average number of access points is particularly evident on the Rudanka-Doboj section and is 35.39 AAP/km. The section with the most part passes through the settlement. Doboj Novi-Doboj (Poljice) has the lowest average number of access points per km, which is 5.63.       Figure 5 shows the overview of arithmetic means of the rise/fall (gradient) of the sections. Since the Obodnik-Klupe section stands out, by its length, it can be classified in the category of mountain sections of two-lane roads with an average gradient of over 1.0%. Some sections, despite the markedly large gradient that exceeds the value of 0.5%, do not have a significant length and cannot be classified as mountain sections. Typical plain sections are Šešlije-Johovac and Johovac-Rudanka.          Figure 7 shows the average number of traffic accidents in the period from 2015 to 2019 by the type of accidents and specified sections. These indicators of the number of accidents per year were analyzed and used in further data synthesis. It is clearly noticeable that DMU4 stands out in terms of the number of traffic accidents by all types.

Calculating the Criteria Weight
The calculation of input and output weights using CRITIC and Fuzzy FUCOM methods is presented below.

Calculation of Input Weights Using the CRITIC Method
Step 1. Development of the decision matrix X is included in Table 3. Step 2. Table 4 shows the normalization of the initial matrix. All inputs belong to a benefit group, therefore, normalization is calculated using Equation (2), for example: Step 3. The symmetric linear correlation matrix m ij is shown in Table 5. Step 4. In this step, objective criteria weights are obtained ( Table 6). In order to integrate these values with the output values, it is necessary to divide the specified values by two and obtain the final input values: According to the obtained results, it can be noticed that the most important input is AADT with the value of 0.143. It can also be noticed that all input values are in a small range, i.e., slightly higher compared to those less significant, which means that there is no input with a significantly higher advantage than others.

Determining the Significance of Outputs Using the Fuzzy FUCOM Method
After defining the criteria, in the second step, their ranking is performed. The criteria are ranked according to the following: Traffic accidents with fatalities (O 1 ) > Traffic accidents with severe injuries (O 2 ) > Traffic accidents with minor injuries (O 3 ) > Traffic accidents with material damage (O 4 ). In the next step (Step 3), based on decision makers' preferences, the significance of the criteria for each expert is defined (Table 7). It is important to note that each of the five experts has at least 15 years of experience in science and practice. All were previously included in the performing professional studies from this field and wrote scientific papers and reports. By applying the fuzzy Bonferroni operator, the averaged values of the criteria are obtained. Therefore, for the O 2 criterion, we obtain the aggregate significance, as follows: q=1 (1, 1.1, 1.35), (1, 1.15, 2.3), (1, 1.05, 1.15), (1.1, 1.25, 1.35 Therefore, we obtain the aggregate significance of the criteria (Table 8). By applying Expression (9), the comparative significance of the criteria is defined: Based on the vector of comparative significance, the limitations of the model are defined (12). By applying Expression (10), we define the first group of constraints: w C1 / w C2 = (1.02, 1.14, 1.465), w C2 / w C3 = (0.743, 1.131, 1.44), and w C3 / w C4 = (1.053, 1.387, 1.836). Based on Expression (11), the constraint arising from the condition of transitivity of relations w C1 / w C3 = (0.757, 1.289, 2.11) and w C2 / w C4 = (0.782, 1.569, 2.644) is defined. Based on the defined constraints, a model (12) is formed to determine the fuzzy optimal values of weight coefficients.
By solving the model, the optimal values of fuzzy weight coefficients

Application of the DEA Model
In the following part of the study, two DEA models are presented: Input and output oriented models in order to define sections of two-lane roads, i.e., DMUs which have in some way a satisfactory level of safety, and then they are included in a model of final ranking using the fuzzy MARCOS method. The formed algorithms of linear programming are solved by applying the Lingo 17 software. The following is an example of an oriented model using Equation (13) for the Doboj Novi-Doboj (Poljice) road section.

Input
Output DEA The following sections of two-lane roads: DMU 2 (3.784), DMU 3 (6.916), DMU 4 (6.348), DMU 5 (2.023) are sections which are further eliminated from the model due to poor performance in terms of traffic safety. It is necessary to make certain corrective measures on them in order to improve traffic safety.
Sections of two-lane roads with a value of 1.000 are included in the next part of the DMU 1 , DMU 6 , DMU 7 , DMU 8 , and DMU 9 models.

Application of Fuzzy MARCOS Method for Final Ranking of Road Sections
Based on the results of the DEA model shown in Table 9, it is determined which DMUs are excluded from further calculation: DMU 2 , DMU 3 , DMU 4 , and DMU 5 , while the remaining DMUs are relatively safe road sections. In order to determine the final ranking of DMUs, the fuzzy MARCOS method is applied. Table 10 shows the data for the remaining DMUs which are included further into the model and on the basis of which the mentioned method is applied. Since only the scale for beneficial criteria was developed in the original fuzzy MARCOS method, the fuzzy linguistic scale for cost criteria has been extended in this paper, which is shown in Table 11.  (7,7,9) (1,1,3) Extremely good EG (7,9,9) (1,1,1) Applying the scale, DMUs are evaluated according to inputs and outputs, which is shown in Table 12. The first part of the table presents linguistic variables, while the lower part of the list shows the evaluation of DMUs based on triangular fuzzy numbers.
After the evaluation using a linguistic scale and transformation into triangular fuzzy numbers by applying Equation (16), the Anti-Ideal Solution (AAI) is defined, i.e., the ideal solution by applying Equation (17), which makes an extended initial decision matrix. After that, the normalization shown in Table 13 is performed.    (27) and (28):  (Table 15). Calculation of the utility function of alternatives f K i is obtained using Equation (29): As the last step of the applied methodology, the final ranking of DMUs according to the descending sequence is performed.
After the application of the overall integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model, it can be noticed that the first section of two-lane roads with a value of 0.832 shows the highest level of traffic safety in terms of certain road network exploitation indicators. The second in a row is the eighth section, DMU 8 , with a value of 0.695, while the other three sections, DMU 7 , DMU 9 , and DMU 6 , have significantly lower values, which means that they are riskier for traffic participants.

Checking the Change in the Weights of the Criteria
In this part of the paper, the influence of the change in the weight values of the criteria is checked. Their change, especially when it comes to the most important parameters, can play a great role in the final ranking [43,44]. Therefore, a total of 24 scenarios were formed in which the values of the four most significant criteria were reduced by applying Equation (30). In the first six scenarios, the values of the most significant criterion (O 1 ) were reduced, and in scenarios S 7 -S 12 , the values of the O 2 output were reduced. Further, in scenarios S 13 -S 18 , the reduced weights of input I 3 were simulated, and in scenarios S 19 -S 24 , of input I 4 . It is important to note that the weight values of the mentioned criteria were reduced in the range of 15%-90%. The values of the criterion weights through new 24 scenarios are shown in Table 16. After forming the scenarios, the calculation starts again with the fuzzy MARCOS method and the new results (ranks) for each scenario are presented in Figure 8. From Figure 8, it can be stated that the model is slightly sensitive when new values of criterion weights are formed. This is the situation in four scenarios out of a total of 24, which is, in percentage, a change of 17%. In scenario S11, the first change in the ranks of DMUs occurs when the first-placed and second-placed alternatives change their positions. In the twelfth scenario, there is the same change with additional mutually changing positions of DMU7 and DMU9. The cause of the changes is the reduction of the second output as a criterion by 75%, i.e., 90%. Two other rank changes occur in scenarios S17 and S18 when DMU7 and DMU9 change their positions, which is a consequence of the From Figure 8, it can be stated that the model is slightly sensitive when new values of criterion weights are formed. This is the situation in four scenarios out of a total of 24, which is, in percentage, a change of 17%. In scenario S 11 , the first change in the ranks of DMUs occurs when the first-placed and second-placed alternatives change their positions. In the twelfth scenario, there is the same change with additional mutually changing positions of DMU 7 and DMU 9 . The cause of the changes is the reduction of the second output as a criterion by 75%, i.e., 90%. Two other rank changes occur in scenarios S 17 and S 18 when DMU 7 and DMU 9 change their positions, which is a consequence of the change of the third input by 75% and 90%. In other scenarios, DMUs retain their initial rank.

Comparison With Other MCDM Methods
The following, i.e., Figure 9, presents the results of the comparative analysis of the proposed integrated model with three methods in a fuzzy form: fuzzy SAW (FSAW), fuzzy WASPAS (FWASPAS), and fuzzy TOPSIS (FTOPSIS). From Figure 8, it can be stated that the model is slightly sensitive when new values of criterion weights are formed. This is the situation in four scenarios out of a total of 24, which is, in percentage, a change of 17%. In scenario S11, the first change in the ranks of DMUs occurs when the first-placed and second-placed alternatives change their positions. In the twelfth scenario, there is the same change with additional mutually changing positions of DMU7 and DMU9. The cause of the changes is the reduction of the second output as a criterion by 75%, i.e., 90%. Two other rank changes occur in scenarios S17 and S18 when DMU7 and DMU9 change their positions, which is a consequence of the change of the third input by 75% and 90%. In other scenarios, DMUs retain their initial rank.

Comparison With Other MCDM Methods
The following, i.e., Figure 9, presents the results of the comparative analysis of the proposed integrated model with three methods in a fuzzy form: fuzzy SAW (FSAW), fuzzy WASPAS (FWASPAS), and fuzzy TOPSIS (FTOPSIS).  Observing Figure 9, it can be seen that there is no change in the ranks compared to the initial rank, regardless of the applied method. Since it is the case, Figure 9 also shows the final values of each DMU after applying the comparative analysis. The values obtained with the fuzzy MARCOS method are in a larger range (0.180-0.832) compared to FSAW (0.357-0.741), FWASPAS (0.224-0.567), and FTOPSIS (0.259-0.786), which is another advantage of the proposed model. Figure 10 shows the results in relation to the change in the size of the initial matrix. Observing the initial rank obtained by the fuzzy MARCOS method, the worst-ranked alternative was eliminated and calculation was performed with a 4 × 8 matrix. In each subsequent set, the worst alternative was eliminated, so the size of the initial matrix decreased.

Influence of Dynamic Initial Matrix Formation
From Figure 10, it can be concluded that changes in the size of the initial matrix do not play any role, since the ranks are the same compared to the initial rank regarding all DMUs that enter the initial matrix. In addition, Figure 10 shows the values of the final results on the left side, so that the constancy of the previously obtained results can be observed.
After the conducted sensitivity analysis, in which results were verified through the three phases: New values of criterion weights through 24 scenarios, a comparative analysis with other methods, and an analysis of the change in the size of the initial matrix, it can be concluded that the model has good performance. Only the first phase showed slight sensitivity when changing the significance of particular criteria by 75% or 90%, while, in the remaining two phases, the stability of the original results was maintained. One of the advantages of the developed integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS approach is the fact that the model eliminated the possibility of inadequate ranking compared to similar approaches, for example, compared to DEA which was proven in the second part of the paper. The obtained results can have great significance in the planning and construction of traffic infrastructure due to the help of DMs to take into account these performances as inputs in such analyzes.
Symmetry 2020, 12, x FOR PEER REVIEW 24 of 28 Observing Figure 9, it can be seen that there is no change in the ranks compared to the initial rank, regardless of the applied method. Since it is the case, Figure 9 also shows the final values of each DMU after applying the comparative analysis. The values obtained with the fuzzy MARCOS method are in a larger range (0.180-0.832) compared to FSAW (0.357-0.741), FWASPAS (0.224-0.567), and FTOPSIS (0.259-0.786), which is another advantage of the proposed model. Figure 10 shows the results in relation to the change in the size of the initial matrix. Observing the initial rank obtained by the fuzzy MARCOS method, the worst-ranked alternative was eliminated and calculation was performed with a 4 × 8 matrix. In each subsequent set, the worst alternative was eliminated, so the size of the initial matrix decreased. From Figure 10, it can be concluded that changes in the size of the initial matrix do not play any role, since the ranks are the same compared to the initial rank regarding all DMUs that enter the initial matrix. In addition, Figure 10 shows the values of the final results on the left side, so that the constancy of the previously obtained results can be observed.

Influence of Dynamic Initial Matrix Formation
After the conducted sensitivity analysis, in which results were verified through the three phases: New values of criterion weights through 24 scenarios, a comparative analysis with other methods, and an analysis of the change in the size of the initial matrix, it can be concluded that the model has good performance. Only the first phase showed slight sensitivity when changing the significance of particular criteria by 75% or 90%, while, in the remaining two phases, the stability of the original results was maintained. One of the advantages of the developed integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS approach is the fact that the model eliminated the possibility of inadequate ranking compared to similar approaches, for example, compared to DEA which was proven in the second part of the paper. The obtained results can have great significance in the planning and construction of traffic infrastructure due to the help of DMs to take into account these performances as inputs in such analyzes.

Conclusions
In this paper, a new integrated model that includes different methods for both the subjective and objective evaluation has been created. In addition, integration with the DEA model and methods in a fuzzy form was performed. The model was created with the aim of determining the level of traffic safety on sections of two-lane roads of the road network in Bosnia and Herzegovina. By determining the level of traffic safety, it is necessary to establish sustainability in the performance of the system. The development of such a model is the main contribution to the field of traffic safety, but also to the problems of decision theory. After the first phase in which four inputs and four outputs were defined, for nine sections of two-lane roads, their quantitative values were determined. Recognizing the complexity of the research area, a set of methods was applied: CRITIC, Fuzzy FUCOM, DEA, Fuzzy MARCOS, forming an integrated model. In addition, the Fuzzy Bonferroni operator was applied to average the values representing the input parameters in the fuzzy FUCOM method. The objective CRITIC method was used to calculate the values of the four inputs, while Fuzzy FUCOM based on the fuzzy BM operator was used to calculate the weight values of the outputs. Then, the DEA model was applied, and it showed that the four sections of the two-lane roads, DMU 2 , DMU 3 , DMU 4 , and DMU 5 , represent the sections which are further eliminated from the model due to their poor characteristics in terms of traffic safety. Finally, Fuzzy MARCOS was applied to perform the final ranking of the five sections that proved to be safe in the DEA model. The results of the overall integrated CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model show that the first section of two-lane roads with a value of 0.832 has the highest level of traffic safety in terms of certain road network exploitation indicators. The second in a row is the eighth section, DMU 8 , with a value of 0.695, while the other three sections, DMU 7 , DMU 9 , and DMU 6 , have significantly lower values, which means that they are riskier for traffic participants. In the developed CRITIC-Fuzzy FUCOM-DEA-Fuzzy MARCOS model, four different approaches were applied, treating each of their advantages equally. In that way, a multi-phase functional model for assessing the condition of road sections from a traffic safety aspect was created. Contrary to that, some limitations of this study can be manifested through the small number of geometric inputs, small total infrastructure length that cover nine road sections, or the difficulty to implement measures of improvements.
Future studies can be conducted on a larger number of road sections with the possibility of implementing the proposed model. That means its practical application as a transfer of the developed technology to society. In addition, they can be modeled by observing different and more geometric input parameters.