Next Article in Journal
Multivalued Common Fixed Points Theorem in Complex b-Metric Spaces
Next Article in Special Issue
Two Approaches to Estimate the Shapley Value for Convex Partially Defined Games
Previous Article in Journal
A Novel Memory Concurrent Editing Model for Large-Scale Video Streams in Edge Computing
Previous Article in Special Issue
ES Structure Based on Soft J-Subset
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

GRA-Based Dynamic Hybrid Multi-Attribute Three-Way Decision-Making for the Performance Evaluation of Elderly-Care Services

1
School of Management Science and Engineering, Shandong University of Finance and Economics, Jinan 250014, China
2
School of Management Engineering, Shandong Jianzhu University, Jinan 250101, China
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(14), 3176; https://doi.org/10.3390/math11143176
Submission received: 15 June 2023 / Revised: 15 July 2023 / Accepted: 18 July 2023 / Published: 19 July 2023
(This article belongs to the Special Issue Mathematical Methods for Decision Making and Optimization)

Abstract

:
As an important branch of modern decision-making theory, multi-attribute decision-making (MADM) plays an important role in various fields. Classic MADM methods can provide a ranking of alternatives, and decision-makers need to evaluate the level subjectively based on the ranking results. Because of the limitation of knowledge, this is likely to lead to potential individual losses. Three-way decision (3WD) theory has good classification ability. Therefore, this paper proposes a dynamic hybrid multi-attribute 3WD (MA3WD) model. First, a new scheme for constructing loss functions is proposed from the perspective of gray relational analysis (GRA), which is an accurate and objective way to describe the relationship between loss functions and attribute values. Then, conditional probabilities are determined by employing the gray relational analysis technique for order preference by similarity to the ideal solution (GRA-TOPSIS). With these discussions, a GRA-based hybrid MA3WD model for a single period is proposed by considering multi-source information. Furthermore, by extending the single-period scenario to a multi-period one, we construct a dynamic hybrid MA3WD model, which can obtain the final three-way decision rules as well as the results of each period and each attribute. Finally, the proposed method is applied to the case of performance evaluation of elderly-care services to demonstrate the effectiveness of the method, and comparative analyses are given to verify the superiority of the proposed method.

1. Introduction

Decision-making is pervasive in our daily life, and multi-attribute decision-making (MADM) is an important part of it. Multi-attribute decision-making has been applied to various fields, including performance evaluation [1,2], risk assessment [3,4], and medical diagnosis [5,6]. With the continuous development of decision theory, classic MADM methods have been progressively proposed, such as TOPSIS [7,8], TODIM [9,10], VIKOR [11,12], etc. However, traditional MADM methods used in most studies are essentially a two-way decision-making procedure [13,14], which only considers two options, namely acceptance and rejection. For instance, when the government makes decisions on the access and exit mechanism of third-party social organizations, the result obtained by MADM methods is to renew a contract or terminate the contract according to the alternative rankings. However, there may be a contradiction where two social organizations with slight differences in evaluation receive opposite decision results, i.e., one is accepted to sign a new contract, while another one is obliged to terminate the contract. To conquer the limitation of traditional MADM methods, this paper introduces three-way decision theory under a multi-attribute decision-making environment. Three-way decision-making, derived from decision-theoretic rough sets (DTRSs), encompasses the ability to objectively classify all alternatives into three regions, i.e., a positive region, a negative region, and a boundary region. Yao provided a reasonable semantic interpretation of these three pairwise disjoint regions from the perspective of cognitive science, corresponding to acceptance, rejection, and non-commitment in the actual decision-making process, respectively [15,16]. Compared with two-way MADM methods, three-way decision-making adds a deferment strategy, which reduces the risk of decision-making. For instance, the two social organizations in the aforementioned example can be classified into a deferment region for further investigation to reduce the loss of misclassification.
In recent years, 3WD has made remarkable achievements in theoretical models and application studies. For theoretical models, many researchers have extended the traditional 3WD model to dynamic environments and proposed improved models, such as sequential three-way decision [17], multi-granulation three-way decision [18], multi-class three-way decision [19], multi-attribute three-way decision [20,21], etc. To solve MADM problems, Jia and Liu [20] proposed a multi-attribute three-way decision model, transforming attribute evaluation values to loss functions in a multi-attribute environment. Zhan et al. [22] introduced the concept of relative utility function, utilized similarity relation to determine conditional probabilities, and proposed a 3WD model in an incomplete fuzzy decision environment. Zhan et al. [23] proposed a MA3WD method based on outranking relations. Huang et al. [24] designed a pair of pre-order relations and further proposed a new method for calculating conditional probabilities. Gao et al. [13] proposed a three-way decision-based target–threat assessment method based on a comprehensive evaluation matrix integrating multi-period information. Subsequently, some researchers incorporated psychological behaviors into the construction of 3WD models. He et al. [25] proposed a behavioral MA3WD model considering the regret-aversion behavior and risk-aversion behavior of human beings. In application studies, 3WD has been widely applied in various practical issues, such as conflict analysis [26,27], medical diagnosis [28,29], investment decisions [30,31], concept analysis [32,33], etc. However, based on the above brief review of 3WD, the existing 3WD methods still have some deficiencies. First, the loss functions and conditional probabilities of 3WD are difficult to determine objectively and scientifically, and lack interpretability; Second, the study of 3WD in a time-dynamic environment is relatively limited, while the time-dynamic environment is the most common type of dynamic decision-making environment in real life.
Due to the complexity of the decision-making environment and the diversity of attributes, exact numbers often struggle to accurately express the fuzzy cognition of decision-makers. Therefore, to describe uncertain information more accurately and reasonably, various fuzzy numbers have been proposed and widely studied, such as interval numbers [34], triangular fuzzy numbers [35], and intuitionistic fuzzy numbers [36]. Moreover, due to the unquantifiable nature of certain attributes, linguistic variables become a powerful tool for describing such information [37]. However, in many practical decision-making problems, it is necessary to use a combination of quantitative and qualitative attribute types. For instance, in the performance evaluation of elderly-care services, attributes such as the number of service programs, professional level, and communication ability are often considered [38,39]. Among them, the number of service programs is a quantitative attribute, which is often expressed by exact numbers; communication ability is a qualitative attribute, often expressed by language variables. Hence, in the decision-making problem with diversity attributes, a single form of information cannot accurately describe all attributes, and it is effective and reliable to adopt a mixed approach using multiple forms of information according to practical needs.
Based on the above analysis of the research subjects, we summarize the motivations and innovations of this paper as follows:
Motivations:
  • In the process of constructing the loss function, many existing methods are unreasonable. For instance, in the classical three-way decision theory [15], the loss function is artificially set by the decision-maker, which is not scientific. In addition, Jia and Liu [20] converted the attribute value into a relative loss function in a multi-attribute decision-making environment. This method is also not rational because it uses 1 and 0 as the maximum and minimum of the attribute values in the construction of the relative loss function.
  • In the existing studies on MA3WD methods, there is limited consideration of the time-dynamic environment. Most methods are based on decision information from a single period [14,21], which leads to an incomplete evaluation of the objects. Gao et al. [13] proposed a 3WD method based on multi-period evaluation information to address the issue of time dynamics, which is an improvement. However, the model constructed by Gao et al. fails to provide continuous evaluations for objects across different periods.
  • Traditional MADM methods have limitations in solving practical problems such as the performance evaluation of elderly-care services. These methods’ essence lies in two-way decision-making, where the outcomes often overlook the necessity of further investigation, leading to potential individual losses.
Innovations:
  • A new scheme for constructing loss functions is proposed from the perspective of GRA, which is an accurate and objective way to describe the relationship between loss functions and attribute values.
  • Conditional probabilities are estimated based on GRA-TOPSIS, which provides comprehensive and objective results for three-way decisions.
  • A GRA-based hybrid MA3WD model considering mixed forms of information is proposed for evaluating objects at a specific period. The model can point out the specific attributes and periods of poor performance of the object so that it accurately improves its shortcomings.
  • By extending the single-period scenario to a multi-period one, we construct a GRA-based dynamic hybrid MA3WD model, extending the study of 3WD in a time-dynamic environment.
  • This paper introduces the 3WD theory into the performance evaluation of elderly-care services, which provides a scientific and reasonable way to solve this issue.
The remainder of this paper is organized as follows. Section 2 lists the basic concepts of multi-attribute three-way decision-making, the definition of a dynamic hybrid multi-attribute information system, and the basic knowledge of gray relational analysis. Section 3 introduces the concept of gray relational analysis into the construction of loss function, proposes a method to determine conditional probability using GRA-TOPSIS, and further proposes a GRA-based hybrid MA3WD model for a single period. In Section 4, a single-period scenario is extended to a multi-period one, and a dynamic hybrid MA3WD model is proposed using information from the multi-period. In Section 5, the proposed model is applied to the performance evaluation of government purchases of elderly-care services, and its effectiveness is verified through comparative analysis. Section 6 summarizes the final concluding remarks in this paper.

2. Preliminaries

2.1. Multi-Attribute Three-Way Decision (MA3WD)

The three-way decision derived from the DTRSs uses a two-state set and a three-action set to describe the decision process [15]. The set of states Ω = C , ¬ C indicates that an object is in C or not in C ; the set of actions A = a P , a B , a N corresponds to three actions x P O S C , x B N D C , x N E G C that are taken when classifying the object x , i.e., acceptance, non-commitment, and rejection. When the object x is in different states, the three actions taken produce corresponding losses or risks, which can be represented by the 3 × 2 loss function matrix [15] in Table 1.
Assume that Pr C | x represents the conditional probability that an object x belongs to a state C , then for any object x , the expected loss after taking all actions can be expressed as:
R a P | x = λ P P Pr C | x + λ P N Pr ¬ C | x
R a B | x = λ B P Pr C | x + λ B N Pr ¬ C | x
R a N x = λ N P Pr C | x + λ N N Pr ¬ C | x
Based on the Bayesian theory that the optimal decision is the one with the least cost, the decision rules are deduced as follows:
( P - 1 )   If   R a P | x R a B | x   and   R a P | x R a N | x ,   decide   x P O S ( C ) ;
( B - 1 )   If   R a B | x R a P | x   and   R a B | x R a N | x ,   decide   x B N D ( C ) ;
( N - 1 )   If   R a N | x R a P | x   and   R a N | x R a B | x ,   decide x N E G ( C ) .
Based on the decision rough sets, the decision rules (P-1)~(N-1) are three-way decisions.
Consider the following relationship between the loss functions:
λ P P λ B P < λ N P , λ N N λ B N < λ P N
Based on the above relationship and Pr ( C | x ) + Pr ( ¬ C | x ) = 1 , the decision rules (P-1)~(N-1) can be further simplified as:
( P - 2 )   If   Pr ( C | x ) α   and   Pr ( C | x ) γ ,   decide   x P O S ( C ) ;
( B - 2 )   If   Pr ( C | x ) α   and   Pr ( C | x ) β ,   decide   x B N D ( C ) ;
( N - 2 )   If   Pr ( C | x ) β   and   Pr ( C | x ) γ ,   decide   x N E G ( C ) .
where the thresholds α , β , and γ are as follows:
α = λ P N λ B N λ P N λ B N + λ B P λ P P
β = λ B N λ N N λ B N λ N N + λ N P λ B P
γ = λ P N λ N N λ P N λ N N + λ N P λ P P
When λ B P λ P P λ B N λ N N < λ P N λ B N λ N P λ B P the thresholds satisfy α > γ > β , then the decision rules (P-2)~(N-2) can be simplified as:
( P - 3 )   If   Pr ( C | x ) α ,   decide   x P O S ( C ) ;
( B - 3 )   If   β < Pr ( C | x ) < α ,   decide   x B N D ( C ) ;
( N - 3 )   If   Pr ( C | x ) β ,   decide   x N E G ( C ) .
Traditional 3WD can only deal with single-attribute decision-making issues. However, most real decisions encompass multi-attribute decision-making. Jia and Liu [20] proposed a new construction method of loss functions based on a multi-attribute environment, which uses attribute values to derive relative loss functions in a data-driven manner. For a MADM problem with m alternatives and n attributes, we denote the evaluation value of object A i on attribute c j as x i j [ N I S j , P I S j ] , where N I S j and P I S j represent the negative ideal solution and positive ideal solution on attribute c j , respectively. The relative loss function derived from x i j is shown in Table 2, where the relative loss of correct classification is 0, i.e., λ P P = λ N N = 0 , d x i j , N I S j and d x i j , P I S j , respectively represent the distance between the evaluation value x i j and the negative and positive ideal solutions under attribute c j , and σ j [ 0 , 0.5 ] represents the risk-avoidance coefficient of c j , reflecting the decision-maker’s preference for vague or accurate outcomes [20].
Then, after converting all attribute values x i j into corresponding loss functions λ x i j , use the attribute weight vector w = w 1 , w 2 , , w n to integrate λ x i j to obtain the comprehensive loss function of each object A i , as shown in Table 3. Next, the decision thresholds are calculated based on the comprehensive loss function, and then the classification result is obtained as the process of traditional 3WD.

2.2. Dynamic Hybrid Multi-Attribute Information System

Let I S = T , U , C , V , f be a dynamic hybrid multi-attribute information system. T = t k | k = 1 , 2 , , p is the set of non-empty finite periods. U = A i | i = 1 , 2 , , m is a non-empty set of finite objects, called the universe. C = C t k | k = 1 , 2 , , p is the set of p attribute sets, where C t k = c j t k | j = 1 , 2 , , n ; k = 1 , 2 , , p is the set of non-empty finite attributes for t k period. In addition, C t k has the following characteristic: N ( C t k ) N ( C t k + 1 ) , which means that the number of attributes in the attribute set in the corresponding period will change dynamically with the increase of k, i.e., it shows an increasing trend. W t k = w j t k | j = 1 , 2 , , n is the set of attribute weights for t k period, w j t k represents the weight of attribute c j for t k period, and satisfying w j k > 0 , j = 1 n w j k = 1 . In addition, the attribute weight has the following features: different attributes have different weights, and even the same attribute has different weights in different t k periods. V = V c j t k | j = 1 , 2 , , n ; k = 1 , 2 , , p is the set of attribute values, and V c j t k is the domain of the attribute c j t k . f : U × C V is a function; for any c j t k C , A i U , there is f A i , c j t k V c j t k , where f A i , c j t k is hybrid data, which may be a real number, interval number, triangular fuzzy number, linguistic term, or intuitionistic fuzzy number.

2.3. Gray Relational Analysis

Gray relational analysis, derived from gray system theory, is a modern statistical method of systematic analysis. The basic idea of this method is to judge the correlation degree between sequences by analyzing the similarity degree of geometric shapes [40,41]. The general steps of this method are as follows:
Suppose there are m alternatives, n attributes, and the decision matrix H = ( h i j ) m × n . Normalize the original decision matrix H to obtain the decision matrix V = ( v i j ) m × n .
(1)
Determine the positive ideal solution (PIS) and the negative ideal solution (NIS)
v j + = m a x i ( v i j ) v j = m i n i ( v i j ) j = 1 , 2 , , n
(2)
Calculate the gray relational coefficient r i j + of alternative A i ( i = 1 , 2 , , m ) from PIS about the attribute c j ( j = 1 , 2 , , n ) .
r i j + = min 1 i m min   1 j n d i j + + ρ max 1 i m max   1 j n d i j + d i j + + ρ max 1 i m max   1 j n d i j +
where d i j + = v j + v i j , ρ is the identification coefficient where 0 < ρ < 1 . There is no specific reference standard for the selection of ρ . The larger the ρ (closer to 1), the smaller the difference of the relational coefficient between different alternatives; the smaller the ρ (closer to 0), the greater the difference of the relational coefficient between different alternatives. Regarding the existing literature [41,42], we take ρ = 0.5 .
By the same token, calculate the gray relational coefficient r i j of alternative A i ( i = 1 , 2 , , m ) from NIS about the attribute c j ( j = 1 , 2 , , n ) .
r i j = min 1 i m min   1 j n d i j + ρ max 1 i m max   1 j n d i j d i j + ρ max 1 i m max   1 j n d i j
where d i j = v j v i j , ρ is the identification coefficient where 0 < ρ < 1 ; here we take ρ = 0.5 .
(3)
Calculate the gray relational degrees r ˜ i + and r ˜ i of alternative A i ( i = 1 , 2 , , m ) corresponding from PIS and NIS.
r ˜ i + = j = 1 n w i j r i j +
r ˜ i = j = 1 n w i j r i j

3. GRA-Based Hybrid MA3WD Model for Single Period

In Section 3, we first construct the relative loss functions from the perspective of gray relational analysis and explore the conditional probability with the GRA-TOPSIS method. Then, a GRA-based hybrid MA3WD model for a single period is established.
As the hybrid MADM problem is the concrete embodiment of a certain period in the dynamic hybrid MADM problem, we should first describe the problem of the latter. A dynamic MADM problem under a hybrid information environment can be explained as follows. Suppose there are a set of alternatives A = A 1 , A 2 , , A m , a set of periods T = t 1 , t 2 , , t p , and ξ = ξ 1 , ξ 2 , , ξ p is the weight vector of time series satisfying ξ k > 0 and k = 1 p ξ k = 1 . The decision attributes set in t k period is C k = c 1 k , c 2 k , , c n k . The attribute weight vector in t k period can be expressed as W k = w 1 k , w 2 k , , w n k , satisfying w j k > 0 and j = 1 n w j k = 1 . The hybrid MADM matrix in the t k period can be expressed as follows:
H k = h i j k   = h 11 k h 12 k h 1 n k h 21 k h 22 k h 2 n k h m 1 k h m 2 k h m n k
where h i j k is hybrid type data, which can be a real number, interval number, triangular fuzzy number, linguistic term, intuitionistic fuzzy number, etc.

3.1. Loss Function Based on Gray Relational Analysis

Jia and Liu used relative loss functions and inverse loss functions to convert attribute values into loss functions, illustrating the correlation between 3WD and MADM [20]. Inspired by Jia and Liu [20], in this section, we first propose a new method for constructing loss functions using GRA. By virtue of its ability to reflect trend differences in data series, GRA is used to represent the relationship between objects in the process of constructing loss functions. In addition, as an effective tool to deal with dynamic decision-making problems, GRA can describe the dynamic process of multi-period evaluation. Then, the comprehensive loss function of each alternative in a single period is constructed based on multiple attributes.
Gray relational coefficients are introduced to construct loss functions of three-way decisions. As the calculating process of gray relational coefficients has been given in Section 2.3, here we extend the traditional method of calculating the gray relational coefficients. When normalizing the decision matrix, we transform the linguistic term sets into triangular fuzzy numbers [43,44]. When determining the PIS and the NIS, the maximum value of the attribute in p periods is taken as the PIS, i.e., v j k + = max i , k v i j k , and the minimum value in p periods is taken as the NIS, i.e., v j k = min i , k v i j k . Since the data in the hybrid decision matrix is not limited to the exact number, d i j + and d i j are calculated according to the distance formula between the corresponding types of data.
In the t k period, let C j k , ¬ C j k present two states, where A i C j k indicates that A i possesses the property c j k and A i ¬ C j k indicates that A i does not possess the property c j k . { a P , a B , a N } presents three actions that indicate decisions of acceptance, non-commitment, and rejection, respectively. The relative loss functions of the alternative A i on the attribute c j k are expressed in Table 4. In Table 4, λ P P i j k = λ N N i j k = 0 indicates that sorting the alternative into the correct region causes no loss. λ N P i j k = 1 r i j k represents the loss incurred by taking rejection action in the state of A i C j k , where r i j k denotes the gray relational coefficient between alternative A i and the N I S j k . The larger the value of r i j k , the closer the association between v i j k and v j k , indicating that the alternative A i has less advantage over the N I S j k . Therefore, decision-makers will suffer the loss 1 r i j k if they reject the object A i belonging to C j k . λ P N i j k = 1 r i j k + represents the loss incurred by taking the acceptance action in the state of A i ¬ C j k , where r i j k + denotes the gray relational coefficient between alternative A i and the P I S j k . The larger the value of r i j k + , the closer the association between v i j k and v j k + , indicating that the alternative A i has less advantage over the P I S j k . Therefore, decision-makers will suffer loss 1 r i j k + if they accept the object A i belonging to ¬ C j k . λ B P i j k and λ B N i j k are related to λ N P i j k and λ P N i j k , which are represented by introducing the parameter σ j k [ 0 , 0.5 ] defined as the risk-avoidance coefficient [20]. Moreover, σ j k reflects the decision-maker’s preference for ambiguous or accurate outcomes. The larger the σ j k , the more accurate the result; the smaller the σ j k , the more ambiguous the result.
To obtain comprehensive three-way decision rules for each object in the tk period, we need to obtain the comprehensive loss function of A i in the t k period. First, the loss function of A i to c j k in Table 4 can be written as
λ v i j k = λ P P i j k λ B P i j k λ N P i j k λ P N i j k λ B N i j k λ N N i j k = 0 σ j k 1 r i j k 1 r i j k 1 r i j k + σ j k 1 r i j k + 0
Then, each attribute is weighted according to its importance to the decision result, and the attribute weight vector W k = w 1 k , w 2 k , , w n k is obtained. Next, the weighted average operator is used to integrate the losses of all attributes in the t k period to acquire the comprehensive loss function of each alternative A i in the tk period, which is expressed as
λ i k = j w j k λ v i j k = 0 j w j k σ j k 1 r i j k j w j k 1 r i j k j w j k 1 r i j k + j w j k σ j k 1 r i j k + 0

3.2. Estimating Conditional Probability by GRA-TOPSIS

Conditional probability is another essential concept of the three-way decision. In many studies, conditional probabilities are given directly according to the experiences of decision-makers; to reduce subjectivity and arbitrariness, Liang adopted TOPSIS to estimate the conditional probability in an objective way [45]. This method describes the differences between each alternative and ideal solution from the positional relationship between the data curves. However, it may cause errors due to data fluctuation. A method to reflect the significance of the alternative by the degree of variation in the data sequence is needed. Therefore, considering the idea in [45], we introduce a GRA-TOPSIS method to estimate the conditional probability. Section 2.3 described the calculation process of the gray relational degrees r ˜ i k + of the alternative Ai from the PIS and r ˜ i k of the alternative Ai from the NIS, and the process of the GRA-TOPSIS method is as follows:
Step 1. Compute the distances of Ai to ideal solutions.
The distance of the alternative Ai from the PIS is
d i k + = j w j k d v i j k , v j k +
The distance of the alternative Ai from the NIS is
d i k = j w j k d v i j k , v j k
where d v i j k , v j k + and d v i j k , v j k represent the distance from v i j k to v j k + and v j k , respectively.
Step 2. Obtain the dimensionless distances d i k + , d i k and dimensionless gray relation degrees r ˜ i k + , r ˜ i k , respectively
D i k + = d i k + max i   d i k + , D i k = d i k max i   d i k
R i k + = r ˜ i k + max r ˜ i k + i , R i k = r ˜ i k max r ˜ i k i
Step 3. Combine the distances D i k + , D i k and gray relational degrees R i k + , R i k as
L i k + = δ 1 D i k + δ 2 R i k + , L i k = δ 1 D i k + + δ 2 R i k
where δ 1 and δ 2 are preference coefficients, which reflect the preferences of decision-makers for the position and shape, and δ 1 + δ 2 = 1 , δ 1 , δ 2 [ 0 , 1 ] . There is no specific standard for the selection of δ i ( i = 1 , 2 ) . When δ 1 is larger than δ 2 , it indicates that decision-makers pay more attention to the influence of distance; when δ 2 is larger than δ 1 , it indicates that decision-makers pay more attention to the influence of shape. Referring to the existing literature [46,47], we take δ 1 = δ 2 = 0.5 , which indicates that decision-makers have no special preference for distance and shape.
Step 4. Calculate the relative closeness
R C i k = L i k + L i k + + L i k
Based on the semantics of Equation (21), R C i k indicates the probability of the object Ai being in the state C, then the conditional probability of the object Ai in the tk period as
Pr k C | [ A i ] = R C i k

3.3. GRA-Based Three-Way Decision Rules for Single Period

In the above two subsections, two basic elements of the proposed 3WD model, namely loss function and conditional probability, have been induced. Subsequently, the decision rules of the GRA-based three-way decisions for a single period are established in this section.
According to the traditional 3WD model, the expected losses R k a c | A i ( = P , B , N ) generated by three actions are calculated as
R k a P | A i = λ P P i k Pr k C | A i + λ P N i k Pr k ¬ C | A i
R k a B | A i = λ B P i k Pr k C | A i + λ B N i k Pr k ¬ C | A i
R k a N | A i = λ N P i k Pr k C | A i + λ N N i k Pr k ¬ C | A i
Since Pr k C | A i + Pr k ¬ C | A i = 1 , we can replace Pr k ¬ C | A i with 1 Pr k C | A i and obtain:
R k a P | A i = 1 Pr k C | A i λ P N i k
R k a B | A i = Pr k C | A i λ B P i k + 1 Pr k C | A i λ B N i k
R k a N | A i = λ N P i k Pr k C | A i
According to the discussion in Section 2.1, the decision rules are deduced as follows:
( P G R A )   If   R k a P | A i R k a B | A i   and   R k a P | A i R k a N | A i ,   decide A i P O S C ;
( B G R A )   If R k a B | A i R k a P | A i   and   R k a B | A i R k a N | A i ,   decide A i B N D C ;
( N G R A )   If   R k a N | A i R k a P | A i   and   R k a N | A i R k a B | A i ,   decide   A i N E G C .
Meanwhile, the loss function satisfies conditions λ P P λ B P < λ N P and λ N N λ B N < λ P N , so ( P G R A )~( N G R A ) can be simplified as:
( P G R A )   If   Pr k C | A i α i k   and   Pr k C | A i γ i k ,   decide   A i P O S C ;
( B G R A )   If   Pr k C | A i α i k   and   Pr k C | A i β i k ,   decide   A i B N D C ;
( N G R A )   If   Pr k C | A i β i k   and   Pr k C | A i γ i k ,   decide   A i N E G C .
The expressions of thresholds α i k , β i k and γ i k are
α i k = j w j k 1 σ j k 1 r i j k + j w j k 1 σ j k 1 r i j k + + j w j k σ j k 1 r i j k
β i k = j w j k σ j k 1 r i j k + j w j k σ j k 1 r i j k + + j w j k 1 σ j k 1 r i j k
γ i k = j w j k 1 r i j k + j w j k 1 r i j k + + j w j k 1 r i j k
Remark 1. 
There are two extreme threshold cases to pay attention to: (1) If  v i j k = v j k  for all j  in the t k  ( k = 1 , 2 , , p ) period, then α i k = β i k = γ i k = 1  will be obtained. This indicates that for the object A i  in the t k  ( k = 1 , 2 , , p ) period only the decision to reject is made; (2) If v i j k = v j k +  for all j  in the t k  ( k = 1 , 2 , , p ) period, then α i k = β i k = γ i k = 0  will be obtained. This indicates that for the object A i  in the t k  ( k = 1 , 2 , , p ) period only the decision to accept is made.

4. Dynamic Hybrid MA3WD Model for Multiple Periods

In this section, the calculation methods of weights for time series and attributes are discussed, respectively. Based on the GRA-based hybrid MA3WD model, a dynamic hybrid MA3WD model is then proposed by considering multi-period information. The specific process of the proposed model is shown in Figure 1.

4.1. Determination of Weights

(1)
Determination of time-series weights
An important feature of dynamic hybrid MA3WD is the time-series information, as data in different time series have different relative importance for final results. In this paper, the entropy-weight method is applied to obtain the weights of the time series. Nonlinear programming based on maximum entropy value is constructed to obtain the time-series weights as follows:
max I = k = 1 p ξ k ln ξ k s t λ = k = 1 p p k p 1 ξ k k = 1 p ξ k = 1 0 ξ k 1 k = 1 , 2 , , p
where ξ k is the weight of the t k period, and p is the number of periods. λ [ 0 , 1 ] means “time degree”, which reflects the importance of time series in the process of operator assembly. A value of λ closer to 0 denotes that the decision-maker places more emphasis on recent data; a value of λ closer to 1 denotes that the decision-maker places more emphasis on distant data. In addition, λ = 0.5 denotes that the decision-maker places equal importance on each period and has no special preference [48].
(2)
Determination of attribute weights
To obtain the weights of attributes, a combination of BWM and entropy-weight methods is utilized. It avoids the excessive subjectivity in BWM, as well as the imbalance of indicator weights in the entropy-weight method. The specific steps of the BWM-entropy-weight method are described as follows.
Step 1. Obtain the subjective weights of attributes by BWM. BWM is a multi-criteria decision-making method proposed by Rezaei [49]. The method is based on the idea of pairwise comparison of indicators. First, the decision-makers select the optimal c B and the worst from the evaluation indicators. Then, c B and c W are compared with other indicators to construct the comparison vectors of the optimal and worst indicators, which are, respectively, expressed as A B = a B 1 , a B 2 , , a B n and A W = a 1 W , a 2 W , , a n W [50]. Next, the optimal weight values are obtained by constructing the following nonlinear programming model.
min ζ s t w B w j a B j ζ w j w W a j W ζ j = 1 n w j = 1 w j 0 , j = 1 , 2 , , n
In this paper, we consider the reality that a set of attributes expands over time. Therefore, when determining the attribute weights in different periods, we first use the BWM method c W to calculate the weight of each attribute in the current period tp that contains the most comprehensive attributes. Then, other periods are regarded as attribute subsets of the tp period. For the acquisition of their attribute weights, we remove the attributes added in the tp period and normalization, thus obtaining the attribute weight vector W ˜ k = w ˜ 1 k , w ˜ 2 k , , w ˜ n k of the tk (k = 1, 2, …, p) period based on the BWM method.
Step 2. Obtain the objective weights of attributes by the entropy-weight method. Entropy is a measure of system disorder in the entropy-weight method. In addition, the lower the entropy value of a set of values, the greater the change in data, and the greater the corresponding entropy weight [51]. The specific calculation steps of this method are
e j = 1 ln m i = 1 m D i j D j ln D i j D j   ( j = 1 , 2 , , n )
where D i j = s = 1 m d v i j , v s j   ( i = 1 , 2 , , m ; j = 1 , 2 , , n ) indicates the deviation of object A i from all other objects, and D j = i = 1 m D i j = i = 1 m s = 1 m d v i j , v s j   ( j = 1 , 2 , , n ) indicates the deviation of all objects from other objects for attribute c j . d v i j , v s j indicates the distance between v i j and v s j . Then, the entropy weight w j of each attribute can be expressed as
w j = 1 e j j = 1 n 1 e j   ( j = 1 , 2 , , n )
and then, we can obtain the attribute weight vector W k = w 1 k , w 2 k , , w n k in the t k   ( k = 1 , 2 , , p ) period based on the entropy-weight method.
Step 3. Calculate the comprehensive weights of attributes. Through BWM and entropy-weight methods, the attribute weight vector of each period can be acquired, respectively. Then we combine the results of the two methods to obtain the attribute weight of the t k ( k = 1 , 2 , , p ) period as
w j k = w ˜ j k + w j k j = 1 n w ˜ j k + w j k
Then, we can obtain the comprehensive weight vector W k = w 1 k , w 2 k , , w n k in t k ( k = 1 , 2 , , p ) .

4.2. Dynamic Hybrid MA3WD for Multiple Periods

In this section, a dynamic hybrid MA3WD model is established. First, a weighted average operator is employed to integrate the comprehensive loss functions of the alternative A i in each period t k shown in Table 5, to obtain the overall loss function of the alternative A i as shown in Table 6. Then, the overall conditional probability of the alternative A i is obtained in the same way. Finally, the three-way decision rules are established.
Table 6 shows that each loss λ ( = P , B , N ; = P , N ) in the overall loss function is obtained by integrating the corresponding losses for p periods. We calculate the thresholds according to the overall loss function, yielding:
α i = k j ξ k w j k 1 σ j k 1 r i j k + k j ξ k w j k 1 σ j k 1 r i j k + + k j ξ k w j k σ j k 1 r i j k
β i = k j ξ k w j k σ j k 1 r i j k + k j ξ k w j k σ j k 1 r i j k + + k j ξ k w j k 1 σ j k 1 r i j k
γ i = k j ξ k w j k 1 r i j k + k j ξ k w j k 1 r i j k + + k j ξ k w j k 1 r i j k
Through the process described in Section 3.2, we can calculate the conditional probability in each independent period Pr k C | [ A i ] . Considering the information from multiple periods, the overall conditional probability Pr C | [ A i ] of A i can be obtained by integrating the corresponding conditional probability Pr k C | [ A i ] for k = 1 , 2 , , p . We denote the overall conditional probability as
Pr C | [ A i ] = k = 1 p ξ k Pr k C | [ A i ]
Then, the three-way decision rules can be described as
( P )   If   k = 1 p ξ k Pr k C | [ A i ] α i ,   decide   A i P O S ( C ) ;
( B )   If   β i < k = 1 p ξ k Pr k C | [ A i ] < α i ,   decide   A i B N D ( C ) ;
( N )   If   k = 1 p ξ k Pr k C | [ A i ] β i ,   decide   A i N E G ( C ) .
where α i and β i are expressed as Equations (37) and (38). In the elderly-service performance evaluation application, if the object A i follows the rule (P), then A i P O S ( C ) , meaning the third-party organization has a good performance and is qualified to obtain a new contract with the government; if the object A i follows the rule (N), then A i N E G ( C ) , meaning that the third-party organization has a poor performance and the government should terminate the contract with it; and if the object A i follows the rule (B), then A i B N D ( C ) , meaning that the third-party organization has a moderate performance and it should take self-inspection and self-rectification for further evaluation.
It should be noted that the third-party social organizations following rule (B) can locate their shortcomings to a certain period and certain attributes according to the proposed model because it also contains the evaluating results for each single period and each single attribute. Therefore, in the final result, the social organizations in the boundary region can determine the certain period in which they have poor performance. Similarly, they can also determine certain attributes with poor performance by the proposed model. Then they can rectify the corresponding shortcomings accurately and strive to establish a cooperative relationship with the government again in the next assessment.

4.3. The Key Steps and Algorithm of Dynamic Hybrid MA3WD Model

To illustrate the dynamic hybrid MA3WD model, we summarize the key steps as well as the Algorithm 1: The specific algorithm for dynamic hybrid MA3WD.of the proposed model.
Step 1. Calculate the relative loss function based on gray relational analysis in the t k ( k = 1 , 2 , , p ) period according to Table 4.
Step 2. Compute the comprehensive loss function of each object A i and the thresholds in t k period by Equations (29)–(31).
Step 3. Calculate the conditional probability Pr k C | [ A i ] through Section 3.2.
Step 4. Obtain the evaluation results of A i through the new GRA-based hybrid MA3WD model.
Step 5. Calculate the overall loss function and overall conditional probability based on the loss function and conditional probability obtained from the GRA-based hybrid MA3WD model.
Step 6. Obtain the final decision result of A i according to the three-way decision rules.
Algorithm 1: The specific algorithm for dynamic hybrid MA3WD.
Input: An information system T , U , C , V , f , the risk-avoidance coefficient vector σ k = σ 1 k , σ 2 k , , σ n k for each period t k , the attribute weight vector W k = w 1 k , w 2 k , , w n k for each period t k , the time-series weight vector ξ = ξ 1 , ξ 2 , , ξ p .
Output: The classification results of objects.
Begin
for k = 1 to p do
 for i = 1 to m and j=1 to n do
 Calculate the relative loss function λ i j k = P , B , N ; = P , N according to Table 4.
 end
 for i = 1 to m and j=1 to n do
 Calculate the comprehensive loss function λ i j k = P , B , N ; = P , N by Equation (15).
 end
 for i = 1 to m do
 Calculate thresholds α i k , β i k and γ i k by Equations (29)–(31).
 end
 for i = 1 to m and j=1 to n do
 Compute conditional probability Pr k C | [ A i ] according to Section 3.2.
 for i = 1 to m do
 Determine the evaluation result of each object in the tk period in light of ( P G R A )~( N G R A ).
 end
end
for i = 1 to m and k=1 to p do
 Calculate the overall loss function λ i = P , B , N ; = P , N according to Table 6.
end
for i = 1 to m do
 Calculate thresholds α i , β i and γ i by Equations (37)–(39).
end
for i = 1 to m and k=1 to p do
 Calculate the overall conditional probability Pr C | [ A i ] by Equation (40).
end
for i = 1 to m do
 Obtain the final decision result of each object according to three-way decision rules.
end

5. Case Illustration

5.1. Example Calculation

The demand for elderly care in China is rising with the acceleration of aging, and traditional elderly-care supply by governments cannot meet the practical requirements. Therefore, the government purchase of elderly-care services has been successfully adopted in China. We take Xinyang, a city in Henan Province, as an example. Given the obvious aging problem in the local area, the Xinyang government has actively invested special funds to purchase elderly-care services. Moreover, the government’s policy is to evaluate the elderly-care service performance of social organizations once every year and decide whether to renew the contract after three years. For this realistic requirement, the GRA-hybrid MA3WD model is used to evaluate the elderly-care service performance each year, and the dynamic hybrid MA3WD model is used to solve the renewal decision after three years.

5.1.1. Evaluation Analysis for Single Period

There are eight social organizations A i ( i = 1 , 2 , , 8 ) that provide elderly-care services purchased by the government. As a period of three years comprises a decision-making cycle, the information of their latest three years (denoted as t1, t2, and t3, respectively, where t3 is the current period) is selected. For the t1 and t2 periods, the following evaluation indexes are used to evaluate the elderly-care services provided by social organizations [38,39].
① Professional training level ( c 1 ). This attribute examines the standardization of training in social organizations and the degree of specialization of the trained personnel.
② Number of personalized service items ( c 2 ). This attribute reflects the level of diversification of services provided by social organizations.
③ Information disclosure ( c 3 ). This attribute examines the ability of social organizations to timely inform the elderly about service adjustment, to ensure the elderly’s right to know about services.
④ Timeliness of service ( c 4 ). Timeliness refers to the swiftness and rapidness of time. The timeliness of providing elderly-care services means that elderly-care services can respond to the needs of service objects in time during the supply process.
⑤ Customer satisfaction ( c 5 ). This attribute is the result of the continuous strengthening of customer satisfaction, reflecting the customer’s recognition, affirmation, and trust in the service.
⑥ Communication ability ( c 6 ). This attribute is used to evaluate the ability of service personnel to communicate effectively with the elderly. With the extension of the service, social organizations provide psychological counseling and treatment to satisfy the needs of the elderly during the period t3. Therefore, in the third year, decision-makers add a new evaluation attribute c7.
⑦ Psychological treatment ability ( c 7 ). Old people often have serious psychological problems such as the lack of company of family members and the fear of death, so social organizations need to provide psychological counseling and treatment.
To verify the feasibility of the method proposed in this paper, after a discussion among experts selected from the field of elderly-care services, the forms of different indexes are determined as shown in Table 7. During the discussion, experts determine the type of each index based on attribute characteristics as well as their own experience and preferences. For example, timeliness of service ( c 4 ) and customer satisfaction ( c 5 ) are best represented on a scale of 1–5 according to expert experience, with 5 indicating strong approval, 1 indicating complete disapproval, and intermediate scores indicating different degrees. Communication ability ( c 6 ) is a qualitative index, so the linguistic variable is the most relevant tool. Furthermore, to demonstrate the applicability of the method to different information, experts were selected to rate based on their different preferences. The evaluation values of different social organizations in the three periods are shown in Table 8, Table 9 and Table 10. In this case, experts use a 7-point linguistic term-set to evaluate c 6 , denoted as
T = { t 0 = e x t r e m l y   p o o r ( E P ) , t 1 = v e r y   p o o r ( V P ) , t 2 = p o o r ( P ) , t 3 = m e d i u m ( M ) , t 4 = g o o d ( G ) , t 5 = v e r y   g o o d ( V G ) , t 6 = e x t r e m l y   g o o d ( E G ) }
Considering that c 3 is not easy to change in the short term, the evaluation value of each object under this attribute is fixed in the three periods.
In the treatment of hybrid information, the distance between the hybrid information and the ideal solution is used as a basis for decision-making in this paper, to avoid information loss due to the interconversion of heterogeneous information [37]. The determination of ideal solution and distance of real numbers, interval numbers, triangular fuzzy numbers, and intuitionistic fuzzy numbers are, respectively, referred to [34,35,36,52], where for these basic concepts we do not elaborate too much. The linguistic term-set T = { t k | k = 1 , 2 , , s } is transformed into triangular fuzzy numbers to participate in the calculation, with the formula t ˜ = ( t k L , t k M , t k U ) = ( max { ( k 1 ) / s , 0 } , k / s , min { ( k + 1 ) / s , 1 } ) , k = 1 , 2 , , s .
Considering the BWM and entropy-weight method, we can calculate the attribute weight vectors as W 1 = ( 0.211 , 0.162 , 0.129 , 0.192 , 0.132 , 0.175 ) , W 2 = ( 0.222 , 0.146 , 0.106 , 0.255 , 0.155 , 0.115 ) , and W 3 = ( 0.297 , 0.088 , 0.096 , 0.146 , 0.136 , 0.104 , 0.132 ) for t 1 , t 2 , and t 3 periods, respectively. According to the uncertain attitude of decision-makers towards different attributes, the risk-avoidance coefficient vector in the t 1 and t 2 periods is σ 1 , 2 = 0.35 , 0.45 , 0.35 , 0.4 , 0.4 , 0.5 . As a new attribute is added in the t 3 period, the risk-avoidance coefficient vector of the period becomes σ 3 = 0.35 , 0.45 , 0.35 , 0.4 , 0.4 , 0.5 , 0.4 . We follow the steps below to evaluate these eight social organizations (the calculation process of the t 2 or t 3 periods is the same as that of the t 1 period; we will not repeat the calculation steps, and only take the t 1 period as an example).
In the evaluation of third-party social organizations providing elderly-care services, there are two states Ω = { C , ¬ C } , where C denotes a good organization and ¬ C denotes a bad organization. Let the set of actions be represented as A = { a P , a B , a N } , where a P , a B , and a N denote good performance, moderate performance, and poor performance, respectively. According to the discussion in Section 3, we first calculate the gray relational coefficients between the eight objects and PIS and NIS in the t1 period by Equations (9) and (10). The relative loss functions of each evaluation value in the t1 period are obtained as shown in Table 11. Then, we use the attribute weight vector W 1 to aggregate the relative loss functions to obtain the comprehensive loss functions of each object as shown in Table 12. Next, we can calculate the conditional probability of each object using the GRA-TOPSIS method mentioned in Section 3.2. Meanwhile, the thresholds of each object can be obtained by Equations (29)–(31). The thresholds and conditional probabilities are listed in Table 13.
For each object, the range of each region and the conditional probability are shown in Figure 2. The region where Pr 1 ( C | [ A i ] ) falls represents the decision of A i in the t 1 period, so we can obtain three-way decision results clearly from Figure 2 as: P O S ( C ) = { A 7 , A 8 } , B N D ( C ) = { A 3 , A 4 , A 5 , A 6 } and N E G ( C ) = { A 1 , A 2 } . From the results, in the t 1 period, social organizations A 7 and A 8 are divided into the positive region, indicating that these two meet the standard of eligibility; A 3 , A 4 , A 5 and A 6 are divided into the boundary region, meaning that they perform moderately and should strive to improve their business in the next period; A 1 and A 2 are divided into the negative region, indicating these two perform poorly and may face the risk of elimination in the final decision if they do not make a change.
As the calculation process of t 2 and t 3 is the same as t 1 , the calculation results of t 2 and t 3 are directly given, which are listed in Table 14.
From the 3WD results of Table 15, we can observe the different results of the t 2 and t 3 periods. In the t 2 period, social organizations A 1 , A 3 , A 6 , A 7 , and A 8 are classified into the positive region, meaning that the five objects perform well; A 4 is classified into the negative region, indicating poor performance during this period; A 2 and A 5 are classified into the boundary region, which means they have mediocre performance and need further improvement to reach the standard of renewal. Similarly, in the t 3 period, social organizations A 3 , A 5 , A 6 , and A 7 perform well; A 4 and A 8 perform poorly; and A 1 and A 2 perform mediocrely, and need to work hard to improve business.

5.1.2. Decision Analysis for Multiple Periods

According to the process of the GRA-based hybrid MA3WD model, the corresponding comprehensive loss function and the conditional probability of each object of the three periods are obtained in Section 5.1.1. Since the decision-maker attaches more importance to the recent data, the value of the time degree λ is set to 0.3, and then the time-series weight vector is obtained according to Equation (32) as follows: ξ = ( 0.154 , 0.292 , 0.554 ) . Next, use the time-series weights to integrate the comprehensive loss functions of the t 1 , t 2 and t 3 periods to obtain the overall loss functions as shown in Table 16. It should be noted that the meaning of the actions set A = { a P , a B , a N } in Table 16 is slightly different from that mentioned in Section 5.1.1. Here, a P , a B and a N denote the renewal of contract, further investigation, and termination of contract, respectively. Like the method used to calculate the overall loss functions, we integrate the corresponding conditional probability of each object in three periods to obtain the overall conditional probability Pr ( C | [ A i ] ) of each object. The thresholds and overall conditional probabilities are listed in Table 17.
From Figure 3, we can visually observe the final 3WD results after comprehensively considering the evaluation information of three periods. Now, according to the classification results, the government can decide on which social organizations to continue to cooperate with. As shown in Figure 3, P O S ( C ) = { A 3 , A 5 , A 6 , A 7 } , B N D ( C ) = { A 1 , A 8 } and N E G ( C ) = { A 2 , A 4 } . Therefore, social organizations A 3 , A 5 , A 6 , and A 7 have good performances and are qualified to obtain new contracts with the government; A 2 and A 4 have poor performances and government should terminate contracts with them; A 1 and A 8 have moderate performances and they should take self-inspection and self-rectification for further evaluation.
As mentioned in Section 4, the model proposed in this paper can precisely locate which attributes are underperforming at which period for the social organization, to rectify and manage the corresponding business. Hence, for the social organizations in the boundary region, they can make targeted rectifications according to the results of the model, and perform outstandingly in the next assessment. For example, in the final decision results, A 8 is divided into the boundary region. According to Table 13 and Table 15, A 8 is divided into the positive region, while in the t 3 period, it is divided into the negative region. Based on this situation, for A 8 , we can explore which attributes underperform at t 3 , i.e., which attributes are divided into the negative region. Considering a single attribute as a 3WD problem, the 3WD rules of A8 under each attribute are obtained based on the 3WD method proposed in this paper.
According to the thresholds and conditional probabilities listed in Table 18, we can judge that A8 is classified into the rejection region under attributes c 1 , c 2 , c 3 , c 6 , and c 7 , which means A 8 performs worse in these five attributes. From this result, if A 8 wants to improve its service performance, it should rectify the business related to these five attributes.

5.2. Comparative Analysis

5.2.1. Comparison between Static and Dynamic Assessment

To illustrate the importance of dynamic evaluation, we compare the results of the 3WD based on the evaluation information in the t 3 (current) period with the results of the 3WD based on the comprehensive information, i.e., when t = t 3 and t = t 1 t 3 . The comparison results are shown in Table 19.
From Table 19, it can be seen that there are differences in the results derived from static and dynamic assessment. However, the result based on dynamic assessment is more reasonable. For example, in the t 3 period, A 8 is divided into the negative region, but if the evaluation information of the t 1 t 3 periods is considered, A 8 is classified in the boundary region. By comparing the decision results of t 1 and t 2 , we can find that A 8 performs well in the first two periods and is classified into the positive region. This suggests that A 8 only performs poorly in the t 3 period, but still has the potential to be renewed and should be given a chance to rectify for the next assessment. Based on the above discussion, integrating multiple periods leads to more credible classification results, which demonstrates the importance of dynamic assessment.

5.2.2. Comparison between the Proposed Method and MADM Methods

In the current section, we compare the proposed method with three traditional MADM methods (the TOPSIS method [53], the VIKOR method [54], and the ELECTRE method [55]). It should be noted that the 3WD method is about the classification of alternatives, while these three MADM methods are only for obtaining the ranking of alternatives and determining an optimal alternative. To compare with these three MADM methods, the ranking of the alternatives is obtained according to the thresholds α , β , and γ , respectively. The smaller the value of α (or β or γ ) of the alternative, the higher the probability of the alternative being accepted [20], so we obtain three sorting results according to the values of α , β , and γ , which are shown in Table 20.
From Table 20, it can be seen that the ranking results of different methods are highly consistent, although there are slight differences. More specifically, there is a difference in the sorting of A 3 , A 5 , and A 6 for all methods, which indicates that the proposed method is valid. In addition, compared with the MADM method, the method proposed in this paper considers the risk preference of the decision-maker and divides all social organizations into three areas, which effectively reduces decision-making risks and helps the decision-maker make the most reasonable decision.

5.2.3. Comparison between the Proposed Method and Existing 3WD Methods

To further illustrate the superiority of the proposed model, we select Jia and Liu’s method [20], Gao et al.’s method [13], and Wang et al.’s method [56] for comparison. These three methods are first applied to the example in this section to obtain the classification results, which are listed in Table 21. Then, we compare and analyze the similarities and differences between the proposed method and the above methods considering four aspects, which are listed in Table 22.
As shown in Table 21, the classification results of the four methods are generally consistent. A 3 , A 5 , A 6 , and A 7 are all classified into the accept region, A 2 and A 4 are both classified into the reject region. There are only differences as to whether A 1 and A 8 are classified in the boundary region or the negative region. Based on the comparison of classification results, the effectiveness of our proposed method can be verified. Next, through the comparison in Table 22, we summarize the advantages of the method proposed in this paper:
(1)
Jia and Liu [20] converted attribute values into loss functions using relative loss and inverse loss functions, which is a great advance on the 3WD model. The determination of conditional probabilities is subjectively given by the decision-maker and lacks interpretability. The method proposed in this paper uses GRA-TOPSIS to estimate conditional probabilities, which overcomes the subjective influence of conditional probabilities given artificially.
(2)
Gao et al.’s method [13] considers the influence of time factors on realistic decision problems and considers the integration of information from multiple periods to make decisions. In some decision-making problems, it is also necessary to evaluate objects in a certain period. Gao et al. [13] lack the evaluation of objects in a single period. In contrast, the model proposed in this paper not only obtains the final decision results but also obtains the results of a certain period and a certain attribute, which facilitates the object to accurately determine which attribute is at a disadvantage for rectification. From this perspective, the proposed model is superior because of its flexibility and universality in the presentation of results.
(3)
Wang et al.’s method [56] introduces regret theory into the 3WD process, which is a great improvement to the 3WD and provides a guiding direction for our future research work. However, both the attribute weights and outcome matrix are decided subjectively by the decision-maker, which lacks transparency and interpretability. The method proposed in this paper uses a combination of BWM and entropy-weight methods to determine the attribute weights, which is more scientific than the method of Wang et al. [56]. At the same time, the proposed model uses GRA to construct the loss functions, which effectively connects the attribute values in the MADM with the loss functions in the 3WD from an objective perspective.

5.2.4. Correlation Analysis

In this section, we introduce the pairwise comparison method Spearman’s correlation coefficient (SCC) [57] to explore the connection of the above methods. There is a consensus that a larger SCC indicates a stronger consistency between the two methods. In general, an SCC greater than 0.8 indicates a strong correlation between the two counterparts. The specific calculation results are shown in Figure 4.
The results in Figure 4 show that all the SCCs between the proposed method and the other selected methods are greater than 0.9; therefore, there is a conclusion that the results derived from the proposed method are valid and credible.

6. Conclusions

In this paper, a GRA-based dynamic hybrid MA3WD model is established, and the main contributions of the method are listed as follows: (1) a new loss function construction method is proposed based on GRA, which adopts a data-driven method to overcome the drawback of excessive subjectivity caused by decision-makers giving loss functions based on experience; (2) with the aid of GRA-TOPSIS method, a new conditional probability determination method is proposed, which is based on the relative closeness reflecting the position and shape similarity; (3) a GRA-based hybrid MA3WD model is proposed for evaluating objects at a specific period. The model can point out the specific attributes and periods of poor performance of the object, which facilitates the object to improve its deficiencies accurately; (4) by extending the single-period environment to a multi-period environment, a dynamic hybrid MA3WD model is proposed, which extends the study of 3WD in a time-dynamic environment.
The proposed model is applied to solve the performance evaluation problem of elderly-care services to illustrate the feasibility and superiority of the proposed method, and the applicability of the proposed method is as follows: (1) the model applies to realistic problems with hybrid information derived from single-period evaluation information as well as multi-period information, which is widely applicable and flexible; (2) the proposed model can be used not only in the performance evaluation of elderly-care services but also in similar decision-making problems in other fields, such as target–threat assessment [13], investment decisions [30], risk analysis [58], medical diagnosis [28], etc.
There are some limitations, namely: (1) some existing 3WD methods have considered the psychological behavior of decision-makers [25], which is not mentioned in this paper, and the proposed method can be combined with regret theory [21] in future work; (2) in the dynamic MA3WD model, we calculate the time-series weights based on the entropy-weight method. Although this method is classical, it can be improved to obtain more reasonable time-series weights. Therefore, we will focus on three aspects in future studies: (1) the proposed method should be combined with group decision-making [59] to make the evaluation information used in decision-making more accurate and reliable; (2) considering the advantages of machine learning, we will introduce the knowledge of machine learning [60] in future studies on 3WD; (3) based on the characteristic that decision-makers are rational, we will consider combining the proposed method with behavioral decision theories such as regret theory [21] and prospect theory [61].

Author Contributions

Conceptualization, F.J.; data curation, Y.W.; formal analysis, Y.S.; funding acquisition, Y.S.; methodology, F.J.; validation, Y.W.; writing—original draft, Y.W.; writing—review and editing, F.J. and Y.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Humanities and Social Sciences Research Project of the Ministry of Education of China (No. 19YJC630059), the Shandong Social Science Planning Fund Program (No. 18DGLJ14), the Natural Science Foundation of Shandong Province (No. ZR2019PG009), and the Higher Educational Social Science Program of Shandong Province (no. J18RA112).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rawat, S.S.; Pant, S.; Kumar, A.; Ram, M.; Sharma, H.K.; Kumar, A. A state-of-the-art survey on analytical hierarchy process applications in sustainable development. Sciences 2022, 7, 883–917. [Google Scholar] [CrossRef]
  2. Rouyendegh, B.D.; Topuz, K.; Dag, A.; Oztekin, A. An AHP-IFT Integrated Model for Performance Evaluation of E-Commerce Web Sites. Inf. Syst. Front. 2019, 21, 1345–1355. [Google Scholar] [CrossRef]
  3. Liu, P.D.; Shen, M.J. Failure Mode and Effects Analysis (FMEA) for Traffic Risk Assessment Based on Unbalanced Double Hierarchy Linguistic Term Set. Int. J. Fuzzy Syst. 2022, 25, 423–450. [Google Scholar] [CrossRef]
  4. Wu, J.; Gong, H.; Liu, F.; Liu, Y. Risk Assessment of Open-Pit Slope Based on Large-Scale Group Decision-Making Method Considering Non-Cooperative Behavior. Int. J. Fuzzy Syst. 2023, 25, 245–263. [Google Scholar] [CrossRef]
  5. Wang, Z.; Xiao, F.Y.; Ding, W.P. Interval-valued intuitionistic fuzzy jenson-shannon divergence and its application in multi-attribute decision making. Appl. Intell. 2022, 52, 16168–16184. [Google Scholar] [CrossRef]
  6. Jiang, H.B.; Zhan, J.M.; Sun, B.Z.; Alcantud, J.C.R. An MADM approach to covering-based variable precision fuzzy rough sets: An application to medical diagnosis. Int. J. Mach. Learn. Cybern. 2020, 11, 2181–2207. [Google Scholar] [CrossRef]
  7. Zhang, L.H.; Zhao, Z.L.; Yang, M.; Li, S.R. A multi-criteria decision method for performance evaluation of public charging service quality. Energy 2020, 195, 116958. [Google Scholar] [CrossRef]
  8. dos Santos, B.M.; Godoy, L.P.; Campos, L.M.S. Performance evaluation of green suppliers using entropy-TOPSIS-F. J. Clean. Prod. 2019, 207, 498–509. [Google Scholar] [CrossRef]
  9. Gong, J.W.; Liu, H.C.; You, X.Y.; Yin, L.S. An integrated multi-criteria decision making approach with linguistic hesitant fuzzy sets for E-learning website evaluation and selection. Appl. Soft Comput. 2021, 102, 107118. [Google Scholar] [CrossRef]
  10. Luo, S.Z.; Liang, W.Z. A hybrid TODIM approach with unknown weight information for the performance evaluation of cleaner production. Comput. Appl. Math. 2021, 40, 23. [Google Scholar] [CrossRef]
  11. Lam, W.S.; Lam, W.H.; Jaaman, S.H.; Liew, K.F. Performance Evaluation of Construction Companies Using Integrated Entropy-Fuzzy VIKOR Model. Entropy 2021, 23, 320. [Google Scholar] [CrossRef] [PubMed]
  12. Lu, M.T.; Hsu, C.C.; Liou, J.J.H.; Lo, H.W. A hybrid MCDM and sustainability-balanced scorecard model to establish sustainable performance evaluation for international airports. J. Air Transp. Manag. 2018, 71, 9–19. [Google Scholar] [CrossRef]
  13. Gao, Y.; Li, D.S.; Zhong, H. A novel target threat assessment method based on three-way decisions under intuitionistic fuzzy multi-attribute decision making environment. Eng. Appl. Artif. Intell. 2020, 87, 103276. [Google Scholar] [CrossRef]
  14. Jiang, H.B.; Hu, B.Q. A decision-theoretic fuzzy rough set in hesitant fuzzy information systems and its application in multi-attribute decision-making. Inf. Sci. 2021, 579, 103–127. [Google Scholar] [CrossRef]
  15. Yao, Y.Y. Three-way decisions with probabilistic rough sets. Inf. Sci. 2010, 180, 341–353. [Google Scholar] [CrossRef] [Green Version]
  16. Yao, Y.Y. The superiority of three-way decisions in probabilistic rough set models. Inf. Sci. 2011, 181, 1080–1096. [Google Scholar] [CrossRef]
  17. Li, H.X.; Zhang, L.B.; Huang, B.; Zhou, X.Z. Sequential three-way decision and granulation for cost-sensitive face recognition. Knowl.-Based Syst. 2016, 91, 241–251. [Google Scholar] [CrossRef]
  18. Mandal, P.; Ranadive, A.S. Multi-granulation fuzzy probabilistic rough sets and their corresponding three-way decisions over two universes. Iran. J. Fuzzy Syst. 2019, 16, 61–76. [Google Scholar] [CrossRef]
  19. Xu, Y.; Tang, J.X.; Wang, X.S. Three sequential multi-class three-way decision models. Inf. Sci. 2020, 537, 62–90. [Google Scholar] [CrossRef]
  20. Jia, F.; Liu, P.D. A novel three-way decision model under multiple-criteria environment. Inf. Sci. 2019, 471, 29–51. [Google Scholar] [CrossRef]
  21. Zhu, J.X.; Ma, X.L.; Zhan, J.M. A regret theory-based three-way decision approach with three strategies. Inf. Sci. 2022, 595, 89–118. [Google Scholar] [CrossRef]
  22. Zhan, J.M.; Ye, J.; Ding, W.P.; Liu, P.D. A Novel Three-Way Decision Model Based on Utility Theory in Incomplete Fuzzy Decision Systems. IEEE Trans. Fuzzy Syst. 2022, 30, 2210–2226. [Google Scholar] [CrossRef]
  23. Zhan, J.M.; Jiang, H.B.; Yao, Y.Y. Three-Way Multiattribute Decision-Making Based on Outranking Relations. IEEE Trans. Fuzzy Syst. 2021, 29, 2844–2858. [Google Scholar] [CrossRef]
  24. Huang, X.F.; Zhan, J.M.; Sun, B.Z. A three-way decision method with pre-order relations. Inf. Sci. 2022, 595, 231–256. [Google Scholar] [CrossRef]
  25. He, S.F.; Wang, Y.M.; Pan, X.H.; Chin, K.S. A novel behavioral three-way decision model with application to the treatment of mild symptoms of COVID-19. Appl. Soft Comput. 2022, 124, 109055. [Google Scholar] [CrossRef]
  26. Sun, B.Z.; Chen, X.T.; Zhang, L.Y.; Ma, W.M. Three-way decision making approach to conflict analysis and resolution using probabilistic rough set over two universes. Inf. Sci. 2020, 507, 809–822. [Google Scholar] [CrossRef]
  27. Yao, Y.Y. Three-way conflict analysis: Reformulations and extensions of the Pawlak model. Knowl.-Based Syst. 2019, 180, 26–37. [Google Scholar] [CrossRef]
  28. Chen, Y.F.; Yue, X.D.; Fujita, H.; Fu, S.Y. Three-way decision support for diagnosis on focal liver lesions. Knowl.-Based Syst. 2017, 127, 85–99. [Google Scholar] [CrossRef]
  29. Li, Z.W.; Xie, N.X.; Huang, D.; Zhang, G.Q. A three-way decision method in a hybrid decision information system and its application in medical diagnosis. Artif. Intell. Rev. 2020, 53, 4707–4736. [Google Scholar] [CrossRef]
  30. Jiang, H.B.; Hu, B.Q. A novel three-way group investment decision model under intuitionistic fuzzy multi-attribute group decision-making environment. Inf. Sci. 2021, 569, 557–581. [Google Scholar] [CrossRef]
  31. Li, X.; Huang, X.J. A Novel Three-Way Investment Decisions Based on Decision-Theoretic Rough Sets with Hesitant Fuzzy Information. Int. J. Fuzzy Syst. 2020, 22, 2708–2719. [Google Scholar] [CrossRef]
  32. Zhi, H.L.; Qi, J.J.; Qian, T.; Wei, L. Three-way dual concept analysis. Int. J. Approx. Reason. 2019, 114, 151–165. [Google Scholar] [CrossRef]
  33. He, X.L.; Wei, L.; She, Y.H. L-fuzzy concept analysis for three-way decisions: Basic definitions and fuzzy inference mechanisms. Int. J. Mach. Learn. Cybern. 2018, 9, 1857–1867. [Google Scholar] [CrossRef]
  34. Wang, Y.J. Combining technique for order preference by similarity to ideal solution with relative preference relation for interval-valued fuzzy multi-criteria decision-making. Soft Comput. 2020, 24, 11347–11364. [Google Scholar] [CrossRef]
  35. Sadabadi, S.A.; Hadi-Vencheh, A.; Jamshidi, A.; Jalali, M. An Improved Fuzzy TOPSIS Method with a New Ranking Index. Int. J. Inf. Technol. Decis. Mak. 2022, 21, 615–641. [Google Scholar] [CrossRef]
  36. Reig-Mullor, J.; Brotons-Martinez, J.M. The evaluation performance for commercial banks by intuitionistic fuzzy numbers: The case of Spain. Soft Comput. 2021, 25, 9061–9075. [Google Scholar] [CrossRef]
  37. Liu, Z.M.; Wang, D.; Wang, W.X.; Liu, P.D. An integrated group decision-making framework for selecting cloud service providers based on regret theory and EVAMIX with hybrid information. Int. J. Intell. Syst. 2022, 37, 3480–3513. [Google Scholar] [CrossRef]
  38. Meng, F.; Chen, B.; Tan, C. Adaptive minimum adjustment consensus model for large-scale group decision making under social networks and its application in Integrated Care of Older People. Appl. Soft Comput. 2023, 132, 109863. [Google Scholar] [CrossRef]
  39. Wang, C.; Ding, S.J. The Quality Evaluation of the Home-care Service for the Elderly Purchased by Government: Multiple Conception, Index Construction and Application Case. Popul. Econ. 2018, 4, 12–20. (In Chinese) [Google Scholar]
  40. Deng, J.L. Introduction to grey system. J. Grey Syst. 1989, 1, 1–24. [Google Scholar]
  41. Kuo, Y.; Yang, T.; Huang, G.-W. The use of grey relational analysis in solving multiple attribute decision-making problems. Comput. Ind. Eng. 2008, 55, 80–93. [Google Scholar] [CrossRef]
  42. Liu, Y.; Du, J.L.; Wang, Y.H. An improved grey group decision-making approach. Appl. Soft Comput. 2019, 76, 78–88. [Google Scholar] [CrossRef]
  43. Sun, B.Z.; Ma, W.M.; Li, B.J.; Li, X.N. Three-way decisions approach to multiple attribute group decision making with linguistic information-based decision-theoretic rough fuzzy set. Int. J. Approx. Reason. 2018, 93, 424–442. [Google Scholar] [CrossRef]
  44. Li, C.C.; Dong, Y.C.; Liang, H.M.; Pedrycz, W.; Herrera, F. Data-driven method to learning personalized individual semantics to support linguistic multi-attribute decision making. Omega-Int. J. Manag. Sci. 2022, 111, 102642. [Google Scholar] [CrossRef]
  45. Liang, D.C.; Xu, Z.S.; Liu, D.; Wu, Y. Method for three-way decisions using ideal TOPSIS solutions at Pythagorean fuzzy information. Inf. Sci. 2018, 435, 282–295. [Google Scholar] [CrossRef]
  46. Liu, D.; Qi, X.C.; Fu, Q.; Li, M.; Zhu, W.F.; Zhang, L.L.; Faiz, M.A.; Khan, M.I.; Li, T.X.; Cui, S. A resilience evaluation method for a combined regional agricultural water and soil resource system based on Weighted Mahalanobis distance and a Gray-TOPSIS model. J. Clean. Prod. 2019, 229, 667–679. [Google Scholar] [CrossRef]
  47. Lu, H.A.; Zhao, Y.M.; Zhou, X.; Wei, Z.K. Selection of Agricultural Machinery Based on Improved CRITIC-Entropy Weight and GRA-TOPSIS Method. Processes 2022, 10, 266. [Google Scholar] [CrossRef]
  48. Yang, Y.Y.; Guo, H.X.; Chen, L.F.; Liu, X.; Gu, M.Y.; Ke, X.L. Regional analysis of the green development level differences in Chinese mineral resource-based cities. Resour. Policy 2019, 61, 261–272. [Google Scholar] [CrossRef]
  49. Rezaei, J. Best-worst multi-criteria decision-making method. Omega 2015, 53, 49–57. [Google Scholar] [CrossRef]
  50. Liu, S.; Hu, Y.N.; Zhang, X.; Li, Y.F.; Liu, L. Blockchain Service Provider Selection Based on an Integrated BWM-Entropy-TOPSIS Method Under an Intuitionistic Fuzzy Environment. IEEE Access 2020, 8, 104148–104164. [Google Scholar] [CrossRef]
  51. Han, Z.S.; Liu, P.D. A fuzzy multi-attribute decision-making method under risk with unknown attribute weights. Technol. Econ. Dev. Econ. 2011, 17, 246–258. [Google Scholar] [CrossRef]
  52. Wang, Y.M.; Liu, P.D.; Yao, Y.Y. BMW-TOPSIS: A generalized TOPSIS model based on three-way decision. Inf. Sci. 2022, 607, 799–818. [Google Scholar] [CrossRef]
  53. Xun, X.L.; Yuan, Y.B. Research on the urban resilience evaluation with hybrid multiple attribute TOPSIS method: An example in China. Nat. Hazards 2020, 103, 557–577. [Google Scholar] [CrossRef] [PubMed]
  54. Liu, H.C.; Mao, L.X.; Zhang, Z.Y.; Li, P. Induced aggregation operators in the VIKOR method and its application in material selection. Appl. Math. Model. 2013, 37, 6325–6338. [Google Scholar] [CrossRef]
  55. Lin, R.J.; Lu, S.F.; Yang, A.; Shen, W.F.; Ren, J.Z. Multi-criteria sustainability assessment and decision-making framework for hydrogen pathways prioritization: An extended ELECTRE method under hybrid information. Int. J. Hydrog. Energy 2021, 46, 13430–13445. [Google Scholar] [CrossRef]
  56. Wang, T.X.; Li, H.X.; Qian, Y.H.; Huang, B.; Zhou, X.Z. A Regret-Based Three-Way Decision Model Under Interval Type-2 Fuzzy Environment. IEEE Trans. Fuzzy Syst. 2022, 30, 175–189. [Google Scholar] [CrossRef]
  57. Pamucar, D.; Petrovic, I.; Cirovic, G. Modification of the Best-Worst and MABAC methods: A novel approach based on interval-valued fuzzy-rough numbers. Expert Syst. Appl. 2018, 91, 89–106. [Google Scholar] [CrossRef]
  58. Duan, C.Y.; Chen, X.Q.; Shi, H.; Liu, H.C. A New Model for Failure Mode and Effects Analysis Based on k-Means Clustering Within Hesitant Linguistic Environment. IEEE Trans. Eng. Manag. 2022, 69, 1837–1847. [Google Scholar] [CrossRef]
  59. Ma, X.Y.; Zhao, M.; Zou, X. Measuring and reaching consensus in group decision making with the linguistic computing model based on discrete fuzzy numbers. Appl. Soft Comput. 2019, 77, 135–154. [Google Scholar] [CrossRef]
  60. Subhashini, L.; Li, Y.F.; Zhang, J.L.; Atukorale, A.S. Integration of fuzzy logic and a convolutional neural network in three-way decision-making. Expert Syst. Appl. 2022, 202, 117103. [Google Scholar] [CrossRef]
  61. Zhan, J.M.; Wang, J.J.; Ding, W.P.; Yao, Y.Y. Three-Way Behavioral Decision Making with Hesitant Fuzzy Information Systems: Survey and Challenges. IEEE-CAA J. Autom. Sin. 2023, 10, 330–350. [Google Scholar] [CrossRef]
Figure 1. Flow chart of the proposed model.
Figure 1. Flow chart of the proposed model.
Mathematics 11 03176 g001
Figure 2. Three-way decision in the t 1 period.
Figure 2. Three-way decision in the t 1 period.
Mathematics 11 03176 g002
Figure 3. Three-way decision of each object.
Figure 3. Three-way decision of each object.
Mathematics 11 03176 g003
Figure 4. The SCC of different decision-making methods.
Figure 4. The SCC of different decision-making methods.
Mathematics 11 03176 g004
Table 1. The loss functions of a three-way decision.
Table 1. The loss functions of a three-way decision.
C ¬ C
a P λ P P λ P N
a B λ B P λ B P
a N λ N P λ N N
Table 2. Relative loss function derived from the evaluation value.
Table 2. Relative loss function derived from the evaluation value.
C j ¬ C j
a P 0 d x i j , P I S j
a B σ j d x i j , N I S j σ j d x i j , P I S j
a N d x i j , N I S j 0
Table 3. Comprehensive loss function.
Table 3. Comprehensive loss function.
C j ¬ C j
a P 0 w j d x i j , P I S j
a B w j σ j d x i j , N I S j w j σ j d x i j , P I S j
a N w j d x i j , N I S j 0
Table 4. Relative loss functions in the t k period.
Table 4. Relative loss functions in the t k period.
A i C j k ¬ C j k
a P 0 1 r i j k +
a B σ j k 1 r i j k σ j k 1 r i j k +
a N 1 r i j k 0
Table 5. Comprehensive loss function of A i in the t k period.
Table 5. Comprehensive loss function of A i in the t k period.
A i C ¬ C
a P 0 j w j k 1 r i j k +
a B j w j k σ j k 1 r i j k j w j k σ j k 1 r i j k +
a N j w j k 1 r i j k + 0
Table 6. Overall loss function of alternative A i .
Table 6. Overall loss function of alternative A i .
A i C ¬ C
a P 0 k j ξ k w j k 1 r i j k +
a B k j ξ k w j k σ j k 1 r i j k k j ξ k w j k σ j k 1 r i j k +
a N k j ξ k w j k 1 r i j k 0
Table 7. Classification of evaluation indexes.
Table 7. Classification of evaluation indexes.
Evaluation IndexesIndex Value FormsIndex Type
c 1 Triangular fuzzy numbersQuantitative (benefit type)
c 2 Real numbersQuantitative (benefit type)
c 3 Intuitionistic fuzzy numbersQuantitative (benefit type)
c 4 Real numbersQuantitative (benefit type)
c 5 Real numbersQuantitative (benefit type)
c 6 Linguistic termsQualitative (benefit type)
c 7 Interval numbersQuantitative (benefit type)
Table 8. Decision matrix for the t 1 period.
Table 8. Decision matrix for the t 1 period.
c 1 c 2 c 3 c 4 c 5 c 6
A 1 (0.64,0.7,0.76)7<0.8,0.1>43G
A 2 (0.6,0.68,0.76)8<0.75,0.15>33P
A 3 (0.7,0.75,0.8)9<0.85,0.1>34VG
A 4 (0.8,0.85,0.9)8<0.8,0.1>43VG
A 5 (0.72,0.8,0.88)6<0.85,0.1>44G
A 6 (0.65,0.7,0.75)7<0.8,0.1>54M
A 7 (0.81,0,85,0.89)10<0.8,0.1>55G
A 8 (0.75,0.85,0.95)7<0.75,0.15>45EG
Table 9. Decision matrix for the t 2 period.
Table 9. Decision matrix for the t 2 period.
c 1 c 2 c 3 c 4 c 5 c 6
A 1 (0.68,0.72,0.76)7<0.8,0.1>54G
A 2 (0.65,0.70,0.75)8<0.75,0.15>45M
A 3 (0.8,0.85,0.9)8<0.85,0.1>44VG
A 4 (0.61,0.66,0.71)8<0.8,0.1>33G
A 5 (0.72,0.77,0.82)6<0.85,0.1>45VG
A 6 (0.75,0.8,0.85)7<0.8,0.1>54G
A 7 (0.82,0.87,0.91)10<0.8,0.1>45EG
A 8 (0.73,0.83,0.93)9<0.75,0.15>45EG
Table 10. Decision matrix for the t 3 period.
Table 10. Decision matrix for the t 3 period.
c 1 c 2 c 3 c 4 c 5 c 6 c 7
A 1 (0.68,0.73,0.78)8<0.8,0.1>54G[77,81]
A 2 (0.7,0.75,0.8)8<0.75,0.15>45G[77,81]
A 3 (0.74,0.79,0.84)10<0.85,0.1>45VG[81,85]
A 4 (0.71,0.76,0.81)9<0.8,0.1>33G[75,79]
A 5 (0.7,0.75,0.8)8<0.85,0.1>55VG[81,85]
A 6 (0.8,0.84,0.9)9<0.8,0.1>54VG[82,86]
A 7 (0.9,0.95,1)10<0.8,0.1>45EG[86,90]
A 8 (0.65,0.7,0.75)8<0.75,0.15>44G[75,79]
Table 11. Relative loss functions in the t 1 period.
Table 11. Relative loss functions in the t 1 period.
c 1 c 2 c 3 c 4 c 5 c 6
C 1 ¬ C 1 C 2 ¬ C 2 C 3 ¬ C 3 C 4 ¬ C 4 C 5 ¬ C 5 C 6 ¬ C 6
A1aP00.44800.49300.10300.39400.56500.484
aB0.0430.1570.1100.2220.0490.0360.1570.15700.2260.2600.242
aN0.12400.24500.14000.3940000.5200
A2aP00.46800.39400.20400.56500.56500.667
aB0.0320.1640.1770.17700.07700.22600.22600.333
aN0.09200.394000000000
A3aP00.39400.2450000.56500.39400.306
aB0.0810.1380.2220.1100.071000.2260.1570.1570.3090.153
aN0.23300.49300.2040000.39400.6190
A4aP00.24500.39400.10300.39400.56500.306
aB0.1350.0860.1770.1770.0490.0360.1570.15700.2260.3090.153
aN0.38600.39400.14000.3940000.6190
A5aP00.33000.5650000.39400.39400.484
aB0.1120.11600.2540.07100.1570.1570.1570.1570.2600.242
aN0.3200000.20400.39400.39400.5200
A6aP00.44800.49300.1030000.39400.594
aB0.0430.1570.1100.2220.0490.0360.22600.1570.1570.1760.297
aN0.12400.24500.14000.56500.39400.3510
A7aP00.2460000103000000.484
aB0.1350.0860.25400.0490.0360.22600.22600.2600.242
aN0.38600.56500.14000.56500.56500.5200
A8aP00.26000.49300.20400.3940000
aB0.1360.0910.1100.22200.0710.1570.1570.22600.3330
aN0.39000.2450000.39400.56500.6670
Table 12. Comprehensive loss functions in the t 1 period.
Table 12. Comprehensive loss functions in the t 1 period.
C ¬ C C ¬ C
A1aP00.423A5aP00.373
aB0.1090.176aB0.1290.159
aN0.2500aN0.3120
A2aP00.489A6aP00.344
aB0.0350.204aB0.1280.146
aN0.0830aN0.3060
A3aP00.337A7aP00.150
aB0.1370.138aB0.1950.065
aN0.3160aN0.4650
A4aP00.333A8aP00.237
aB0.1480.138aB0.1650.095
aN0.3470aN0.3890
Table 13. Decision thresholds and conditional probabilities in the t 1 period.
Table 13. Decision thresholds and conditional probabilities in the t 1 period.
A1A2A3A4A5A6A7A8
α i 1 0.6930.8890.5920.5680.6240.6060.3030.463
β i 1 0.5550.8810.4360.4100.4650.4520.1940.297
γ i 1 0.6280.8550.5160.4890.5450.5290.2440.378
Pr 1 C | A i 0.4170.2680.5230.5450.4820.4880.7170.632
Table 14. Decision thresholds and conditional probabilities in the t 2 and t 3 periods.
Table 14. Decision thresholds and conditional probabilities in the t 2 and t 3 periods.
A1A2A3A4A5A6A7A8
t 2 α i 2 0.5310.6000.5170.8000.5500.4910.3320.426
β i 2 0.3600.4270.3440.6880.3710.3250.1790.253
γ i 2 0.4440.5140.4280.7500.4590.4040.2500.336
Pr 2 C | A i 0.5430.4830.5670.3110.5360.5750.6940.643
t 3 α i 3 0.5930.6010.4520.7520.4480.4370.1840.699
β i 3 0.4050.4130.2680.5920.2690.2610.0900.529
γ i 3 0.5000.5080.3560.6790.3550.3440.1300.620
Pr 3 C | A i 0.4610.4580.5770.3300.5750.5940.7090.371
Table 15. Three-way decision results in t 2 and t 3 periods.
Table 15. Three-way decision results in t 2 and t 3 periods.
POS(C)BND(C)NEG(C)
t 2 { A 1 , A 3 , A 6 , A 7 , A 8 } { A 2 , A 5 } { A 4 }
t 3 { A 3 , A 5 , A 6 , A 7 } { A 1 , A 2 } { A 4 , A 8 }
Table 16. Overall loss functions.
Table 16. Overall loss functions.
C ¬ C C ¬ C
A1aP00.319A5aP00.256
aB0.1300.130aB0.1460.105
aN0.3140aN0.3550
A2aP00.348A6aP00.246
aB0.1140.142aB0.1560.101
aN0.2740aN0.3820
A3aP00.257A7aP00.105
aB0.1560.103aB0.1860.042
aN0.3780aN0.4540
A4aP00.403A8aP00.313
aB0.0860.163aB0.1320.124
aN0.2010aN0.3140
Table 17. Decision thresholds and overall conditional probabilities.
Table 17. Decision thresholds and overall conditional probabilities.
A1A2A3A4A5A6A7A8
α i 0.5920.6450.4960.7370.5110.4820.2540.588
β i 0.4140.4680.3180.5850.3330.3100.1350.405
γ i 0.5040.5590.4050.6670.4190.3920.1880.499
Pr C | A i 0.4780.4360.5660.3570.5490.5720.7050.491
Table 18. Calculation results of A8 under each attribute.
Table 18. Calculation results of A8 under each attribute.
c 1 c 2 c 3 c 4 c 5 c 6 c 7
A 8 α 0.8700.55010.60.60.4821
β 0.6610.45010.40.40.4821
Pr ( C j | [ A 8 ] ) 0.2590.3560.2850.50.50.3460.264
Table 19. Comparison at t = t3 and t = t1 ~ t3.
Table 19. Comparison at t = t3 and t = t1 ~ t3.
POS(C)BND(C)NEG(C)
t = t 3 { A 3 , A 5 , A 6 , A 7 } { A 1 , A 2 } { A 4 , A 8 }
t = t 1 t 3 { A 3 , A 5 , A 6 , A 7 } { A 1 , A 8 } { A 2 , A 4 }
Table 20. Ranking comparison of alternatives obtained by different methods.
Table 20. Ranking comparison of alternatives obtained by different methods.
Proposed Method
(Based on α )
Proposed Method
(Based on β )
Proposed Method
(Based on γ )
TOPSIS [53]VIKOR [54]ELECTRE [55]
A1666666
A2777787
A3333434
A4888878
A5444342
A6222223
A7111111
A8555555
Table 21. Comparison of classification results.
Table 21. Comparison of classification results.
POS(C)BND(C)NEG(C)
Jia and Liu’s method [20] A 3 , A 5 , A 6 , A 7 A 1 , A 8 A 2 , A 4
Gao et al.’s method [13] A 3 , A 5 , A 6 , A 7 A 1 , A 2 , A 4 , A 8
Wang et al.’s method [56] A 3 , A 5 , A 6 , A 7 A 1 , A 8 A 2 , A 4
Proposed method A 3 , A 5 , A 6 , A 7 A 1 , A 8 A 2 , A 4
Table 22. Comparative analysis among different methods.
Table 22. Comparative analysis among different methods.
Different MethodsAttribute WeightsOutcome or Loss FunctionsConditional
Probability
Dynamic Decision-Making
Jia and Liu’s method [20]SubjectiveObjectiveSubjective×
Gao et al.’s method [13]ObjectiveObjectiveObjective
Wang et al.’s method [56]SubjectiveSubjectiveObjective×
Proposed methodObjectiveObjectiveObjective
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Jia, F.; Wang, Y.; Su, Y. GRA-Based Dynamic Hybrid Multi-Attribute Three-Way Decision-Making for the Performance Evaluation of Elderly-Care Services. Mathematics 2023, 11, 3176. https://doi.org/10.3390/math11143176

AMA Style

Jia F, Wang Y, Su Y. GRA-Based Dynamic Hybrid Multi-Attribute Three-Way Decision-Making for the Performance Evaluation of Elderly-Care Services. Mathematics. 2023; 11(14):3176. https://doi.org/10.3390/math11143176

Chicago/Turabian Style

Jia, Fan, Yujie Wang, and Yiting Su. 2023. "GRA-Based Dynamic Hybrid Multi-Attribute Three-Way Decision-Making for the Performance Evaluation of Elderly-Care Services" Mathematics 11, no. 14: 3176. https://doi.org/10.3390/math11143176

APA Style

Jia, F., Wang, Y., & Su, Y. (2023). GRA-Based Dynamic Hybrid Multi-Attribute Three-Way Decision-Making for the Performance Evaluation of Elderly-Care Services. Mathematics, 11(14), 3176. https://doi.org/10.3390/math11143176

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop