Identifying Key Performance Indicators for Holistic Hospital Management with a Modified DEMATEL Approach

Performance analysis is an important way for hospitals to achieve higher efficiency and effectiveness in providing services to their customers. The performance of the healthcare system can be measured by many indicators, but it is difficult to improve them simultaneously due to the limited resources. A feasible way is to identify the central and influential indicators to improve healthcare performance in a stepwise manner. In this paper, we propose a hybrid multiple criteria decision making (MCDM) approach to identify key performance indicators (KPIs) for holistic hospital management. First, through integrating evidential reasoning approach and interval 2-tuple linguistic variables, various assessments of performance indicators provided by healthcare experts are modeled. Then, the decision making trial and evaluation laboratory (DEMATEL) technique is adopted to build an interactive network and visualize the causal relationships between the performance indicators. Finally, an empirical case study is provided to demonstrate the proposed approach for improving the efficiency of healthcare management. The results show that “accidents/adverse events”, “nosocomial infection”, ‘‘incidents/errors”, “number of operations/procedures” are significant influential indicators. Also, the indicators of “length of stay”, “bed occupancy” and “financial measures” play important roles in performance evaluation of the healthcare organization. The proposed decision making approach could be considered as a reference for healthcare administrators to enhance the performance of their healthcare institutions.


Introduction
In the past decades, healthcare systems have been involved in many changes ranging from technological to normative. But various issues have been raised in the healthcare sector and the operation of hospital management needs great improvement. As evidence for how well organizational objectives or goals are achieved, performance measurements have emerged in healthcare organizations [1,2]. Evaluating clinical and service performance is significant to drive quality excellence and achieve more effective performance. Thus, many national indicator projects [3,4] as well as international projects [5,6] have been initiated for measuring the quality and performance of healthcare services. These existing performance assessment frameworks include a large number of indicators [7], but it is difficult to improve them simultaneously since the resources of hospital management are constrained. Furthermore, the optimization of a specific process deviates from the overall welfare of a whole system, but few of the current frameworks can analyze and enhance healthcare performance systematically. Therefore, it is of great importance to determine a limited number of essential performance indicators based on their influences and relationships to address the management issues of healthcare organizations from a holistic point of view.
In previous studies, key performance indicators (KPIs) for healthcare management are commonly identified through case studies or interviews with healthcare experts [4,8,9]. However, these studies mainly focus on a generic framework of screening effective indicators for performance measurement; almost none has examined the creation of a mechanism that distinguishes causal relationships between KPIs. In other words, there is a lack of articulation of the cause and effect relationships between performance indicators when implemented for performance improvement. Consequently, many hospitals dilute the efforts of their performance measurement systems because of basic mistakes in mapping [10]. The decision making trial and evaluation laboratory (DEMATEL) method [11], initiated by the Geneva Research Centre of the Battelle Memorial Institute, is pragmatic to visualize the structure of complicated causal relationships and clarify the essentials of a system. It is a method based on digraph and matrix theory and has the ability to divide multiple factors into cause and effect groups [12]. Through analyzing direct and indirect relations between system elements, an ideal way to solve the intertwined problems can be obtained by the DEMATEL [13][14][15]. Although there is an abundance of literature using DEMATEL, this effective structural modeling tool has not yet been applied to the healthcare performance management. Essentially, for a large number of indicators impacting each other, healthcare performance measurement can be considered as a complex system. Hence, it is an interesting research topic for applying the DEMATEL technique to capture the network of interdependencies among performance indicators and identify KPIs.
In many real-life situations, experts' judgments are vague and it is difficult to estimate direct effect between elements with an exact numerical value. Instead, domain experts feel more comfortable providing their knowledge with linguistic expressions due to the complexity of the considered objects. Moreover, they may demonstrate different opinions from each other and produce different types of assessment information of performance indicators, some of which may be precise or imprecise, certain or uncertain, and complete or incomplete. These various types of linguistic assessments are very difficult to incorporate into the classical DEMATEL approach. The evidential reasoning (ER) approach was developed based on Dempster-Shaffer theory and distributed assessment framework for analyzing decision making problems with different types of uncertainties [16,17]. It can handle various kinds of human judgments such as uncertainty, fuzziness and ignorance and allow decision makers to express opinions in a flexible manner. Due to its characteristics and benefits, the ER approach has been widely applied to address decision-making problems in different areas [18,19]. To manage linguistic information without loss of information, the 2-tuple linguistic computational model was proposed by Herrera, Martínez [20] for computing with words. A well-known extension of 2-tuple linguistic method is the interval 2-tuple linguistic model [21], which employs uncertain linguistic variables to represent linguistic evaluation information. Recently, many of studies have reported decision-making models under the interval 2-tuple linguistic setting [22,23]. Consequently, the normal DEMATEL method can be extended with the ER and interval 2-tuples for more accurately constructing the network structure of interdependent healthcare performance indicators.
Based on the above analyses, this paper aims to develop a hybrid multiple criteria decision-making (MCDM) approach, called linguistic evidential DEMATEL, to analyze the total relations of performance indicators and identify KPIs for healthcare performance improvement. For the convenience of expressing uncertainties, linguistic assessments provided by domain experts are modeled via the ER approach. The interval 2-tuple linguistic model is utilized to aggregate individual opinions of experts into group assessments. Next, the DEMATEL technique is applied on the group assessment matrix, and the indicators that are of more essential importance for the whole system can be recognized. For doing so, the remainder of this paper is organized as follows: Section 1 reviews the literature on interval 2-tuples and the DEMATEL method. Section 2 introduces some basic concepts and definitions which will be utilized in this study. In Section 3, we present an extended DEMATEL approach to identify KPIs for holistic hospital management. In Section 4, a case study is conducted to illustrate the application of the proposed hybrid approach. Finally, conclusions and future research directions are given in the Section 5.

Applications of Interval 2-Tuples
Since its introduction, the interval 2-tuple linguistic method has received great attention from researchers and has been applied to many fields. For example, Xue et al. [24] combined hesitant fuzzy linguistic term sets with interval 2-tuples to deal with the hesitancy and uncertainty of decision makers and developed a QUALIFLEX-based linguistic method for robot evaluation and selection. Singh et al. [23] addressed the energy-planning problem in a new power plant set up by extending the PROMETHEE II approach to the interval 2-tuple linguistic environment. Lu et al. [25] solved the healthcare waste treatment technology selection problem with a hybrid decision-making model using interval 2-tuple induced distance operators and technique for order preference by similarity to an ideal solution (TOPSIS) method. You et al. [26] proposed an interval 2-tuple linguistic VIKOR approach for group multi-criteria supplier selection within incomplete and uncertain information environment, and Liu et al. [27] presented an attitudinal-based interval 2-tuple linguistic VIKOR method to select the best disposal site in municipal solid waste management. Liu et al. [28] suggested a failure mode and effect analysis (FMEA) method using interval 2-tuple linguistic variables and grey relational analysis (GRA), and Liu et al. [29] assessed the risk of healthcare failure modes based on interval 2-tuple hybrid weighted distance measure. In addition, Meng et al. [30] defined some generalized interval 2-tuple linguistic correlated aggregation operators based on Choquet integral and generalized Shapley function for solving the multiple criteria decision-making (MCDM) problems where the elements in a set are interdependent. Wang et al. [31] integrated interval 2-tuple linguistic method with Choquet integral aggregation operators to deal with the linguistic MCDM problems with correlated criteria. Lin et al. [32] put forward the generalized interval 2-tuple linguistic Shapley chi-square averaging operator to cope with interactive phenomenon among experts (or attributes) in MCDM. Liu et al. [33] proposed some interval 2-tuple linguistic Bonferroni mean operators for MCDM considering the interrelationships among input arguments, and Liu et al. [34] defined some dependent interval 2-tuple linguistic aggregation operators, in which the associated weights only depend on the input arguments and can relieve the effect of biased arguments on the decision results.

Applications of the DEMATEL Method
In the literature, a lot of researchers have adopted the decision making trial and evaluation laboratory (DEMATEL) method for the identification of key factors through analyzing the interaction relationships among them. For instance, Luthra et al. [35] evaluated the recognized key enablers in solar power initiatives and developments in India's current scenario based on fuzzy DEMATEL methodology. Liang et al. [36] identified the critical success factors influencing the sustainable development of biofuel industry in China using grey DEMATEL approach and suggested appropriate strategic measures to promote the sustainability of biofuel industry in China. Xia et al. [37] employed grey DEMATEL method to analyze the significant internal barriers for automotive parts remanufacturers in China and help them increase the chance of practicing more productive remanufacturing. Wu et al. [38] used fuzzy DEMATEL method to investigate the decisive factors for green supply chain practices implementation in the Vietnamese automobile manufacturing industry. Liu et al. [39] determined critical failure modes in system FMEA by combining fuzzy weighted average with fuzzy DEMATEL model, and Guo et al. [40] evaluated green corporate social responsibility indicators and the degrees of influence and causal relationships among them via a novel DEMATEL technique. Akyuz, Celik [41] adopted fuzzy DEMATEL approach to evaluate critical operational hazards during gas freeing process in crude oil tankers, and Wu [42] utilized the DEMATEL and balanced scorecard (BSC) methods to identify critical central and influential KPIs for improving banking performance. Cheng et al. [43] investigated the concepts and relationships of KPIs through mind map and extract KPIs of new product development by analyzing their cause and effect relationships. Additionally, Govindan et al. [44] developed an intuitionistic fuzzy DEMATEL method for handling the important and causal relationships between green supply chain management practices and performances. Quader et al. [45] used a 2-tuple DEMATEL technique to evaluate the top influencing factors for CO 2 capture and storage in the iron and steel industry. Li et al. [46] developed an evidential DEMATEL method based on intuitionistic fuzzy numbers and Dempster-Shaffer theory to identify critical success factors in emergency management.

Interval 2-Tuple Linguistic Variables
The interval 2-tuple linguistic model [21] is an extension of 2-tuple linguistic approach [20,47] for better representing uncertain linguistic decision information. The main advantages of this formalism are that decision makers can express their judgments by the use of different linguistic term sets, and various uncertainties in the assessments of decision makers can be well reflected and modeled [48][49][50]. Definition 1. Let S = s 0 , s 1 , ..., s g be a linguistic term set. An interval 2-tuple linguistic variable is composed of two 2-tuples, denoted by [(s k , α 1 ), (s l , α 2 )], where (s k , α 1 ) ≤ (s l , α 2 ), s k (s l ) and α 1 (α 2 ) represent the linguistic label of the linguistic term set S and symbolic translation, respectively. The interval 2-tuple that expresses the equivalent information to an interval value [β 1 , is computed by the following function [21]: On the other hand, there is always a function ∆ −1 so that an interval 2-tuple can be converted into an interval value [β 1 , It is worth highlighting that if (s k = s l ) and (α 1 , α 2 ), the interval 2-tuple linguistic variable reduces to a 2-tuple linguistic variable.
then the main operations of interval 2-tuples are expressed as follows [48]: Definition 3. Let S = s 0 , s 1 , ..., s g be a linguistic term set and a = [(s k , α 1 ), (s l , α 2 )] is an interval 2-tuple, then the interval 2-tuple linguistic variable can be transferred into a crisp value bŷ where S( a) ∈ [0, 1] is the score function of the interval 2-tuple a [21].

Evidential Reasoning Approach
The ER approach is well suited for dealing with multiple criteria decision analysis problems considering the quantitative and qualitative measurements assessed using subjective judgments with uncertainties. Details of the ER approach can be easily found in current literature, see, e.g., [16,17]. Here, we only briefly review some basic concepts which will be used in our proposed decision making approach.
. For a frame of discernmentĤ = H ij , i, j = 1, ..., n , a belief structure is expressed by where E(C) stands for the performance evaluation in terms of a particular criterion, H ii for i = 1, . . . ,n are individual assessment grades defined to classify that criterion, H ij for i = 1, . . . .,n−1 and j = i + 1 to n are the interval assessment grades between H ii and H jj , and β ij are the belief degrees attached to the corresponding assessment grades. A belief degree generally represents the strength to which an answer is believed to be true and it must be equal to or less than 1. If ∑ n i,j=1 β ij = 1, the evaluation is complete, otherwise the evaluation is considered as incomplete.

Definition 5. (Rule of Combination). Let
H ij , β k ij , i, j = 1, ..., n; k = 1, ..., K be the belief structure of the kth piece of evidence on the frame of discernmentĤ, that has an associated weighing vector w = (w 1 , w 2 , ..., w K ) T , with w k ∈ [0, 1] and Σ K k=1 w k = 1, the combination of the K pieces of evidence, called collective belief structure, is denoted as where β K ij is referred to as collective belief degree and determined by

The DEMATEL Method
The decision making trial and evaluation laboratory (DEMATEL) methodology is an appropriate approach for extracting the interdependent relationships and the intensity of interdependence between complex components of a system [13,14]. The main steps involved in the DEMATEL can be explained as follows: Step 1: Obtain the initial direct-influence matrix First, a committee of l experts E = {E 1 , E 2 , ..., E l } is formed to evaluate the direct influence between different factors based on a numeric scale from 0 to 4, representing "no influence (0)", "low influence (1)", "medium influence (2)", "high influence (3)", and "very high influence (4)", respectively.
As the result of evaluation, the initial direct-influence matrix Z k = z k ij n×n can be obtained for each expert, where all principal diagonal elements are equal to zero and z k ij represents the judgment provided by the kth expert on the degree to which factor i affects factor j. By aggregating all the experts' opinions, the group direct-influence matrix Z = z ij n×n is computed by Step 2: Calculate the normalized direct-influence matrix The normalized direct-influence matrix X = x ij n×n is calculated through Equations (8)- (9).
where s denotes the biggest value among the sums of each row and each column.
Step 3: Derive the total-influence matrix Once the normalized direct-influence matrix X is obtained, the total-influence matrix T = t ij n×n can be derived by summing all direct and indirect effects as shown in Equation (10).
where I is the identity matrix.
Step 4: Build the influential relation diagram (IRD) Using the values of R + C and R -C, the IRD can be mapped, where R and C are the sum of the rows and the sum of the columns in the total-influence matrix T, respectively.
Note that r i denotes the sum of all direct and indirect influences dispatched from factor i to other factors and is called the degree of influential impact. In contrast, c j represents both direct and indirect impacts that factor j receives and is called the degree of influenced impact. Let i = j and i, j ∈ {1, 2, ..., n}; the value of r i + c j named "prominence" indicates the total effects given and received by the ith factor. That is, it represents both the ith factor's impact on the whole system and other factors' impact on the ith factor. So, r i + c j is a measure of importance degree of the ith factor in the entire system. Similarly, the difference r i − c j called "relation" shows the net effect that contributed by the ith factor to the system. Specifically, if r i − c j > 0, the ith factor is a net cause and belongs to the cause group, and if r i − c j < 0, the ith factor is a net receiver and belongs to the effect group. Finally, an IRD can be plotted using r i + c j as the horizontal axis and r i − c j as the vertical axis. Moreover, the entries in the matrix T can be mapped into the diagram by using solid lines and broken lines for two-way and one-way significant relationships, respectively.
Above, an IRD is constructed based on the matrix T to explain the structure relations of factors. However, in some situations, the IRD may be too complex to show valuable information for decision making if all the relations are considered. Therefore, a threshold value θ is often set to filter out negligible effects to diminish this complexity of the whole system. That is, only the elements of the matrix T, whose influence level is greater than the threshold value, are selected and converted into the IRD. If the threshold value is too low, many factors are included and the IRM will be too complex to comprehend. On the contrary, some important factors may be excluded if the threshold value is too high. In this study, the threshold value is determined by computing the average of the matrix T [12].

The Proposed Linguistic Evidential DEMATEL Approach
Quantifying the interrelationships of indicators and identifying key performance indicators (KPIs) for healthcare management can be considered as a group decision behavior. In many complex systems, subjective evaluations of experts are always expressed in linguistic terms rather than crisp values. Furthermore, considering their different expertise and backgrounds, decision makers may provide different types of assessment information for the interdependent relationships between performance indicators, some of which may be precise or imprecise, certain or uncertain, and complete or incomplete. Under such circumstances, the classical decision making trial and evaluation laboratory (DEMATEL) method cannot be applied directly to figure out the KPIs in relation to the success of healthcare performance measurement. In this section, we propose a linguistic evidential DEMATEL approach to address the uncertainties of human judgments to identify the importance and classification of performance indicators. The flowchart in Figure 1 shows the general framework of the proposed approach. In this framework, belief structures of the evidential reasoning (ER) approach are implemented to express ambiguous concepts associated with decision makers' subjective judgments. Then, the interval 2-tuple linguistic model is utilized to aggregate individual evaluations of direct influence of performance indicators into group opinions. At last, the DEMATEL is applied on the acquired results to identify a minimum set of KPIs for holistic hospital management.

The Proposed Linguistic Evidential DEMATEL Approach
Quantifying the interrelationships of indicators and identifying key performance indicators (KPIs) for healthcare management can be considered as a group decision behavior. In many complex systems, subjective evaluations of experts are always expressed in linguistic terms rather than crisp values. Furthermore, considering their different expertise and backgrounds, decision makers may provide different types of assessment information for the interdependent relationships between performance indicators, some of which may be precise or imprecise, certain or uncertain, and complete or incomplete. Under such circumstances, the classical decision making trial and evaluation laboratory (DEMATEL) method cannot be applied directly to figure out the KPIs in relation to the success of healthcare performance measurement. In this section, we propose a linguistic evidential DEMATEL approach to address the uncertainties of human judgments to identify the importance and classification of performance indicators. The flowchart in Figure 1 shows the general framework of the proposed approach. In this framework, belief structures of the evidential reasoning (ER) approach are implemented to express ambiguous concepts associated with decision makers' subjective judgments. Then, the interval 2-tuple linguistic model is utilized to aggregate individual evaluations of direct influence of performance indicators into group opinions. At last, the DEMATEL is applied on the acquired results to identify a minimum set of KPIs for holistic hospital management. The procedure of the proposed linguistic evidential DEMATEL is summarized as follows: Step 1. Determine possible performance indicators influencing system objective A systematic review of literature is required to search and collect relevant information at this stage. Additionally, it is necessary to establish a committee of domain experts, who provide group knowledge for related issues. In this step, possible performance indicators used for measuring and managing healthcare systems should be determined based on collected information and experts' opinions.
Step 2. Invite experts to assess direct influence between performance indicators After the agreement about initial performance indicators is formed, questionnaire survey of the healthcare experts should be conducted to investigate the interaction between each pair of indicators. As a result, the linguistic assessments regarding the direct influence which indicators have impacts on each other can be obtained. The procedure of the proposed linguistic evidential DEMATEL is summarized as follows: Step 1. Determine possible performance indicators influencing system objective A systematic review of literature is required to search and collect relevant information at this stage. Additionally, it is necessary to establish a committee of domain experts, who provide group knowledge for related issues. In this step, possible performance indicators used for measuring and managing healthcare systems should be determined based on collected information and experts' opinions.

Step 2. Invite experts to assess direct influence between performance indicators
After the agreement about initial performance indicators is formed, questionnaire survey of the healthcare experts should be conducted to investigate the interaction between each pair of indicators. As a result, the linguistic assessments regarding the direct influence which indicators have impacts on each other can be obtained.
As mentioned before, different types of uncertainties are inevitably involved in the judgments the expert group provided on the interrelationships among indicators due to the subjectivity and incompleteness. Such a problem can be solved using the ER approach, which can capture information of different forms and accommodate uncertainties of different types under a unified framework. Suppose the individual evaluation grade set H is defined as then experts can provide their subjective judgments in the following flexible ways: No judgment, which means the expert is not willing to or cannot provide an assessment due to a lack of evidence or data. In other words, the grade by this expert could be anywhere between Very Low and Very High and can be represented as {(H 15 , 1.0)}.
Step 3. Aggregate experts' assessments to compute the initial direct-influence matrix Aggregate the linguistic evaluations of experts into group assessments by applying the combination rule of Equations (5)-(6) to every element of the linguistic belief decision matrixes. Then convert the obtained collective belief structures into interval 2-tuples using the weighted sum method and turn them into crisp scores by Equation (3). In this way, we can obtain z ij , which represents the direct impact of indicator i on indicator j. Thereby, an initial direct-influence matrix Z = z ij n×n can be constructed, where Z is a nonnegative matrix.

Step 4. Draw IRD and establish the structural model of performance indicators
Based on the initial direct-influence matrix Z, the total-influence matrix T can be easily calculated by applying the DEMATEL approach in Equations (8)- (10). Then the importance degree and net effect degree are computed for each performance indicator to construct the IRD; thus the structural relationship and importance of performance indicators can be visualized.
Step 5. Analyze the structure of performance indicators to identify KPIs Analyze each performance indicator in comprehensive consideration of the indexes r i , c i , r i + i , r i − c i and the visualized IRD. According to the position of each performance indicator in the whole system, the ones which have great effect on other indicators or have complicated relationship with other indicators can be identified. Clearly, performance indicators of this kind are KPIs and improvement of them would bring improvement to the overall healthcare system.

Application
Based on a systematic literature review on healthcare performance assessment, the performance indicators frequently used in previous studies were identified from the perspectives of patient, employee and management [1]. The literature survey searched PubMed and PubMed Central databases for articles applying performance indicators for the measurement of healthcare performance. As a result, 24 relevant papers were elicited according to the predefined screening criteria, from which a total of 428 indicators were extracted. Among them, 15 indicators appeared in five or more of the selected studies are considered as "frequently used" indicators (cf. Table 1). Next, these performance indicators are analyzed and structured following the procedure of our proposed linguistic evidential decision making trial and evaluation laboratory (DEMATEL) approach.
To measure interrelationships among the 15 performance indicators, surveys and interviews were conducted in a university hospital to make assessments in terms of influence and relationship. For doing so, an expert committee of five decision makers was built to complete the questionnaire. These experts from different departments or institutions included a head nurse and a medical doctor who is responsible for an internal medicine department, a hospital leader and a deputy hospital leader who are also physicians, and a professor of industrial engineering who is researcher on healthcare operation management. All these decision makers possess professional knowledge of healthcare performance measurement and have worked in related fields for more than three years. To reflect their differences in decision-making, the five experts were assigned the following weights: 0.15, 0.20, 0.25, 0.25 and 0.15, respectively. In addition, due to the difficulty in precisely assessing the interdependency of indicators, the experts agreed to evaluate them using the linguistic term set S: The five decision makers gave their linguistic assessments about the direct influence between each pair of performance indicators. Because it was different from one to another and included incomplete information, the experts' assessments were expressed by using linguistic belief structures. For example, the initial direct-influence matrix given by the first expert Z 1 is presented in Table 2. Next, the linguistic belief structures are aggregated into collective belief structures using Equations (5) and (6). Table 3 displays the group assessments of the five experts on the patient related indicators. Then Equation (3) is applied to transform the interval 2-tuple linguistic matrix into crisp values to construct the initial direct-relation matrix Z = z ij 15×15 as in Table 4.
Next, the data of the interrelationships between performance indicators are were analyzed with the DEMATEL approach. The total-influence matrix T is computed by Equations (8)-(10) and given in Table 5. Then, the indexes of R + C and R -C of each performance indicator are computed via Equations (11)- (12). The computed results are summarized in Table 6. Finally, an IRD of the 15 indicators is acquired as shown in Figure 2. Note that the arithmetic mean of the matrix T (i.e., 0.0951) is used as threshold value to filter out weaker influential relations among the performance indicators.
As seen in Figure 2, the 15 indicators are visually divided into two groups according to whether their values of r i − c i are positive or negative. As a result, the cause group with positive r i − c i values includes P4, P5, P3, M2, M3, E2, and other indicators including E3, M1, M4, P6, E1, M5, M6, P2, P1 are in the effect group. There are many other valuable clues that can be acquired from Figure 2 to facilitate making healthcare improvement decisions. In the following, the IRD will be discussed in detail, and the KPIs for holistic healthcare management are identified taking into account the scores in Table 6.

Discussions
Among the indicators of the cause group, "accidents/adverse events" (P4) has the highest value of r i − c i , which implies that P4 exerts more effect on the whole system than that it receives from other indicators. In addition, as shown in Table 6, the degree of influential impact of P4 is 2.293, which ranks first among all the causal indicators. It means that P4 has great effect on other indicators and that improvement of P4 can lead to improvement of the overall system. To enhance the performance of healthcare system, a low rate of accidents/adverse events should be first achieved. Therefore, P4 is a KPI that should be paid more attention to in the healthcare performance measurement system. Meanwhile, "nosocomial infection" (P5) and "incidents/errors" (P3) can be clustered as KPIs for similar reasons.
With respect to "number of operations/procedures" (M2), its net effect value r i − c i is in the fourth place among the cause indicators and its influential impact index r i is relatively high. The impact given by M2 to the whole system is great, and optimization of M2 will definitely improve the effectiveness and efficiency of the healthcare system. Hence, it is reasonable to identify "number of operations/procedures" (M2) as one of the KPIs.
Next, the characteristics of effect indicators are analyzed to figure out possible KPIs. Among the performance indicators, "financial measures" (M6) has the highest value of r i + c i , indicating that it is quite important for the healthcare management system. However, as can be seen in Table 6, the r i − c i value of M6 is −0.829, a value less than zero showing M6 is a net effect indicator. To further elucidate this scenario, its r i and c i values (i.e., 1.499 and 2.328) are relatively high in all the indicators. This denotes that although M6 is a net receiver, it plays a significant role in promoting and improving healthcare management activities. Thus, "financial measures" (M6) is a KPI considering its complex relationships with other indicators. Following the similar logic, "length of stay" (M4) and "bed occupancy" (M5) can be considered as KPIs in the healthcare system.
Based on the aforementioned discussions, seven KPIs for hospital performance management from the patient, employment and management perspectives are identified. These KPIs include four-cause indicators and three-effect indicators which are "accidents/adverse events", "nosocomial infection", "incidents/errors", "number of operations/ procedures", "length of stay", "bed occupancy", and "financial measures". Note that depending on the particular hospital applied, the assessments of performance indicators may be dissimilar and thus different KPIs may be derived. However, as for the considered case study, it is a favorable way to improve the healthcare performance through optimizing the seven KPIs. Additionally, the identified KPIs can be considered as crucial criteria to improve the efficiency of healthcare management during the planning phase.
To verify the performance of the proposed hybrid MCDM approach, we gathered medical doctors and managers to check the obtained KPIs. According to the domain experts, the proposed linguistic evidential DEMATEL model is highly suitable for the considered application to analyze the cause-effect relationships among indicators and identify KPIs which are of more fundamental importance for the whole healthcare system. Accordingly, by monitoring and managing these KPIs, the effects of prevention, continual improvements, and innovations can be achieved to shape core competitiveness of the healthcare organization.      Based on the aforementioned discussions, seven KPIs for hospital performance management from the patient, employment and management perspectives are identified. These KPIs include four-cause indicators and three-effect indicators which are "accidents/adverse events", "nosocomial infection", "incidents/errors", "number of operations/ procedures", "length of stay", "bed occupancy", and "financial measures". Note that depending on the particular hospital applied, the assessments of performance indicators may be dissimilar and thus different KPIs may be derived. However, as for the considered case study, it is a favorable way to improve the healthcare performance through optimizing the seven KPIs. Additionally, the identified KPIs can be considered as crucial criteria to improve the efficiency of healthcare management during the planning phase.
To verify the performance of the proposed hybrid MCDM approach, we gathered medical doctors and managers to check the obtained KPIs. According to the domain experts, the proposed linguistic evidential DEMATEL model is highly suitable for the considered application to analyze the cause-effect relationships among indicators and identify KPIs which are of more fundamental

Conclusions
In this paper we introduced a novel performance indicator identification and assessment framework, named as linguistic evidential decision making trial and evaluation laboratory (DEMATEL), for holistic hospital management. This model extends the traditional DEMATEL method to capture decision information of different forms and to accommodate uncertainties of different types by applying the evidential reasoning (ER) approach and interval 2-tuple linguistic variables. It successfully determined the importance of each indicator and acquired the causal relationships between performance indicators. As a result, a set of complex indicators were separated into cause and effect groups and a visible influential relation diagram (IRD) was formed. Based on the results of the case study, seven key performance indicators (KPIs) were defined for the healthcare performance management. Therefore, hospital managers can apply a phased quality control method to improve the healthcare performance gradually with the limited resources.
Despite its contributions, this study has several limitations, which may be addressed by future research. First, the weights of domain experts in the proposed approach are determined through a direct assignment method. In many practical situations, however, the information of expert weights is completely unknown or partially known. Therefore, future research should be conducted in developing an optimal model to derive expert weights objectively. Second, in the case study, the data are collected from a small number of healthcare experts, which may limit the generalization of its findings to other hospitals. Thus, in the future, more case studies with accepted number of participants can be carried out to improve or further validate the findings of this study. Additionally, the proposed model of identifying KPIs for healthcare management can be easily applied to diverse real applications, such as determining critical success factors of emergency management and analyzing major internal barriers in electric vehicle industrialization.