{'fraction_non_alphanumeric': 0.037306835460942746, 'fraction_numerical': 0.018741822553838536, 'mean_word_length': 5.0764933390631715, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 1, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 1, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Integrating fuzzy analytic hierarchy process with PROMETHEE method for total quality management consultant selection

Evaluation of the proper and appropriate consultants can play an important role in successful total quality management (TQM) program implementation and helps the manufacturing organizations to attain competitive advantage. In general, many conflicting factors affect the appropriate consultant selection problem which adheres to uncertain and imprecise data. In this paper, a simple, systematic and logical scientific approach is structured to evaluate TQM consultant through integrating Fuzzy Analytical Hierarchy Process with the Preference Ranking Organization Method for Enrichment Evaluations. The proposed decision-making approach takes advantage of the synergy between these two well-known multi-criteria decision-making methods in which their strong and weak points are detected and a ranking is provided which facilitates the final selection for the decision-maker. To accredit the proposed model, it is implemented in a furniture industry in Bangladesh. The results indicate that technical/administrative is the most significant criteria whereas work experience in related field is the most important sub-criteria.


Introduction
Due to rapid globalization, manufacturing organizations that can compete in today's competitive market will be those that can win in global competition of both domestic and international markets. Due to the dynamic business environment and changing customers' preferences, total quality management (TQM) becomes an essential management philosophy for manufacturing organizations to remain competitive (Khan, 2011). Today, a number of manufacturing organizations implement TQM to generate a competitive advantage (Shenawy, Baker, & Lemak, 2007), improve organizational effectiveness (Macinati, 2008), productivity (Kumar, Choisne, Grosbois, & Kumar, 2009), achieve customer satisfaction (Hur, 2009), improve the organization's overall business process (Zehir & Sadikoglu, 2010), generate improved products and services (Munhurrun & Munhurrun, 2008), and overall reputation (Jung, Wang, & Wu, 2009;Psychogios & Priporas, 2007).
On the contrary, Preference Ranking Organization METHod of Enrichment Evaluation (PROMETHEE), one of the most efficient outranking methods (Herva & Roca, 2013;Kalogeras, Baourakis, Zopounidis, & Dijk, 2005) is closely coinciding with the human perspectives and it determines the preferences among multiple decisions (Brans & Mareschal, 2005;Kabir, Sadiq, & Solomon Tesfamariam, 2013). For an integrated assessment analysis, PROMETHEE is a suitable approach as its flexible algorithm can enable tailor-made enhancements to meet specific requirements (Kabir et al., 2013). PROMETHEE also has some limitations like it does not yet provide any formal guidelines on how weights can be elicited (Gervasio & da Silva, 2012;Turcksin, Bernardinia, & Macharisa, 2011). It assumes that the decision-maker is able to weigh the criteria appropriately, no specific guidance is provided for its determination (Kabir & Sumi, 2014). Therefore, to eliminate these limitations and to develop an improved and systematic decision support tool for TQM consultant selection, the proposed decision-making approach takes advantage of the synergy between these two well-known multi-criteria methods where FAHP will apply to determine the weights of the criteria and sub-criteria, and PROMETHEE will be used for the aggregation of the criteria and ranking of the alternatives.
The remainder of this paper is organized as follows. In the next section, a brief literature review is presented. In the following section, the proposed methodology has been described with an overview and the concepts of the Fuzzy AHP, and PROMETHEE methods. In the next section, the proposed methodology is applied to evaluate and select the appropriate consultant for a furniture manufacturing company in Bangladesh. Finally, the last section presents the conclusion and discusses the limitations and scope for future research.
2. Literature review Different MCDM methods have been proposed for the selection or evaluation of consultant for the quality improvement. Cheung, Kuen, and Skitmore (2002) implemented analytic hierarchy process (AHP) for the architectural consultant selection for the design and construction of the projects in Hong Kong. Tsai, Shen, Lee, and Kuo (2007) also applied AHP to set priority weights for consultant alternatives in order to solve the enterprise resource planning (ERP) consultant selection problems. In another study, Tsai et al. (2009) applied ANOVA and regression analysis to examine the impact of ERP consultant selection factors on ERP project management and the impact of ERP project management on the performance improvement on ERP systems. Ling (2003) presented a conceptual model for the selection of architects by project managers using multi-attribute value technique (MAVT). Ng and Chow (2004) also used MAVT for engineering consultant's performance evaluation. The fuzzy gap analysis model is proposed to select and evaluate performance of engineering consultants by Chow and Ng (2007). Saremi et al. (2009) proposed nominal group technique (NGT) with technique for order preference by similarity to ideal solution (TOPSIS) under fuzzy environment for the consultant selection of an Iranian auto part manufacturer. NGT is used for selecting the criteria and fuzzy TOPSIS is used to rank the consultant firms considering qualitative and quantitative criteria. But the authors considered all the criteria as qualitative including implementation cost. Omar et al. (2009) presented design science research method (DSRM) and hierarchical fuzzy TOPSIS method for consultant selection for Malaysian Department of Irrigation and Drainage infrastructure projects. Alencar and Almeida (2010) applied PROMETHEE VI to select consultants for a construction project considering three subjective (culture, design, and quality) and three objective (cost, time, and experience) criteria. Chu (2011) suggested centroid ranking approach-based fuzzy MCDM method to evaluate consulting firm. The ranking approach of Euclidean distance based on centroid points was applied to determine the ranking order of the alternatives. El-Santawy and El-Dean (2012) combined multi-objective optimization on the basis of ratio analysis (MOORA) method with standard deviation weight method to select a consulting firm for a multi-national manufacturing company to determine the price for its new product. The authors considered three benefit criteria (company size, potential profit, and expected growth) and two cost criteria (cost of the initial consultation and the monthly current installments paid to the firm as its monthly fees). In another study, El-Santawy and Ahmed (2012) presented VIKOR to select the consulting firm considering three benefit criteria (company size, potential profit, and expected growth) and one cost criterion (cost of the consulting). In both studies, the authors did not consider any qualitative criteria. FAHP was applied to select the best Turkish quality consultant for the implementation of ISO 9000 projects considering 3 main attributes and 11 sub-attributes by Cebeci and Ruan (2007). But the authors considered only the qualitative criteria in their analysis. Moreover, due to limitation of Chang (1996) extent analysis, the weights of willingness to help clients, communication and interpersonal skills of the consultant, professional knowledge, expertise of consultants in client's sector, and problem-solving ability are considered zero. According to Wang, Luo, and Hua (2007), Chang (1996) extent analysis method might assign an irrational zero weight to some useful decision criteria and sub-criteria, leading to them not to be considered in decision analysis. In this study, modified Chang extent analysis is proposed to resolve the limitations of Chang extent analysis. Vayvay et al. (2012) compared AHP, Fuzzy AHP and analytic network process (ANP) for ERP consultant's selection. However, the authors considered transportation cost, consultancy cost, companies employed, and projects completed as qualitative criteria during their evaluation. Kabir (2014) proposed VIKOR method with fuzzy set theory for consultant selection for automotive battery manufacturing company. But implementation cost and experience in relevant field are also evaluated as qualitative criteria by the decision-makers (DMs).

An overview of the Fuzzy AHP and PROMETHEE approaches
In this study, FAHP will apply to structure the decision-making problem and to determine the weights of the criteria and sub-criteria, and PROMETHEE will be used for the aggregation of the criteria, ranking of the alternatives, and sensitivity analysis (Figure 1). FAHP and PROMETHEE are detailed theoretically in the following sub-sections.

Fuzzy AHP
Fuzziness and vagueness are common characteristics in many real-life decision-making problems. Therefore, an effective decision-making model needs to consider vagueness or ambiguity (Kabir & Hasin, 2012). The transformation of qualitative preferences to point estimates may not be sensible as DMs often provide uncertain answers rather than precise values (Kabir & Sumi, 2014). Conventional AHP that requires the selection of subjective values in pairwise comparison may not be sufficient and uncertainty should be considered in some or all pairwise comparison values (Hofmann & Knébel, 2013;Kabir & Hasin, 2012;Kahraman et al., 2003). The application of Saaty's AHP has some shortcomings as follows: (1) the AHP method is mainly used in nearly crisp decision applications, (2) the AHP method does not take into account the uncertainty associated with the mapping of one's judgment to a number, (3) the subjective judgment, selection, and preference of DMs have great influence on the AHP results, (4) ranking of the AHP method is rather imprecise, (5) the AHP method creates and deals with a very unbalanced scale of judgment. Therefore, conventional AHP seems inadequate to capture DMs requirements explicitly (Kabir & Hasin, 2012;Kabir & Sumi, 2014;Padhi & Mohapatra, 2009). FAHP should be more appropriate and effective than conventional AHP as it can take the pessimism/optimism rating attitude of DMs into account, linguistic values, whose membership functions are usually characterized by triangular or trapezoidal fuzzy numbers (Kabir & Hasin, 2012;Kutlu & Ekmekçioglu, 2012;Sen & Çinar, 2010).
Researchers proposed several methods like geometric mean method, direct fuzzification method, synthetic extent analysis, fuzzy least square method, fuzzy modification of the logarithmic least squares method, fuzzy preference programming, and two-stage logarithmic programming to attain the priorities in FAHP (Kabir & Hasin, 2012). Modified Chang's extent analysis is utilized in this research to evaluate the focusing problem. Chang (1996) develops a new approach for handling Fuzzy AHP, using triangular fuzzy numbers (TFN) for the pairwise comparison scale of Fuzzy AHP, and using the extent analysis method for the synthetic extent values of the pairwise comparisons. A TFN denoted asM = (l,m,u) where l < m < u, has the following triangular type membership function: In this methodology, decision-makers are asked to indicate the relative importance of two evaluation criteria in the same level. Linguistic pairwise comparisons are transformed to corresponding TFNs illustrated in Table 1.
Geometric mean operations have been used to combine pairwise comparisons of DMs or aggregating group decisions (Davies, 1994;Kahraman et al., 2003): It is necessary to check the consistency ratio of the comparison after aggregating group decisions. In this study, graded mean integration approach is utilized for defuzzifying the matrix (Kabir & Sumi, 2014). According to the graded mean integration approach, a fuzzy numberM = (l,m,u) can be transformed into a crisp number by employing Equation (3) (Kutlu & Ekmekçioglu, 2012): After the defuzzification of each value in the matrix, 'Consistency Ratio' (CR) of the matrix has been calculated and checked whether CR is smaller than .10 or not (Hofmann & Knébel, 2013;Kabir & Hasin, 2012;Kahraman et al., 2003). In this study, modified Chang extent analysis method is used to derive the importance weights of criteria from pairwise comparisons (Kabir & Sumi, 2014). Let W be the normalized weight vector of triangular fuzzy comparison matrix A, which includes the importance weights of criteria in the crisp form. The steps for calculating this vector are as follows.
The first step is to calculate the fuzzy synthetic extent value of each pairwise comparison. DenoteC ij as the TFN related to the pairwise comparison of criterion i over criterion j, which is represented as (l ij , m ij , u ij ). According to Chang (1996), the value of Table 1. Linguistic variables and fuzzy scales (Kabir & Hasin, 2012;Kabir & Sumi, 2014 fuzzy synthetic extent with respect to the criterion i, denoted asS i = (l i , m i , u i ), can be obtained via Equation (4): where n is the size of the fuzzy judgment matrix A.
To obtain P n j¼1C ij , fuzzy addition operation can be performed as and to obtain P n i¼1 P n j¼1C ij " # À1 , fuzzy addition operation ofC ij values can be performed as The inverse vector of the Equation (6) can be determined using Equation (7): Wang et al. (2007) corrected the normalization formula (Equation (7)) as Equation (8) X The second step is to derive fuzzy ranking value ofS i . In this step,S i is compared to other synthetic extent values of A,S i = (l j , m j , u j ). According to Chang (1996), the degree of possibility ofS i !S j is obtained by the following equation: where d is the ordinate of the highest intersection point between lS i and lS j . To compareS i andS j , both the values of V (S j !S i ) and V (S i !S j ) are needed.
The degree of possibility defined by the extent analysis method is an index for comparing two TFNs rather than an index for calculating their relative importance (Kabir & Hasin, 2012). Therefore, normalized degrees of possibility can only show to what degree a TFN is greater than all the others, but cannot be used to represent their relative importance (Kabir & Sumi, 2014;. To resolve this problem, Liou and Wang (1992) proposed total integral value with index of optimism, which derives the priorities of the synthetic extent values of A by the following equation: where α is the index of optimism which represents degree of optimism for decisionmakers. If α approaches 1 in [0, 1], the decision makers are more optimistic and otherwise they are more pessimistic (Kabir & Hasin, 2012;Kabir & Sumi, 2014;Sen & Çinar, 2010). Finally, the normalized importance weight vector W = (w 1 , w 2 , … , w n ) T of the fuzzy judgment matrix A is determined using following equation: where w i is a non-fuzzy number.

Theoretical basis
In PROMETHEE, the preference structure for each criterion is based on pairwise comparisons (Brans & Mareschal, 2005). The deviation between the evaluations of two alternatives on a particular criterion is considered. For larger deviations, the decision-maker allocates a larger preference to the best alternative or possibly no preference if the deviation is negligible. The larger the deviation, the larger the preference will be (Kabir & Sumi, 2014). These preferences may vary between 0 and 1.
Thus, function for each criterion will be: For which, 0 P j ðx; yÞ 1 Equation (12) indicates the preference of x over y for observed deviations between their evaluations on criterion g j (.), in case of a criterion to be maximized. For criteria to be minimized, the preference function should be reversed: A generalized criterion {g j (.), P j (x,y)} has to be defined for each criterion. Table 2 shows six types of particular preference functions which are used in this study (Brans & Mareschal, 2005;Gervasio & da Silva, 2012;Kabir & Sumi, 2014). However, other types of generalized criteria may also be considered.

PROMETHEE I methodology
For the application of PROMETHEE I, two major actions like the calculation of aggregated preference indices and outranking flows are needed. The aggregate preference indices π(x,y) expresses with which degree x is preferred to y over all the criteria, while π(y,x) indicates how y is preferred to x. Assuming x; y 2 A: pðx; yÞ ¼ P k j¼1 P j ðx; yÞw j ; pðy; xÞ ¼ P k j¼1 P j ðy; xÞw j ; π(x,y) ≈ 1, implies a strong global preference of x over y, while, π(x,y) ≈ 0, implies a weak global preference of x over y (Brans & Mareschal, 2005). Each alternative x faces (n − 1) other alternatives in A. Two outranking flows can be defined: (i) The positive outranking flow: expressing how an alternative x is outranking all the others. The higher [ þ ðxÞ, the better the alternative (Kabir & Sumi, 2014).
(ii) The negative outranking flow: Type III -Criterion with linear preference (V-shape) Type IV -Level criterion representing how an alternative x is outranked by all the others. The lower [ À ðaÞ, the better the alternative (Kabir & Sumi, 2014). The PROMETHEE I partial ranking is obtained from the positive and negative outranking flows: where PR I , ID I , and IC I represent the preference, indifference, and incomparability, respectively. In the case of xIC I y, a higher power of one alternative is associated to a lower weakness of the other. In this case, the information provided by both flows is not consistent and both alternatives are considered incomparable (Gervasio & da Silva, 2012;Kabir & Sumi, 2014).

PROMETHEE II methodology
PROMETHEE II shows a complete ranking of the alternatives. In this case, the net outranking flow is defined as the balance between the positive and the negative outranking flows: Type VI -Gaussian criterion In PROMETHEE II, all the alternatives are comparable; no incomparability's remain (Brans & Mareschal, 2005). When Ø(x) > Ø(y), x outranks y on all the criteria (Kabir & Sumi, 2014).

GAIA plane
Geometrical analysis for interactive aid (GAIA) provides a visual guidance for the principal criteria that are used for ranking of the alternatives based on the PROMETHEE II method (Kabir & Sumi, 2014). GAIA matrix is constructed from a decomposition of the net outranking flows Ø(x). The matrix data is then processed by a Principal component analysis (PCA) algorithm and displayed on GAIA biplot (Keller, Massart, & Brans, 1991). PCA is used to reduce the number of dimensions while minimizing the loss of information. This transformation of a multi-criteria problem to a two-dimensional space and geometrical representation of relations between criteria and alternatives provides a new perspective to the problem with the inevitable loss of some relation characteristics (Kabir & Sumi, 2014). It measures the amount of information preserved in the GAIA plane (values larger than 70% correspond to reliable GAIA planes) to control its quality (Herva & Roca, 2013;Kabir & Sumi, 2014;Visual PROMETHEE, 2013).

Application of the proposed methodology
The proposed method has been implemented on a furniture manufacturing company situated in the eastern part of Bangladesh and having more than one decade of successful operations. This study was unable to mention the real name of the organization under investigation as a condition of obtaining access for data collection. The company, which began manufacturing, retail, and export in the early 1990s, was one of the largest furniture manufacturing plants in Bangladesh. The area of the factory is more than 350,000 square feet and total number of employees is almost 2500. The company has more than 60 outlets and showrooms for sales. Currently, the company has joint venture with Bangladesh-Denmark and is exporting to Denmark, Netherland, India, Ukraine, Sweden, Germany, Switzerland, Singapore, Thailand, and Japan. The organization has decided to increase the overall production and customer satisfaction from 2012. For this, Board of Directors wants to apply TQM for getting better competitive advantages. That is why the company was willing to select appropriate TQM consultants.
To select the most appropriate TQM consultant, the company adopted the proposed method and set up a consultant selection decision panel. Twelve specialists with average work experiences of 13 years were chosen to form the decision team. Table 3 shows the profiles of the experts.
Initially, a literature review was conducted to identify various criteria and subcriteria considered by the researchers and decision-makers (Cebeci & Ruan, 2007;Cheung et al., 2002;Chow & Ng, 2007;Chu, 2011;Kabir, 2014;Ling, 2003;Ng & Chow, 2004;Omar et al., 2009;Saremi et al., 2009;Tsai et al., 2007Tsai et al., , 2009Vayvay et al., 2012). Thereafter, these criteria were placed in front of the expert panel to elicit the most appropriate and widely considered criteria from their experience and judgment for selection of TQM consultant. The expert panel selected six TQM consultants (further to be named as Con A to Con F) for further evaluation. The expert panel generate a decision hierarchy model which is shown in Figure 2.
After the decision hierarchy model was defined, the proposed approach in the earlier section was applied. To determine the importance weights of the criteria and sub-criteria, the expert panel were asked to make pairwise comparisons for the main and sub-criteria using fuzzy scales shown in Table 1. The geometric mean of the individual assignments has been determined to obtain overall or group results using Equation (2). Table 4 indicates the comparison matrix for the criteria with respect to the overall goal or objective.
The consistency ratio of the defuzzified comparison matrix (Table 4) is .032. As CR < .1, the level of inconsistency present in the information stored in comparison matrix is   The integral values of these values were calculated as below: I TA = .366; I ER = .345; I EC = .142; I KR = .171. After obtaining the total integral values using Equation (10), the importance weights of the main criteria were obtained via normalization (Equation (11)): W = (.357, .337, .139, .167).
Similarly, comparison matrix for the sub-criteria was determined and their weights were calculated, which are given in Table 5.
The input details (i.e. it has to be minimized or maximized, quantitative and qualitative data, and unit for each criterion) of six consultants are presented in Table 6. Among all the criteria or attributes, average project completion time (PC, months), work experience in related field (WE, number of years), and implementation cost (IC, Tk. × 1000) are considered as quantitative criteria and rest of the attributes are considered as qualitative criteria. For this purpose, a 7-point qualitative scale ranging from very low, low, below average, average, above average, high to very high has been applied in this study. Among them average project completion time and implementation cost are considered as non-beneficial attributes and lower values are desirable.
In this study, seven scenarios have been considered. For quantitative and qualitative assessments, the PROMETHEE guidelines advise to apply a linear preference function and Usual or Level preference function, respectively (Turcksin et al., 2011). Therefore, for the first scenario, Linear and Usual preference functions has been applied for quantitative and qualitative assessment, respectively. For the remaining six scenarios, Usual, U-shape, V-shape, Level, Linear, and Guassian preference function have been used for both qualitative and quantitative criteria. After the determination of the evaluation matrix and the preference functions, the scenarios are evaluated and ranked by means of the PROMETHEE decision-making software, Visual PROMETHEE (2013). The positive (Ø + ) and negative (Ø − ) flows (PROMETHEE I) and the net flow (Ø)  Figure 3 indicate that Con E is the right choice for the given application under the given conditions. The next choice is Con A and the last choice is Con F. Figure 4 indicates the GAIA plan of the decision problem, where nine-dimensional space of criteria is projected on a two-dimensional plane by means of PCA. The weighing of the criteria and the PROMETHEE II ranking are represented by the decision axis. Alternatives are represented by points and criteria by axes in this plan. Thus, alternatives scoring high on a particular criterion are represented by points located in the direction of the corresponding criterion axis (Brans & Mareschal, 1994). To differentiate the   directions (Kabir & Sumi, 2014). The quality level for visualization is 92.7%, which means that very few information got lost by the projection. Figure 4 indicates that technical skills, management skills, and responsible & trustworthy are close to each other. Similarly, work experience in related field, reputation, and average project completion time are close to each other. This means that consultants like Con E exhibiting a higher technical skills, management skills and responsible & trustworthy whereas Con A exhibiting a higher work experience in related field, reputation and lower average project completion time. It is thus possible to find consultants that are good on these criteria at the same time. Con F and Con B perform well with respect to communication & interpersonal skills but perform worse for rest of the criteria.
The sensitivity analysis of the proposed model is performed as PROMETHEE I and II rankings are greatly influenced by the weights allocated to the criteria. To modify the weights and to observe the resulting modifications of the PROMETHEE II ranking, a special feature of the Visual PROMETHEE software, called The Walking Weights is applied in this study ( Figure 5). The lower bar chart is showing the weights of the criteria, whereas the upper bar chart indicate PROMETHEE II complete ranking. The slider allows to change the weight of the selected criteria and to see the impact (Kabir & Sumi, 2014). For the following preference function, it can be easily observed that Con E still dominates the other ones followed by Con A.

Conclusions
Meditating entities such as consultants can play an important role for successfully implementing TQM program and help the manufacturing organizations to attain competitive advantage and improve organizational effectiveness. However, the TQM consultant selection problem is often influenced by uncertainty in practice and in real-life decision-making process as the decision-makers are unable (or unwilling) to express his/her preferences precisely in numerical values and the evaluations are very often expressed in linguistic terms.
To deal with this kind of situation, fuzzy set theory is an appropriate tool. In this study, FAHP has been proposed to determine the weights of the criteria and sub-criteria as PROMETHEE does not provide any formal guidelines for weighing up to now. On the other hand, PROMETHEE enriches FAHP by associating a preference function to each criterion. For this purpose, an integrated approach of FAHP and PROMETHEE has been proposed, in which the strengths of both methodologies are combined in a single MCDM tool to select the most appropriate TQM consultant.
The proposed decision support tool is very flexible and can handle any number of quantitative and qualitative attributes simultaneously. The proposed method not only enables us to determine the outranking order of the consultants, but also assess and rate them. It is also capable to identify groups of criteria expressing similar preferences and to better understand the conflicts that have to be solved in order to make a decision. The GAIA plane enables a graphical representation of the alternatives and criteria and helps to explore the strong and weak points of the alternatives. As research limitations, to get more effective decision, detail firm evaluation criteria, and sub-criteria can be incorporated. For further research, the results of this study can be integrated and/or compared with that of other MCDM methods like TOPSIS, ELECTRE, VIKOR, and MOORA. Other network-based techniques like ANP, Cognitive mapping can be also implemented to handle to the interdependencies between the criteria.