Enriched preference modeling and robustness analysis for the ELECTRE Tri-B method

We consider multiple criteria sorting problems with preference-ordered classes delimited by a set of boundary profiles. While significantly extending the ELECTRE Tri-B method, we present an integrated framework for modeling indirect preference information and conducting robustness analysis. We allow the Decision Maker (DM) to provide the following three types of holistic judgments: assignment examples, assignment-based pairwise comparisons, and desired class cardinalities. A diversity of recommendation that can be obtained given the plurality of outranking-based sorting models compatible with the DM’s preferences is quantified by means of six types of results. These include possible assignments, class acceptability indices, necessary assignment-based preference relation, assignment-based outranking indices, extreme class cardinalities, and class cardinality indices. We discuss the impact of preference information on the derived outcomes, the interrelations between the exact results computed with mathematical programming and stochastic indices estimated with the Monte Carlo simulations, and new measures for quantifying the robustness of results. The practical usefulness of the approach is illustrated on data from the Financial Times concerning MBA programs.


Introduction
Multiple criteria sorting problems involve an assignment of alternatives to pre-defined and ordered classes in the presence of numerous pertinent viewpoints. The definition of decision classes is related to how the alternatives placed in a given class should be processed or treated. For example, in portfolio decision analysis, the stocks can be assigned to three categories: attractive, to be studied further, and non-attractive . When evaluating software producing technology companies, the latter ones can be sorted into predefined efficiency classes (Bilich and da Silva 2008). Furthermore, the assessment of natural gas pipelines involves assigning different accident scenarios into risk categories (Brito et al. 2010). In nanotechnology, one can predict a precaution level while handling nanomaterials in certain conditions (Kadziński et al. 2020). As far as medical diagnosis is concerned, patients can be judged in terms of the risk levels of developing different illnesses (e.g., diabetic retinopathy Saleh et al. 2018). In the context of a comprehensive assessment of insulating materials, we can consider preference-ordered sustainability classes .
The field of Multiple Criteria Decision Aiding (MCDA) offers a plethora of methods that support Decision Makers (DMs) in providing recommendations regarding sorting problems. These approaches differ in terms of the underlying assumptions, contexts of application, required preference information, incorporated models, and applied sorting rules. Nonetheless, all of them expect the DM to interact in the specification of a consistent family of criteria, a set of alternatives and their performances, as well as a set of parameters that represent his/her preferences.
One of the most prevailing sorting methods is ELECTRE Tri-B (Yu 1992), which incorporates an outranking-based preference model (Figueira et al. 2005) and defines decision classes by means of limiting profiles. Such profiles are composed of performances that can be seen as natural boundaries between a pair of successive classes. The intuitiveness of ELECTRE Tri-B was found appealing in numerous real-world case studies in such various application areas as climate change (Diakoulaki and Hontou 2003), economy and finance (Dimitras et al. 1995), land-use suitability assessment (Joerin et al. 2001), skills accreditation system (Siskos et al. 2007), water resources management (Raju et al. 2000), or zoning risk analysis (Merad et al. 2004).
The basic variant of ELECTRE Tri-B has also been advanced in numerous ways with the aim of extending its applicability. In particular, it was extended to account for both a hierarchical structure of criteria and interaction effects (Corrente et al. 2016). Moreover, it was generalized to admit both a set of limiting class profiles for better characterization of the boundaries between categories (Fernandez et al. 2017) and interval data (Fernandez et al. 2019). Also, the requirement of defining class boundaries was alleviated by admitting the use of characteristic (i.e., the most typical) profiles (Almeida-  or defining the class boundaries implicitly by means of example assignment examples (Köksalan et al. 2009). Nevertheless, the greatest number of methodological advancements of ELECTRE Tri-B was related to the elicitation of parameters compatible with the DM's preferences. Since the DMs may have some difficulties in specifying directly a coherent set of such parameters (involving criteria weights, credibility level, comparison thresholds, and limiting profiles), a preference disaggregation approach was envisaged to deal with the DM's holistic preference information and indirect inference of parameters.
The pioneering preference disaggregation approach for ELECTRE Tri-B was provided in Mousseau and Słowiński (1998). The authors proposed an interactive non-linear optimization model to infer all parameters of the sorting method from the desired assignments for a subset of reference alternatives relatively well-known to the DM. Then, the linear counterpart of this model was introduced to derive only the criteria weights and credibility level (Mousseau et al. 2001). On the contrary, Mixed Integer Linear Programming was used in Ngo The and  to learn the category limits. Furthermore, the definition of the credibility of an outranking relation was revised in Dias and Mousseau (2006) for allowing easier indirect inference of veto-related parameters. A group decision preference disaggregation setting was considered in Damart et al. (2007) with the aim of preserving consistency of sorting examples at both individual and collective levels. Moreover, when considering assignment examples of multiple DMs, (Cailloux et al. 2012) discussed how to infer collective class boundaries. While all aforementioned works considered the pessimistic sorting rule of ELECTRE Tri-B, (Zheng et al. 2014) studied the efficiency of inference procedures in the context of an optimistic procedure. Finally, the algorithms for dealing with inconsistent holistic preference judgments were proposed in Mousseau et al. (2003) and Mousseau et al. (2006).
Apart from inferring a set of compatible parameters from the DM's assignment examples, an inherent part of the preference disaggregation paradigm consists in using such parameters through a sorting model to derive well-justified decisions. Traditionally, the analysis of the DM's holistic decisions resulted in constructing a single set of compatible parameter values, whose employment led to the univocal class assignments of decision alternatives. However, when using indirect and/or imprecise preference information, there usually exist infinitely many compatible sets of parameters. Although their application on the set of reference alternatives leads to reproducing the DM's desired assignments, the sorting recommendation derived for the non-reference alternatives which were not judged by the DM may differ from one set of parameters to another.
In this perspective, (Dias et al. 2002) combined a preference disaggregation approach with robustness analysis providing for each alternative a set of classes confirmed by at least one compatible set of parameters. In this case, the variables of proposed mathematical programming models were the criteria weights and credibility level. Furthermore, while admitting imprecision in the specification of admissible parameter values related to the formulation of sorting model (including the boundary profiles), (Tervonen et al. 2007) and (Tervonen et al. 2009) adapted Stochastic Multicriteria Acceptability Analysis (SMAA) to ELECTRE Tri-B. Specifically, the SMAA-TRI method incorporated the simulation techniques to estimate a proportion of feasible sets of parameters that assign each alternative to a given class.
This paper advances ELECTRE Tri-B by proposing an enriched framework for preference modeling and robustness analysis. When it comes to preference learning, we account for three types of indirect and imprecise preference information. They refer to a small part of a desired final recommendation or requirements to be imposed on the derived assignments. Specifically, apart from the assignment examples which have been traditionally incorporated into preference disaggregation approaches for ELECTRE Tri-B (see Mousseau and Słowiński 1998;Dias et al. 2002), we consider assignment-based pairwise comparisons (Kadziński et al. 2015) and desired class cardinalities . The former refer to the relative class assignments desired for pairs of alternatives without specifying any concrete classes (e.g., alternative a should be assigned to a class better than alternative b), whereas the latter specify the bounds on the number of alternatives that can be assigned to a particular class (e.g., at least 5 alternatives should be assigned to the most preferred class, and at most half of the alternatives can be sorted into the worst class). While pairwise comparisons are strictly related to the DM's preferences, the constraints on the number of alternatives that can be assigned to each class are rather implied by a particular decision aiding context. For all types of preference information, we discuss the mathematical programming models that contribute to the definition of all compatible sets of parameters. Accounting for such diverse forms of preference is beneficial in terms of offering greater flexibility to the DM, reducing the space of feasible parameter sets, and making the proposed framework suitable for various sorting contexts.
As far as robustness analysis is concerned, we incorporate a number of robust results, which are derived from the exploitation of a set of parameter sets compatible with the DM's preferences. The diversity of these outcomes comes from using different tools for assessing the robustness of the same part of the recommendation and from considering various per-spectives on the stability of recommendation for a given problem. Specifically, we jointly use linear programming-for constructing the necessary (certain), possible, and extreme outcomes in an exact manner-and the Monte Carlo simulation-for estimating the probability of results, which are possible, though not certain. Indeed, some recent case studies indicate that a joint consideration of the robust and stochastic results can be even more beneficial for real-world decision making, as these results offer complementary perspectives on the stability of recommendation (see Dias et al. 2018;. The constructed results differ in terms of an adopted perspective while representing multidimensional and interrelated conclusions (Kadziński and Ciomek 2016). For the individual alternatives, we consider possible assignments and class acceptability indices, for pairs of alternatives-necessary assignment-based preference relation  and assignment-based outranking indices, whereas for decision classes-extreme class cardinalities and class cardinality indices . Note that the so-far existing robustness analysis methods for ELECTRE Tri-B focussed on the stability of class assignments. Accounting for various types of results, including assignment-based preference relation and class cardinalities, is beneficial for understanding one's own preferences, informing decision making, enhancing the DM's learning on the problem, stimulating his/her reaction to interactively enrich preference information, and arriving at a more credible recommendation.
When considering different types of preference information, the space of compatible sets of parameters is non-convex. Thus, to derive the exact robust results, we propose some Mixed-Integer Linear Programming models, whereas to estimate the acceptability indices, we combine Hit-And-Run (HAR)  with the rejection sampling technique. The additional contribution of the paper comes from discussing an impact of the provided preference information on the truth or falsity of robust results as well as the interrelations between the exact and stochastic outcomes obtained with mathematical programming and the Monte Carlo simulations. The practical usefulness of the approach is illustrated on data from the Financial Times concerning MBA programs. Overall, the novelty of the paper derives from the following developments: -introducing original models based on mathematical programming for defining a set of outranking models compatible with three types of indirect and imprecise preference information concerning alternatives, pairwise comparisons, and decision classes; in particular, we propose two alternative models for representing assignment-based pairwise comparisons, which have not been used before in the context of ELECTRE Tri-B; -introducing novel procedures based on mathematical programming and the Monte Carlo simulation for exploiting a set of compatible outranking models and deriving six types of complementary results; the exact and stochastic outcomes concerning assignment-based preference relations and class cardinalities are considered for the first time in view of ELECTRE Tri-B; -implementing an original idea for combining Hit-And-Run and the rejection sampling techniques for efficient exploitation of a non-convex space of compatible outranking models; -comprehensive discussion of the theoretical properties of the results and interrelations between inputs and outputs which enhance interactive elicitation of preferences; -providing the richest framework for preference modeling and robustness analysis, which can be generalized to other sorting methods; -introducing novel measures for quantifying the robustness of sorting results.
The remainder of the paper is organized as follows. In Sect. 2, we remind the basic variant of ELECTRE Tri-B. Section 3 introduces an integrated framework for preference modeling and robustness analysis with ELECTRE Tri-B. In Sect. 4, we discuss the properties of the model outcomes, which are relevant for real-world applications. Section 5 demonstrates the method's applicability on a realistic data set. The last section concludes and provides avenues for future research.

Reminder on the ELECTRE Tri-B method
We use the following notation throughout the paper: . . , a n }-a finite set of n alternatives; -A R = {a * 1 , a * 2 . . .}-a set of reference alternatives on which the DM accepts to express preferences; A R ⊆ A; . . , b p }-a set of limiting class profiles such that b h−1 and b h are, respectively, the lower and upper boundaries of class C h , h = 1, . . . , p; -G = {g 1 , . . . , g j , . . . , g m }-a set of m evaluation criteria, g j : a → R for all j ∈ J = {1, 2, . . . , m}; performances of a i ∈ A or b h ∈ B on g j are denoted by, respectively, g j (a i ) and g j (b h ); if different is not stated explicitly, we assume that all criteria are of gain type (i.e., greater performances are more preferred); q j (a i ), p j (a i ), v j (a i ) for j = 1, 2, . . . , m and a i ∈ A-indifference, preference, and veto thresholds for criterion g j and alternative a i such that v j (a i ) > p j (a) ≥ q j (a) ≥ 0; -w j for j = 1, 2, . . . , m-weight of criterion g j , expressing its relative power in set G; the weights are non-negative w j ≥ 0, for j = 1, . . . , m, and normalized to sum up to one, i.e., m j=1 w j = 1; -λ ∈ [0.5, 1]-a credibility level indicating a minimal value of outranking credibility σ (a, b h ) implying the truth of a crisp outranking relation S of alternative a i over profile b h (a i Sb h ) (since λ ≥ 0.5, the support of a weighted majority of all criteria is always required to instantiate a i Sb h ).
Construction of outranking relation ELECTRE Tri-B supports the assignment of each alternative a i ∈ A to one of decision classes in C. For this purpose, the method constructs an outranking relation S, which verifies whether a i is at least as good as boundary profile b h , h = 0, . . . , p. The truth of a i Sb h is captured by means of the concordance and discordance tests. The former involves computation of a comprehensive concordance index C(a i , b h ), which quantifies the strength of a coalition of criteria confirming that a i is at least as good as b h , i.e.: where ϕ j (a i , b h ) is the marginal concordance index indicating a degree to which criterion g j supports the hypothesis a i Sb h : In particular, if a i is not worse than b h by more than indifference threshold q j (a i ) on all criteria, then C(a i , b h ) = 1, which means that the hypothesis a i Sb h is fully supported. On the contrary, if a i is worse than b h by at least preference threshold p j (a i ) for j = 1, . . . , m, then C(a i , b h ) = 0, indicating an absolute lack of support for validating a i Sb h . The discordance index measures the strength of opposition against the truth of outranking on each criterion g j , j = 1, . . . , m, in the following way: Specifically, when a i is worse than b h on g j by at least veto threshold v j (a i ), criterion g j strongly opposes to a i Sb h , hence invalidating this relation. The credibility of outranking aggregates the concordance and discordance indices to measure an overall support given to a i Sb h by all criteria: We apply a revised credibility index proposed in Dias and Mousseau (2006) so that to account for the veto-related effects within the preference disaggregation procedure. Note that when using σ (a i , b h ) originally coupled with ELECTRE Tri-B, the inference models would become non-linear. The truth of a i Sb h is confirmed iff σ (a i , b h ) is not lesser than credibility level λ, i.e.: Otherwise, i.e., if σ (a i , b h ) < λ, a i cannot be judged comprehensively at least as good as b h (a i S c b h ).

Exploitation of outranking relation
The assignment procedures of ELECTRE Tri-B exploit the outranking relation for all pairs of (a i , b h ) ∈ A × B. Specifically, a pair of disjoint sorting rules can be applied. The optimistic one assigns a i to class C h in case b h+1 is the least preferred profile such that a i S c b h+1 and b h+1 Sa i . In this case, b h+1 is strictly preferred to a i , thus preventing its assignment to a class better than C h . The pessimistic rule assigns a i to class C h in case b h is the most preferred profile such that a i Sb h . Then, a i is at least as good as b h , which hampers its assignment to a class worse than C h . Throughout the paper, we will use the pessimistic rule.

Integrated framework for preference modeling and robustness analysis with ELECTRE Tri-B
In this section, we extend the ELECTRE Tri-B method by providing the mathematical programming models that allow the DM to provide three types of indirect and imprecise preference information. We also discuss the algorithms for deriving six types of results, which contribute to the multi-dimensional robustness analysis.

Preference modeling
ELECTRE Tri-B employs an outranking-based sorting model. We expect the DM to provide the indifference q j , preference p j , and veto v j thresholds for each criterion. These thresholds can be defined as constants or affine functions. The latter ones make the threshold values dependent on the performance of particular alternatives. Moreover, the DM is asked to specify directly the boundary profiles B = {b 1 . . . , b h , . . . , b p−1 }, which serve as frontiers between classes. The extreme profiles b 0 and b p are composed of, respectively, the worst and the best performances on all criteria. Overall, we assume the variables in our model are criteria weights w j , j = 1, . . . , m, and credibility level λ. Hence, the basic set of constraints can be formulated as follows: The admissible values of weights and credibility level are delimited by the indirect and imprecise preference information. The use of such statements reduces the cognitive effort on the part of DM. In what follows, we present the mathematical programming models that allow to incorporate the following three types of preferences: assignment examples, assignment-based pairwise comparisons, and desired class cardinalities. These models involve binary variables, which are subsequently used to minimize the number of preference information pieces that need to be removed to ensure consistency between the DM's preference information and an assumed outranking-based sorting model of ELECTRE Tri-B.
Assignment examples Assignment example consists of a reference alternative a * i ∈ A R ⊆ A and its desired assignment being, respectively, the least and the most preferred classes admissible for a * i . The example statements of this type are as follows: "a * 1 needs to be assigned to at least class C 3 and at most class C 5 " or "a * 2 should be assigned to C 1 ". The following constraint set allows to model in terms of the parameters incorporated by the pessimistic rule of ELECTRE Tri-B: where M and ε are, respectively, arbitrarily large and small positive constants (this interpretation is valid also for the constraints sets that follow). In case u(a * i ) = 0, constraint [AE 1 ] stands for assigning a * i to class at least C L DM (a * i ) , whereas constraint [AE 2 ] implies that a * will be assigned to class at most C U DM (a * i ) . In case u(a * i ) = 1, the set of constraints is always satisfied, whichever the values of all remaining parameters.

Assignment-based pairwise comparisons Assignment-based pairwise comparison involves specification of a desired difference between classes to which a pair of reference alternatives
We distinguish the following two sub-types of such comparisons: -Alternative a * i is better than alternative a * j by at least k ≥ 0 classes (a * i → ≥k,DM a * j ), e.g., "a * 1 is better than a * 3 by at least 2 classes". The following constraint set translates a * i → ≥k,DM a * j to the parameters of an assumed sorting model: is equal to one, implying the occurrence of the above mentioned scenario for some h.
-Alternative a * i is better than alternative a * j by at most l ≥ 0 classes, (a * i → ≤l,DM a * j ), e.g., "a * 1 is by at most 1 class better than a * 3 ". The following constraint set translates a * i → ≤l,DM a * j in terms of the parameters of an assumed sorting model: In case u(a * i , a * j ) = 1, all constraints are satisfied, irrespective of the remaining parameter values.
In case a * i → ≥0,DM a * j and a * i → ≤0,DM a * j , alternatives a * i and a * j should be assigned to the same class.

Desired class cardinalities
In classical sorting problems, each alternative is comprehensively judged in terms of its intrinsic value through its comparison with some norms or references (in the case of ELECTRE Tri-B, these are defined by a set of boundary profiles separating the classes) (see Almeida-Dias 2010 andDias et al. 2003). Under such a scenario, sorting incorporates an absolute judgment of each alternative to be assigned. As a result, depending on the performances of alternatives on particular criteria, some classes may remain empty, whereas other classes may accommodate the vast majority of alternatives. However, the context of some decision problems may suggest that the number of alternatives assigned to some or all classes should be bounded. The example problems involving this type of constraints have been discussed in, e.g., (Almeida-Dias 2010; Dias et al. 2003;Kadziński and Słowiński 2013;. They concern credit risk assessment, assignment of bonus packages, the evaluation of the performance of retails or R&D projects, accreditation of qualification and skills, establishing national priorities in the energy sector, or student admission programs. Such requirements may be linked to both problem's nature and the DM's preferences. A variety of these examples confirms the need for methods that tolerate the lower and/or upper bounds on the category sizes.
In this regard, desired class cardinalities indicate extreme numbers of alternatives C min h,DM and C max h,DM which can be assigned to a given class C h . The example preference information pieces of this kind are as follows: "at least 3 and at most 7 alternatives can be assigned to C 1 " or "class C 2 needs to accommodate at least 5 alternatives". The respective mathematical programming model can be formulated as follows: for all a i ∈ A : Note that the use of desired class cardinalities makes sense when the applied sorting model parameters are not defined precisely. Otherwise, i.e., if ELECTRE Tri-B is parameterized with precise parameter values, the method would assign each alternative to some class in a deterministic way. Then, the class cardinalities would simply correspond to the number of alternatives accommodated in each class. Whether they would meet the pre-defined constraints cannot be controlled by the method. On the contrary, when the model parameters are defined through indirect or imprecise preference information (such as assignment examples or assignment-based pairwise comparisons), the inclusion of the requirements on the category size has an impact on which parameter sets can be used for deriving a sorting recommendation and which should be neglected, because of not being compatible with the DM's requirements of all types.
As a result, in constrained sorting problems, an assignment does not depend solely on the intrinsic value of a given alternative, but also on the assignment of remaining alternatives (Almeida-Dias 2010). Thus, a partial dependence in the considered set of alternatives occurs, incorporating relative comparisons, typical for choice and ranking, into sorting. Indeed, even if the assignment of some alternative to a given class would be possible without cardinality constraints, their incorporation may exclude such possibility in favor of the other possible assignments. This is, however, consistent with the specificity of thus defined sorting problems.
Once the desired class cardinalities are specified and the binary variables u(a i , h) are incorporated into the model, the assignment-based preference relations can be modeled in a different, more transparent way. Specifically, a * i → ≥k,DM a * j can be translated into the parameters of an assumed sorting model in the following way: In turn, a * i → ≤l,DM a * j can be modeled as follows:

Consistency verification
A set of parameter sets compatible with the DM's preference information is denoted by To ensure that it is non-empty, the following optimization problem needs to be solved: Minimi ze It minimizes the number of preference information pieces (i.e., assignment examples for a * i ∈ A R , assignment-based pairwise comparisons for (a * i , a * j ) ∈ P R , and desired class cardinalities for C h , h ∈ H ) that need to be removed to restore consistency. The instantiation of each statement is controlled by the respective binary variable, i.e., u(a * i ) for assignment examples, u(a * i , a * j ) for assignment-based pairwise comparisons, and u(C h ) for desired class cardinalities. In case all preference information pieces can be reproduced together, the optimal value of the objective function would be equal to zero. Otherwise, the binary variables equal to one would indicate the DM's statements that need to be removed to reinstate compatibility with an assumed outranking-based sorting model.

Robustness analysis
In case a set of indirect and imprecise preference information pieces is not contradictory, infinitely many sorting models can reproduce it. Their exploitation can be conducted in two ways. On the one hand, a representative set of parameters can be selected using some arbitrary procedure. The application of thus constructed sorting model on the set of alternatives implies univocal recommendation. On the other hand, all compatible sets of parameters can be exploited to conduct robustness analysis. In the proposed approach, we implement the latter idea, while extending it to provide multiple complementary results quantifying the stability of sorting recommendation given the plurality of outranking-based models compatible with the DM's incomplete preferences. Since we treat criteria weights and credibility levels as variables, all compatible sets of parameters defined by E(A R ) are denoted by S(w, λ).

Robustness analysis with mathematical programming
In this section, we focus on deriving exact robust results by means of mathematical programming. These outcomes include: -Possible assignment C P (a i ) indicating a set of classes to which alternative a i ∈ A is assigned by at least one compatible set of parameters (a i → P C h ), i.e. C P (a i ) = {C h : ∃(w, λ) ∈ S(w, λ) for which a i → C h }; -Necessary assignment-based preference relation →,N which holds for pairs of alternatives (a i , a j ) ∈ A × A, such that a i is assigned to a class at least as good as a j for all compatible sets of parameters, i.e., a i →,N a j ⇔ ∀(w, λ) ∈ S(w, λ) a i → C h and a j → C k where h ≥ k; -Extreme class cardinalities indicating the minimal C min h and maximal C max h numbers of alternatives assigned to class C h for some compatible set of parameters, i.e., C min Constraints [P A 1 ] and [P A 2 ] are responsible for assigning alternative a i to class C h by ensuring that a i outranks the lower profile and does not outrank the upper profile of C h . Let us denote the indices of the least and the most preferred classes to which a i can be possibly assigned by, respectively, L(a i ) and R(a i ).
Necessary assignment-based preference relation The purpose of the assignment-based preference relation is to enable the pairwise comparisons of alternatives in a manner compatible with the sorting method. This is not possible when referring to the assignment-or cardinalityoriented results. The analysis for pairs may be particularly useful when the DM is interested in the recommendation obtained for some particular alternatives or when there are numerous alternatives possibly assigned to the same class range. By comparing the alternatives in terms of the robust assignment-based relations, one can arrive at conclusions such as: -"all compatible outranking models assign a to a class at least as good as b" (in case of analyzing necessary assignment-based preference relation), or -"for 50% of compatible outranking models a is assigned to a class better than b, for 20% of feasible models the order of classes is inverse, and for the remaining 30%-a and b are placed in the same class" (in case of analyzing Assignment-based Outranking Indices).
These conclusions are compatible with the sorting problem, and cannot be obtained through analyzing the stability of an outranking or preference relation considered in the context of ranking problems.
When it comes to the necessary assignment-based preference relation →,N , for a i , a j ∈ A, a i →,N a j holds if ∀h ∈ H , E h,C (a i →,N a j ) given below is not feasible or, in case of feasibility, ε * = max ε subject to E h,C (a i →,N a j ) is not greater than 0.
Constraints [P N 1 ] and [P N 2 ] ensure that a j is assigned to class at least C h+1 , whereas a i is assigned to class at most C h . Hence, to verify if a i is assigned to a class at least as good as a j for all compatible sets of parameters, we check if a j can be assigned to class better than a i for at least one feasible parameter set. Minimi ze/Maximi ze :

Extreme class cardinalities
where u(a i , h) is equal to one in case the conditions justifying an assignment of a i to C h are satisfied. Hence, by minimizing (maximizing) the sum of u(a i , h) for all a i ∈ A, we identify an outranking model in the set of all compatible models defined by E(A R ) that minimizes (maximizes) the number of alternatives for which the conditions supporting an assignment to C h are met at the same time.
The possible assignments, necessary assignment-based preference relation, and extreme class cardinalities offer complementary perspectives on the robustness of results. None of them can be fully captured based on the analysis of the remaining ones. This justifies their joint employment. For example, when alternative a i is possibly assigned to classes between C 3 and C 5 and alternative a j is possibly assigned to classes between C 2 and C 4 , a i is not guaranteed to be always assigned to a class at least as good as a j , so it may only be verified by means of the necessary assignment-based preference relation. Furthermore, there may exist numerous alternatives that are possibly assigned to class C h , but this may happen for different compatible sets of parameters. Thus, only the analysis of extreme class cardinalities captures how many alternatives can be placed in the same class for some feasible set of parameters.

Robustness analysis with the Monte Carlo simulations
The exploitation of a set of compatible sorting models using mathematical programming techniques leads to the exact robust results. They reveal what happens for some, all, and extreme compatible sets of parameters. Nonetheless, the possible assignment often includes multiple classes, the necessary relation leaves many pairs of alternatives, whereas the range of observable class cardinalities can be wide. In this perspective, it would be useful to compute the results which quantify the probability of attaining different outcomes given the multiplicity of compatible sorting models. Specifically, we consider the following results, which are defined as multi-dimensional integrals over the space of uniformly distributed weights and credibility levels compatible with the DM's preference information: -Class Acceptability Index C AI (a i , h) indicating the proportion of compatible sets of parameters assigning alternative a i ∈ A to class C h , i.e.
where m(w, λ, i, h) is the class membership function: -Assignment-Based Outranking Index AO I (a i , a j ) indicating the share of compatible sets of parameters for which alternative a i ∈ A is assigned to a class at least as good as alternative a j ∈ A, i.e.: where r (w, λ, i, j) is the assignment-based outranking relation confirmation function: r (w, λ, i, j) = 1, if a i → C h and a j → C k and h ≥ k for (w, λ), 0, otherwise.
-Class Cardinality Index CC I (h, r ) indicating the share of compatible sets of parameters for which the size of class C h (i.e., the number of alternatives assigned to it) is equal to r , i.e.
where c(w, λ, h, r ) is the class cardinality function: Since S(w, λ) is multi-dimensional and non-convex, it is not possible to compute the stochastic acceptability indices accurately. Computing the exact value of these indices would require us to (1) calculate the volume of the full space exactly, (2) for each alternative, pair, or class to partition the space in regions that grant a particular result, and (3) to compute the volumes of these regions. Exact computation of the volume of a polytope is extremely hard from the computational viewpoint (Kadziński and Tervonen 2013). In turn, to estimate C AI s, AO I s, and CC I s, we use the Monte Carlo simulation to sample from the space of uniformly distributed criteria weights and credibility levels. Since there do not exist efficient algorithms for sampling the compatible outranking models from the non-convex space, we combine the HAR algorithm ) with a naive rejection technique. Specifically, we first sample from a convex space of compatible parameter sets defined by basic normalization and non-negativity constraints as well as the DM's assignment examples, i.e., E B ASE ∪ E ASS−E X . Note that these constraints introduce the linear cuts in the space of feasible weights and credibility levels, which makes HAR suitable for efficiently sampling from such a space (Lovász 1999). However, the constraints introduced by the assignment-based pairwise comparisons and desired class cardinalities involve the binary variables, thus making the space of compatible parameters sets non-convex. In this regard, to find the samples which are compatible with all types of DM's preference information, we reject these samples returned by HAR that do not align with E PC L ∪ E PCU ∪ E C ARD . The latter is justifiable because a rejection sampling is feasible in low dimensionality problems, and acceptable error limits for the stochastic indices can always be achieved, given sufficient computation time, with a certain number of Monte Carlo iterations (Kadziński and Tervonen 2013).

Properties of the model outcomes
In this section, we present the properties of robust sorting results. We first discuss the impact of preference information on the outcomes. Then, we present how these results evolve with an incremental specification of preference information. We exhibit the interrelations between the exact results derived with mathematical programming and the stochastic acceptability indices estimated with the Monte Carlo simulations. All these properties are of high practical relevance. However, since their proofs are obvious, we omit them. Finally, we discuss some measures for quantifying comprehensive robustness of results.

Relations between provided preference information and model outcomes
The approach proposed in this paper has been designed so that to ensure a correspondence between the preference information and the robust results. The relations between the provided input and the derived outputs exist on three levels concerning individual alternatives, pairs of alternatives, and decision classes. They can be summarized as follows: -for each reference alternative a * i ∈ A R its possible assignment is a subset of the desired assignment, i.e. C P (a ; that is, the set of classes confirmed by at least one compatible parameter set can be narrower than the range specified by the DM, but never wider; -for each pair of reference alternatives a * i , a * j ∈ A R such that a * i → ≥k,DM a * j with k ≥ 0, the necessary assignment-based preference relation holds, i.e. a * i → ≥k,DM a * j , k ≥ 0 ⇒ a * i →,N a * j ; -for each class C h with specified desired class cardinalities, the extreme class cardinalities are within the bounds provided by the DM, i.e., C min h,DM ≤ C min h and C max h ≤ C max h,DM . Evolution of robust results with incremental specification of preference information The analysis of the robust result stimulates the DM to enrich his/her preference information. For example, the analysis of C P (a i ) and C AI s may lead to new or more precise assignment examples, →,N and AO I s may enhance the DM to provide additional assignment-based pairwise comparisons, whereas C min h , C max h , and CC I s may stimulate more constrained desired class cardinalities. With each new additional piece of preference information, the set of compatible parameter sets becomes smaller. Thus, when considering the outcomes obtained in two subsequent iterations t and t + 1 with preference information in (t + 1)th iteration enriching the DM's requirements in tth iteration, the following properties concerning the evolution of results hold: -for each alternative a i ∈ A its possible assignment can become narrower, but not wider, i.e., C t+1 P (a i ) ⊆ C t P (a i ); -the assignment-based necessary preference relation is enriched, i.e., →,N t ⊆ →,N t+1 ; -the difference between extreme class cardinalities can become lesser, but not greater, i.e., C min h,t ≤ C min h,t+1 and C max h,t+1 ≤ C max h,t .

Interrelations between the exact outcomes and estimates of stochastic acceptability indices
The results computed with mathematical programming are exact. Specifically, the necessary outcomes are confirmed by all compatible parameter sets, the possible results-by at least one feasible set of parameters, whereas the extreme ones indicate the recommendation obtained in the most and the least advantageous cases. On the contrary, the results derived from the Monte Carlo simulation are only estimates of the exact acceptability indices. Thus, although they can be estimated up to an arbitrarily selected accuracy threshold by analyzing a required number of samples, they would never be exact. In what follows, we list the interrelations between the necessary, possible, and extreme outcomes and the estimates of stochastic acceptability indices. Let us start with the results concerning the class assignments of individual alternatives: When it comes to the results observed for pairs of alternatives, the following properties hold: a i →,N a j ⇒ AO I (a i , a j ) = 1; -AO I (a i , a j ) < 1 ⇒ ¬(a i →,N a j ).
As far as the class cardinalities are concerned, the following interrelations can be considered: The above properties imply that the DM's preference information has an impact also on the stochastic acceptability indices. As a result, the latter ones evolve with an incremental specification of preferences. However, it is not possible to formulate some general conclusions on how a particular acceptability index would be modified in case the respective exact result would be either possible, though not necessary, or it would be contained within the extreme, non-precise bounds.
One needs to bear in mind that C AI s, AO I s, and CC I s are not sufficient to derive the exact results. Since we consider the estimates of acceptability indices, C AI (a i , h) = 0 and . Moreover, C AI (a i , h) = 1 and AO I (a i , a j ) = 1 do not imply C P (a i ) is a singleton or a i →,N a j , respectively.

Robustness measures for multiple criteria sorting Let us introduce a few measures quantifying comprehensive robustness of results in view of six types of outcomes obtained from the analysis of S(w, λ):
-an average difference between indices of the most and the least preferred possible classes, built on the possible assignments C P (a i ) of all alternatives a i ∈ A: -a mean class assignment entropy, built on the class acceptability indices for all alternatives a i ∈ A and all classes C h , h = 1, . . . , p: -the share of ordered pairs of alternatives for which the order of classes is not fixed, which means that the necessary assignment-based preference relation is not true, but the possible relation holds (which is approximated with AP O I (a i , a j ) > 0): where: -a mean assignment-based outranking entropy, built on the assignment-based outranking indices for all pairs of alternatives (a i , a j ) ∈ A × A, i = j: -a mean class cardinality entropy, built on the class cardinality indices for all classes C h , h = 1, . . . , p, and all possible cardinalities r = 1, . . . , n: The above measures can be analyzed by the DM, who incrementally provides preference information. All measures have been designed in such a way that lesser values indicate greater robustness of results. In this perspective, the DM can judge, after each iteration, whether the interactive process should be continued, in case the robustness is satisfying, or stopped, otherwise.

Illustrative study: multiple criteria evaluation of MBA programs
In this section, we present an illustrative study to demonstrate the applicability and practical usefulness of the proposed framework for preference modeling and robustness analysis with ELECTRE Tri-B. The actual data concerns evaluation of 30 MBA programs (Köksalan et al. 2009) (see Table 1). The considered gain-type criteria are defined on a 0-100 scale, and involve alumni career progress (g 1 ), diversity (g 2 ) and idea generation (g 3 ) (the performances are provided in Table 2).  Table 3 Performances of the boundary class profiles and indifference q j , preference p j , and veto v j thresholds for all criteria g j Boundary profiles We aim at assigning the MBA programs to five quality classes C 1 , C 2 , C 3 , C 4 , C 5 , with C 1 and C 5 indicating the least and the most preferred categories, respectively. The performances of boundary profiles b 0 , b 1 …, b 5 as well as the indifference, preference, and veto thresholds are provided in Table 3.
We will illustrate the use of the proposed methodological framework in its most general form, i.e., by referring to all types of indirect and imprecise preference information and robust outcomes. However, there is no obligation to use all kinds of inputs and outputs in case some of them are not found appealing by the DM (e.g., in case of assignment-based pairwise comparisons) or they are not imposed by the problem's nature or DM's preferences (e.g., in case of desired class cardinalities). The versatile robustness analysis can be conducted for a set of outranking models compatible with any preference information.
Also, the proposed framework for robustness analysis may support the elicitation of indirect preference information already in the first iteration. In particular, it would be possible to determine the results when only the boundary profiles and comparison thresholds are provided. Subsequently specified holistic judgments (e.g., assignment examples) or requirements (e.g., desired class cardinalities) should be consistent with the outcomes derived from such an analysis (e.g., the possible assignments or extreme class cardinalities, respectively). In this way, they can be linked to a prior specification of an outranking model through profiles and thresholds.

Robustness analysis in the first iteration
We start with different statements representative for the three allowed types of preference information. In particular, we consider five assignment examples, one for each class: a 9 → C 1 , a 23 → C 2 , a 3 → C 3 , a 21 → C 4 , and a 1 → C 5 . Moreover, we account for three assignment-based pairwise comparisons: a 3 is better than a 2 by at most one class, a 5 is better than a 17 by at least two classes, and a 27 is worse than alternative a 23 by at most two classes. Finally, we tolerate the extreme class cardinalities provided in Table 4. They impose that each class accommodates at least a few alternatives and prevent too many MBA programs to be classified in the same way. The provided preference information is consistent with an assumed sorting model. Thus, the set of compatible weights and credibility levels is non-empty and can be exploited in the spirit of robustness analysis. The possible, necessary, and extreme results are computed with mathematical programming techniques, whereas the values of acceptability indices are estimated with the Monte Carlo simulation. In Fig. 1, we present the samples considered in the first iteration. Each dot corresponds to a compatible set of parameters composed of three weights (w 1 , w 2 , w 3 ) corresponding to the coordinates according to the three axes and credibility level λ, whose value is represented with a color on a pre-defined scale. For clarity of presentation, some randomly selected dots are described with precise parameter values. All dots, including small and large ones, are counterparts of 10,000 sets of parameters compatible with the assumptions of the outranking-based model as well as assignment examples. They were obtained with HAR. The analysis of their distribution in the space of all feasible parameters indicates that the space of compatible weights and credibility levels is not constrained much by the assignment examples (see dots in Fig. 2). The large dots in Fig. 1 represent the sets of parameters compatible with all pieces of preference information, including assignment-based pairwise comparisons and desired class cardinalities (see Fig. 3). Hence, the small dots correspond to the samples that were neglected due to the application of a rejection sampling technique, whereas the large ones are used to estimate the stochastic results. As indicated by Fig. 1, assignment-based pairwise comparisons and desired class cardinalities significantly reduced the space of compatible weights and credibility levels. In Table 5, we provide the possible assignments for all alternatives. They are precise for 24 alternatives. For example, 3 alternatives (a 1 , a 5 , a 10 ) are univocally assigned to the best  Table 5 Possible assignments for all alternatives in the first iteration a 10 C 5 a 20 C 3 , C 4 a 30 C 2 class (C 5 ), whereas 6 alternatives (a 9 , a 13 , a 15 , a 19 , a 27 , a 29 ) are assigned to the worst class (C 1 ) with all compatible sets of parameters. All DM's assignment examples are reproduced in the possible assignments (e.g., since a 23 was assigned by the DM to class C 2 , we have C P (a 23 ) = {C 2 }). For the remaining 6 alternatives (a 8 , a 12 , a 20 , a 24 , a 26 , and a 28 ), the possible assignments are imprecise. Specifically, these alternatives are sorted into one of two consecutive classes depending on the selected compatible set of parameters.
When it comes to C AI s, for the alternatives with precise possible assignments, they are equal to one for the sole admissible class and zero for all remaining classes. The analysis of C AI s is the most informative for the alternatives with imprecise possible assignments. Indeed, for some of them, the probability of being assigned to either of two admissible classes is similar. For example, the shares of compatible sets of parameters assigning a 24 to C 3 and C 4 are equal to, respectively, 51.09% and 48.91%. On the contrary, for some other alternatives, the vast majority of feasible weights and credibility levels assign them to a single class. In particular, for a 11 , the probability of being assigned to C 4 is eight times greater than the respective probability for C 3 .
The other perspective on the robustness of sorting recommendation concerns the relative comparison of assignments for pairs of alternatives. A Hasse diagram indicating a partial pre-order imposed on the set of alternatives by the necessary assignment-based preference relation is presented in Fig. 4. It respects the DM's assignment-based pairwise comparisons. For example, since a 23 and a 5 were required to be assigned to a class better than, respectively, a 27 and a 17 , →,N is true for these pairs.
To support understanding of the interrelations between →,N and C P (a i ), the nodes representing alternatives with the same possible assignments are marked with the same color. In this regard, some relations follow the possible assignments directly. For example: -alternatives assigned to C 5 with all compatible sets of parameters (a 1 , a 5 , a 10 ) are preferred over all remaining alternatives, whereas alternatives possibly assigned only to C 1 (a 9 , a 13 , a 15 , a 19 , a 27 , a 29 ) are worse than all other alternatives; -alternatives assigned to C 3 (a 3 , a 12 , a 17 ) are preferred over those possibly assigned to either C 3 or C 2 (a 8 and a 26 ); -alternatives precisely assigned to the same class are indifferent in terms of the necessary assignment-based preference relation, and hence they are grouped in the same node in Fig. 4.
However, some relations cannot be derived solely from the analysis of possible assignments. For example, when considering alternatives possibly assigned to class C 3 or C 4 , a 11 is incomparable with a 20 . This means that for some feasible sets of parameters, a 11 is assigned to a class better than a 20 , whereas for some other compatible weights and credibility levelsthe order of classes for this pair is inverse. Furthermore, although the possible assignments for a 20 and a 24 are the same, a 20 is always is assigned to a class at least as good as a 24 , whereas the opposite is not true. This suggests that for some feasible sorting models, a 20 is assigned to a class strictly better than a 24 . The analysis of the necessary assignment-based preference relation can be enriched with consideration of Assignment-based Outranking Indices. We report them in the form of a heatmap, which translates the index values between 0 and 1 on a pre-defined color scale (see Fig. 5). As indicated by the associated color scale, the darker is a specific cell in a matrix, the greater is the proportion of compatible sets of parameters for which an alternative from the row is assigned to a class at least as good as an alternative from the column.
For pairs of alternatives related by the necessary assignment-based preference relations, AO I s are equal to one. In particular, AO I (a 1 , a 5 ) = 1 and AO I (a 5 , a 1 ) = 1. Furthermore, for pairs with disjoint sets of possible assignments, AP O I is equal to zero in case one of them is always assigned to a class strictly worse than the other (e.g., AO I (a 12 , a 14 ) = 0). Let us emphasize that even if an alternative is not necessarily preferred to some other alternative, it can be assigned to a class at least as good for some feasible sets of parameters. In particular, a 3 is assigned to a class not worse than a 11 , a 20 , and a 24 by, respectively, 11%, 41%, and 51% of compatible weights and credibility levels. Finally, the analysis of AO I s is the most informative for pairs, which are incomparable in terms of →,N . For example, the proportion of compatible sets of parameters for which a 11 is assigned to a class not worse than a 20 is greater than the share of feasible models for which the inverse relation holds. The exact values of Assignment-based Outranking Indices are provided for all pairs in the e-Appendix (supplementary material available online).
The last perspective concerns class cardinalities. In Table 6, we present the extreme numbers of alternatives that can be assigned to each class by some compatible sets of parameters. For all classes, they are more precise than the cardinalities desired by the DM. For example,  the minimal number of alternatives that were allowed to be assigned to C 5 was 5, whereas the least observable cardinality for this class is 8. In the same spirit, the maximal number of alternatives that could be placed in C 3 was 8, whereas the greatest size of C 3 that is observed for some feasible set of parameters is 6. Furthermore, for all feasible sets of weights and credibility levels, only 3 alternatives are assigned to C 5 . This is understandable as they were the only alternatives that could possibly be assigned to the most preferred class. For the remaining classes, the minimal and maximal sizes differ. For class C 1 these extreme numbers (C min h = 6 and C max h = 7) agree with the analysis of possible assignments as they were 7 alternatives which could be possibly assigned to C 1 , but only six of them were assigned to this class for all compatible sets of parameters. However, for some other classes, these extreme cardinalities cannot be derived from the analysis of other results. For example, the minimal and maximal numbers of alternatives assigned to C 3 are equal to, respectively, 4 and 6. However, the numbers of alternatives that are always or sometimes assigned to C 3 are equal to 3 and 8, respectively. These sizes are not observed for any feasible sets of parameters.
The analysis of extreme class cardinalities can be supported with consideration of class cardinality indices. Table 7 represents a heatmap of CC I s. The cells without a number correspond to the cardinalities, which are not attained for any compatible set of parameters (i.e., CC I (h, r ) = 0). Clearly, for C 5 -the size of 3 is confirmed by all considered samples (CC I (5, 3) = 1). For the remaining classes, the analysis of CC I s offers rich insights on how many alternatives they accommodate for different compatible sets of parameters. In particular, the vast majority of feasible sets of weight and credibility levels indicate that there are 6 and 10 alternatives assigned to, respectively, C 1 and C 2 (i.e., CC I (1, 6) = 0.74 and CC I (2, 10) = 0.76). Furthermore, for probabilities of 4, 5, and 6 alternatives being assigned to C 3 are very similar. Finally, for C 4 -the size of 4 is the least probable, being confirmed only by 1% of compatible sets of parameters.

Robustness analysis in the second iteration
Let us assume that the analysis of robust results in the first iteration stimulated the DM to provide additional preference information of different types: -alternative a 11 should be assigned to class C 3 implied by the analysis of its imprecise possible assignment and Class Acceptability Indices; -alternative a 24 should be assigned to a class better than a 26 based on the analysis of Assignment-based Outranking Indices; Table 7 Class Acceptability Indices in the first iteration The compatible sets of parameters derived with the Monte Carlo simulation for thus enriched preference information are presented in Fig. 7. The sets of parameters compatible solely with assignment examples (see Fig. 8) as well as with the remaining types of preference information (see Fig. 9) have been vastly constrained when compared with the first iteration. This is understandable in view of how the preferences were formulated in the second iteration, i.e., they excluded some results previously deemed as possible.
The outcomes that were confirmed by all compatible sets of parameters in the previous iteration are also validated in the next iteration (e.g., a univocal possible assignment of a 14 to C 3 ; a necessary assignment-based preference of a 2 over a 28 , or a precise cardinality of class C 5 ). Thus, in what follows, we focus only on the evolution of results that were not validated for all feasible sets of weights and credibility levels in the first iteration.
When it comes to the possible assignments (see Table 8), they have become precise for another five alternatives. For example, a 8 and a 26 are assigned to C 2 for all compatible sets of parameters, whereas previously, they could also be placed in C 3 . In fact, there is only one alternative (a 20 ) for which the possible assignment is still imprecise. The analysis of C AI s indicated that 53.17% compatible models assign it to C 3 , whereas 46.83% sort a 20 to C 4 . Thus, class C 3 is slightly more probable for a 20 given indeterminacy of the DM's preference model due to the incompleteness of his/her preference information, whereas in the previous iteration-class C 4 was indicated by almost 60% of feasible sorting models.
The necessary assignment-based preference relation is presented in the form of a Hasse diagram in Fig. 10. The subsets of alternatives which are assigned to the same class by all feasible sets of parameters have now become more numerous. For example, 10 alternatives are indifferent in terms of →,N due to being precisely assigned to C 2 . Moreover, there is no single pair of alternatives that would be incomparable in terms of the necessary relation. For example, a 20 is now assigned to a class at least as good as a 11 for all compatible weights and credibility levels. When it comes to the Assignment-based Outranking Indices, the only values which are neither 0 nor 1 can be observed for the comparison of a 20 with other alternatives. Specifically, it is assigned to the same class (i.e., C 4 ) as a 14 , a 16 , a 21 , and a 25 for less than half of feasible sorting models, whereas a 3 , a 11 , a 12 , a 17 , and a 24 are assigned to the same class (i.e., C 3 ) for more than 50% of compatible parameter sets (Fig. 11).
The extreme class cardinalities indicate that the numbers of alternatives assigned to each class are now more precise for all classes. On the one hand, there are always 7, 10, and 3 alternatives which are assigned to, respectively, C 1 , C 2 , and C 5 . On the other hand, C 3 accommodates from 5 to 6 alternatives, whereas the cardinality of C 4 is between 4 and 5. The respective CC I s are given in Fig. 12. They confirm that it is slightly more probable that 4 and 6 alternatives are assigned to, respectively, C 4 and C 3 , which is consistent with the previously discussed results related to the assignment of a 20 .
The DM should analyze the obtained class assignments and cardinalities as well as assignment-based preference relation after each iteration. (S)he should judge whether the suggested recommendation is convincing and decisive enough, and the correspondence between the output of the model and the preferences (s)he has at the moment is satisfying. This decision can be supported by the analysis of robustness measures (see Table 9). For example, the average difference between indices of the most and the least preferred possible classes has decreased from 1.2 to 1.033 (see f C PW ). In fact, the possible assignment is now imprecise Fig. 9 Sets of parameters derived with the Monte Carlo simulation compatible with all types of preference information provided in the second iteration Table 8 Possible assignments for all alternatives in the second iteration only for 1 out of 30 alternatives as compared to 6 alternatives in the first iteration. In the same spirit, the share of ordered pairs of alternatives f N AP O I for which the order of classes is not fixed has decreased from 0.0793 to 0.0103 (i.e., from 69 to 9 out of 870 ordered pairs), whereas an average difference f CCW between the greatest and the least class cardinalities has dropped from 2.6 to 1.4. We can also observe significant improvement of the entropy-based If the DM judges the robustness of results as satisfying, the interactive process stops. We terminate the illustration after the second iteration because the recommendation is already very robust. The only uncertainty is related to the assignment of a 20 . Once the DM would assign it to either C 3 or C 4 , all compatible models would already recommend the same sorting decision.
In case the robustness of sorting outcomes would not be found sufficient by the DM, (s)he should pursue the exchange of preference information. If (s)he changed her/his mind or discovered that the expressed judgments were inconsistent with some previous judgments that (s)he considers more important, (s)he may backtrack to one of the earlier iterations and continue from this point. In this way, the process of preference construction is either continued or restarted (Corrente et al. 2013). programming models for incorporating three types of indirect and imprecise preferences. These include assignment examples (e.g., "alternative a 1 should be assigned to C 2 " or "alternative a 2 should be placed in either C 3 or C 4 "), assignment-based pairwise comparisons (e.g., "alternative a 3 should be assigned to a class better than alternative a 4 " or "alternatives a 5 and a 6 need to be assigned to the same class"), and desired class cardinalities (e.g., "the number of alternatives assigned to class C 1 should be between 3 and 6"). On the other hand, we provided the algorithms for computing multiple robust results. The exact outcomes, including the possible assignments, necessary assignment-based preference relation, and extreme  class cardinalities, can be derived with mathematical programming techniques. Furthermore, the stochastic acceptability indices quantifying the probability of a particular assignment, the truth of assignment-based relation, or some specific number of alternatives assigned to a given class are estimated with the Monte Carlo simulation. For this purpose, we used an original combination of the Hit-And-Run algorithm and the rejection sampling technique. Overall, we provide the richest MCDA framework for robustness analysis. Indeed, (Kadziński and Ciomek 2016) and (Kadziński et al. 2015) considered only exact robust outcomes in the context of outranking-and value-based sorting procedures. Let us emphasize that when applying the method to real-world decision problems, the DMs can use only preference information that they find relevant for a particular context and refer to the results that they consider useful for supporting decision making. The usefulness of multiple accounted inputs and outputs was discussed theoretically based on several properties and demonstrated on a realistic example concerning the evaluation of MBA programs. Firstly, we revealed the impact of provided preference relation on the results. It is enhanced by the correspondence between different perspectives concerning individual alternatives, pairs of alternatives, and decision class. Secondly, we showed how the results evolve with an incremental specification of preference information. In this perspective, diverse robust outcomes play a crucial role in stimulating the DM to supply additional preference information pieces, whereas the contraction of the set of compatible parameter sets implies that the results become more and more robust. Thirdly, we focussed on why it is useful to analyze the outcomes concerning assignments, assignment-based relation, and class sizes. Although they are interrelated, none of them can be fully derived from the analysis of the remaining ones. Fourthly, we discussed the benefits of using exact robust outcomes and stochastic acceptability indices. Specifically, we demonstrated how these two types of outcomes complement each other with exact results indicating which recommendation is certain, possible, or extreme and stochastic indices filling a gap, in a probabilistic way, between what is possible and certain. Finally, we proposed some measures for quantifying the comprehensive robustness of sorting results. Their analysis may stimulate the decision on terminating an interaction with the method or pursuing an incremental specification of preference information.
The idea of processing diverse human-readable preference information and providing sophisticated, robust results at the method's output could be transferred to other MCDA approaches. The most appealing extension refers to the ELECTRE Tri-C method (Almeida- , which -unlike ELECTRE Tri-B-requires the DM to provide characteristic rather than boundary class profiles and employs two conjoint rules to work out an assignment for each alternative.
In what follows, we discuss the limitations of our approach and the related directions for future research. First, in our sorting model, we considered only criteria weights and credibility levels as variables whose values need to be inferred from indirect and imprecise preference information. However, an outranking-based sorting model requires specification of other parameters such as boundary profiles or indifference, preference, and veto thresholds. As demonstrated by a rich literature on preference disaggregation in the context of ELECTRE Tri-B, it is feasible to learn these parameters from the DM's holistic judgments. However, it is not possible to infer them all at once with linear programming techniques. Thus, it would require the development of dedicated non-linear optimization algorithms.
Second, our method has been designed for an interactive specification of relatively small subsets of indirect and imprecise statements. However, when larger sets of preferences are involved, inconsistencies occur more frequently. Moreover, the number of binary variables involved in specifying compatible sorting models is proportional to the number of alternatives and classes. When considering the insufficient efficiency of existing solvers, it would be neither possible to identify inconsistencies nor to model the preferences using Mixed-Integer Linear Programming in case the number of binary variables would go into thousands. Hence we envisage an extension of the method for dealing with large sets of preference information pieces and large decision problems by developing some heuristic preference learning approaches.
Third, we left the DM freedom to select the alternatives, pairs of alternatives, or classes for which (s)he specified some desired outcomes or constraints that should be respected by the recommended decision. As demonstrated by the recent research on active learning (see, e.g., Ciomek et al. 2017), it is possible to design some algorithms to indicate the objects that should be judged by the DM so that to increase both information gain from the DM's answer and robustness of the results. These approaches should incorporate the measures for quantifying the robustness of sorting outcomes proposed in this paper.
Fourth, when discussing the results of a case study, we visualized the compatible sets of parameters. This is very useful for understanding how the DM's indirect and imprecise statements constrain the space of an outranking-based sorting model's parameters. In the considered study, the number of model variables was four. Thus, we used three-dimensional coordinate systems to illustrate the weight values and color for reflecting the values of a credibility level. When the fifth variable would be involved, it would be possible to show its values using points' sizes. However, with an even greater number of model parameters, it would be already challenging to ensure the figure's readability.
Finally, we demonstrated the method's applicability on a study concerning the evaluation of MBA programs. We simulated all types of preference information pieces and analyzed all types of robust results. In this way, we could illustrate the capabilities of the proposed approach. The ultimate verification of its practical usefulness will involve applying it in case studies with real-world DMs.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.