Evaluation of Criteria in Fruit Production Using the Interval Fuzzy Rough PIPRECIA Method
Abstract:
This study investigates the application of Multi-Criteria Decision-Making (MCDM) techniques in fruit production, specifically focusing on the use of the interval fuzzy rough pivot pairwise relative criteria importance assessment (PIPRECIA) method for criteria evaluation. A total of 11 criteria were evaluated to rank various combinations of plum varieties and rootstocks. The criteria selected represent key aspects of plum production, including phenology, yield, physical fruit characteristics, and the chemical composition and quality of the fruit. Data for the study were collected through surveys of 17 experts and plum producers. The results indicated that the criteria related to overall yield and fruit weight were deemed the most significant, while those concerning the chemical composition and fruit quality were considered of lesser importance. The findings highlight the potential of the interval fuzzy rough PIPRECIA method in addressing both research and managerial challenges in fruit production. It is suggested that future research expand the application of this method to other geographical regions and agricultural sectors. Additionally, the development of accessible software tools featuring user-friendly interfaces could facilitate broader adoption of MCDM techniques in agricultural decision-making.1. Introduction
Among fruit species, plums rank ninth globally in terms of production and fourth among temperate fruits, following apples, pears, and peaches. China is the leading producer of plums worldwide, with a production of 6.8 million tons in 2023, followed by Romania (645,090 tons) and Serbia (362,713 tons). The average global plum production between 2021 and 2023 was 12.5 million tons [1]. The widespread cultivation of plums can be attributed to the versatile use of its fruits, modest growth requirements, high yields, and excellent fruit quality. The variety and rootstock are crucial factors for successful fruit cultivation, including plum production. However, the choice of rootstock is equally significant, as it affects not only vigor and yield [2], [3] but also fruit quality [4], [5]. Greater attention must be paid to rootstock selection compared to the variety, as no universal rootstock exists that is suitable for varying agroecological conditions, cultivation systems, and intensive production [6]. Successful plum production in high-density orchards with high yields and quality fruit depends on both the variety and the rootstock [7].
A well-balanced ratio between generative and vegetative buds is essential for balancing yield and fruit quality. Research over the past two decades highlights the influence of vegetative rootstocks in reducing tree vigor and increasing plum yield per unit area [8], [9]. Through interaction with the grafted rootstock, the variety can significantly influence rootstock performance regarding tree vigor [10]. Fertility is one of the most important breeding goals for European plum varieties [11].
In fruit research, comparisons of varieties and/or variety/rootstock combinations, as well as rootstocks, often involve numerous alternatives—sometimes exceeding ten. Additionally, numerous characteristics (criteria) have been used, which in some studies exceed 20. Ali et al. [12] conducted research on 13 apple varieties with 21 traits each. Some alternatives excel in certain traits, while others perform better in different ones. Objective ranking of these varieties is not feasible without using MCDM techniques. Comparative research with five plum varieties and ten traits was conducted by Milošević et al. [13]. The study found that varieties exhibited different values across various criteria. For instance, the "Katinka" variety had the highest yield per tree, while the "Boranka" and "Čačanska Rana" varieties had the largest fruit size and weight.
Radović et al. [14] focused on the comparative influence of different rootstocks on the chemical composition of plum fruits. The study involved four varieties and three rootstocks, yielding 12 combinations (alternatives). The content of phenolic compounds, sugars, and mineral composition was assessed using ten individual traits (criteria). Significant variation was observed among different variety/rootstock combinations regarding values across individual criteria. These examples illustrate the complexity of ranking and selecting the best alternative (variety or variety/rootstock combination) in fruit production, highlighting the need for MCDM techniques to achieve greater objectivity.
Recently, some studies have incorporated MCDM techniques for ranking alternatives in decision-making processes in fruit production. The fuzzy compromise ranking of alternatives from distance to ideal solution (CRADIS) and Criteria Importance Through Intercriteria Correlation (CRITIC) methods were used for market evaluation of six pear varieties based on ten criteria [15]. In the process of ranking and selecting plum varieties for establishing new orchards, Nedeljković et al. [16] employed several MCDM techniques. They ranked six plum varieties using five criteria and 20 sub-criteria. The weights of the selected criteria were determined using the Full Consistency Method (FUCOM) and the CRITIC, while the fuzzy Measurement Alternatives and Ranking according to Compromise Solution (MARCOS) method was used for ranking alternatives.
2. Methodology
An interval fuzzy rough number (IFRN) was marked with "$A$".
where, $a_{j q}^L=\underline{\operatorname{Lim}}\left(I *\left(a_j\right)_{l q}\right)$ and $a_{j q}^U=\overline{\operatorname{Lim}}\left(I *\left(a_j\right)_{u q}\right)$; $(j=1,2,3 ; 1 \leq q \leq k)$.
The expressions for addition, subtraction, multiplication and division are as follows:
Addition:
Subtraction:
Multiplication:
Division:
Fuzzy PIPRECIA method was created by Stević et al. [17] and has been applied many times in various studies [18], [19], [20]. IFRN PIPRECIA version was developed in 2024 by Chen et al. [21]. The procedure of the IFR PIPRECIA is shown below.
Step 1: Defining the set of $n$ criteria.
Step 2: Forming the group of $t$ experts who must use two scales shown in Figure 1 and Figure 2.
Step 3: Linguistic scale 1-2 is applied by the experts if the $C_j$ criterion has higher importance compared to $C_{j-1}$. Abbreviations from Figure 1 and Figure 2 were explained by Chen et al. [21].

A linguistic scale of 0-1 should be applied when $C_j$ has lesser importance compared to $C_{j-1}$.

Step 4: In this step, Eq. (6) should be applied for each expert.
where, Fuzzy$\left[s_j^r\right]$ is the assessment of criteria by each expert $r$.
Step 5: Defining the group IFRN matrix $x_j$.
Step 6: Computation of the coefficient $\text{IFRN}\left[k_j\right]$.
Step 7: Computation of $\text{IFRN}\left[q_j\right]$.
Step 8: Computation of $\text{IFRN}\left[w_j\right]$.
Step 9: In this step, the experts evaluate the criteria again, starting from the penultimate one.
where, Fuzzy$\left[s_j^{r^{\prime}}\right]$ is the assessment of criteria by each expert $r^{\prime}$.
Step 10: Computation of the coefficient $\text{IFRN}\left[k_j{ }^{\prime}\right]$.
Step 11: Computation of $\text{IFRN}\left[q_j{ }^{\prime}\right]$.
Step 12: Computation of $\text{IFRN}\left[w_j^{\prime}\right]$.
Step 13: Final weights can be obtained:
3. Case Study
This study is part of the project related to determining criteria weights and creating a web application in fruit production to support producers and decision-makers. Totally, 17 experts were involved in this project from various structures. For the purpose of this study, independent producers were selected and the calculation process of criteria weights based on the preference of three experts from this field was shown.
The IFRN PIPRECIA method was used for calculating criteria weights. The following 11 criteria integrating plum cultivars/rootstocks were observed: C1 - the number of flower buds per 1 m of fruiting twig; C2 - the number of flower buds per 1 m of 2-year-old twig; C3 - germination of pollen; C4 - final fruit set; C5 - trunk cross-sectional area; C6 - cumulative yield efficiency; C7 - fruit weight; C8 - flesh ratio; C9 - soluble solids; C10 - total phenolic content in skin of fruit; C11 - radical-scavenging activity. The first activity in applying this method was assessing the mutual significance of criteria by experts using the linguistic scales, as shown in Table 1.
DM1 | DM2 | DM3 | ||||
IFRNP | IIFRNP | IFRNP | IIFRNP | IFRNP | IIFRNP | |
C1 | - | - | - | - | - | - |
C2 | AE | WL | SM | MLS | AE | WL |
C3 | AE | WL | AE | WL | SM | MLS |
C4 | MMS | L | SM | WL | SM | MLS |
C5 | ML | MMS | L | MMS | RL | M |
C6 | M | AL | MMS | RL | M | RL |
C7 | AE | WL | WL | AE | AE | WL |
C8 | ML | M | WL | AE | AE | WL |
C9 | AE | MLS | MLS | SM | L | MMS |
C10 | DL | M | WL | AE | WL | AE |
C11 | WL | AE | WL | AE | WL | AE |
The next step was transforming linguistic values into Triangular Fuzzy Numbers (TFNs), as shown in Table 2.
The procedure of converting TFNs into interval fuzzy rough is as follows:
A rough matrix for C5 is shown below.
Based on the expert assessment, as shown in Table 2, three classes of objects $l, m$ and $u$ were selected: $l=(0.29;0 .40 ; 0.33), m=(0.33 ; 0.50 ; 0.40)$ and $u=(0.40 ; 0.67 ; 0.50)$.
\[\begin{aligned} l: \quad \quad & \underline{\operatorname{Lim}}(0.29)=0.29, \overline{\operatorname{Lim}}(0.29)=\frac{1}{3}(0.29+0.40+0.33)=0.34 \\ & \underline{\operatorname{Lim}}(0.40)=\frac{1}{3}(0.29+0.40+0.33)=0.34, \overline{\operatorname{Lim}}(0.40)=0.40 \\ & \underline{\operatorname{Lim}}(0.33)=\frac{1}{2}(0.29+0.33)=0.31, \overline{\operatorname{Lim}}(0.33)=\frac{1}{2}(0.40+0.33)=0.37\end{aligned}\]
\[\begin{aligned} m: \quad \quad & \underline{\operatorname{Lim}}(0.33)=0.33, \overline{\operatorname{Lim}}(0.33)=\frac{1}{3}(0.33+0.50+0.40)=0.41 \\ & \underline{\operatorname{Lim}}(0.50)=\frac{1}{3}(0.33+0.50+0.40)=0.41, \overline{\operatorname{Lim}}(0.50)=0.50 \\ & \underline{\operatorname{Lim}}(0.40)=\frac{1}{2}(0.33+0.40)=0.37, \overline{\operatorname{Lim}}(0.40)=\frac{1}{2}(0.50+0.40)=0.45\end{aligned}\]
\[\begin{aligned} u: \quad \quad & \underline{\operatorname{Lim}}(0.40)=0.40, \overline{\operatorname{Lim}}(0.40)=\frac{1}{3}(0.40+0.67+0.50)=0.52 \\& \underline{\operatorname{Lim}}(0.67)=\frac{1}{3}(0.40+0.67+0.50)=0.52, \overline{\operatorname{Lim}}(0.67)=0.67 \\& \underline{\operatorname{Lim}}(0.50)=\frac{1}{2}(0.40+0.50)=0.45, \overline{\operatorname{Lim}}(0.50)=\frac{1}{2}(0.67+0.50)=0.58\end{aligned}\]
DM1 | DM2 | DM3 | ||||
IFRNP | IIFRNP | IFRNP | IIFRNP | IFRNP | IIFRNP | |
C1 | - | - | - | - | - | - |
C2 | (1,1,1.05) | (0.67,1,1) | (1.1,1.15,1.2) | (0.5,0.67,1) | (1,1,1.05) | (0.67,1,1) |
C3 | (1,1,1.05) | (0.67,1,1) | (1,1,1.05) | (0.67,1,1) | (1.1,1.15,1.2) | (0.5,0.67,1) |
C4 | (1.2,1.3,1.35) | (0.4,0.5,0.67) | (1.1,1.15,1.2) | (0.67,1,1) | (1.1,1.15,1.2) | (0.5,0.67,1) |
C5 | (0.29,0.33,0.4) | (1.2,1.3,1.35) | (0.4,0.5,0.67) | (1.2,1.3,1.35) | (0.33,0.4,0.5) | (1.3,1.45,1.5) |
C6 | (1.3,1.45,1.5) | (0.22,0.25,0.29) | (1.2,1.3,1.35) | (0.33,0.4,0.5) | (1.3,1.45,1.5) | (0.33,0.4,0.5) |
C7 | (1,1,1.05) | (0.67,1,1) | (0.67,1,1) | (1,1,1.05) | (1,1,1.05) | (0.67,1,1) |
C8 | (0.29,0.33,0.4) | (1.3,1.45,1.5) | (0.67,1,1) | (1,1,1.05) | (1,1,1.05) | (0.67,1,1) |
C9 | (1,1,1.05) | (0.5,0.67,1) | (0.5,0.67,1) | (1.1,1.15,1.2) | (0.4,0.5,0.67) | (1.2,1.3,1.35) |
C10 | (0.25,0.29,0.33) | (1.3,1.45,1.5) | (0.67,1,1) | (1,1,1.05) | (0.67,1,1) | (1,1,1.05) |
C11 | (0.67,1,1) | (1,1,1.05) | (0.67,1,1) | (1,1,1.05) | (0.67,1,1) | (1,1,1.05) |
C12 | (1,1,1.05) | (0.67,1,1) | (1.1,1.15,1.2) | (0.5,0.67,1) | (1,1,1.05) | (0.67,1,1) |
In this way, IFRNs were computed as follows:
$ \begin{aligned} & \text{IFRN}\left(E_1\right)=[(0.29,0.34),(0.33,0.41),(0.40,0.52)] \\ & \text{IFRN}\left(E_2\right)=[(0.34,0.40),(0.41,0.50),(0.52,0.67)] \\ & \text{IFRN}\left(E_3\right)=[(0.31,0.37),(0.37,0.45),(0.45,0.58)] \end{aligned} $
Using the averaging equation, the final IFRN for C5 was calculated as follows:
$ \text{IFRN}\left(C_5\right)=[(0.31,0.37),(0.37,0.45),(0.46,0.59)] $
Then $\text{IFRN}\left(X_j\right)$ was calculated (Table 3). For the inverse IFRN PIPRECIA method, the same methodology was applied.
| IFRN PIPRECIA | Inverse IFRN PIPRECIA |
| IFRN (Xj) | |
C1 | [(0,0),(0,0),(0,0)] | [(0.573,0.648),(0.816,0.963),(1,1)] |
C2 | [(1.01,1.017),(1.017,1.067),(1.067,1.133)] | [(0.573,0.648),(0.816,0.963),(1,1)] |
C3 | [(1.01,1.017),(1.017,1.067),(1.067,1.133)] | [(0.457,0.582),(0.6,0.796),(0.816,0.963)] |
C4 | [(1.11,1.153),(1.167,1.217),(1.217,1.283)] | [(1.21,1.253),(1.317,1.367),(1.367,1.433)] |
C5 | [(0.312,0.368),(0.371,0.45),(0.457,0.589)] | [(0.274,0.306),(0.317,0.362),(0.382,0.477)] |
C6 | [(1.247,1.29),(1.367,1.417),(1.417,1.483)] | [(0.704,0.853),(1,1),(1.007,1.03)] |
C7 | [(0.816,0.963),(1,1),(1.02,1.043)] | [(0.829,1.047),(1.05,1.07),(1.07,1.32)] |
C8 | [(0.472,0.604),(0.604,0.64),(0.64,0.967)] | [(0.743,0.872),(0.872,1.093),(1.093,1.27)] |
C9 | [(0.493,0.6),(0.6,0.802),(0.802,0.997)] | [(1.033,1.05),(1.05,1.1),(1.1,1.3)] |
C10 | [(0.437,0.54),(0.602,0.631),(0.631,0.927)] | [(1,1),(1,1),(1.05,1.05)] |
C11 | [(0.667,0.667),(1,1),(1,1)] | [(0,0),(0,0),(0,0)] |
In the next step, $\text{IFRN}\left(k_j\right)$ was calculated.
The first part of the matrix $\text{IFRN}\left(k_j\right)$ is equal to IFRN of one, because that is the rule defined in the applied methodology.
$\begin{aligned}\text{IFRN}\left(k_j\right)=\left(\begin{array}{c}{[(1.000,1.000),(1.000,1.000),(1.000,1.000)]} \\ {[(0.867,0.933),(0.933,0.983),(0.983,0.990)]} \\ {[(0.867,0.933),(0.933,0.983),(0.983,0.990)]} \\ {[(0.717,0.783),(0.783,0.833),(0.847,0.890)]} \\ {[(1.411,1.543),(1.550,1.629),(1.632,1.688)]} \\ {[(0.517,0.583),(0.583,0.633),(0.710,0.753)]} \\ {[(0.957,0.980),(1.000,1.000),(1.037,1.184)]} \\ {[(1.033,1.360),(1.360,1.396),(1.396,1.528)]} \\ {[(1.003,1.198),(1.198,1.400),(1.400,1.507)]} \\ {[(1.073,1.369),(1.369,1.398),(1.460,1.563)]} \\ {[(1.000,1.000),(1.000,1.000),(1.333,1.333)]}\end{array}\right)\end{aligned}$
An example of calculating the elements of matrix $\text{IFRN}\left(k_j\right)$ is as follows:
$\begin{aligned}\text{IFRN}\left[k_2\right] & =\left[\begin{array}{l}\left(2-x_j^{U 3}\right),\left(2-x_j^{L 3}\right),\left(2-x_j^{U 2}\right),\left(2-x_j^{L 2}\right),\left(2-x_j^{U 1}\right),\left(2-x_j^{L 1}\right) \\ (2-1.333),(2-1.067),(2-1.067),(2-1.017),(2-1.017),(2-1.010)\end{array}\right] \\ & =[(0.867,0.933),(0.933,0.983),(0.983,0.990)]\end{aligned}$
The next activity represents the calculation of $\text{IFRN}\left[q_j\right]$ as follows:
$\begin{aligned}\text{IFRN}\left[q_j\right]=\left(\begin{array}{c}{[(1.000,1.000),(1.000,1.000),(1.000,1.000)]} \\{[(1.010,1.017),(1.017,1.071),(1.071,1.154)]} \\{[(1.020,1.034),(1.034,1.148),(1.148,1.331)]} \\{[(1.146,1.221),(1.241,1.465),(1.465,1.858)]} \\{[(0.679,0.748),(0.762,0.945),(0.950,1.316)]} \\{[(0.901,1.054),(1.203,1.621),(1.628,2.548)]} \\{[(0.761,1.017),(1.203,1.621),(1.661,2.663)]} \\{[(0.498,0.729),(0.862,1.192),(1.221,2.578)]} \\{[(0.331,0.520),(0.616,0.995),(1.020,2.569)]} \\{[(0.211,0.356),(0.440,0.727),(0.745,2.393)]} \\{[(0.159,0.267),(0.440,0.727),(0.745,2.393)]}\end{array}\right)\end{aligned}$
Again, the first element is equal to IRFN of one. In order to calculate the shown matrix, the following procedure should be applied:
$\begin{aligned}&\text{IFRN}\left[q_3\right]=\left[\left(\frac{q_2^{L 1}}{k_3^{U 3}}\right),\left(\frac{q_2^{U 1}}{k_3^{L 3}}\right),\left(\frac{q_2^{L 2}}{k_3^{U 2}}\right),\left(\frac{q_2^{U 2}}{k_3^{L 2}}\right),\left(\frac{q_2^{L 3}}{k_2^{U 1}}\right),\left(\frac{q_2^{U 3}}{k_3^{L 1}}\right)\right] \\& \text{IFRN}\left[q_3\right]=\left[\left(\frac{1.010}{0.990}\right),\left(\frac{1.017}{0.983}\right),\left(\frac{1.017}{0.983}\right),\left(\frac{1.071}{0.933}\right),\left(\frac{1.017}{0.933}\right),\left(\frac{1.154}{0.867}\right)\right] \\& \text{IFRN}\left[q_3\right]=[(1.020,1.034),(1.034,1.148),(1.148,1.331)]\end{aligned}$
The matrix $\text{IFRN}\left[w_j\right]$ was obtained after summing the previous matrix and calculating
\[\left[\delta_j\right]=[(7.717,8.964),(9.818,12.512),(12.654,21.804)]:\]
$\begin{align}\text{IFRN}\left[w_j\right]=\left(\begin{array}{c}{[(0.046, 0.079), (0.080, 0.102), (0.112, 0.130)]} \\{[(0.046, 0.080), (0.081, 0.109), (0.120, 0.150)]} \\{[(0.047, 0.082), (0.083, 0.117), (0.128, 0.173)]} \\{[(0.053, 0.097), (0.099, 0.149), (0.163, 0.241)]} \\{[(0.031, 0.059), (0.061, 0.096), (0.106, 0.171)]} \\{[(0.041, 0.083), (0.096, 0.165), (0.182, 0.330)]} \\{[(0.035, 0.080), (0.096, 0.165), (0.185, 0.345)]} \\{[(0.023, 0.058), (0.069, 0.121), (0.136, 0.334)]} \\{[(0.015, 0.041), (0.049, 0.101), (0.114, 0.333)]} \\{[(0.010, 0.028), (0.035, 0.074), (0.083, 0.310)]} \\{[(0.007, 0.021), (0.035, 0.074), (0.083, 0.310)]}\end{array}\right)\end{align}$
\[\text{IFRN}\left[w_1\right]=\left[\left(\frac{1.000}{21.804}\right),\left(\frac{1.000}{12.654}\right),\left(\frac{1.000}{12.512}\right),\left(\frac{1.000}{9.818}\right),\left(\frac{1.000}{8.964}\right),\left(\frac{1.000}{7.717}\right)\right]\]
Using equations for the inverse IFRN PIPRECIA, part of the methodology criteria weights is as follows:
$\begin{align*}\text{IFRN}\left[w_j{ }^{\prime}\right]=\left(\begin{array}{c}{[(0.005,0.020),(0.036,0.110),(0.142,0.547)]} \\{[(0.006,0.027),(0.043,0.114),(0.142,0.547)]}\\{[(0.009,0.037),(0.050,0.118),(0.142,0.547)]}\\{[(0.014,0.052),(0.071,0.142),(0.169,0.567)]}\\{[(0.011,0.039),(0.048,0.090),(0.107,0.321)]}\\{[(0.019,0.066),(0.081,0.147),(0.173,0.489)]}\\{[(0.025,0.075),(0.081,0.147),(0.172,0.475)]}\\{[(0.029,0.072),(0.077,0.137),(0.160,0.323)]}\\{[(0.037,0.081),(0.087,0.124),(0.145,0.236)]}\\{[(0.036,0.077),(0.083,0.112),(0.130,0.165)]}\\{[(0.036,0.077),(0.083,0.112),(0.124,0.157)]}\end{array}\right)\end{align*}$
The final criteria weights are shown in Table 4.
The obtained results show that the order of criteria is as follows: C7$>$C6$>$C4$>$C8$>$C3$>$C2$>$C1$>$C9$>$C10$>$C11.
$W_j$ | |
C1 | {[(0.025,0.050),(0.058,0.106),(0.127,0.338)]} |
C2 | {[(0.026,0.054),(0.062,0.111),(0.131,0.348)]} |
C3 | {[(0.028,0.059),(0.067,0.117),(0.135,0.360)]} |
C4 | {[(0.033,0.074),(0.085,0.146),(0.166,0.404)]} |
C5 | {[(0.021,0.049),(0.055,0.093),(0.106,0.246)]} |
C6 | {[(0.030,0.075),(0.089,0.156),(0.177,0.410)]} |
C7 | {[(0.030,0.078),(0.089,0.156),(0.178,0.410)]} |
C8 | {[(0.026,0.065),(0.073,0.129),(0.148,0.328)]} |
C9 | {[(0.026,0.061),(0.068,0.113),(0.129,0.284)]} |
C10 | {[(0.023,0.053),(0.059,0.093),(0.107,0.238)]} |
C11 | {[(0.022,0.049),(0.059,0.093),(0.103,0.233)]} |
4. Conclusions
The research was conducted on 11 different criteria relevant to the selection of plum variety/rootstock combinations, involving 17 domain experts and fruit producers. The interval fuzzy rough PIPRECIA method proved to be applicable for addressing complex decision-making problems in fruit production. Based on the independent evaluations of experts and the application of the IFRN PIPRECIA method, the results highlight the prioritization of criteria related to yield and physical fruit characteristics. Specifically, the criteria of fruit weight (C7), cumulative yield efficiency (C6), and final fruit set (C4) were determined as the most important. At the same time, criteria related to the chemical composition and fruit quality were identified as the least important (C9 - soluble solids; C10 - total phenolic content in skin of fruit; and C11 - radical-scavenging activity).
The experience gained from this research demonstrates the significant potential of using the interval fuzzy rough PIPRECIA method in addressing research and managerial problems in fruit production. Future studies should focus on applying this and similar methods to decision-making problems in other agricultural sectors and developing digital tools to facilitate the application of MCDM techniques for a broader range of users.
Conceptualization, M.R. and R.B.; methodology, Ž.S.; software, Ž.S.; validation, D.M., and A.G.L.; formal analysis, G.V.; investigation, M.R.; resources, M.R.; data curation, G.V.; writing—original draft preparation, Ž.S.; writing—review and editing, G.V.; visualization, Ž.S.; supervision, D.M.; project administration, R.B.; funding acquisition, M.R. All authors have read and agreed to the published version of the manuscript.
The data supporting our research results are included within the article or supplementary material.
The authors declare no conflict of interest.
