Evaluation of Human Resources in Transportation Companies Using Multi-Criteria Model for Ranking Alternatives by Deﬁning Relations between Ideal and Anti-Ideal Alternative (RADERIA)

: Multi-criteria decision-making methods (MCDM) represent a very powerful tool for making decisions in different areas. Making a rational and reliable decision, while respecting different factors, is a challenging and difﬁcult task; MCDM models have a great impact on achieving this goal. In this paper, a new MCDM technique is presented—ranking alternatives by deﬁning relations between the ideal and anti-ideal alternative (RADERIA), which was tested for the evaluation of human resources (HR) in a transportation company. The RADERIA model has three key advantages that recommend it for future use: (1) the RADERIA model has a new approach for data normalization that enables deﬁning the normalization interval according to the judgments of a decision-maker; (2) an adaptive model for data normalization of the RADERIA model allows tough conversion into various forms of decreasing functions (linear, quadratic equation, etc.); and (3) the resistance of the RADERIA model to the rank reversal problem. Furthermore, in many simulations, the RADERIA method has shown stability when processing a larger number of datasets. This was also conﬁrmed by a case study with 36 alternatives, as considered in this paper. The results and veriﬁcation of the proposed new method were acquired through a comprehensive veriﬁcation of the complexity of the results. The complexity of the results was executed through (1) comparison with four other multi-criteria methods, (2) checking the resistance of the RADERIA model to the rank reversal problem, and (3) the analysis of the impact of changes in the measurement scale on the ranking results. enough to change the rank of the alternatives. This experiment shows the stability of the results of the RADERIA model, recommending it for solving real-world problems.


Introduction
Managing a successful business in any field depends primarily on the ability to maintain the needs and meet the expectations of all participants in the examined chain. It is necessary to define the goals and strive to achieve them. Logistics plays an important role in this, giving the answers to the question of how to rationalize, considering primarily the performances that have the most impact on rationalization. Lately, especially in the past few years, even companies in developing countries have recognized the importance of logistics in the rationalization of their entire system. According to [1], rationalization of logistics activities and processes is very important in doing business and achieving efficiency in a company. This is also confirmed by Koskinen and Hilmola [2]; according to them, the rationalization of the basic logistics subsystems is the key factor in achieving efficiency and effectiveness in companies. On the one hand, in the structure of the basic logistics subsystems, transportation presents the subsystem with the largest structure of logistics costs expressed on a daily basis, especially in large transportation companies. This is confirmed by Guasch [3] who emphasizes that the transportation costs represent the most important part of all logistics costs-percentage participation determined by Rantasila and Ojala [4]. They calculated that more than a third, or more precisely around 40%, of all logistics costs are made by transportation subsystems. On the other hand, this logistics subsystem enables achieving the purpose and goals of logistics, thanks to its importance and irreplaceability. Borzacchiello et al. [5] also confirm this, emphasizing that in today's modern economic system, transportation represents dynamic part, especially because it achieves great economic benefits and increases economic growth. In such companies, the rationalization of their activities plays a significant role in providing business efficiency [6]. Taking all of the above into consideration, it is necessary to constantly monitor the indicators that are responsible for making transportation costs, then focus on these performances and strive to improve them because they have an enormous impact on increasing the efficiency of transportation activities, as well as reducing time and costs.
To rationalize costs acquired through different transportation activities and processes, it is necessary to obtain adequate management and define adequate strategies [7]. According to Stević et al. [1], this can be seen especially in big and medium transportation companies, which is the case in this research. As already stated, in recent years, companies have realized the significance of logistics and they have started managing logistics, i.e., transportation costs. Furthermore, it is necessary to pay attention to the rolling stock, its structure, drivers, and management; according to Stojić et al. [8], the rolling stock represents the essential importance of the effective functioning of transportation companies. Higher rate mobility, which has become our everyday reality, causes a higher frequency of commodity flows, which contributes to the increasing importance of freight transport, and, according to Ghiani et al. [9], it plays a key role in today's economies.
Taking into consideration all of the above-mentioned, it is necessary to create great management and carry out adequate managing decisions in transportation companies. Basically, transportation companies that are working as logistics providers have to provide the required quality service for a minimal price. The outsourcing approach is, according to Liu et al. [10], widely present in all logistics aspects of business, especially in the field of transportation, characterized by significant and direct participation in total logistics costs. Transportation companies that are working as logistics outsourcing are largely dependent on the quality of the rolling stock, i.e., transportation vehicles [11], and therefore on HR, too, that is, the drivers.
If we observe the entire transportation market with a large number of transportation (logistics) companies, reducing the company's costs is a mandatory factor in order to achieve greater competitiveness. In transportation service companies, the largest share of costs is fuel consumption that is, consciously or unconsciously, influenced by the drivers themselves. One of the ways to decrease fuel consumption is the implementation of eco training for the drivers [12], or implementation of a reward system for drivers who achieve higher performance, i.e., the lowest fuel consumption. As a result of the training, reduction in fuel consumption can be higher than 10%; for example, in the study [13], 13.6% on average was achieved. In addition, the integration of these two approaches could have an adequate impact on reaching the target function. We should take into account the fact that, according to the performed study in [14], truck drivers were found to be happier during off-job activities. Therefore, managers of transportation companies should constantly work on increasing the satisfaction of their drivers. According to study [15], the high demands of their jobs can lead them to overlook their health, which can negatively impact every participant in supply chains.
In this paper, research was carried out in a company that is primarily engaged in international transportation and works on the EU territory. It was found that there are large differences in fuel consumption depending on the driver, even when it comes to the same transport routes. Therefore, an evaluation of the drivers was conducted in this paper, taking five criteria into consideration, where gas consumption is the most significant. A multi-criteria decision-making model (MCDM) was formed, and for these purposes, a new MCDM method called RADERIA was developed, which defines relations between an ideal and anti-ideal alternative and in that way, it respects their correlations, which can have a high impact on forming the final model, or on reaching a final decision. Solving the problem this way represents a new direction in the field of multi-criteria decision-making, which can have a positive impact on the entire field.
Based on everything outlined in this study, several important goals can be identified. The first goal is to form an adequate model for managing HR in transportation companies in order to achieve effectiveness and efficiency. This includes evaluating the work of the drivers on a monthly basis using five criteria, as well as additional money rewards for the best ranking drivers. The second goal represents the development of the new RADERIA method in order to get the MCDM method that also respects different variability factors and eliminates the rank reversal problem, which is present in many MCDM models. The third goal of this paper represents the contribution of Prospective Multiple Attribute Decision Making (PMADM) as defined in [16].
The present paper contains seven sections. The first section introduces the motivation of the research, objectives, and importance of the field. The second section shows the literature review in which similar problems in managing HR in transportation and logistics companies are analyzed. After that, in the third section, the methodology developed in this paper is presented. It implies defining all the steps for carrying out the research, implementing the LBWA method for determining the importance of data as well as the new RADERIA method. The fourth section consists of the study case where the advantages of the developed methodology are represented. In the fifth section, the validity analysis is conducted by comparing the developed methodology and other methodologies using different scenarios. The sixth section consists of a discussion on the obtained results and reviewing the advantages of the new method. Contributions and future research are described in the seventh section.

Literature Review
Micro, small and medium companies are popular and suitable for considering the implementation of the concept of business process management [17]. Based on the above, it is possible to see a necessity for creating a model for evaluating human performance [18], especially because, according to Gürbüz and Albayrak [19], human performance assessment is one of the key areas for analysis of the continuity of a company. This was confirmed in a study [20] where the authors noticed that companies can form and keep a competitive advantage by leveraging the expertise of HR development. One of the benefits of this study is the creation of a model for evaluating the work performance of drivers periodically based on multiple parameters, which has not been observed in the literature so far. Transportation, as a part of a supply chain, often represents an area of interest for various authors from different aspects, but not from the mentioned aspect.
According to Dubey and Gunasekaran, [21] can be a remark that drivers do not care about the environment or society, as manifested in their behavior. Therefore, they should have skills that make them more ecologically aware. Essentially, the goal of studying such problems and optimizing them leads to the achievement of environmental, economic, and social optimization aspects. Adequate choice and evaluation of staff in logistics represents an extremely significant factor of optimization because Klumpp and Abidi [22] noticed that one of the most significant tasks in logistics in the future is the evaluation and choice of workers in accordance with their competences and knowledge levels. Wu and Hou [23] carried out research on evaluating workers in a distributive center in Taiwan. The goal of their study was the development of employee performance estimation (EPE), which included three different models: direct performance determination (DPD), indirect performance determination (IPD), and performance score analysis (PSA). Workers on different levels were evaluated: 36 distinct first-line operators, 4 office-level managers, and 2 division-level managers. Chang [24] evaluated six employees in logistics companies according to four criteria: potential for future, corporate business achievement, organizational commitment, and working ability. He applied an integrated MCDM model, which consists of the Analytic Hierarchy Process (AHP) and technique for order preference by similarity to an ideal solution (TOPSIS). Paper [25] shows the modeling and evaluating of five different strategies in HR areas in an Indian logistics company using system dynamics modeling. This study also explored the impact of employees on the logistics performance index. After the study was carried out, the results were presented as a percentage for each investing strategy. Kampf and Ližbetinová [26] in their study used the AHP method in order to identify talents in logistics as one way to improve market competitiveness. Their goal was, after identifying such HR, to bring new values so that the company could optimize their business. Similar research was done by Qu et al. [27] and Badi et al. [28] who used the AHP method and triangle-definite weighted functions. The received results can be used as a reference for the companies in making the right and best decision on talents and provide long-term development strategies.
Xin et al. [29] state that talent is the most precious resource in the modern world, which was also confirmed in a study by Li et al. [30]. They carried out research whose aim was the application of fuzzy comprehensive appraisal in the human capital evaluation of logistics companies.
The purpose of carrying out such evaluations, according to Yue [31], is to understand the compatibility level between the demands of the workplace and the employees. In this research, according to the same author who applied the AHP method, it was given the possibility of acquiring a clearer perspective of the condition in the HR area in a logistics company by applying adequate measures and strategies. Samad [32] carried out research in the logistics sector on 292 employees. This study presents the findings of a study that examined the relationship between organizational justice and employees' work outcomes. In [33], research was performed in the area of logistics distribution and it was explained that HR, and investing in these intangible assets, can contribute to improving company performance and increasing competitiveness in the market. In a research conducted by Gürbüz and Albayrak [19], a model for evaluating employees in a pharmaceutical company marketing sector was presented. Analytical Network Process (ANP) and Choquet Integral (CI) were applied simultaneously to perform an adequate evaluation.

Methodology
In this part of the paper, a methodology for evaluating drivers in a logistics organization is presented. The methodology for evaluating drivers consists of three phases, as shown in Figure 1.
In the first phase, an expert team group has to be formed, and the set of alternatives is formed, as are criteria for evaluating alternatives. In the second phase, weighted coefficients are calculated using the LBWA (Level Based Weight Assessment) method. Weighted coefficients obtained in this phase are used further on in the RADERIA model for evaluating drivers. A detailed review of the steps in the LBWA-RADERIA multi-criteria model is shown in the following part.

The LBWA Model for Determining Weighted Coefficient of the Criteria
The LBWA model is based on a pairwise comparison of the criteria through forming a non-decreasing array of the criteria importance levels. Methods that subjectively determine the weighted coefficients of the criteria focus on the information obtained based on the judgments of the decision-makers [34]. This model implies the use of an original algorithm for grouping criteria based on the level of importance in order to eliminate the need for predefining an ordinary scale for pairwise comparison of the criteria. After grouping levels, the judgments of the decision-makers are defined [35] because modeling expert judgments in decision-making problems [36] is an interesting field of research. This way, it is possible to obtain relations between criteria that satisfy the condition of maximal consistency and also represent realistic relations, which are defined through expert judgments. The proposed model eliminates inconsistencies in expert judgments which can be seen in other famous models (Best-Worst Method-BWM and AHP) because the LBWA results are always highly consistent. We emphasize this because when increasing the level of consistency in the AHP and the BWM, expert judgments are distorted as a result of the violation of the consistency of the model. The mentioned problem is often due to the use of the nine-point scale, which has limited possibilities for expressing expert judgments [37]. The LBWA model has the following advantages: (1) enables rationally expressing the judgments of a decision-maker with a minimal number of comparisons, n − 1; (2) eliminates the shortcomings of nine-point that are described in many studies (Asadabadi et al.,

The LBWA Model for Determining Weighted Coefficient of the Criteria
The LBWA model is based on a pairwise comparison of the criteria through forming a non-decreasing array of the criteria importance levels. Methods that subjectively determine the weighted coefficients of the criteria focus on the information obtained based on the judgments of the decision-makers [34]. This model implies the use of an original algorithm for grouping criteria based on the level of importance in order to eliminate the need for predefining an ordinary scale for pairwise comparison of the criteria. After grouping levels, the judgments of the decision-makers are defined [35] because modeling expert judgments in decision-making problems [36] is an interesting field of research. This way, it is possible to obtain relations between criteria that satisfy the condition of maximal consistency and also represent realistic relations, which are defined through expert judgments. The proposed model eliminates inconsistencies in expert judgments which can be seen in other famous models (Best-Worst Method-BWM and AHP) because the LBWA results are always highly consistent. We emphasize this because when increasing the level of consistency in the AHP and the BWM, expert judgments are distorted as a result of the violation of the consistency of the model. The mentioned problem is often due to the use of the nine-point scale, which has limited possibilities for expressing expert judgments [37]. The LBWA model has the following advantages: (1) enables rationally expressing the judgments of a decision-maker with a minimal number of comparisons, n − 1; (2) eliminates the shortcomings of nine-point that are described in many studies (Asadabadi et al., 2019); and (3) the shown model has great flexibility while expressing expert judgments and it eliminates discrepancies which can occur with models that use the pairwise comparison technique. The LBWA algorithm is shown in Figure 2.
Mathematics 2021, 9, x FOR PEER REVIEW 6 of 27 2019); and (3) the shown model has great flexibility while expressing expert judgments and it eliminates discrepancies which can occur with models that use the pairwise comparison technique. The LBWA algorithm is shown in Figure 2.
Step 1: Determining the most important criteria from the criteria set   1 2 , , , n S C C C  Step 2: Ranking criteria in accordance with experts' judgments: 1 2 ...
Step 3: Grouping criteria according to the levels of importance: Level L 1 : Step 4: Defining importance criteria function Step 5: Calculating optimal values of the weighted coefficients of the criteria ( ) Detailed explanation of the steps that are used in the LBWA algorithm can be seen in the paper written by Zizovic and Pamucar [35].

Ranking Alternatives by Defining the Relationship between Ideal and Anti-Ideal Alternatives (RADERIA)
A multi-criteria model with n criteria 12 , ,..., n C C C is given. Criteria are aligned according to their importance based on subjective judgments of a decision maker, while two or more criteria can be of the same importance. We shall assume that there is m alternative which is denoted by 12 ij a value is added and the initial decision-making matrix ij mn Xa     is formed. The criteria can be of the maximization or minimization type. We shall assume that all criteria are of the maximization type (which does not minimize the generality). Based on the defined settings of the multi-criteria model, the following steps of the RADERIA model are presented.
Step 1. For each criterion  Detailed explanation of the steps that are used in the LBWA algorithm can be seen in the paper written by Zizovic and Pamucar [35].

Ranking Alternatives by Defining the Relationship between Ideal and Anti-Ideal Alternatives (RADERIA)
A multi-criteria model with n criteria C 1 , C 2 , . . . , C n is given. Criteria are aligned according to their importance C 1 > C 2 > C 3 > . . . > C n based on subjective judgments of a decision maker, while two or more criteria can be of the same importance. We shall assume that there is m alternative which is denoted by A 1 , A 2 , . . . , A m . To each of these alternatives A i (i = 1, 2, . . . , m) on the criterion C j (j = 1, 2, .., n), a ij value is added and the initial decision-making matrix X = a ij m×n is formed. The criteria can be of the maximization or minimization type. We shall assume that all criteria are of the maximization type (which does not minimize the generality). Based on the defined settings of the multi-criteria model, the following steps of the RADERIA model are presented.
Step 1. For each criterion C j (j = 1, 2, . . . , n), a decision-maker determines the ideal value for the alternative I j and anti-ideal value for the W j , where the following condition exists W j < I j (j = 1, 2, . . . , n). Anti-ideal value (W j ) represents the least acceptable value for the decision-maker on the given criterion, while the ideal value represents the most acceptable value on the given criterion. The values I j and W j can be defined using the set of the alternatives A i on the criterion C j , but mainly they are given the values which are not in the set of the considered alternatives.
Step 2. Defining the advantage of I j compared to W j ; A decision-maker defines how many times the ideal alternative I j is better than the anti-ideal alternative W j with respect to the criterion C j (j = 1, 2, . . . , n). Here, we should emphasize that the term "better alternative" is not the same as "bigger alternative". We stress this because these two terms are constantly used as synonyms. For example, if we observe two alternatives (A 1 and A 2 ) with respect to the criterion C 1 , A 1 = 5 and A 2 = 15 (assume that the values are defined based on a ratio scale), we can say that the alternative A 2 is bigger than the alternative A 1 by three times, exactly. At the same time, we can estimate that the alternative A 2 is better than the alternative A 1 by five times.
Let us have α j (α j > 1) that is a value which determines how many times I j is better than W j . The limit values α for every criterion C j are defined in the same interval. Simultaneously, the decision-maker defines k j points (k j ≥ 1) from the interval W j , I j , e 1,j , e 2,j , . . . , e k j ,j where e k j +1,j = W j < e k j ,j < e k j −1,j < . . . < e 2,j < e 1,j < I j = e 0,j , that satisfies the following: j times better than e 1,j I j is α (2) j times better than e 2,j ...
Step 3. Forming the advantage matrix I j compared to the alternative A i , X = α ij m×n : In this step, we define how many times the ideal alternative (I j ) is better than the alternative A i . Let us have the value for the alternative A i with respect to the criterion C j defined in the initial decision-making matrix X = a ij m×n . Then, if the value is a ij where e s+1,j < a ij ≤ e s,j (s = 0, 1, . . . , k j ), then it will be as follows: If a ij = e s,j , then the alternative I j is better than A i by exactly α (s) j times. If a ij < e s,j , then the alternative I j is α If a ij = W j , then α (i) j = α j , or I j is better than A i exactly α j times. This happens only when the alternative A i uses a ij for its limiting value from the intervals W j , e k j ,j .
The relation between I j and the alternatives A i and W j is shown in Figure 3. In Figure 3, we can see that the value () i j  can be defined using decreasing functions. In Figure 3a, the values () i j  are shown using a linear function, which is represented with Expression (1). However, values () i j  can be defined with other decreasing functions, In Figure 3, we can see that the value α (i) j can be defined using decreasing functions. In Figure 3a, the values α (i) j are shown using a linear function, which is represented with Expression (1). However, values α (i) j can be defined with other decreasing functions, such as with the quadratic function, etc. (Figure 3b).
After defining values α (i) j for each alternative A i with respect to the criterion C j (j = 1, 2, . . . , n), we get the advantage matrix of the alternative I j compared to the alternative A i (i = 1, 2, . . . , m) (Expression (2)).
Step 4. Defining weighting matrix U = u ij m×n . Let us have weighted coefficients for all criteria, where n ∑ j=1 w j = 1, w j > 0 (j = 1, 2, . . . , n) and satisfy the condition w 1 ≥ w 2 ≥ , . . . , ≥ w n . Multiplying elements of the matrix X = α ij m×n with the weighted coefficients of the criteria w j , we get the weighted advantage matrix of the ideal alternative compared to the alternative A i : For values u ij we get u ij = α ij · w j , (i = 1, 2, . . . , m; j = 1, 2, . . . , n).
Step 5. Determining the final value of the criteria function: Adding elements of the matrix U = u ij m×n by columns, we acquire final values of the criteria alternative functions (Expression (4)): It is obvious that, if for all values A i ∈ (I 1 , I 2 ,..., I n ) for each j (j = 1, 2, . . . , n), then In the case that A i for all values a ij has anti-ideal values, that is A i ∈ (W 1 , W 2 ,..., W n ), then: Then we can state that show how many times the alternative A i is weaker than the ideal alternative.
Step 6. Ranking alternatives: Ranking alternatives is done based on the values of the criteria function-Q(A i ). In order for the alternative to be as highly ranked as possible, it is preferable to have the lowest value Q(A i ): (1) Let us have A x and A y , two arbitrary alternatives, and let us have the following values Q(A x ) and Q(A y ); then, if Q(A x ) > Q(A y ), it is true that A x → A y . Sign "→" represents domination of the alternative A x compared to the alternative A y .
(2) If the condition Q(A x ) = Q(A y ) is satisfied, then it will be as follows: • From the matrix U = u ij m×n we eliminate the values of all alternatives based on the weakest criteria, and from such a reduced matrix we can observe the new values • If it is still Q(A x ) = Q(A y ) even after eliminating the weakest criterion, then the next weakest criterion is eliminated. We repeat this until we get If the equality Q(A x ) = Q(A y ) still exists even after eliminating all the criteria, then we analyze if any of the alternatives, A x and A y , on the criterion C j acquire values which are better than the value which was defined as the ideal alternative on the criterion C j . In the case that there is such an alternative, then that alternative is chosen as the better one.

•
If the equality Q(A x ) = Q(A y ) still exists even after eliminating all the criteria, then we can say that the alternatives A x and A y are equal.

Case Study and Experimental Results
The research was carried out in a transportation company for international transportation. For each subsystem in the company, elements for improvement were determined. In order to create an action plan for improving the performance of the company, it was concluded that HR represents the critical factor for the business that can be significantly improved. Even though there exists an adequate motivational system for the employees, i.e., truck drivers, HR represents the focus for improvement because the adequacy of the transportation depends mainly on them. It was concluded that there is a need for a periodical evaluation of the drivers, rewarding two to three of the best drivers, and also a need for creating such a multi-criteria model.
In order to find an adequate solution, a model for evaluating drivers is formed based on five quantitative and qualitative criteria. The first criterion (C1) represents fuel consumption expressed in liters per 100 km. This criterion is quantitative and of the costs type. The company manages very meticulous paperwork on every kilometer driven and fuel consumption, so this data was taken from the database. The second criterion (C2) is also quantitative and of the costs type, and it refers to the damage per km, expressed in money value damage on the vehicle, caused by driver's carelessness. This criterion includes all damages, from the small ones to the large ones, which are very rare. The other three criteria are of a qualitative character and of the useful type. The third criterion (C3) implies the adequate maintenance of the vehicle by the drivers, periodically examined by the managers of the company.
The fourth criterion (C4) represents the driver's ability to give adequate information on time. This mainly refers to the driver's ability to manage adequate communication with the dispatchers, especially in unexpected situations, which are very common in the transportation process. The fifth criterion (C5) represents loyalty that can be seen through the execution of all duties without complaint and also the driver's flexibility in executing their tasks.

Determining Criteria Weights Using the LBWA Model
For determining criteria weight, the LBWA model [35] was used, which allows reaching credible and reliable values of the weighted coefficients, forming a rational judgment and, at the same time, receiving credible results while making decisions. In the process of determining weighted coefficients, four experts participated, who are marked as E1, E2, E3, and E4, and who represent the dispatchers in the examined company. Given the fact that the company is very stable and has been very successful on the market for over 20 years, the dispatchers have at least five years of experience in dealing with this kind of task. Since it is a specific problem that refers to evaluating drivers, it is logical that the superiors who are in constant contact with them evaluate them, and they together resolve problems.
Using the LBWA methodology, the following weighted coefficients for each expert are acquired, which can be seen in the following steps.
Step 2. Grouping criteria based on levels of significance: According to the decisionmakers' judgments, the criteria are grouped by the following subsets/levels (see Table 1). The criteria are grouped by levels according to the experts' judgments, based on the rules determined in the LBWA algorithm. The criteria grouping, which is shown in Table 1, can also be presented in the following way: Step 3. Based on the condition λ = max{|S 1 |, |S 2 |, . . . , |S k |}, it is as defined the maximum value of the scale for the comparison of the criteria (λ).
From the conditions above, we can conclude that the scale for the comparison of the criteria for the experts E1 and E3 is I i p ∈ [0, 1] while for the experts E2 and E4, is In each preference criteria level, the comparison of the criteria is conducted using the values I i p ∈ [0, 1] and I i p ∈ [0, 2]. Experts' criteria comparisons by levels of significance are presented in Table 2. Step 4. According to the maximal value of the scale for the comparison of the criteria (λ(E1) = λ(E3) = 1 and λ(E2) = λ(E4) = 2), the elasticity coefficient (γ) is defined for each expert, where γ > λ, that is γ > 1 and γ > 2.
Step 5. Defining influence criteria function: We know that γ > 1 for the experts E1 and E3 the value γ = 2 is determined arbitrarily. Because for the experts E2 and E4, the value of the elasticity coefficient satisfies condition γ > 2, the value γ = 3 is arbitrarily defined. Using expression f (C j ) = γ/(i · γ + I j ) the influence criteria functions are calculated for all criteria (Table 3). Table 3. Influence criteria functions.
The example for calculating criteria functions for the expert E1 is shown in the following part: Similarly, the values of the criteria functions for the experts E2, E3, and E4 are calculated.
Step 6. Calculating the optimal values of the criteria weighted coefficients: Using expressions w 1 = (1 + f (C 2 ) + · · · + f (C n )) −1 and w j = f (C j ) · w 1 (where w 1 is the most favorable criteria), we acquired optimal values of the weighted coefficients of the criteria for each expert (Table 4). Calculation of the weighted coefficients of the criteria for expert E1 is shown in the following part. Using Expression (4), the value of the weighted coefficient of the most favorable criterion C 1 is achieved.
Using Expression (5), we get the value of the weighted coefficients of the remaining criteria: Finally, we get a vector of the weighted coefficient w E1 j = (0.457, 0.222, 0.130, 0.102, 0.088) T .

Evaluation of the Alternatives Using the RADERIA Model
The initial decision-making matrix consists of 36 alternatives, or drivers, who are evaluated. The alternatives are marked as A i (i = 1, 2, . . . , 36). Criteria C1 and C2 are of the minimizing type, while the criteria C3, C4, and C5 are of the maximizing type. Elements of the initial matrix that refer to criteria C1 and C2 are received through measuring real indicators, while the criteria C3, C4, and C5 are evaluated by the decision-maker. For evaluating, a five-step scale is used: 1-very bad; 3-bad; 5-average; 7-good; 9-very good. The total number of evaluated drivers is 36. After forming the starting decisionmaking matrix, in Table 5 the RADERIA method algorithm is used. Step 1.
For each criterion C j (j = 1, 2, . . . , 5) the ideal value of the alternative I j is defined, as is the anti-ideal value W j (Table 6). The values I j (22, 0, 10, 10, 10) and W j (40, 2, 1, 1, 1) are defined according to the decisionmaker's judgments based on the ideal values of the evaluating criteria.
In the following part, for each evaluating criterion, an advantage I j is defined compared to W j . The advantage I j compared to W j is presented with α Step 3: In the next part, the procedure for the development of the matrix X = α ij m×n is shown, in which elements α ij represent the advantage of the alternative I j compared to the alternative A i . In order to define the value α ij , for each criterion C j (j = 1, 2, . . . , 5), a unique point e j is defined (Table 7). j (i = 1, 2, . . . , 36; j = 1, 2, . . . , 5) for all elements of the initial decisionmaking matrix. Using Expression (1), the elements α ij of the matrix X are calculated. In the following part, a procedure for calculating elements of the initial matrix for the criterion C5 and alternatives A1, A2, A3, and A4 is shown: The shown relations also apply to the criteria C3 and C4 because the values of the alternatives on the mentioned criteria are at the same range of [1,9]. Based on the received values () i j  , the advantage matrix of the ideal alternative compared to the alternative i A is formed.
Step 4. The weighted advantage matrix of the ideal alternative compared to other alternatives is defined. Multiplying the matrix elements X (Table 8) with the weighted coefficients of the criterion, j w we get the weighted advantage matrix of the ideal alternative relative to the alternative i A (Table 9).   Figure 4.
The shown relations also apply to the criteria C3 and C4 because the values of the alternatives on the mentioned criteria are at the same range of [1,9]. Based on the received values α (i) j , the advantage matrix of the ideal alternative compared to the alternative A i is formed.
Step 4. The weighted advantage matrix of the ideal alternative compared to other alternatives is defined. Multiplying the matrix elements X (Table 8) with the weighted coefficients of the criterion, w j we get the weighted advantage matrix of the ideal alternative relative to the alternative A i (Table 9).
Steps 5 and 6: The final value of the criteria function and ranking of the alternatives is determined. Using Expression (4), we get the final values of the alternative's criteria functions Q(A i ) ( Table 9). Ranking of the alternatives is done based on the values of the criteria function Q(A i ). The value of the function Q(A i ) shows how much further the alternative A i is from the ideal alternative. For the alternative to be ranked as high as possible, it is desirable to have a lower value Q(A i ), that is, to be as close as possible to the ideal alternative.

Validity of the Results
In this section of the paper, the validity of the results is shown using the RADERIA model. The sensitivity analysis is conducted through three phases: (1) In the first phase, the validity of the results is done through comparison with other multi-criteria methods; (2) the second phase implies checking the stability of the RADERIA model in dynamic conditions, that is, the resistance of the model to the rank reversal problem; and (3) the third phase represents analyzing the influence of the change of the measuring scale on ranking results.

Comparison with the Other Multi-Criteria Methods
In this part of the paper, the results' validity is shown by comparing the new RADERIA method with other traditional multi-criteria techniques. The comparison was made with the following multi-criteria models: TOPSIS [38], VIKOR (VIsekriterijumsko KOmpromisno Rangiranje) [39], MULTIMOORA (Multi-Objective Optimization by a Ratio Analysis plus the Full Multiplicative Form) [40], and COPRAS (Complex Proportional ASsessment) [41]. The same set of alternatives and the same values of the weighted criteria were used in the TOPSIS, VIKOR, MULTIMOORA, and COPRAS methods. The results of using the above-mentioned methods can be seen in Figure 5.
The ranking results show that all multi-criteria models confirm that the alternative A12 ranked first. In the VIKOR model, as well as at the RADERIA model, the alternative A1 is second-ranked. Simultaneously, the alternative A1 is fourth-ranked in the COPRAS model and sixth-ranked in the TOPIS model, while it is seventh-ranked in the MULTI-MOORA model. Simultaneously, all the models suggest the set {A9, A18} as the set of the worst alternatives. Comparing the RADERIA model ranks with the other multi-criteria techniques, it can be observed that several models suggest no unique rank. The reason for having the differences in rankings is having many alternatives with not many differences between them.
To determine the relationships among the results acquired through different MCDM techniques, Spearman's rank correlation coefficient (SCC) is used. SCC is one of the useful and significant measures for determining statistical correlation when we have ordinal variables or ranking variables. In this paper, SCC is used in order to determine the statistical importance of the difference among the rankings which were acquired through the RADERIA model and other approaches, Table 10. The ranking results show that all multi-criteria models confirm that the alternative A12 ranked first. In the VIKOR model, as well as at the RADERIA model, the alternative A1 is second-ranked. Simultaneously, the alternative A1 is fourth-ranked in the COPRAS model and sixth-ranked in the TOPIS model, while it is seventh-ranked in the MULTI-MOORA model. Simultaneously, all the models suggest the set {A9, A18} as the set of the worst alternatives. Comparing the RADERIA model ranks with the other multi-criteria techniques, it can be observed that several models suggest no unique rank. The reason for having the differences in rankings is having many alternatives with not many differences between them.
To determine the relationships among the results acquired through different MCDM techniques, Spearman's rank correlation coefficient (SCC) is used. SCC is one of the useful and significant measures for determining statistical correlation when we have ordinal variables or ranking variables. In this paper, SCC is used in order to determine the statistical importance of the difference among the rankings which were acquired through th e RADERIA model and other approaches, Table 10. Despite the fact that the COPRAS, TOPSIS, and MULTIMOORA methods apply different ways of data normalization, a high correlation was achieved. Correlation with the COPRAS, TOPSIS, and MULTIMOORA models is in the interval [0.930, 0.968], while the correlation with the VIKOR model is 0.832. There are numerous studies in the literature [42][43][44][45] that show that the application of different normalization techniques can lead to  Despite the fact that the COPRAS, TOPSIS, and MULTIMOORA methods apply different ways of data normalization, a high correlation was achieved. Correlation with the COPRAS, TOPSIS, and MULTIMOORA models is in the interval [0.930, 0.968], while the correlation with the VIKOR model is 0.832. There are numerous studies in the literature [42][43][44][45] that show that the application of different normalization techniques can lead to rank differences in the applied MCDM techniques. This is the reason for the small correction between the ranks of the VIKOR and MULTIMOORA methods and the TOPSIS method. The VIKOR and MULTIMOORA methods apply linear normalization, while the TOPSIS method applies vector normalization; therefore, the small correlation is a consequence of the application of different normalization techniques.
The high correlation of the RADERIA method results with the results of other MCDM techniques proves the effectiveness and credibility of the RADERIA methodology. The mean SCC value for all of the models is 0.917. A lower correlation between the VKOR and RADERIA models is expected due to the specific algorithm of the VIKOR method. Namely, the VIKOR methodology takes into account the optimistic and pessimistic index when ranking alternatives. Thus, the final working alternative depends on a pessimistic and optimistic scenario, which is not the case with the RADERIA method. If the pessimistic index is eliminated with the VIKOR methodology, the results are obtained (SCC = 0.973), resulting in a significantly higher correlation with the RADERIA method. The results show that there is a high correlation of the rankings between the RADERIA model and other used techniques. According to some research [41], all the SCC values higher than 0.8 show a significantly high correlation. Because we have all the values for SCC significantly higher than 0.8 in this paper and the average value is 0.917, we can conclude that there is a high correlation between the suggested approach and other tested MCDM techniques. Based on the achieved results, we can conclude that the suggested rank is confirmed and credible.

Rank Reversal Problem
One way to verify the validity of the decision-making model is to construct the dynamic matrixes and then evaluate the results given by the model in newly formed conditions. Suppose the results show some logical inconsistencies, which can be seen through unwanted alternative rankings. In that case, there might be a concern that there is a problem with the applied method's mathematical apparatus.
To observe this problem more clearly, we shall examine the example of three candidates applying for a working position. We shall denote the candidates as A, B, and C, and we shall assume that one of the multi-criteria methods for candidate rankings was used. The applied multi-criteria model suggested the following ranking A > B > C, according to which the first candidate is the most suitable for the working position, and then candidates B and C, respectively. Then, let us assume that candidate B, who is in the second position, is replaced with a worse candidate, D. At the same time, candidates A and C are considered for the position with the same characteristics as before. It is necessary that the new set of alternatives (A, D, and C) be ranked, assuming that the criteria have exactly the same weight as before. The question is whether candidate A, after the ranking, would still be the best candidate.
The rank reversal problem was noticed for the first time by Belton and Gear [46] and they examined the use of the AHP model for ranking alternatives. They conducted simple research considering three alternatives and two criteria. They introduced the copy of a non-optimal alternative, which resulted in increasing the set of alternatives to four. After they evaluated the new set of alternatives, assuming that all the criteria were of the same weight as before, they noticed that the indication of the best alternative could change. They proved the possibility of a rank reversal in the AHP model. A few years later, Triantaphyllou and Mann [47] noticed that the rank reversal problem can occur even when a non-optimal alternative is replaced with the worse one, and not with a copy, which was the case in Belton and Gear's experiment. Later on, many authors have shown that this problem occurs in many multi-criteria models [48][49][50][51][52][53].
Taking all the above into consideration, testing the resistance of the RADERIA model to the rank reversal problem represents a logical step to validate the results of the model, and also the recommendation of the new multi-criteria model for solving real-world problems. In this section, an experiment was carried out to evalutate the resistance of the RADERIA model to the rank reversal problem. Thirty-five scenarios were formed in the experiment, where the change of elements in the decision-making matrix was simulated. For each scenario, the number of alternatives was changed, eliminating the worst alternative from further consideration. After defining the new set of alternatives, in the newly formed conditions, the ranking of the remaining alternative was done using the RADERIA model. This experiment has two goals: (1) observing the complexity of the received result in uncertain conditions and (2) analyzing the performance of the RADERIA model in dynamic initial decision-making matrix conditions.
In the initial scenario, the worst alternative, A18, was eliminated from further consideration. The RADERIA model was applied to the new set of 35 alternatives and the new ranking was achieved ( Figure 6). The new result confirms that A12 is still the best alternative, while A9 is the worst alternative. Further implementation of the presented approach allows acquiring the preferences of the remaining scenarios ( Figure 6). This experiment has two goals: (1) observing the complexity of the received result in uncertain conditions and (2) analyzing the performance of the RADERIA model in dynamic initial decision-making matrix conditions.
At the same time, the other MCDM techniques used in this paper were tested in the same conditions-TOPSIS, VIKOR, MULTIMOORA and COPRAS. In Figure 7, the results of using the TOPSIS, VIKOR, MULTIMOORA, and COPRAS models in a dynamic environment are shown.  Figure 7, the results of using the TOPSIS, VIKOR, MULTIMOORA, and COPRAS models in a dynamic environment are shown.
In Figure 6, we can see that the rank reversal problem does not appear in the RADERIA model, while in TOPSIS, VIKOR, MULTIMOORA, and COPRAS models ( Figure 7) the problem of rank reversal appears in many scenarios. Through all 36 scenarios in TOPSIS, VIKOR, MULTIMOORA, and COPRAS, the change in function criteria appears. That further leads to the change in the ranking of the alternatives. The rank of the best alternative (A12) does not change in the TOPSIS, MULTIMOORA, and COPRAS models. However, in the VIKOR model, after the first scenario, the alternative A12 falls to second place, after the seventh scenario, to third place, then finally in S32, it falls to fifth place. We have a similar situation with the change in the ranks of the remaining alternatives in the VIKOR model, as well as the TOPSIS, MULTIMOORA, and COPRAS models.
In the RADERIA model, during the experiment, the values of the alternatives (Q(A i )) remain unchanged, so rank reversal does not happen regardless of the number of alternative changes. In the RADERIA method, the alternative A12 ranked best through all 36 scenarios which confirms the complexity and stability of the RADERIA model in a dynamic environment. Besides that, it is shown that the RADERIA model is resistant to the rank reversal of the alternatives.

Changes of the Range in Measurement Scale
In multi-criteria models, decision-makers often express their judgments using different ranges in measurement scales. Pamucar and Cirovic [42] showed that the use of different ranges in measurement scales, in certain multi-criteria models, can result in rank reversal of the alternatives. In the following part, the analysis was conducted, showing how the change of the ranges in measurement scales affects the ranking results in the RADERIA model. As stated in the previous part of the paper, for presenting values of the criteria C3, C4, and C5, a 5-point scale was used in the interval [1,9]. From the interval [1,9], only the odd values were used. The transformation of the range in the measurement scales was done using the following function y = x/2 + 1/2; that is how the 5-point scale was achieved in the interval [1,5].

Discussion of Results
To improve the complexity of the model of multi-criteria decision making [54][55][56][57][58][59], the RADERIA model is proposed based on the relation between the alternatives. Relations between the alternatives are defined according to the reference point (ideal and anti-ideal alternatives) which enables comprehensive information on alternatives and making an objective decision. The RADERIA model, as a multi-criteria model for making decisions, represents a powerful and complex tool for optimizing multiple goals. Compared to other methods, this method is simple and effective, enables easy sorting and optimization of the analyzed goals. In comparison with other multi-criteria methods, this method has an effective mathematical tool that eliminates the rank reversal problem.
The RADERIA model improves the MCDM field, introducing an algorithm that eliminates the rank reversal problem. The RADERIA method integrates three starting points for making a complex decision: (1) defining reference points (ideal and anti-ideal values); (2) determining relations between alternatives and ideal/anti-ideal values; and (3) aggregation of the alternative values for expressing the estimated distance of the alternatives compared to the ideal solution. Moreover, the results achieved through the RADERIA method are more rational because of the new approach which enables aggregation of the estimated distance alternatives from reference points. In the following part, the advantages of the RADERIA model and the reason it stands out compared to other traditional multi-criteria models are outlined:

•
The adaptivity of the RADERIA model through the model transformation for data normalization: The model for normalization of the data which is applied in the RADERIA model represents a linear model for normalization (Figure 3a), which is based on a linear function. Adaptivity of the model can be seen in its transformation into other forms of decreasing functions, for example, quadratic functions ( Figure 3b). This characteristic allows the real presentation of expert judgments and making rational decisions.

•
The resistance of the RADERIA model to rank reversal problem: One of the most significant disadvantages of many MCDM methods is the rank reversal problem. The question of rank reversal was noticed by Belton and Gear [46] for the first time.
To this day, this issue has received great attention. Moreover, there are many studies [48,50,51] that analyze this phenomenon and show its importance in the process of making decisions. This phenomenon especially comes to expression in dynamic conditions of making decisions where the number of alternatives changes during the process of making a decision. This phenomenon can be noticed in many traditional MCDM methods. For example, the phenomenon of rank reversal is confirmed for AHP, ELECTRE (ELimination and Choice Expressing REality), TOPSIS, and PROMETHEE (Preference Ranking Organization METhod for Enrichment of Evaluations). Besides the mentioned methods, the rank reversal problem is also confirmed in this paper in MULTIMOORA, COPRAS, and VIKOR methods. However, in the same experiment, the RADERIA model shows resistance to rank reversal. Consequently, the RADERIA method shows significant stability and reliability of the results in a dynamic environment. It is also important to mention that in numerous simulations, the RADERIA method shows stability while processing a larger number of datasets. The study case also confirmed this with 36 alternatives, as examined in this paper.
In further analysis of the results, we compared the RADERIA method with other methods in the MCDM. The results showed mild deviations of the rankings compared to the examined methods (Table 10). The highest closeness of ranks is indicated with the COPRAS method. Moreover, a high correlation is shown with the TOPSIS and the MULTIMOORA models. A slightly lower correlation among the ranks can be seen in the VIKOR model with SCC = 0.832. These results were expected because the VIKOR model possesses a specific mathematical algorithm that examines the relationship between the optimistic and pessimistic results and finds a compromise between the two. The specificity of the VIKOR algorithm was also confirmed when determining correlation with other models: COPRAS-SCC = 0.712, TOPSIS-SCC = 0.644 and MULTIMOORA-SCC = 0.620. When observing an optimistic solution in VIKOR, then the correlation among the ranks increases from SCC = 0.832 to SCC = 0.966.
From the above analysis, we can conclude that the RADERIA method represents a new approach for rational decision-making that processes information on preferences in a new and original way.
The new data normalization model contributes to a realistic and stable assessment of decision-making problems. This further improves the accuracy and reliability of the results, which is especially important for solving real-world problems in various fields. As MCDM seeks to implement a reliable and rational decision-making model, the RADERIA model is the logical choice for future multi-criteria solving of problems.

Conclusions
The results from the conducted research and study make a double contribution, which can be seen through the possibility of acquiring experts and scientific optimal approximation. The expert contribution represents the possibility of rationalization of the costs that appear in the transportation process alone. Creating a motivational model where every driver will try to increase their engagement, on the one hand, and decrease fuel costs, on the other hand, affects the entire organization and achieves a more efficient business. This is possible if all the factors are satisfied by the drivers. Companies in developing countries, as is the case in Bosnia and Herzegovina, do not adequately understand the role and the importance of HR because they mainly invest in assets and equipment, and much less into human resources. This research should serve as one of the most important in creating evaluating and self-evaluating models for companies in Bosnia and Herzegovina.
In this paper, two innovative multi-criteria models are suggested (the LBWA model and the RADERIA model) that together form a hybrid LBWA-RADERIA multi-criteria model. The LBWA model was used for determining weighted coefficients of criteria, while the RADERIA model was used for evaluating drivers in the analyzed transportation company. The original approach for grouping criteria by importance levels in the LBWA model enables easier defining relations among criteria and maximal consistency of the results. The LBWA model is used for expressing real relations among criteria, without approximations of the expert judgments. Moreover, the LBWA model eliminates deviations of expert judgments that appear after using scales with a small range of values [37]. By eliminating this kind of problem, we also eliminate the approximation of expert judgments. This further has an impact on the distortion of the consistency of the model, i.e., distancing weighted coefficients from the optimal values. Besides the LBWA model, a new multi-criteria model for evaluating alternatives is suggested-the RADERIA model. The RADERIA model has shown a high level of flexibility and complexity of the results, which makes it recommended for solving other real problems in different areas. The proposed model is based on expert knowledge, which enables expressing expert judgments in accordance with the specificity of the observed problem. The original mathematical tool of the RADERIA model helps to convert data from the initial decision matrix into any interval, which is suitable for Mathematics 2021, 9, 976 23 of 25 rational decision making. This adaptability can help managers use this methodology for making decisions in different, real situations.
In accordance with the foregoing, we can outline the importance of this paper: (1) developing the new multi-criteria methodology, the LBWA-RADERIA model, for solving real problems in the business world; (2) reviewing the new RADERIA method, which is based on the rational defining relations between the ideal and anti-ideal alternative; (3) developing the new multi-criteria model which eliminates the rank reversal problem and which offers reliable results for achieving rational decisions; and (4) developing the original universal problem for evaluating HR, which enables rational managing with human potentials in a company.
The new multi-criteria methodology, which is suggested through the RADERIA model, is validated through three phases: (1) in the first phase, comparison with the TOP-SIS, VIKOR, MULTIMOORA, and COPRAS models was conducted. The results showed a high level of correlation between the results of the RADERIA model and other tested models; (2) in the second phase, testing of the RADERIA, TOPSIS, VIKOR, MULTIMOORA, and COPRAS models in dynamic conditions was conducted, that is, the resistance of the model to the rank reversal problem was tested. The testing showed that the RADERIA model eliminates the rank reversal problem, while other models (TOPSIS, VIKOR, MUL-TIMOORA, and COPRAS) do not give reliable results in dynamic conditions; and (3) in the third phase, the validity of the results was executed in the case when there was a change in the measurement scale of the qualitative criteria. The RADERIA method showed complexity in the case of different measurement scales. It showed that a slight change occurs in the values of the criteria functions, but the rankings remain the same. The stability of the RADERIA method is confirmed, making it recommendable for use in future research to address multi-criteria problems.
The directions of future research should be directed toward the implementation of the RADERIA method for solving other real problems. Additionally, one of the directions of future research should be the extension of the RADERIA method in a fuzzy and rough environment [60], where the application of linguistic variables would allow rational expression of human judgments. In addition, the RADERIA method can be used with other methods for determining the weight of the criteria. It is possible to generate random data and perform a correlation study in order to prove the robustness of the RADERIA model additionally. This could be a direction for future research, also.