TY - JOUR AU1 - Ye, Qiwen AB - 1 Introduction The engineering [1, 2], economic and social systems [3, 4] are all complex and dynamic systems. Systematic elements interact and profoundly change the world via casual relationships, such as the effect of greenhouse gases on global climate; the dynamic process of resource allocation occurring in the transportation system [5]. Some of these casual relationships among elements are conceivably straight forward whereas the others are complicated to identify and may even impose perils on our survival. Facing such situations, decision-makers are growingly apprehensive that the traditional tools not only prove ineffective in addressing their enduring challenges but might, in fact, exacerbate them. Regrettably, earnest attempts to resolve pressing issues frequently give rise to unforeseen and adverse consequences. Therefore, the analysis of causal relationships among systematic elements raises an urgent challenge. The Decision-Making Trial and Evaluation Laboratory (DEMATEL) technique, originally developed by the Battelle Memorial Institute in collaboration with the Geneva Research Centre [6], is designed to elucidate intricate interconnections within dynamic systems [7]. DEMATEL leverages graph theory and matrix tools, representing a widely employed approach for comprehensive element assessments. By scrutinizing the logical interdependencies and direct impacts among system elements, DEMATEL effectively discerns the significance of these elements and classifies them as either effect factors or cause factors [8]. The significance of an element pertains to its role and the degree of control it necessitates. This categorization facilitates the analysis of intricate relationships between system elements. Essentially, DEMATEL excels in examining bidirectional and transitive causal associations among numerous system elements in complex, dynamic systems, employing causal and effect diagrams to visualize system causality. In recent decades, the DEMATEL method becomes prevailing in different areas such as operation management [9], innovation policy portfolio reconfiguration [10], human resource management [11], knowledge management [12], transportation management [13], and environmental economics decision [14, 15]. In addition to researching DEMATEL applications, scholars are continuously improving and expanding the method to increase its accuracy and applicability. This line of research can be categorized into two main streams. One stream involves combining DEMATEL with other decision-making methods, primarily to better adapt to specific decision-making scenarios and stratify different decision-making requirements. For example, Ortiz-Barrios et al. (2020), Karasan et al. (2022), WH El-Garaihy (2021) integrated the analytic hierarchy process (AHP) or data envelopment analysis (DEA) with DEMATEL [16–18]. The factor hierarchy analysis of AHP and DEA can be complemented with DEMATEL’s factor association analysis to comprehensively analyze the logical associations of decision-making factors. Similarly, research such as Chen (2021), Liang et al. (2022) combined interpretive structural modelling (ISM) with DEMATEL to enrich the analysis of factors’ relations [19, 20]. Furthermore, extending the DEMATEL method with the technique for order reference by similarity to ideal solution (TOPSIS) aids in optimizing sorting schemes [21, 22]. Scholars also extended DEMATEL with partial least squares (PLS) and Vise Kriterijumska Optimizacija Kompromisno Resenje (VIKOR), in order to reduce the difficulty of obtaining direct correlation matrices and solve decision-making problems with conflicting attributes [23–25]. The other research stream of extended DEMATEL focuses on the quantitative representation of expert preference information and is dominated by the gray DEMATEL and fuzzy DEMATEL. Grey DEMATEL combines gray system theory to extend the DEMATEL method for the incompleteness of decision-making information [14, 26]. Fuzzy DEMATEL, on the other hand, addresses the fuzziness and uncertainty of expert decision-making information, as decision makers prefer to apply linguistic terms to express their thoughts over crisp values due to fuzzy cognitions of human beings. In order to better handle the mentioned issue, different formats of fuzzy number are involved in previous scholarly efforts. Triangular fuzzy numbers, trapezoidal fuzzy numbers, Z-numbers, and 2-tuple linguistics are the most common fuzzy numbers integrated into DEMATEL, which can more scientifically portray the ambiguity and uncertainty of the experts in the DEMATEL judgment process and capture the fuzzy relationships between the factors. For example, Mahmoudi et al. (2019), Abdullah and Goh (2019), Çelik and Arslankaya (2023) tailored the fuzzy DEMATEL method with triangular fuzzy numbers for group decision making regarding the health care, waste management and quality management [27–29]. Compared with triangular fuzzy numbers, trapezoidal fuzzy numbers can be more flexible in dealing with complex fuzzy variables. Stephen and Felix (2023) applied DEMATEL in complex modifiable risk factors identification of cardiovascular disease with triangular fuzzy numbers [30]. Triangular fuzzy numbers are also employed in integrated multi-criteria decision-making methods, which involve complicated process [31, 32]. Similarly, Gaussian fuzzy numbers and Pythagorean fuzzy numbers are also commonly employed in extended DEMATEL methods [33–35]. The former accurately characterizes the distribution of fuzzy variables, thereby suitable for delineating continuous fuzzy variables. Conversely, the latter is capable of representing the fuzzy range of variables rather than fixed fuzzy values, thereby affording a degree of flexibility. However, it is unable to capture the distribution of variables within the fuzzy range. Additionally, more complex fuzzy numbers have been integrated into the DEMATEL method, such as the commonly encountered 2-tuple linguistics, Z-numbers, and T-spherical fuzzy sets. In comparison to traditional single-value fuzzy numbers, 2-tuple linguistics, through the combination of membership degree and non-membership degree, offers a more precise depiction of the degree and scope of fuzziness. Suo et al. (2019) introduced a DEMATEL approach that incorporates 2-tuple linguistics, utilizing the 2-tuple linguistic representation model, serving the purpose of simultaneously determining the significance and categorization of infrastructures risk factors [36]. In addition, Zhang et al. (2020) used 2-tuple linguistic data to create an interval-value intuitionistic fuzzy sets and applied the DEMATEL method to analyze the dimensions and standards of youth unemployment risk factors [37]. Z-numbers can comprehensively consider the deterministic components and uncertainty range of numerical values, but the calculation process is relatively complex and depends heavily on parameter settings. Wang et al. (2021) extended the DEMATEL method with Z-numbers and proposed the Z-numbers power weighted average operator for evaluating the human error probability [38]. In the same vein, Z-numbers are also integrated into DEMATEL for handling fuzzy information in different decision-making scenarios, such as the value propositions evaluation of smart product-service systems and the indicator identification of hospital performance [39, 40]. The membership degree of T-spherical fuzzy sets gradually decreases as the distance between elements and the center of the fuzzy set increases, allowing it to better describe the fuzzy or uncertain relationship between elements and the fuzzy set. However, this characteristic also limits its ability to describe fuzziness, which may lead to poor performance in addressing certain complex fuzzy problems. Therefore, T-spherical fuzzy sets usually applied in complex integrated DEMATEL methods for handling partial fuzzy information in the decision problems. For example, this fuzzy set is common in the integrated multi-criteria decision-making methods, such as the TOPSIS-DEMATEL methods proposed by Eti et al. (2023), Özdemirci et al. (2023) [41, 42]. In addition, T-spherical fuzzy sets also employed in the DEMAEL methods combining with other fuzzy numbers, such as 2-tuple linguistic [43], Pythagorean fuzzy numbers [44]. The aforementioned fuzzy numbers are typically applied in engineering, decision analysis, control systems, and similar domains, primarily derived using mathematical methods or statistical analyses. They are suitable for problems that can be described through quantifiable data or models, rather than relying on subjective judgments from experts. However, in practical situations, decision-makers usually come from different areas, bringing varying levels of knowledge and expertise in their respective domains [45–48]. Consequently, decision-makers may not have enough expertise to precisely articulate their judgment regarding the objects or candidates, resulting in fuzziness and uncertainty in the process of group decision-making. Essentially, the mentioned fuzzy numbers are primarily defined by their membership degrees, making it challenging to quantify non-membership and hesitation degrees. This inherent limitation leads to initial information loss and subsequent decrease in precision. Under this circumstance, it is advisable to represent decision-maker evaluations through the application of intuitionistic fuzzy numbers [49, 50]. This approach is believed to effectively address the limitations mentioned above. Thus, using intuitionistic fuzzy numbers to convey the judgment of decision-makers is more suitable than relying on exact numerical values or conventional fuzzy variables [51]. Emerging literature has largely focused on the intuitionistic fuzzy sets (IFSs) theory and applied it to such areas as cluster analysis [52], medical diagnosis [53], decision making [54, 55], and pattern recognition [56–58]. Notwithstanding the pervasive use, there is a lack of an appropriate method to tackle the dependent factor analysis problem using intuitionistic fuzzy information. In several practical problems of dependent factor analysis such as risk factor analysis, adoption factor analysis and pricing factor analysis, decision-makers may need additional comprehension to express their judgment containing the information affirmation, negation and hesitation. Given this situation, it is conceivable to notice that the intuitionistic fuzzy information is more effective and appropriate for addressing the problems identified. However, the intuitionistic fuzzy information is not involved in the process of dependent factors analysis in existing literature. In addition, the prior fuzzy DEMATEL methods are incapable of dealing with the decision problems with the intuitionistic fuzzy information. In an effort to advance this line of research, this paper proposes a novel method to more efficiently and effectively cope with the pressing problem. To demonstrate the proposed method, this paper applies the extended DEMATEL method in factor identification of electric vehicle (EV). EV is the innovative vehicle for transportation sustainability and smart logistics, and its decision-making is under an intuitionistic fuzzy environment with uncertain technical factors and fuzzy judgement information. The rest of this paper is structured as follows: Section 2 provides a concise review of the basic concepts and definitions of IFSs and the DEMATEL method. Subsequently, in Section 3, an expanded DEMATEL approach is introduced for analyzing the interrelationships among factors within the intuitionistic fuzzy context. The potential of the proposed method is illustrated through a real case in Section 4. Finally, Section 5 serves to summarize and highlight the characteristics and contribution of the proposed method. 2 Preliminaries In this section, the author provides the illustration of the fundamental concepts and definitions related to the DEMATEL method and the IFSs. 2.1 The intuitionistic fuzzy set The IFS proposed by Atanassov (1986) offers a viable approach to address vagueness. It is built on the foundation of classical fuzzy set theory [50]. Definition 1 Let set Y is fixed. An IFS A in Y is defined as (1) where uA(y),vA(y):Y→[0,1] are the membership function and non-membership function, respectively, with the condition 0≤uA(y)+vA(y)≤1. For each IFS in Y, let’s call (2) the intuitionistic index y of in A. It represents the hesitancy degree of y to A [59]. The author observed that (3) For each fuzzy set F in Y, there is . Hence, fuzzy sets are considered as the particular cases of IFSs. For convenience, let’s denote an intuitionistic fuzzy number (IFN) [50, 60] by a = (ua,va) where ua∈[0,1], va∈[0,1], ua+va≤1, and let Θ be the set of all IFNs. The score of S(a) can be assessed using the score function S shown as [61] (4) where S(a) = [0,1]. For any three IFNs a = (ua,va), and the following operational laws hold true. (1) , (2) , (3) , (4) . Definition 2 Let be a collection of IFNs, and let the intuitionistic fuzzy weighted average (IFWA) [62], IFWA:Θn→Θ, if (5) then IFWA is the intuitionistic fuzzy weighted averaging value, where δ = (δ1,δ2,…,δn)T is the weight vector of aj(j = 1,2,…,n) with δj∈[0,1] and . 2.2 The DEMATEL method This section will present the principal of the DEMATEL method. Suppose that the set of factors is F = {F1,F2,…,Fn} the procedure steps of the method are as follow [63–65]: Steps 1: Establish the initial direct-relation matrix. Let zij denotes the indicated degree of the dependence between factors Fi and Fj. Particularly, there does not exist dependence between Fi and itself. Then the initial direct-relation matrix Z = [zij]n×n is constructed. Step 2: Calculate the normalized initial direct-relation matrix. Normalize the initial direct-relation matrix Z with the following method and obtain the normalized initial direct-relation matrix Y = [yij]n×n. (6) In (6), by eliminating all rows and columns related to the absorbing states, the sub-stochastic matrix Y is obtained from an absorbing Markov chain matrix, and equipped with the following two properties [66]: (1) , where O is the null matrix; (2) , where I is the identity matrix. Steps 3: Calculate the over-relation matrix. Let T = [tij]n×n be over-relation matrix, and it can be obtained with (7) where tij denotes the overall degree of correlation between factors Fi and Fj. Steps 4: Determine the factor prominence and relations. Let ci represent the overall degree that factor Fi influences others and it can be expressed as (8) Let hi denote the overall degree to which factor Fi is influenced by others. hi can be obtained as follow (9) Let Di be the prominence of factor Fi that determine the importance of factor Fi and it is (10) The larger the value of Di, the more significant factor Fi is. Let Ri denote the relation of factor Fi an indicator for assessing the role of factor Fi and it is defined as (11) If Ri>0, then Fi is a cause factor; if Ri<0 then Fi is an effect factor. Steps 5: Build the causal and effect diagram. A causal and effect diagram with the horizontal axis D and vertical axis R is constructed to visualize the importance and classification of all factors based on prominence Di and relation Ri. 2.1 The intuitionistic fuzzy set The IFS proposed by Atanassov (1986) offers a viable approach to address vagueness. It is built on the foundation of classical fuzzy set theory [50]. Definition 1 Let set Y is fixed. An IFS A in Y is defined as (1) where uA(y),vA(y):Y→[0,1] are the membership function and non-membership function, respectively, with the condition 0≤uA(y)+vA(y)≤1. For each IFS in Y, let’s call (2) the intuitionistic index y of in A. It represents the hesitancy degree of y to A [59]. The author observed that (3) For each fuzzy set F in Y, there is . Hence, fuzzy sets are considered as the particular cases of IFSs. For convenience, let’s denote an intuitionistic fuzzy number (IFN) [50, 60] by a = (ua,va) where ua∈[0,1], va∈[0,1], ua+va≤1, and let Θ be the set of all IFNs. The score of S(a) can be assessed using the score function S shown as [61] (4) where S(a) = [0,1]. For any three IFNs a = (ua,va), and the following operational laws hold true. (1) , (2) , (3) , (4) . Definition 2 Let be a collection of IFNs, and let the intuitionistic fuzzy weighted average (IFWA) [62], IFWA:Θn→Θ, if (5) then IFWA is the intuitionistic fuzzy weighted averaging value, where δ = (δ1,δ2,…,δn)T is the weight vector of aj(j = 1,2,…,n) with δj∈[0,1] and . 2.2 The DEMATEL method This section will present the principal of the DEMATEL method. Suppose that the set of factors is F = {F1,F2,…,Fn} the procedure steps of the method are as follow [63–65]: Steps 1: Establish the initial direct-relation matrix. Let zij denotes the indicated degree of the dependence between factors Fi and Fj. Particularly, there does not exist dependence between Fi and itself. Then the initial direct-relation matrix Z = [zij]n×n is constructed. Step 2: Calculate the normalized initial direct-relation matrix. Normalize the initial direct-relation matrix Z with the following method and obtain the normalized initial direct-relation matrix Y = [yij]n×n. (6) In (6), by eliminating all rows and columns related to the absorbing states, the sub-stochastic matrix Y is obtained from an absorbing Markov chain matrix, and equipped with the following two properties [66]: (1) , where O is the null matrix; (2) , where I is the identity matrix. Steps 3: Calculate the over-relation matrix. Let T = [tij]n×n be over-relation matrix, and it can be obtained with (7) where tij denotes the overall degree of correlation between factors Fi and Fj. Steps 4: Determine the factor prominence and relations. Let ci represent the overall degree that factor Fi influences others and it can be expressed as (8) Let hi denote the overall degree to which factor Fi is influenced by others. hi can be obtained as follow (9) Let Di be the prominence of factor Fi that determine the importance of factor Fi and it is (10) The larger the value of Di, the more significant factor Fi is. Let Ri denote the relation of factor Fi an indicator for assessing the role of factor Fi and it is defined as (11) If Ri>0, then Fi is a cause factor; if Ri<0 then Fi is an effect factor. Steps 5: Build the causal and effect diagram. A causal and effect diagram with the horizontal axis D and vertical axis R is constructed to visualize the importance and classification of all factors based on prominence Di and relation Ri. 3 The proposed method In this section, an extended DEMATEL method with the IFSs is proposed to better identify the importance and classification of factors. A formal procedure of the proposed method is provided as illustrated in Fig 1. Firstly, the initial intuitionistic fuzzy direct-relation matrices are aggregated into a group intuitionistic fuzzy direct-relation matrix using IFWA. Then, normalized intuitionistic fuzzy direct-relation matrix and intuitionistic fuzzy over-relation matrix are constructed to derive the prominence and relation of factors. Furthermore, establish the importance of factors with respect to the prominences, and the classification of factors based on the relations. Finally, a causal and effect diagram is constructed based on the prominences and relations. Download: PPT PowerPoint slide PNG larger image TIFF original image Fig 1. Resolution process of the presented problem. https://doi.org/10.1371/journal.pone.0314650.g001 Let E = {E1,E2,…,Em} be a finite set of decision-makers and Ek (k = 1,2,…,m) represent the kth decision-maker. The author makes the assumption that the decision makers share equal importance. F = {F1,F2,…,Fn} is a finite set of factors and Fi (i = 1,2,…,n) denotes the ith factor. Let be the judgment of decision maker Ek on the correlation degree between factors Fi and Fj. Particularly, indicates that there is no correlation between Fi and itself. Then the initial intuitionistic fuzzy direct-relation matrix Zk provided by decision-maker Ek can be built as follows. (12) Then, the author aggregates all the individual decision opinions into a group opinion, namely the group intuitionistic fuzzy decision matrix, which can be derived with Eq (5) and denoted as Z, i.e., (13) where . The normalized intuitionistic fuzzy direct-relation matrix of the group intuitionistic fuzzy direct-relation matrix Z, denoted as Y is given by (14) where . Two crisp values matrices from Y are obtained as follows: (15) The following proposition enables the computation of Y to be achieved by the multiplication of crisp matrices. Proposition 1 Let Yτ = [(yij)τ]n×n, where . Two matrices are further defined, (16) Proof The proof is skipped due to its straightforwardness by using matrix multiplication. Proposition 2 , and . Proof The augmented matrix U′ as below is obtained by adding a row and a column to the matrix U. , where u1,n+1,u2,n+1,…, un,n+1 make U′ be a stochastic matrix. Since , there is at least one u1,n+1,u2,n+1,…, un,n+1 is positive. Therefore, U′ is a stochastic matrix of an absorbing Markov chain and matrix U is the sub-stochastic matrix of U′. Let ϱ(U) be the spectral radius of matrix U, and the sufficient and necessary condition of ϱ(U)<1 is [65, 66]. Hence, there is . can be proved with the analogous procedures. The author defined the intuitionistic fuzzy over-relation matrix T as Eq (17), following the classical DEMATEL method [67] (17) Proposition 3 Let T = [tij]n×n,i,j = 1,2,…,n, where then , and . Proof By (17) and Propositions 1 and 2, there is and then there will be . This completes the proof of Proposition 3. Then the author will prove the elements of the T also are IFNs. Proposition 4 Let then Proof By (13) and Proposition 1 there are 0≤uij,vij≤1 and O≤U,V≤I. Based on the matrix operation properties [68], there are (I−U)−1