A New Integrated Fuzzy Approach to Selecting the Best Solution for Business Balance of Passenger Rail Operator: Fuzzy PIPRECIA-Fuzzy EDAS Model

: The analysis of operations of the passenger tra ﬃ c operator in the Republic of Srpska (RS) showed that the volume of passenger transport has, for the last ﬁfteen years, been in constant decline. It is of particular importance that the operator has, year after year, recorded a negative balance of business. The way out of the current unfavorable situation in the sector of passenger tra ﬃ c is based on the application of Public Service Obligation (PSO) based on the Regulation 1370 / 2007. In order to solve the problems, seven realistically possible variants have been identiﬁed. This paper deﬁnes the criteria for selecting the best variant, as well as a new integrated fuzzy model for the selection of the best variant that will enable the operator to make a proﬁt. To deﬁne the weights of criteria in this paper, we have used the fuzzy PIvot Pairwise RElative Criteria Importance Assessment (F-PIPRECIA) method, while for ranking and selection of the best variant, we have used the Fuzzy Evaluation based on Distance from Average Solution (F-EDAS) method. Results show that the seventh variant: “Increase in revenue from ticket sales and PSO services and reduction in costs“ is the best solution in current conditions. Validation tests are performed with di ﬀ erent scenarios and approaches and show that the model is stable. A validity test was created consisting of variations in the signiﬁcance of model input parameters, testing of reverse rank, applying the fuzzy Measurement Alternatives and Ranking according to the COmpromise Solution (F-MARCOS), fuzzy Simple Additive Weighing (F-SAW) method, and fuzzy Technique for Order of Preference by Similarity to Ideal Solution (F-TOPSIS). As a part of the validation tests, Spearman’s coe ﬃ cient of correlation (SCC) in some scenarios is performed and weights of the criteria have been obtained using the Fuzzy Analytic Hierarchy Process (F-AHP) and Full Consistency Method (FUCOM). S.M.;


Introduction
One of the most important factors for the functioning and development of cities and regions according to Stojić et al. [1] is the public transport of passengers. The transport policy developed by the

Brief Background
In the paper about public suburban transport in Germany, Beck [4] analyzes the state of the so-called commercial and non-commercial transport. In doing so, he notes that, after a decade of stagnation due to non-commercial transport, in the performance of the PSO, there is a positive change and the intensification of competition. The methodology for assessment of the future transport needs in PPT by Rojo et al. [5] is upgraded by the inclusion of the subjective value of time and readiness users pay for the improvement of services in order to determine the optimal concept of PSO. The system is optimized in two ways: With and without considering the economic business of the company in the function of the goal. Vesković et al. [6] used fuzzy logic for the assessment of the liberalization of rail passenger traffic on the example of Serbia, and one of the criteria in the model for assessment is PSO. Nash et al. [7] used quantitative and qualitative methods to investigate the impact on the cost of the vertical separation of railways in cases of a radical approach to restructuring. They are suspicious that reforming the railways through vertical and horizontal separation leads to cost savings. They state that precisely determining the methods and control of the distribution "of state money," subventions (PSO, maintenance and infrastructure development) have primarily led to cost reductions. In order to achieve the aims defined by the overall transport policy according to Ibarra-Rojas and Rios-Solis [8], cities and municipalities choose to subsidize PPT. These aims are different and range from providing transport options to all social categories to increase mobility for all residents. As a special advantage of these systems, Tirachini and Hensher [9] and Kim and Schonfeld [10] point out that the implementation of such a transport policy reduces the need for the use of personal vehicles. This, in turn, offers the opportunity to better manage urban space and transform the environment for the sustainable development of urban communities. In his paper, Van Reeven [11] developed a model aiming to demonstrate that the costs on the principle of consumer spending time do not provide justification for public transport subventions. PPT subventions are common in developing countries and are often justified by the availability of traffic accessibility, but not efficiency. In view of this justification, it is of interest to know how to use and distribute transport subventions.
To understand the idea behind public transport subsidies, Vuchic [12] and Hanson and Giuliano [13] emphasize that cities and municipalities do not subsidize operators, but the actual public transport service offered to citizens. In the absence of a subsidy, carriers are forced to charge the full cost of transportation for passengers through the price of tickets, which would lead to a significant reduction in transport demand and thus a decrease in traffic supply. Such a transport strategy implies, on the one hand, reduced mobility and, on the other, increased citizen dissatisfaction.

Proposed Methodology
The multi-criteria decision-making (MCDM) methods are widely used for the facilitation of the decision-making process in various fields [14][15][16]. The original developed MCDM methodology shown in Figure 1 was applied for selection of the best solution for the business balance of the passenger rail operator. As part of the first phase of the research, data were collected. After that, an adequate database on transport policy was created in order to obtain and to analyze their effects on the business operators. Based on the collected data and created base, the forming of the MCDM model represents the second phase of the proposed methodology. Five most important criteria, explained in detail in the further text, were considered, while seven different variants were identified. Based on such parameters, an initial fuzzy decision matrix was formed. The third part of the methodology represents the most important part of the research and consists of two steps that are causally linked both to each other and to the elements of the following phase. These steps represent the development of an original integrated fuzzy MCDM model. First, the significance of the criteria was determined using the F-PIPRECIA method [17] according to the assessment of three decision-makers. Evaluation of various variants for selecting the best solution for the business balance of the passenger rail operator was performed using the F-EDAS method [18]. The fourth phase of the methodology involves the validation and sensitivity analysis of the proposed model. It is implemented throughout a few steps, where the first step relates to variations in the significance of the criteria. All individual approaches are individually included in the calculation of the F-EDAS method and a comparative analysis is given with respect to the proposed model. Testing the influence of dynamic factors-of the reverse rank and calculation of the criteria weights using the Fuzzy Analytic Hierarchy Process (F-AHP) [19] and Full Consistency Method (FUCOM) [20] methods-is also a part of the validity test. The next step includes the comparison of the developed model with three other fuzzy MCDM methods: fuzzy Measurement Alternatives and Ranking according to the COmpromise Solution (F-MARCOS) [21], fuzzy Simple Additive Weighing (F-SAW) [22], and fuzzy Technique for Order of Preference by Similarity to Ideal Solution (F-TOPSIS) [23]. As part of the first phase of the research, data were collected. After that, an adequate database on transport policy was created in order to obtain and to analyze their effects on the business operators. Based on the collected data and created base, the forming of the MCDM model represents the second phase of the proposed methodology. Five most important criteria, explained in detail in the further text, were considered, while seven different variants were identified. Based on such parameters, an initial fuzzy decision matrix was formed. The third part of the methodology represents the most important part of the research and consists of two steps that are causally linked both to each other and to the elements of the following phase. These steps represent the development of an original integrated fuzzy MCDM model. First, the significance of the criteria was determined using the F-PIPRECIA method [17] according to the assessment of three decision-makers. Evaluation of various variants for selecting the best solution for the business balance of the passenger rail operator was performed using the F-EDAS method [18]. The fourth phase of the methodology involves the validation and sensitivity analysis of the proposed model. It is implemented throughout a few steps, where the first step relates to variations in the significance of the criteria. All individual approaches are individually included in the calculation of the F-EDAS method and a comparative analysis is given with respect to the proposed model. Testing the influence of dynamic factors-of the reverse rank and calculation of the criteria weights using the Fuzzy Analytic Hierarchy Process (F-AHP) [19] and Full Consistency Method (FUCOM) [20] methods-is also a part of the validity test. The next step includes the comparison of the developed model with three other fuzzy MCDM methods: fuzzy Measurement Alternatives and Ranking according to the COmpromise Solution (F-MARCOS) [21], fuzzy Simple Additive Weighing (F-SAW) [22], and fuzzy Technique for Order of Preference by Similarity to Ideal Solution (F-TOPSIS) [23].

The Position of Public Transport Services (PTS) for Passengers by Rail in the Transportation System of the Republic of Srspka
The Railways of the Republic of Srpska (RRS) have been established as a public transport company, and it is important to emphasize that by "under the railway traffic of interest for the Republic of Srpska," we mean "railway public passenger transport." Irrespective of commercial interest, RRS must have at their disposal adequate capacity (material and human) and organizational conditions for the provision of public transport services (PTS) for passengers. Therefore, the authorities of the Republic of Srpska exert significant impact on the results of operations and the balance sheet of the company. RRS, and the segment of the company that deals with passenger transport (passenger transport operations) in particular, must establish an original system of determining results (revenues, expenses, profit-loss) on the grounds of the public transport of passengers. The opening of the railways to competition in the market of transport services brought about a separation of management and accounts (balance sheet) of infrastructure and transport. Consistent realization of this process means that "RRS shall-through a special type of bookkeeping-present to its founder the state and the railway infrastructure costs compared to the costs of operators." Separate reporting of costs is aimed at expressing the impact of business segments upon the operating results, which are determined by the balance sheet. Therefore, it is necessary to separately determine the balance of infrastructure and the balance of transport (assets, debts, obligations, liabilities, equity, revenues, expenses, results), as well as the consolidated balance of the corporation. In its efforts to provide for the traffic of interest for the Republic of Srpska, the government participates through partial financing. This means that RRS provide funding for a part of the public transport system that is of interest for the Republic of Srpska. According to Gangwar and Raghuram [35], one of the options is structuring public private partnerships. The volume of passenger transport is in a constant downward trend, and the largest volume of transport was recorded in 1996, amounting to 1,648,000 of transported passengers, while in 2009, the RRS transported no more than 368,289 passengers. The negative trend has continued in the years to come, so in the last two years, the annual number was at the level of about 150,000 passengers. Financial results regarding passenger traffic have been made according to the planning documents: Annual report of RRS for 2014, and business plan for the period 2012 to 2014.
The revenue and expenditure plan in passenger traffic is projected at the level of the financial loss of over −19,000,000 KM for each considered year, which is why the plan of inflows and outflows of funds remains at the level of loss of −26,869,280 KM in 2012 to −38,616,384 KM in 2014. The increase in expenditure in 2012 was by 67% higher compared to 2011, amounting to 10,371,589 KM, and the revenues were lower by about 30%, i.e., by 2,559,281 KM. The increase in expenditure in the said amount was the result of an increase in the following: Cost of fees for access and use of railway infrastructure in the amount of KM + 5,458,404 (+52%); cost of wages, salaries, and other employee benefits in the amount of KM +2,996,647 (+29%), cost of materials for the work in the amount of KM + 1,020,676 (+10%), and costs of production services in the amount of + 966,014 (9%) (Figure 2).
The problem with the above, i.e., the problem with the operations in the reported period with a huge financial loss, lies in reduced business competitiveness of RRS as the operator at a future liberalized market of transport services, and, therefore, its uncertain business future. In this sense, the "experts" of the Railways of the Republic of Srpska have reduced passenger traffic for the 2008/2009 timetable by eliminating 22 passenger trains that were, by internal calculations, within the area of unprofitable business. Reactions of passengers to this move were completely understandable, so the reduction in the number of trains by 25% (from 76 to 54 trains a day) led to a reduction in the number of passengers for close to 50%, or to be more exact, by 45% in that same (first) year when they implemented the reduction in the number of trains (from 635,000 annually to 368,000). This trend of reducing the number of passengers due to an unsatisfactory timetable and reduced frequency of trains was carried out on almost all routes. The authorities of the Republic of Srpska noted that by this move, they achieved a reduction in operating costs of about 2,000,000 KM but failed to note the loss and reduction in income due to a drastic reduction in the number of passengers.
The downward trend in train numbers has led to an increase in the company's financial losses. It is true that the cost of doing business has been somewhat reduced ( Figure 2), but revenue has fallen significantly, leading to greater financial losses ( Figure 3). The financial loss of the company in 2011 amounted to 6,969,205 KM. In 2012, expenditures of KM 10,441,741 were still high, although they were somewhat reduced, but revenues were significantly lower and decreased by over KM 2.5 million so that the company's negative balance in fiscal 2012 increased to 19,970,227 million KM.  Reactions of passengers to this move were completely understandable, so the reduction in the number of trains by 25% (from 76 to 54 trains a day) led to a reduction in the number of passengers for close to 50%, or to be more exact, by 45% in that same (first) year when they implemented the reduction in the number of trains (from 635,000 annually to 368,000). This trend of reducing the number of passengers due to an unsatisfactory timetable and reduced frequency of trains was carried out on almost all routes. The authorities of the Republic of Srpska noted that by this move, they achieved a reduction in operating costs of about 2,000,000 KM but failed to note the loss and reduction in income due to a drastic reduction in the number of passengers.
The downward trend in train numbers has led to an increase in the company's financial losses. It is true that the cost of doing business has been somewhat reduced ( Figure 2), but revenue has fallen significantly, leading to greater financial losses ( Figure 3). The financial loss of the company in 2011 amounted to 6,969,205 KM. In 2012, expenditures of KM 10,441,741 were still high, although they were somewhat reduced, but revenues were significantly lower and decreased by over KM 2.5 million so that the company's negative balance in fiscal 2012 increased to 19,970,227 million KM. Reactions of passengers to this move were completely understandable, so the reduction in the number of trains by 25% (from 76 to 54 trains a day) led to a reduction in the number of passengers for close to 50%, or to be more exact, by 45% in that same (first) year when they implemented the reduction in the number of trains (from 635,000 annually to 368,000). This trend of reducing the number of passengers due to an unsatisfactory timetable and reduced frequency of trains was carried out on almost all routes. The authorities of the Republic of Srpska noted that by this move, they achieved a reduction in operating costs of about 2,000,000 KM but failed to note the loss and reduction in income due to a drastic reduction in the number of passengers.
The downward trend in train numbers has led to an increase in the company's financial losses. It is true that the cost of doing business has been somewhat reduced ( Figure 2), but revenue has fallen significantly, leading to greater financial losses ( Figure 3). The financial loss of the company in 2011 amounted to 6,969,205 KM. In 2012, expenditures of KM 10,441,741 were still high, although they were somewhat reduced, but revenues were significantly lower and decreased by over KM 2.5 million so that the company's negative balance in fiscal 2012 increased to 19,970,227 million KM.

Problem Identification and Solving Methodology
The analysis of the passenger traffic subsystem showed the following characteristics: From 1996 until today, the number of passengers has been in constant decline, the number of passenger trains has decreased in domestic traffic by 28 trains, in inter-entity transport by 18 trains, and by 8 trains in international transport. The railway fee for the infrastructure in domestic services amounts to 4,176,295 KM, in inter-entity transport amounts to 601,836 KM, and in international transport amounts to 1,603,022 KM. The other elements of the RRS business operations are shown in [36]. The problem-solving methodology is based on: increase in revenue from direct ticket sales and increase in revenue from agreements on PSO, 2.
The expected result of the mentioned activities according to the given methodology should be sustainable business operations. Figure 4 shows the cost realization plan for the next fiscal year, as well as the perception of cost coverage by government revenues and government subsidies (PSOs), and a model of the long-term business stabilization goal ( Figure 5).

Problem Identification and Solving Methodology
The analysis of the passenger traffic subsystem showed the following characteristics: From 1996 until today, the number of passengers has been in constant decline, the number of passenger trains has decreased in domestic traffic by 28 trains, in inter-entity transport by 18 trains, and by 8 trains in international transport. The railway fee for the infrastructure in domestic services amounts to 4,176,295 KM, in inter-entity transport amounts to 601,836 KM, and in international transport amounts to 1,603,022 KM. The other elements of the RRS business operations are shown in [36]. The problem-solving methodology is based on: 1. increase in revenue from direct ticket sales and increase in revenue from agreements on PSO, 2. reduction in expensing, i.e., operating costs.
The expected result of the mentioned activities according to the given methodology should be sustainable business operations. Figure 4 shows the cost realization plan for the next fiscal year, as well as the perception of cost coverage by government revenues and government subsidies (PSOs), and a model of the long-term business stabilization goal ( Figure 5). The aim of the problem solving is to bring RRS as the operator in passenger traffic into the domain of positive operations, thus providing the necessary conditions for successful operation at a liberalized transport market. In [6], details about the model for liberalization in Serbia can be found.

Problem Identification and Solving Methodology
The analysis of the passenger traffic subsystem showed the following characteristics: From 1996 until today, the number of passengers has been in constant decline, the number of passenger trains has decreased in domestic traffic by 28 trains, in inter-entity transport by 18 trains, and by 8 trains in international transport. The railway fee for the infrastructure in domestic services amounts to 4,176,295 KM, in inter-entity transport amounts to 601,836 KM, and in international transport amounts to 1,603,022 KM. The other elements of the RRS business operations are shown in [36]. The problem-solving methodology is based on: 1. increase in revenue from direct ticket sales and increase in revenue from agreements on PSO, 2. reduction in expensing, i.e., operating costs.
The expected result of the mentioned activities according to the given methodology should be sustainable business operations. Figure 4 shows the cost realization plan for the next fiscal year, as well as the perception of cost coverage by government revenues and government subsidies (PSOs), and a model of the long-term business stabilization goal ( Figure 5). The aim of the problem solving is to bring RRS as the operator in passenger traffic into the domain of positive operations, thus providing the necessary conditions for successful operation at a liberalized transport market. In [6], details about the model for liberalization in Serbia can be found.  The aim of the problem solving is to bring RRS as the operator in passenger traffic into the domain of positive operations, thus providing the necessary conditions for successful operation at a liberalized transport market. In [6], details about the model for liberalization in Serbia can be found.
In this paper, a new integrated F-PIPRECIA-F-EDAS model is created for solving problems. Multi-criteria methods for decision-making are used to resolve a large number of problems in all spheres of business, and they represent an area that is developing rapidly, primarily due to a large number of methods that have been developed, particularly within the last decade. The combination of these methods with fuzzy logic gives excellent results because classical methods cannot, with such precision, perform the required quantification, and this is where fuzzy logic shows all its advantages [37,38].

Possible Solutions
In order to resolve the problem, seven realistically possible variants (V) have been identified. All variants are described in [36]: V1: Reduction in operating costs ( Figure 6 left): where T denotes costs.
Symmetry 2020, 12, x FOR PEER REVIEW 8 of 20 In this paper, a new integrated F-PIPRECIA-F-EDAS model is created for solving problems. Multi-criteria methods for decision-making are used to resolve a large number of problems in all spheres of business, and they represent an area that is developing rapidly, primarily due to a large number of methods that have been developed, particularly within the last decade. The combination of these methods with fuzzy logic gives excellent results because classical methods cannot, with such precision, perform the required quantification, and this is where fuzzy logic shows all its advantages [37,38].

Possible Solutions
In order to resolve the problem, seven realistically possible variants (V) have been identified. All variants are described in [36]: o V1: Reduction in operating costs ( Figure 6 left): where T denotes costs.

Identification of Evaluation Criteria
Identification and quantification of the criteria for evaluating the manner of implementing the principles and concluding a PSO contract were carried out over four steps: Defining the required level or volume of service, reduction in business costs, increase in revenues from ticket sales, and increase in revenues under the PSO contract.
The selection of an optimal variant depends on many factors. Therefore, in the proposed methodology, as the criteria, the following values have been adopted: The reality of the feasibility of

Identification of Evaluation Criteria
Identification and quantification of the criteria for evaluating the manner of implementing the principles and concluding a PSO contract were carried out over four steps: Defining the required level or volume of service, reduction in business costs, increase in revenues from ticket sales, and increase in revenues under the PSO contract.
The selection of an optimal variant depends on many factors. Therefore, in the proposed methodology, as the criteria, the following values have been adopted: The reality of the feasibility of the proposed variant, means available to the public authority -budget, the ability of the operator, the effect of realization, and the period of realization (Table 1). Table 1. Criteria for evaluation of identified variants.

Mark
Title of Criterion Type C 1 variant feasibility profit-type C 2 public authorities-the budget size profit-type C 3 operators' ability profit-type C 4 the effect of realization profit-type C 5 Period of realization loss-type

Results
As for obtaining the weight value of criteria, we have used the F-PIPRECIA method; after identifying the criteria on which the ranking of potential variants will be made, it is required that we compare the criteria by using the scale presented in [17]. In order to derive the relative importance of the criteria, a team of three experts had been established; for many years, they have been performing managerial functions in the field of railway transport. As this is an already exploited method, detailed procedures for calculating the values of criteria will not be shown, but rather the summed results by each step (Table 2). Where s j represents the group matrix obtained by expert's assessment, starting from the second criterion, and k j is the coefficient obtained when s j is subtracted from number 2, except for s 1 . q j is the fuzzy weight, w j is the relative weight of the criterion, and DF is the defuzzified value. Based on the aggregation of the values wj shown in Table 2, the final criterion values are obtained:w 1 = 0.423; w 2 = 0.230; w 3 = 0.141; w 4 = 0.242; w 5 = 0.098. After calculating the weight value of criteria, we then begin the selection of the optimal variant by using the F-EDAS method. On the basis of the linguistic scale, the experts evaluate variants according to each criterion individually (Table 3).   Next, we need to calculate positive distances (PDA) and negative distances (NDA) from the average solutions depending on the criteria type. In this case, only the fifth criterion is useless, while the others are useful criteria. First, we obtain the values of the positive distance (PDA) and the values of the negative distance from the average solution. In order to obtain the values shown in Table 5, it is necessary to first apply step 5 of the F-EDAS method, and this represents the sum of the weighted matrix for positive sp i and negative distance sn i for all variants. Further, it is necessary to normalize previous values in order to obtain nsp i and nsn i . Finally, it is necessary to calculate the assessment of the results, the appraisal score ( as i ), and make the defuzzification appraisal score ( as i ) ( Table 5). Based on the performed analysis, and in accordance with the task, implementation of Variant A 7 is recommended as the most acceptable solution. As good enough solutions, we might accept variants A 6 and A 5 ; Variant A 3 could possibly represent a satisfactory solution. Thus, it is evident that the most acceptable variant is essentially the scenario in which the positive result stems from joint "efforts" of the operator (decreased costs and increased revenues from the ticket sales) and public authorities through increased subsidies for PSO. Another acceptable variant is a scenario where, because of the limitations of the market (low flow and low purchasing power of the population-passengers), there Symmetry 2020, 12, 743 13 of 20 lacks any significant increase in revenue from ticket sales; the solution is then sought through reduction in costs and increase in PSO subsidies. The variants where the problem is solved only by increased PSO subsidies by the public authorities and the combined approach based on the increase in revenues from ticket sales and operator's cost reduction are not favorable.

Changing the Significance of Criteria
In this phase of validation test, the impact of changing the three most important criteria C 1 , C 2 , and C 4 on the ranking results was analyzed. Using Equation (8), a total of 18 scenarios were formed.
In scenarios S 1 -S 6 , the first criterion was changed, criterion C 2 was changed in scenarios S 7 −S 12 , and criterion C 4 was changed in scenarios S 13 -S 18 . In Equation (8), W nβ represents the new value of criteria C 2 -C 5 for scenarios S 1 −S 6 ; then, C 1 , C 3 -C 5 for scenarios S 7 -S 12 , i.e., C 1 -C 3 , and C 5 for scenarios S 13 -S 18 . W nα represents the corrected value of criteria C 1 , C 2 , and C 3 respectively by groups of scenarios, W β represents the original value of the criterion considered, and W n represents the original value of the criterion whose value is reduced, in this case, C 1 , C 2 , and C 4 .
In all scenarios, the value of criteria was reduced by 15%, while the values of the remaining criteria were proportionally corrected by applying Equation (8). After forming 18 new vectors of the weight coefficients of the criteria (Table 6), new model results were obtained, as presented in Figure 10. In most scenarios, there is no change in initial rank, as shown in Figure 10. However, it is important to emphasize that the model is very sensitive to the change in the most important criterion, and in scenarios S1-S6, significant changes occur. With a slight decrease in the value of the first criterion, the ranks slightly change; for example, variants V 1 and V 4 change their positions in the second scenario. As the value of the first criterion decreases drastically, the ranks also change drastically. In the fourth scenario, V 7 loses the first position, while in the sixth scenario, it comes in last place. Practically, the most important role is played by the first criterion in the set decision conditions. In accordance with the rank changes in the mentioned scenarios, a statistical check of the rank correlation was performed using Spearman's correlation coefficient, as shown in Figure 11.  Figure 10. Comparison of obtained results by F-PIPRECIA-Fuzzy Evaluation based on Distance from Average Solution (F-EDAS) model with all formed scenarios S1-S18.
In most scenarios, there is no change in initial rank, as shown in Figure 10. However, it is important to emphasize that the model is very sensitive to the change in the most important criterion, and in scenarios S1-S6, significant changes occur. With a slight decrease in the value of the first criterion, the ranks slightly change; for example, variants V1 and V4 change their positions in the second scenario. As the value of the first criterion decreases drastically, the ranks also change drastically. In the fourth scenario, V7 loses the first position, while in the sixth scenario, it comes in last place. Practically, the most important role is played by the first criterion in the set decision conditions. In accordance with the rank changes in the mentioned scenarios, a statistical check of the rank correlation was performed using Spearman's correlation coefficient, as shown in Figure 11. The calculated Spearman's correlation coefficient (Figure 11), despite significant deviations in some scenarios, shows a high correlation of ranks in total, 0.821. Generally, in 13 out of 18 scenarios, variants have a full correlation. The correlation between the initial results obtained by the F-PIPRECIA-F-EDAS model and the S2 and S3 scenarios is 0.964, while in the S4 scenario, it is 0.786. The biggest deviation in the rankings is in the fifth and sixth scenarios when the negative correlations are −0.071 and -0.857, respectively.

Impact of Reverse Rank Matrices
One of the ways to test the validity of the obtained results is to construct dynamic matrices that analyze the solutions that the model provides under new conditions. A change in the number of variants is made for each scenario, eliminating the worst variant from further consideration. In the test, six scenarios are formed in which the change in elements of the decision matrix is simulated.
As can be seen in Figure 12, there is no change in ranks for any variant. That means that the proposed F-PIPRECIA-F-EDAS model is stable and gives good results. The calculated Spearman's correlation coefficient (Figure 11), despite significant deviations in some scenarios, shows a high correlation of ranks in total, 0.821. Generally, in 13 out of 18 scenarios, variants have a full correlation. The correlation between the initial results obtained by the F-PIPRECIA-F-EDAS model and the S 2 and S 3 scenarios is 0.964, while in the S 4 scenario, it is 0.786. The biggest deviation in the rankings is in the fifth and sixth scenarios when the negative correlations are −0.071 and -0.857, respectively.

Impact of Reverse Rank Matrices
One of the ways to test the validity of the obtained results is to construct dynamic matrices that analyze the solutions that the model provides under new conditions. A change in the number of variants is made for each scenario, eliminating the worst variant from further consideration. In the test, six scenarios are formed in which the change in elements of the decision matrix is simulated.
As can be seen in Figure 12, there is no change in ranks for any variant. That means that the proposed F-PIPRECIA-F-EDAS model is stable and gives good results.

Comparison with Other fuzzy MCDM Methods
In this part, a validation test is performed, including comparison with three other fuzzy methods: F-MARCOS, F-SAW, and the F-TOPSIS method. Obtained results are presented in Figure  13. As can be seen in Figure 13, there is no change in ranks for any variant. In Figure 13, in addition to the rankings of variants, values for each variant are given so that a cross-sectional comparison can be made.

Determining Criteria Weights with F-AHP and FUCOM Methods
In this part of the paper, the criteria weights were re-determined using the F-AHP and FUCOM methods, and the results compared to the original F-PIPRECIA-F-EDAS model are shown in Figure  14.

Comparison with other Fuzzy MCDM Methods
In this part, a validation test is performed, including comparison with three other fuzzy methods: F-MARCOS, F-SAW, and the F-TOPSIS method. Obtained results are presented in Figure 13.

Comparison with Other fuzzy MCDM Methods
In this part, a validation test is performed, including comparison with three other fuzzy methods: F-MARCOS, F-SAW, and the F-TOPSIS method. Obtained results are presented in Figure  13. As can be seen in Figure 13, there is no change in ranks for any variant. In Figure 13, in addition to the rankings of variants, values for each variant are given so that a cross-sectional comparison can be made.

Determining Criteria Weights with F-AHP and FUCOM Methods
In this part of the paper, the criteria weights were re-determined using the F-AHP and FUCOM methods, and the results compared to the original F-PIPRECIA-F-EDAS model are shown in Figure  14. As can be seen in Figure 13, there is no change in ranks for any variant. In Figure 13, in addition to the rankings of variants, values for each variant are given so that a cross-sectional comparison can be made.

Determining Criteria Weights with F-AHP and FUCOM Methods
In this part of the paper, the criteria weights were re-determined using the F-AHP and FUCOM methods, and the results compared to the original F-PIPRECIA-F-EDAS model are shown in Figure 14. Applying the above methods for determining the significance of the criteria and including them into the F-EDAS method yield the results presented in Figure 14. In addition to the ranks shown on the left, values of variants on the right are defined. It can be observed that F-PIPRECIA and FUCOM give identical ranks, while, applying F-AHP, there are deviations in the ranks of the first and fourth variants.

Additional Correction of Criteria Weights Obtained Using F-AHP
After presenting the previous results, the stability of the model is additionally determined as changing the significance of particular criteria. Therefore, a sensitivity analysis has been performed, which is presented throughout two parts in this subsection. Figure 15 shows the ranking of variants in all ten scenarios, while Figure 16 shows Spearman's coefficient of correlation for the ranking of variants. In the first set, the three most important criteria reduced the values by 10%, while the others increased by 15%. In the second set, the two most important criteria reduced the values by 15%, while the others increased by 10%. In the third set, the first criterion reduced by 20%, while the others increased by 5%. In next set, the second criterion reduced by 20%, while the others increased by 5%. In the fifth set, the fourth criterion reduced by 20%, while the others increased by 5%. In the next set,  Applying the above methods for determining the significance of the criteria and including them into the F-EDAS method yield the results presented in Figure 14. In addition to the ranks shown on the left, values of variants on the right are defined. It can be observed that F-PIPRECIA and FUCOM give identical ranks, while, applying F-AHP, there are deviations in the ranks of the first and fourth variants.

Additional Correction of Criteria Weights Obtained Using F-AHP
After presenting the previous results, the stability of the model is additionally determined as changing the significance of particular criteria. Therefore, a sensitivity analysis has been performed, which is presented throughout two parts in this subsection. Figure 15 shows the ranking of variants in all ten scenarios, while Figure 16 shows Spearman's coefficient of correlation for the ranking of variants. In the first set, the three most important criteria reduced the values by 10%, while the others increased by 15%. In the second set, the two most important criteria reduced the values by 15%, while the others increased by 10%. In the third set, the first criterion reduced by 20%, while the others increased by 5%. In next set, the second criterion reduced by 20%, while the others increased by 5%. In the fifth set, the fourth criterion reduced by 20%, while the others increased by 5%. In the next set, the first three criteria have values of 0.25, the fourth has 0.15, and the fifth has 0.1. In the seventh set, the criteria have values as follows: C 1 = C 2 = C 4 = 0.30, C 3 = 0.10, and the last criterion has a value of zero. In set 8: C 1 = C 2 = C 4 = 0.30, C 5 = 0.10, and the third criterion has a value of zero. In set 9: C 1 = 0.34, C 2 = 0.27, C 3 = 0.20, C 4 = 0.13, C 5 = 0.06. In the last set, C 1 = 0.30, C 2 = 0.20, C 3 = 0.15, C 4 = 0.20, C 5 = 0.15.
As it can be seen in Figure 15, the seventh variant in seven, from ten formed sets, represents the best solution, while in the other scenarios, the best solution is variant six. The fifth variant is stable in all formed scenarios and has a third position. Variant three and two are also very stable and, only in the first and sixth sets, changing the position. Variant three has position five in the first set, while variant two has position six in the first and sixth sets. The ranking of the first variant varies from the fourth to seventh position in different scenarios, while the fourth variant varies from the fifth to seventh position. We can conclude that with the decrease, the three most important criteria by the 10% results and ranking of variants are very sensitive. As it can be seen in Figure 15, the seventh variant in seven, from ten formed sets, represents the best solution, while in the other scenarios, the best solution is variant six. The fifth variant is stable in all formed scenarios and has a third position. Variant three and two are also very stable and, only in the first and sixth sets, changing the position. Variant three has position five in the first set, while variant two has position six in the first and sixth sets. The ranking of the first variant varies from the fourth to seventh position in different scenarios, while the fourth variant varies from the fifth to seventh position. We can conclude that with the decrease, the three most important criteria by the 10% results and ranking of variants are very sensitive. Figure 16 shows the SCC throughout all scenarios. From Figure 15, it can be seen that the model is sensitive to changes in the weight of the criteria and that each criterion can play an important role in the variant ranking. Spearman's coefficient of correlation has the range of 0.786-1.00, which represents a high degree of correlation, and the results obtained using the integrated fuzzy model are considered stable. The average SCC value for all ten formed scenarios in relation to the initial rank is 0.948. Figure 16. SCC through ten formed scenarios of F-AHP-F-EDAS model.

Discussion and Conclusion
In certain cases, there is risk of insufficient financial resources for the execution of the PSO. The costs for the realization of PSO by the operators that are in a state or local ownership may also affect the possibility of implementing the model. Within the framework of the realization of this model,  As it can be seen in Figure 15, the seventh variant in seven, from ten formed sets, represents the best solution, while in the other scenarios, the best solution is variant six. The fifth variant is stable in all formed scenarios and has a third position. Variant three and two are also very stable and, only in the first and sixth sets, changing the position. Variant three has position five in the first set, while variant two has position six in the first and sixth sets. The ranking of the first variant varies from the fourth to seventh position in different scenarios, while the fourth variant varies from the fifth to seventh position. We can conclude that with the decrease, the three most important criteria by the 10% results and ranking of variants are very sensitive. Figure 16 shows the SCC throughout all scenarios. From Figure 15, it can be seen that the model is sensitive to changes in the weight of the criteria and that each criterion can play an important role in the variant ranking. Spearman's coefficient of correlation has the range of 0.786-1.00, which represents a high degree of correlation, and the results obtained using the integrated fuzzy model are considered stable. The average SCC value for all ten formed scenarios in relation to the initial rank is 0.948. Figure 16. SCC through ten formed scenarios of F-AHP-F-EDAS model.

Discussion and Conclusion
In certain cases, there is risk of insufficient financial resources for the execution of the PSO. The costs for the realization of PSO by the operators that are in a state or local ownership may also affect the possibility of implementing the model. Within the framework of the realization of this model, Figure 16. SCC through ten formed scenarios of F-AHP-F-EDAS model. Figure 16 shows the SCC throughout all scenarios. From Figure 15, it can be seen that the model is sensitive to changes in the weight of the criteria and that each criterion can play an important role in the variant ranking. Spearman's coefficient of correlation has the range of 0.786-1.00, which represents a high degree of correlation, and the results obtained using the integrated fuzzy model are considered stable. The average SCC value for all ten formed scenarios in relation to the initial rank is 0.948.

Discussion and Conclusions
In certain cases, there is risk of insufficient financial resources for the execution of the PSO. The costs for the realization of PSO by the operators that are in a state or local ownership may also affect the possibility of implementing the model. Within the framework of the realization of this model, there are several possible sensitive situations that can appear from the moment of planning to the realization: Poor implementation of "business cost reduction" activities, especially with operators owned by government authorities; the lack of interest in "increasing revenue" in the gross contract, especially with operators owned by the authorities, regardless of whether they are revenue from the sale of tickets or other effects; incomplete and untimely realization of the fee for the execution of the PSO; lack of sufficient financial resources from the authorities to increase the fee for the execution of the PSO; the weakness of the state operators in the realization of other effects that can be realized on the basis of the granted right to perform PSO.
In this paper, a dynamic model for optimal application of the PSO system in the PPT process is proposed, which can contribute to the development of appropriate systems for the implementation of PPT services. In addition, it contributes to raising the service quality with the achievement of minimal costs of the functioning of these systems from the aspect of state and local government. By applying this model, it is possible to achieve a large number of effects (increase in passenger transport volume, higher and more stable quality of transport services, reduction in travel costs, better and more efficient cost control, etc.) and achieve significant savings in the functioning of the PPT system. Optimization of the PPT system has an indirect influence on the optimization of transport capacities and improvement in the quality of the transport service with economic quantification and cost savings.
The model was tested in the case of the organization of passenger traffic in the RRS (B&H). Based on the performed analysis, and in accordance with the task, the implementation of Variant V 7 is recommended as the most acceptable solution. As good enough solutions, we might accept variants V 6 and V 5 . The contribution of this research represents the possibility for rationalization of the PTT system in RRS. The new F-PIPRECIA-F-EDAS model developed in this research uses the strengths of fuzzy logic and multicriteria decision-making methods. One of the reasons for the F-PIPRECIA method application is its ability to equally handle quantitative and qualitative criteria. One of the reasons for using the F-EDAS method is a mathematical apparatus that assumes the evaluation of variants on the basis of positive and negative deviations from the average solution. The development of the new F-PIPRECIA-F-EDAS model based on TFNs represents the main scientific novelty of this paper. Future research related to this paper should be the implementation of the best variant and post-analysis of PPT systems.