Multiple Criteria Assessment of Insulating Materials with a Group Decision Framework Incorporating Outranking Preference Model and Characteristic Class Proﬁles

We present a group decision making framework for evaluating sustainability of the insulating materials. We tested thirteen materials on a model that was applied to retroﬁt a traditional rural building through roof’s insulation. To evaluate the materials from the socio-economic and environmental viewpoints, we combined life cycle costing and assessment with an adaptive comfort evaluation. In this way, the performances of each coating material were measured in terms of an incurred reduction of costs and consumption of resources, maintenance of the cultural and historic signif-icance of buildings, and a guaranteed indoor thermal comfort. The comprehensive assessment of the materials involved their assignment to one of the three preference-ordered sustainability classes. For this purpose, we used a multiple criteria decision analysis approach that accounted for preferences of a few tens of rural buildings’ owners. The proposed methodological framework incorporated an outranking-based preference model to compare the insulating materials with the characteristic class proﬁles while using the weights derived from the revised Simos procedure. The initial sorting recommendation for each material was validated against the outcomes of robustness analysis that combined the preferences of individual stakeholders either at the output or at the input level. The analysis revealed that the most favorable materials in terms of their overall sustainability were glass wool, hemp ﬁbres, kenaf ﬁbres, polystyrene foam, polyurethane, and rock wool.


Introduction
This paper presents a group decision framework for evaluating sustainability of the insulating materials to retrofit traditional rural buildings. The importance of this research derives from the previous studies on both retrofitting solutions tailored to traditional rural buildings as well as judging an overall desirability of coating materials (see, e.g., Krarti 2015;Fabbri et al. 2012;Ma et al. 2012;Yung and Chan 2012;Martínez-Molina et al. 2016). These studies prove that energy efficiency and thermal comfort are crucial for the maintenance of historic buildings.
The context of the study is that of a typical farmhouse in central Italy. The incorporated building model derives from the analysis of over 800 farmhouses surveyed by the census of the scattered rural buildings of the municipality of Perugia (Umbria region). The high landscape values of traditional buildings and the legislation about their preservation prevent external alterations (Mazzarella 2015). Therefore, the most viable solutions are to intervene on the roof of these structures, increasing their thermal inertia with coating materials (Verbeeck and Hens 2005;Kumar and Suman 2013;Taylor et al. 2000).
We comprehensively evaluate the materials for the roof insulation by considering economic, social, and environmental viewpoints. For this purpose, we incorporate a life cycle costing (LCC) approach, a life cycle assessment (LCA), and a dynamic thermal simulation for the evaluation of energy savings and thermal comfort. As such, we aim at identifying the materials that guarantee the indoor thermal comfort, at the same time reducing the consumption of resources in their entire life cycle as well as maintaining cultural and historic significance of the buildings. In this perspective, we differentiate from the vast majority of previous studies concerning coating materials which incorporate a mono-disciplinary approach (Copiello 2017).
To provide an overall sustainability assessment of coating materials, we incorporate Multiple Criteria Decision Analysis (MCDA). MCDA offers a diversity of approaches designed for providing the decision makers (DMs) with a recommendation concerning a set of alternatives evaluated in terms of multiple conflicting points of view. Few applications of MCDA methods for the evaluation of building materials, which are reported in the literature (Ginevicius et al. 2008) deal mainly with the environmental sustainability of materials (Papadopoulos and Giama 2007;Khoshnava et al. 2016). Some combinations of LCA and MCDA were considered by Santos et al. (2017) and Piombo et al. (2016). Applications which included both LCC and LCA for the definition of criteria to be used in MCDA are still rare (Piombo et al. 2016). Decision analysis methods used in the above-mentioned studies involved different variants of AHP (Motuziene et al. 2016;Khoshnava et al. 2016), PROMETHEE II (Kumar et al. 2017), Weighted Sum, TOPSIS (Čuláková et al. 2013), VIKOR, and COPRAS (Ginevicius et al. 2008).
From the viewpoint of MCDA, our study differs from the aforementioned ones in terms of the following major aspects: • We formulate the considered problem in terms of multiple criteria sorting, thus aiming at assigning the materials to a set of pre-defined and ordered sustainability classes (categories) rather than at ordering them from the best to the worst; • We assess the insulating materials while taking into account preferences of multiple DMs (owners of rural houses), thus incorporating group decision making tools into the evaluation framework; • The adopted assignment procedure builds upon outranking-based comparison of the insulating materials with the characteristic profiles composed of the per-class most representative performances on all criteria (Kadziński et al. 2015b); • The research results are validated against the outcomes of robustness analysis that takes into account all sets of weights compatible with either the ranking of criteria provided by each DM within the revised Simos (SRF) procedure (Figueira and Roy 2002) or a group compromise ranking of criteria that is constructed with an original procedure proposed in this paper.
The remainder of the paper is organized in the following way. In the next section, we review the existing group decision making methods for multiple criteria sorting. Section 3 describes a three-stage decision aiding method that has been used to evaluate the insulating materials while taking into account preferences of a group of stakeholders. Section 4 exhibits comprehensive results of multiple criteria assessment of the insulating materials. The last section concludes.

Review of Multiple Criteria Sorting Group Decision Methods
The objective of the case study presented in this paper is to give an easily interpretable comprehensive assessment of the insulating materials' sustainability. This is achieved by assigning them to a set of pre-defined and ordered decision classes based on their performances on multiple criteria (Kadziński et al. 2015b). While computing the sorting recommendation, we account for the preferences of a group of experts and stakeholders. This requires implementation of a group decision making framework. As real-world situations often involve multiple stakeholders, some methods have been proposed to support groups in making collective sorting decisions (Daher and Almeida 2010). These approaches can be distinguished at different levels. In particular, they differ in terms of a preference model employed to represent preferences of the DMs. Furthermore, an underlying classification rule may involve analysis of a single preference model instance or all sets of parameters compatible with the DMs' preference information. Moreover, sorting methods can be divided with respect to the level on which individual viewpoints are aggregated (Dias and Climaco 2000). Finally, some approaches account for the importance degrees of the involved DMs, while other methods assume that all DMs play the same role in the committee.
Among multiple criteria sorting group decision methods, outranking-based approaches are prevailing. Most decision support systems in this stream incorporate Electre TRI-B (Yu 1992;Roy 1996). For example, Dias and Climaco (2000) proposed an approach that admits each DM to specify imprecise constraints on the parameters of an outranking model, then exploits a set of compatible parameters using robust assignment rule, and finally aggregates individual perspectives in a disjunctive or conjunctive manner (thus, not accounting for the DMs' powers). The former accepts an assignment if it is justified by at least one DM, whereas the latter confirms some classification only if it is consistent with the preferences of all DMs. In this way, a group may agree on some result even if its members do not share the same model parameters. This idea was extended by Damart et al. (2007) to an interactive preference disaggregation approach that accepts assignment examples provided by different DMs. The method incorporates robustness analysis by deriving for each DM the possible class assignments (confirmed by at least one compatible preference model instance) and guides the group on sorting exemplary alternatives by exhibiting the levels of consensus between the DMs. Analogously, Shen et al. (2016) developed an adaptive approach under intuitionistic fuzzy environment that allows to reach a classification with an acceptable individual and group consensus levels. Moreover, de Morais Bezerra et al. (2017) enriched Electre TRI-B with the tools for visualizing the comparison of individual results and procedures for guiding the changes of model parameters for deriving a better consensus.
Furthermore, Jabeur and Martel (2007) proposed a framework, which derives a collective sorting decision at the output level from the individual non-robust classifications by additionally accounting for the relative importance of group members. Then, Morais et al. (2014) used a stochastic variant of Electre TRI-B, called SMAA-TRI, to consider uncertainty in criteria weights and to derive for each DM the shares of the relevant parameter vectors that assign a given alternative to a certain category. An overview of thus obtained individual results leads to a collective recommendation. Conversely, Cailloux et al. (2012) employed assignment examples provided by multiple DMs for reaching an agreement at the input level. In particular, they proposed some linear programming models for deriving a joint set of boundary class profiles and veto thresholds.
As far as outranking-based sorting approaches incorporating a model typical for PROMETHEE are concerned, Nemery (2008) extended the FlowSort method to group decision making. His proposal derives an assignment for each alternative from its relative comparison (strength and weakness) against the boundary or central class profiles specified by each individual DM. A similar idea was implemented by Lolli et al. (2015) in FlowSort-GDSS. The underlying procedure derives class assignments by comparing comprehensive (global) net flows of alternatives and reference profiles. The proposed sorting rules distinguish between scenarios in which analysis of the individual assignments leads to either univocal or non-unanimous recommendation. Although the viewpoints of different DMs are aggregated at the output level, the method defines some consistency conditions on the preference information (in particular, reference profiles) provided by the individual DMs.
The majority of existing value-based approaches derive a sorting recommendation incorporating robustness analysis and not differentiating between the roles played by the DMs. In particular, the UTADIS GMS -GROUP method (Greco et al. 2012) accounts for the assignment examples provided by each DM and derives collective results that concern two levels of certainty. The first level refers to the necessary and possible consequences of individual preference information, which is typical for Robust Ordinal Regression (ROR) (Greco et al. 2010;Kadziński et al. 2015b). The other level is related to the necessity or possibility of a support that a particular assignment is given in the set of DMs. This method was further adapted by Liu et al. (2015) to account for the uncertain evaluations represented with the evidential reasoning approach, to provide some measures on the agreement between the DMs, and to derive a collective univocal assignment.
Conversely,  aimed at a joint representation of assignment examples provided by all DMs by a set of additive value functions and investigating the necessary and possible consequences of applying the latter on the set of alternatives. When there is no value function compatible with preferences of all DMs, some linear programming techniques can be used to remove a minimal subset of inconsistent assignment examples. A similar approach was proposed by Cai et al. (2012), though additionally accounting for the DMs' priorities. The latter ones intervene in the selection of a representative value function and in resolving inconsistency in the provided assignment examples. These priorities are updated with the progressive preference elicitation process to reflect the preciseness, quantity and consistency of the example decisions supplied by each DM.
Finally, when it comes to using "if …then …" decision rules for representing preferences of the DMs, one proposed various extensions of the Dominance-based Rough Set Approach (DRSA) (Greco et al. 2001). These accept preference information in form of individual assignment examples. First, Greco et al. (2006) introduced some concepts (e.g., multi-union and mega-union) related to dominance with respect to minimal profiles of evaluations provided by different DMs. Then, Chen et al. (2012) proposed to aggregate the recommendations suggested by individual linguistic decision rules into an overall assignment be means of a Dempster-Shafer Theory. The crucial concepts incorporated in the DRSA sorting method proposed by Sun and Ma (2015) are a dominance relation on the set of multiple sorting decisions (each provided by an individual DM) and a multi-agent conflict analysis framework. Furthermore, Chakhar and Saad (2012) and Chakhar et al. (2016) illustrated how to combine individual approximations of class unions and derive collective decision rules that permit classification of all alternatives in a way consistent with the judgments of all DMs. These approaches measure the contribution of each expert to the collective assignment in terms of the individual quality of classification. Finally, Kadziński et al. (2016) adapted the principle of ROR to a group decision framework with DRSA, thus considering all sets of rules compatible with the individual assignment examples and combining their indications only at the output level.
In this paper, we propose an outranking-based group decision approach that incorporates Electre TRI-rC. Thus, it derives the assignments by comparing alternatives with the characteristic class profiles rather than with the boundary profiles as in Electre TRI-B. The basic procedure we use takes into account a single preference model instance (incorporating criteria weights derived from the SRF procedure) for each DM and aggregates the individual viewpoints at the output level. While still aggregating the preferences at the output level, we extend the basic framework to offer results of robustness analysis with multiple sets of parameters compatible with the DMs' value systems. Additionally, we propose a new algorithm for constructing a group compromise ranking of criteria, hence offering aggregation of the individual viewpoints also at the input level. At all stages, we assume that the involved stakeholders have the same importance degrees. Moreover, instead of providing precise assignments, our framework offers acceptability indices indicating the support that is given to the assignment of each alternative to various classes by different DMs and/or preference model instances compatible with their preferences.

Multiple Criteria Decision Analysis Method for the Assessment of Insulating Materials
This section describes a three-stage multiple criteria decision analysis method that has been used to evaluate the insulating materials while taking into account preferences of a group of stakeholders. Firstly, we discuss the Electre TRI-rC method (Kadziński et al. 2015b) that has been employed to assign the materials to a set of pre-defined and ordered classes. It incorporates the SRF procedure to compute the criteria weights (Figueira and Roy 2002). The method has been extended to a group decision setting to derive for each material some group class acceptability indices, which indicate the proportion of stakeholders that accept an assignment of the material to a given class. Secondly, we have adapted Stochastic Multi-criteria Acceptability Analysis (SMAA; Lahdelma and Salminen 2001;Tervonen and Figueira 2008;Tervonen et al. 2007) to the context of Electre TRI-rC and SRF procedure. It has been used to conduct robustness analysis (Roy 2010) for the results obtained in the first part, i.e., to validate their certainty while avoiding the arbitrary choice of criteria weights, which is conducted by the SRF procedure. Thirdly, we have proposed an algorithm for constructing a group compromise ranking of criteria based on the orders provided by the individual DMs. This ranking of criteria has been used as an input for SMAA to offer yet another view on the stability of computed results. Let us use the following notation (Kadziński et al. 2015a): • A = {a 1 , a 2 , . . . , a n } is a set of alternatives (insulating materials); . . , g m } is a family of evaluation criteria that represent relevant points of view on the quality of assessed alternatives; • g j (a) is the performance of alternative a with respect to criterion g j , j = 1, . . . , m (when presenting the method, without loss of generality, we assume that all criteria are of gain type, i.e., the greater the performance, the better); • C 1 , C 2 , . . . , C p are the preference ordered classes to which alternatives should be assigned; we assume that C h is preferred to C h−1 for h = 2, . . . , p.

Assessment of Insulating Materials Within a Group Decision Framework Incorporating Electre TRI-rC and the SRF Procedure
In this section, we present the Electre TRI-rC method (Kadziński et al. 2015b) that is used to assign the materials to a set of pre-defined and ordered classes. The method derives for each material a possibly imprecise assignment by constructing and exploiting an outranking relation S (Figueira et al. 2013). This relation quantifies an outcome of the comparison between the materials and a set of characteristic class profiles (Rezaei et al. 2017). In what follows, we discuss the main steps of the incorporated approach.
Step 1 For each class C h , provide the most typical (representative) performances on all criteria g j , j = 1, . . . , m, thus specifying the characteristic profiles b h , h = 1, . . . , p (Almeida Dias et al. 2010). Defining such profiles was found intuitive and manageable by the involved experts, which was the main reason for incorporating Electre TRI-rC in the study. The set of all characteristic profiles is denoted by B.
Steps 2-7 are conducted separately for each Decision Maker ( Step 2 Determine the weight w k j of each criterion g j , j = 1, . . . , m, using the SRF procedure (Figueira and Roy 2002). This method expects DM k to: • Assign some importance rank l k ( j) to each criterion g j ; this is attained by ordering the cards with criteria names from the least to the most important (the greater l k ( j), the greater w k j ; some criteria can be assigned the same rank, thus being judged indifferent); • Quantify a difference between importance coefficients of the successive groups of criteria judged as indifferent, L k s and L k s+1 , by inserting e k s white (empty) cards between them (the greater e k s , the greater the difference between the weights assigned to the criteria contained in L k s+1 and L k s ); • Specify ratio Z k between the importances of the most and the least significant criteria denoted by L k v(k) and L k 1 .
These inputs are used to derive the criteria weights as follows (Figueira and Roy 2002;Corrente et al. 2016): Steps 3-6 are conducted for each pair consisting of alternative a and profile b h .
Step 3 For each criterion g j compute a marginal concordance index c k j (a, b h ) defined as follows: The index quantifies a degree to which a is at least as good as b h on g j . Let us remark that in our study the experts defined the performances of characteristic profiles on all criteria by selecting them from the performances of the considered materials. This facilitated the preference elicitation process when dealing with a set of criteria with heterogeneous performance scales. In this perspective, when comparing the alternatives with the characteristic class profiles, we decided to exploit only the ordinal character of criteria and not use the discrimination (indifference and preference) thresholds, which can be, in general, employed in Electre. That is, in our application, the outranking of alternative a over profile b h on g j means that g j (a) is at least as good as the most typical (representative) performance for class C h on g j of some considered material.
Step 4 Compute a comprehensive concordance index σ k (a, b h ) defined in the following way: The index quantifies a joint strength of a subset of criteria supporting the hypothesis about a outranking b h (aS k b h ). Note that in our study, no criterion was judged strong enough to be attributed a power to veto against the outranking relation. Thus, no discordance effect has been considered.
Step 5 Specify the cutting level λ k (also called majority threshold), and compare σ k (a, b h ) with λ k to verify the truth of a crisp outranking relation aS k b h in the following way: The truth of relation b h S k a can be verified analogously.
Step 6 Use information on the truth or falsity of aS k b h and b h S k a to check the validity of: Step 7 For alternative a determine its desired class interval C k (a) = C k L (a) , C k R (a) by applying the assignment rules of ELECTRE TRI-rC (Kadziński et al. 2015b). To compute the worst class C k L (a), compare a successively to b h , for h = p − 1, . . . , 1, seeking the first (i.e., the best) characteristic profile b h such that: and select C k L (a) = C h+1 . When no such a profile is found, C k L (a) = C 1 . To compute the best class C k R (a), compare a successively to b h , for h = 2, . . . , p, seeking the first (i.e., the worst) characteristic profile b h such that: In case no such a profile is found, C k R (a) = C p .
Step 8 Combine the individual class assignments for all DMs into group class acceptability indices E ∂ (a, h) (Damart et al. 2007;Kadziński et al. 2016). These are defined as the proportion of DMs (stakeholders) that accept an assignment of alternative a to class C h , i.e.: where for k = 1, . . . , K : This measure indicates a cumulative support given to the assignment of a to C h by all group members.

Stochastic Multi-criteria Acceptability Analysis with Electre TRI-rC
The SRF procedure derives the precise weight values from the ranking of criteria, intensities of preference, and ratio between the most and the least important criteria provided by DM k applying some arbitrary rule (Figueira and Roy 2002). However, there exist multiple weight vectors compatible with such incomplete preference information. Recently, many researchers have raised the robustness concern in view of the SRF procedure to quantify the impact of uncertainty in the selection of an arbitrary weight vector on the stability of computed recommendation. In particular, Siskos and Tsotsolas (2015) proposed a set of robustness rules for the SRF procedure to obtain tangible and adequately supported results. Then, Govindan et al. (2017) suggested to exploit the whole set of compatible weight vectors to construct the necessary and possible results being confirmed by, respectively, all or at least one compatible vector. Further, Corrente et al. (2017) adapted the stochastic analysis of recommendation with the SRF procedure to the context of Electre III. We follow the latter research direction and integrate Stochastic Multi-criteria Acceptability Analysis (Lahdelma and Salminen 2001;Tervonen et al. 2007) to handle possibly imprecise weight values compatible with the ranking of criteria and to derive robust recommendation with Electre TRI-rC. SMAA applies the Monte Carlo simulation to provide each DM with the acceptability indices which measure the variety of different preferences (in particular, weight vectors) that confirm the validity of particular elements of the recommendation. In our case, the space w k (S R F) of weight vectors compatible with preferences of DM k is defined by the following constraint set E k (S R F) : where the interpretation of different constraints is as follows: • [O1] ensures that criteria ranked better by DM k will be assigned greater weight; • [O2] guarantees that criteria deemed indifferent by DM k will be assigned equal weights; • [O3] sets the ratio Z between weights of the most and the least significant criteria; • [O4] respects the intensities of preference for different pairs of criteria that have been quantified with the number of inserted empty cards; • [O5] normalizes the weights.
These constraints also ensure that all weights are positive. For each DM k , each weight vector w ∈ w k (S R F) and each alternative a ∈ A, we compute the resulting class assignment C k w (a) = C k w,L (a) , C k w,R (a) with Electre TRI-rC. We define the class range stochastic acceptability index C RS AI k (a, [L , R]) ) on a range of classes C k L (a) , . . . , C k R (a) with L ≤ R as the proportion of compatible weights w ∈ w k (S R F) that assign alternative a precisely to the range of classes C k L (a) , . . . , C k R (a) . Formally, the index is computed as follows: i.e., the proportion of weights that either precisely or imprecisely assign a to C h (Kadziński and Tervonen 2013;Kadziński et al. 2014). Let us define such a cumulative class stochastic acceptability index CuC S AI k (a, h) as: We estimate C RS AI s with acceptable error bounds by sampling the space w k (S R F) with the Hit-And-Run (HAR) algorithm . Overall, C RS AI k (a, [h L , h R ]) and CuC S AI k (a, h) can be interpreted as a support given by DM k to the assignment of a to, respectively, To measure a cumulative support given to the assignment of a to C h by all DMs in ∂ K , we consider a cumulative group class stochastic acceptability index CuC S AI ∂ K (a, h), defined as follows (Kadziński et al. 2016(Kadziński et al. , 2018:

Selection of a Group Compromise Ranking of Criteria
In this section, we introduce a procedure for deriving a compromise complete ranking of criteria based on the rankings provided individually by each DM k within the SRF procedure. The procedure builds on the algorithm that was introduced by Govindan et al. (2017) for constructing a utilitarian ranking of alternatives. Hence, we adopt an idea of minimizing a sum of of distances between the compromise ranking and all individual rankings. When considering a complete ranking of criteria for DM k , for each pair (g j , g l ) one of the three relations holds: g j is preferred to g l (g j k g l ), or g j is indifferent with g l (g j ∼ k g l ), or g l is preferred to g j (g j ≺ k g l ). Let R jl k and R jl k denote the relations holding between g j and g l in the rankings provided by, respectively, DM k and DM k (e.g., R  Table 1 (for a detailed justification of these values, see Roy and Słowiński 1993). A distance between two rankings of criteria provided by DM k and DM k involving all ordered pairs of criteria (g j , g l ) is defined as follows: In what follows, we present a Binary Linear Program (BLP) for constructing a compromise ranking of criteria for group ∂ K involving K DMs. Following Govindan et al. (2017), for each pair of criteria (g j , g l ), we introduce two binary variables p jl ∂ and i jl ∂ (see constraint [R1] in E ∂ (S F R)) with the following interpretation: • p jl ∂ represents a weak preference of g j over g l in the compromise ranking (i.e., in case p jl ∂ = 1, then g j ∂ g l or g j ∼ ∂ g l ); note that p jl ∂ and p l j ∂ can be used to instantiate one of the three relations Since we impose completeness and transitivity on a weak preference relation, we require that p jl ∂ = 1 or p l j ∂ = 1 (see [R2]) and that p jr ∂ = 1 and p rl ∂ = 1 imply p jl ∂ = 1 (see [R4]). When constructing a utilitarian complete ranking of criteria, we aim at minimizing a comprehensive distance between relations ( ∂ , ≺ ∂ , or ∼ ∂ ) instantiated for all pairs of criteria in the compromise ranking and relations observed for these pairs in the individual DMs' rankings (for DM k , the relation between g j and g l ( j < l) is denoted by R jl k ): ∂ to a value of the objective function (for a detailed explanation, see Govindan et al. 2017).
Once a group compromise ranking of criteria is constructed, we conduct robustness analysis with SMAA in the same way as described in the previous section for an individual DM. This leads us to deriving cumulative group compromise class stochastic acceptability indices CuCC S AI ∂ K (a, h).

Decision Aiding with the Proposed Approach
Multiple criteria sorting decisions can be aided with the proposed group decision making framework through the process illustrated in Fig. 1. It starts with specifying the sets of alternatives, criteria, and ordered classes as well as the alternatives' evaluations (performances) on the criteria.
Then, the preference information is elicited from the involved experts and/or stakeholders. Each stakeholder is required to provide a cutting level as well as a ranking of criteria that incorporates the intensities of preference and the ratio between the importance coefficients of the most and the least significant criteria, as required by the SRF procedure. Moreover, the experts are expected to define a characteristic profile for each class. In our study, the profiles were agreed by multiple experts, but, in general, the methodological framework admits that each stakeholder provides his/her individual set of profiles.
Further, the method derives three types of results. These indicate a support that is given to the assignment of considered alternatives to different classes via the application of Electre TRI-rC for different sets of weights and cutting levels compatible with the preferences of the involved experts. In two cases, the preferences of the individual stakeholders are aggregated only at the output level. Depending on whether these individual preferences are processed using the SRF procedure or the Monte Carlo simulation, the method computes, respectively, group class acceptability indices or

Fig. 1
Decision aiding process with the proposed group decision methodological framework cumulative group class stochastic acceptability indices. In the third case, the preferences are aggregated at the input level by constructing a group compromise ranking of criteria. Then, the method applies SMAA to derive cumulative group compromise class stochastic acceptability indices. Finally, these three types of outcomes should be analyzed and combined into the recommended assignments. This is straightforward in case the support given to the assignment of alternatives to decision classes by different results is similar. In case of ambiguous indications by different procedures, the inconsistency needs to be raised by a decision analyst.
Obviously, it is not required to use all three types of procedures and respective results for each study. This may be useful when offering different viewpoints on the robustness of sorting recommendation is desired. Otherwise, one can employ just a single procedure for processing the experts' preferences depending on whether they should be aggregated at the input or output level and whether the robustness analysis should be incorporated into a particular study.

Results of Multiple Criteria Assessment of Insulating Materials with the Outranking Preference Model and Characteristic Class Profiles
The study aims at evaluating overall sustainability of coating materials used in buildings retrofitting. We consider 13 materials listed in Table 2 (they are denoted by A = {a 1 , a 2 , . . . , a 13 }). All materials having a thickness of 15cm were placed internally on the roof of a model building typical for central Italy, and evaluated from the socio-economic and environmental viewpoints. The six relevant criteria which have been used to assess the materials are: hour of discomfort (g 1 ; DH),CO 2 avoidance (g 2 ); Net Present Value (g 3 ; NPV), human health (g 4 ); ecosystem quality (g 5 ), and consumed resources (g 6 ).
In what follows, we explain their meaning. Discomfort degree Hour (g 1 ; the less, the better) evaluates a thermal performance of a building on an annual basis (CEN 2007) in accordance with the EN 15251 standard. Thus defined, it serves as a measure of comfort. The performance on g 1 is quantified as an overall time during which the temperature falls outside the second comfort category that was considered in the study (Carlucci and Pagliano 2012), and then weighing it by how much the limit has been exceeded. For this purpose, we have used the following equation: where CC 2 is the lower or upper limit of the assumed comfort category, OT i is the operative temperature at hour ì, and the multiplier 10 60 refers to an employed time step of 10 minutes. CO 2 avoidance (g 2 ; the more, the better) measures the energy saved during the building life by using a particular insulating material when compared to the case of no insulation in the following way: where E S is the estimated Energy Saved in GJ at time t with a time horizon of 25 years, 277.78 is a conversion factor to GJ in kWH, while 406.31 is the conversion factor for Italy from kWH to kg of CO 2 per year (EIA, 2015). Therefore, the CO 2 avoided refers only to the use phase, which is not considered in the LCA study. Net present value (g 3 ; the more, the better) is the difference between the present values of cash outflows and inflows. On one hand, the outflows involve Primary Energy Input (P E I ) cost, installation cost I at time t = 0, and the dismissing cost E L T after the lifespan T of the investment (25 years). On the other hand, the inflows refer to the Cost of Energy Saved E S t in different time periods t. Overall, we have computed N PV as follows: where i is the discount rate. For a detailed justification of this measure, see Menconi and Grohmann (2014). Thus defined, N PV can be seen as an outcome of Life Cycle Costing, which is an economic methodology for assessing the profitability of using different alternatives by taking into account the costs they incur at different stages of a life cycle (e.g., construction, operations, and maintenance).
For the assessment of environmental impacts, we have used the Eco-indicator 99 method (Goedkoop and Spriensma 2001) implemented in the SimaPro software (Product Ecology Consultants 1990). The method aggregates the results of Life Cycle Assessment into a set of parameters that can be interpreted as damage categories. In general, LCA is useful for identifying the environmental implications of a given alternative through the quantification of consumed resources (e.g., energy, raw materials, water) and related emissions (e.g., emissions into the air, water and soil, waste and co-products) (Paolotti et al. 2017). We used the following three environmental Eco-indicators expressed on a dedicated point scale: • Human health (g 4 ; the less, the better) which is derived from the analysis of the following normalized impact categories: carcinogens, respiratory organics and inorganics, climate change, radiation, and ozone layer; • Ecosystem quality (g 5 ; the less, the better) which is made up by the following three normalized impact categories: ecotoxicity, acidification/eutrophication, and land use; • Resources (g 6 ; the less, the better) which aggregates two normalized impact categories: minerals and fossil fuels.
The LCA focused on the production phase, starting from the production of a raw material to the obtaining of its complete version. We omitted the use and disposal phases, hence implementing an LCA "from cradle to gate" (Paolotti et al. 2016). All the impacts were calculated considering a functional unit of 1 m 3 of insulating material. The performances of 13 insulating materials with respect to 6 criteria are provided in Table 2. For all materials but hemp fibres, Ecoinvent Database (Ecoinvent 2010) was used as a source of foreground and background data related to both production and assembly processes as well as to the transport, electricity and fuel consumption. Instead, for the hemp processes the underlying data was derived from Zampori et al. (2013).
The objective of the case study is to give an easily interpretable comprehensive assessment of the materials' sustainability. This is achieved by assigning them to a set of three pre-defined and ordered classes: C 1 (low sustainability), C 2 (medium sustainability), and C 3 (high sustainability).
The study involved elicitation of preferences from the two groups of stakeholders. On one hand, a characteristic profile b h for each class C h , h = 1, 2, 3, has been collectively specified by the experts from the university-based engineering team specialized in the materials and retrofitting of rural buildings. On the other hand, the preferences on the importance of individual criteria have been elicited individually from multiple stakeholders who were owners of rural buildings interested in a renovation of their houses for improving the energetic performance. Thus, they can be perceived as potential consumers of the insulating materials.
When it comes to the characteristic profiles, the experts decided to define them by indicating one of the performances observed in the set of materials. The consensus between the experts on the most typical performance levels for each class has been reached during an interactive focus group. These levels are summarized in Table 3.

Results of Multiple Criteria Assessment of the Insulating Materials Within a Group Decision Framework Incorporating Electre TRI-rC and the SRF Procedure
The weights representing the importance of individual criteria have been elicited from the rural buildings' owner. In what follows, we call them stakeholders. Overall, we approached 63 owners by explaining them the characteristics of different materials, the interpretation of all criteria and their relation to different phases of the materials' life cycle. Among them, 38 stakeholders (let us denote them by ∂ K = {DM 1 , DM 2 , . . . , DM 38 }) claimed to understand the meaning and role of different criteria, and expressing their willingness to provide preferences on the criteria importance.
In Table 4, we present the incomplete preference information required by the SRF procedure, which was provided by three selected stakeholders. We also report the computed weights w k j and cutting level λ k . All stakeholders agreed that λ k should be equal to the sum of weights of the three most important criteria. The complete data for all group members is provided in the supplementary material available as an e-Appendix (the same remark applies to the results discussed in the following sections).
The results of a comprehensive comparison between 13 materials and 3 characteristic profiles are quantified with the comprehensive concordance indices. In Table 5, we present such indices for four exemplary materials for DM 1 . Table 5 exhibits also the justification of delivered assignment for the exemplary materials. For instance, a precise assignment of a 6 to C 1 can be explained with b 2 being preferred to a 6 and there existing sufficiently strong support in favor of b 1 outranking a 6 (σ 1 (a 6 , b 2 ) = 0.000 < σ 1 (b 1 , a 6 ) = 0.524).
In Table 6, we report the assignments obtained for all materials for different DMs. In particular, for DM 1 there are 6 materials assigned to the best class ( a 5 , a 7 , a 8 , a 11 , a 12 , a 13 ), 3 materials whose quality is evaluated as medium (a 1 , a 2 , a 3 ), and 4 materials judged as bad (a 4 , a 6 , a 9 , a 10 ). The assignments for DM 5 are the same except for a 4 being imprecisely assigned to [C 1 , C 2 ].
The spaces of consensus and disagreement with respect to the assignments obtained for all DMs are quantified with the group class acceptability indices E ∂ (a, h) (see Table 7). For example, for a 1 none stakeholder confirmed its assignment to the worst class C 1 , 36 out of 38 stakeholders supported its assignment to the medium class C 2 , and 3 stakeholders suggested the assignment of a 1 to the best class C 3 . These numbers have been translated to the following group acceptability indices: E ∂ (a 1 , 1) = 0 38 = 0, E ∂ (a 1 , 2) = 36 38 = 0.95, and E ∂ (a 1 , 3) = 3 38 = 0.08. On the contrary, for a 2 all stakeholders agreed with respect to its assignment to C 2 (E ∂ (a 2 , 2) = 38 38 = 1.0), while the results obtained for 6 of them additionally indicated hesitation in terms of its assignment to C 1 (E ∂ (a 2 , 1) = 6 38 = 0.16). The analysis of E ∂ (a, h) leads to indicating the assignments which are necessary in terms of the support they are provided in the group of stakeholders. Additionally, these results clearly indicate the most and the least probable assignments. In particular, for each material we are able to indicate the class with the greatest support among all stakeholders. It is C 1 for a 6 , a 9 and a 10 , C 2 for a 1 , a 2 , a 3 and a 4 , or C 3 for a 5 , a 7 , a 8 ,  a 11 , a 12 , and a 13 . The support which is given to the assignment of the materials to other classes is significantly smaller. For clarity of presentation, in all tables exhibiting stochastic acceptability indices (Tables 7,8,9 and 11), the text in bold indicates the class with the greatest support for a given material.

Table 4
The order of cards with criteria names (ranks l k ( j)), white cards e k s , and ratio Z k provided by the three selected DMs in the SRF procedure, the weights w k j derived from the SRF procedure, and the cutting level  Table 5 Credibility indices and class assignments obtained with ELECTRE TRI-rC for four exemplary materials for DM 1 (cutting level λ 1 = 0.714) a 1 , b h ) 1.000 0.799 0.238 σ 1 (a 6 , b h ) 0.585 0.000 0.000 a 6 ) 0.524 1.000 1.000 σ 1 (a 11 , b h ) 1.000 1.000 0.476 σ 1 (a 12 , b h ) 1.000 1.000 0.524 a 11 ) 0.000 0.000 0.524 σ 1 (b h , a 12 ) 0.000 0.476 0.738

Results of Stochastic Multi-criteria Acceptability Analysis with Electre TRI-rC
To validate the recommendation for insulating materials against the arbitrary choice of weights conducted with the SRF procedure, we applied SMAA. For each stakeholder, we considered a sample of 10000 uniformly distributed weight vectors compatible with the ranking of criteria (s)he provided within the SRF procedure.   (2) 19 (16) 28 (1) 21 (10) 6 (31) vectors compatible with preferences of all stakeholders. Thus, such a recommendation needs to be treated with certainty. Conversely, CuC S AI ∂ K (a, h) = 0 (e.g., a 2 to C 3 , a 3 to C 1 , or a 9 to C 3 ) indicates the no classification model of any stakeholder confirmed the respective assignment. This makes it excluded from the potential recommendation.

Results of Stochastic Multi-criteria Acceptability Analysis for a Group Compromise Ranking of Criteria
The results presented in the previous sections were derived by aggregating the outcomes obtained individually for each stakeholder. In this section, we offer another perspective on the stability of results by searching for a compromise between different stakeholders already at the stage of provided preferences. In Table 10, we report the numbers of DMs indicating preference or indifference for all pairs of criteria in the ranking they provided for the purpose of applying the SRF procedure. For example, 14 out of 38 stakeholders preferred g 1 to g 2 , 22 stakeholders opted for an inverse preference, and only 2 stakeholders judged this pair indifferent. Conversely, when comparing g 5 to g 6 , 31 experts opted for an indifference, and only one claimed that g 5 was more important than g 6 . The information from the DMs' individual rankings has been used as an input for the algorithm constructing a compromise utilitarian ranking of criteria, i.e., the one which is on average the closest to 38 individual rankings. In this way, the following group compromise order of criteria has been constructed: Thus, the greatest importance has been attributed to ecosystem quality (g 5 ) and resources (g 6 ), while the least important criteria are NPV (g 3 ) and hour of discomfort (g 1 ). The relation instantiated for different pairs of criteria is consistent with the opin- Table 11 Cumulative group compromise class stochastic acceptability indices CuCC S AI ∂ K (a, h) for all materials h\a a 1 a 2 a 3 a 4 a 5 a 6 a 7 a 8 a 9 a 10 a 11 a 12 a 13 1 0.000 0.419 0.000 0.533 0.000 1.000 0.000 0.000 1.000 1.000 0.000 0.000 0.000 2 0.581 1.000 1.000 0.467 0.000 0.000 0.000 0.419 0.000 0.000 0.000 0.000 0.000 3 0.419 0.000 0.000 0.000 1.000 0.000 1.000 1.000 0.000 0.000 1.000 1.000 1.000 ion expressed by the significant number of stakeholders. For example, 24 stakeholders ranked g 5 and g 6 as the two most important criteria, while 19 of them ranked this pair tied for the first place. Furthermore, 21 stakeholders judged g 3 as the least important criterion.
Obviously, one needs to bear in mind that the compromise ranking of criteria minimizes the sum of distances between relations observed for all pairs of criteria in all individual rankings. In this perspective, it may not be considered representative by all individuals (see, e.g., DM 7 , DM 9 , DM 12 , DM 17 , DM 19 , DM 20 , or DM 36 ) whose preferences are represented in the compromise ranking to a marginal degree (i.e., an overall distance between their ranking and the compromise one is substantial).
Such a compromise ranking of criteria has been used to simulate DMs' joint preferences within SMAA. Consistently with the previous sections, the cutting level λ was assumed to be equal to the sum of weights of the three most significant criteria. The results of robustness analysis are materialized with the cumulative group compromise class stochastic acceptability indices CuCC S AI ∂ K (a, h) (see Table 11).
For most materials, the variability of results is lesser than in case of deriving the recommendation by aggregating the individual viewpoints. Indeed, for 11 out of 13 materials there is some class which is recommended with certainty (then, CuCC S AI ∂ K (a, h) = 1). Also, for all materials but a 4 the class assignments with the greatest support have not changed with respect to those reported in the previous sections. The main differences concern a lesser support for the assignment of a 1 , a 4 and a 9 to C 2 in favor of judging the quality of a 1 as high (C 3 ) and the quality of a 4 or a 9 as low (C 1 ). Finally, although the assignments of a 2 and a 8 to, respectively, C 2 and C 3 are robust, the acceptability for their assignment to some worse classes (C 1 and C 2 , respectively) has increased to 0.419.
The probability of other assignments was often non-negligible though significantly lower than for the above indicated classes. Nevertheless, the results obtained from the stochastic analysis allowed to nullify the risk of a false declaration that some material was assigned to a class which was not confirmed by any compatible set of weights for any expert. For each insulating material, the recommended decision can be justified by comparing its performances on different criteria with those of the characteristic class profiles. In Table 12, we indicate the subsets of criteria on which the materials outrank (i.e., are at least as good as) the characteristic profiles. In this regard, let us explicitly explain the most likely assignments suggested for some materials: • a 10 is worse than b 1 on all criteria, thus being assigned to the worst class C 1 ; in the same spirit, a 6 is worse than b 1 on g 2 , g 4 , and g 5 (thus, on 3 out of 4 considered environmental criteria), and not better than b 2 on any criterion, which makes C 1 its most desired class; • a 3 is better than b 1 and worse than b 3 on all criteria, which makes its performance vector typical for C 2 ; • a 12 and a 13 are at least as good as b 2 on all criteria and better than b 3 on four criteria (g 1 , g 2 , g 3 , g 5 or g 1 , g 2 , g 3 , g 6 , respectively (note that both scenarios include two accounted socio-economic criteria, g 1 and g 3 )), which makes their assignment to C 3 the most justified.

Conclusions
We considered a multiple criteria problem of sustainability assessment of insulating materials. We combined Life Cycle Costing, Life Cycle Assessment, and adaptive comfort evaluation to derive performances of these materials on six socio-economic and environmental criteria. The comprehensive assessment of the materials involved their assignment to three preference-ordered sustainability classes. The classification was performed with a group decision counterpart of the Electre TRI-rC method that compares alternatives with the characteristic class profiles defined by the experts.
To derive a recommendation that would reflect viewpoints of a wide spectrum of potential customers, we accounted for the preference information of a few tens of rural buildings' owners being interested in the roof's insulation. The initial recommendation was derived by computing the proportion of stakeholders who accepted an assignment of a particular material to a given class. These results were subsequently validated against the outcomes of a two-fold robustness analysis realized with the Monte Carlo simulation. The latter exploited the space of all criteria weights compatible with either each stakeholder's preference information provided in the SRF procedure or collective ranking of criteria that was derived with an original algorithm proposed in this paper.
The three-stage analysis revealed that the most sustainable materials were glass wool, hemp fibres, kenaf fibres, polystyrene foam, polyurethane, and rock wool. This was mainly due to their favorable performances quantified with the Net Present Value and Eco-indicators. On the contrary, gypsum fibreboard, mineralized wood and plywood were assessed as the least sustainable materials. This can be justified in terms of their poor performances on thermal comfort, human health, and ecosystem quality. Overall, the proposed method provided greater clarity for decision making and guaranteed credibility in the eyes of the traditional rural houses' owners. Moreover, all research results-concerning both materials' performances on the individual criteria and comprehensive sorting recommendation-were well perceived by the experts on insulating materials in Italy.
The proposed framework can be applied to other decision contexts than that of a typical farmhouse in central Italy. This would require, however, accounting for a comfort model as well as warm and cold periods suitable to a particular geographical context, specification of a relevant lifespan for the investment, and adapting life cycle assessment to the reality of a particular study.
From the methodological viewpoint, we envisage the following future developments. Firstly, we plan to extend the SRF procedure to a group decision context so that it tolerates intensities of preference for different pairs of criteria and accepts information on different roles (weights) of the decision makers. Secondly, we aim at extending the proposed group decision framework to methods dealing with choice and ranking problems. This would require elaboration of the algorithms for deriving a compromise recommendation that would appropriately combine results of robustness analysis computed individually for each stakeholder.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.