An Entropy-Based Method for Probabilistic Linguistic Group Decision Making and its Application of Selecting Car Sharing Platforms

As the tourism and mobile internet develop, car sharing is becoming more and more popular. How to select an appropriate car sharing platform is an important issue to tourists. The car sharing platform selection can be regarded as a kind of multi-attribute group decision making (MAGDM) problems. The probabilistic linguistic term set (PLTS) is a powerful tool to express tourists’ evaluations in the car sharing platform selection. This paper develops a probabilistic linguistic group decision making method for selecting a suitable car sharing platform. First, two aggregation operators of PLTSs are proposed. Subsequently, a fuzzy entropy and a hesitancy entropy of a PLTS are developed to measure the fuzziness and hesitancy of a PLTS, respectively. Combining the fuzzy entropy and hesitancy entropy, a total entropy of a PLTS is generated. Furthermore, a cross entropy between PLTSs is proposed as well. Using the total entropy and cross entropy, DMs’ weights and attribute weights are determined, respectively. By defining preference functions with PLTSs, an improved PL-PROMETHEE approach is developed to rank alternatives. Thereby, a novel method is proposed for solving MAGDM with PLTSs. A car sharing platform selection is examined at length to show the application and superiority of the proposed method.


Introduction
With the development of mobile internet and "Internet +", sharing economy, such as traffic-sharing, health care-sharing and food-sharing, has boomed in recent years. Carsharing, a new mode of traffic-sharing, is becoming more and more popular along with the development of tourism. People usually complete the car rental and provide their feedbacks on consumer experiences on sharing platforms. Therefore, how to select a suitable car sharing platform is important for tourists. Since the car sharing platforms are often evaluated in safety, convenience, the brand of car and so on, the selection of car sharing platform can be described as a kind of multi-attribute decision making problems. Ordinarily, several tourists travel together. Hence, more than one tourist (decision maker, shorted by DM) decides which car sharing platform is selected. Thereby, the car sharing platform selection can be considered as a kind of multi-attribute group decision making (MAGDM) Wan et al., 2015;Dong et al., 2018;Kou et al., 2020;.
Owing to the complexity of problems and vagueness of human thinking, it is difficult for DMs to describe the decision information as crisp numbers. DMs are more apt to use linguistic variables (Zadeh, 1975) to evaluate alternatives. However, the linguistic variable allows DMs to express decision information with only one linguistic term (LT). Sometimes, DMs cannot adequately describe their evaluations using one exact LT and have hesitancy among several LTs. To overcome this drawback of the linguistic variable, Rodríguez et al. (2012) introduced the hesitant fuzzy linguistic term set (HFLTS) which permits DMs to express their preferences on alternatives with several possible LTs. Nevertheless, all possible LTs in a HFLTS have the same importance. In fact, DMs may prefer a LT to others. Thus, some useful information may be lost (Pang et al., 2016). To make up for this limitation, Pang et al. (2016) generalized HFLTS and presented the probabilistic linguistic term set (PLTS) in which each possible LT is assigned a probability (weight). Hence, the PLTS contains more useful information compared with the HFLTS. In recent years, the PLTS has gained more attention and some research achievements have been obtained. These achievements can be roughly divided into four categories.
(1) Operations and aggregation operators of PLTSs. Operations of PLTSs and aggregation operators are important to aggregate decision information in the PLTS circumstance. Pang et al. (2016) firstly defined some operational laws of PLTSs, and then put forward PLWA and PLWG operators of PLTSs, respectively. Subsequently, Zhang et al. (2017) pointed that the operation result by the laws (Pang et al., 2016) is a linguistic term rather than a PLTS. To remedy this defect, Zhang et al. (2017) defined new operations of PLTSs. Nevertheless, the operation results by these new operations in Zhang et al. (2017) may exceed the bounds of LTSs. To overcome this limitation, Gou and Xu (2016) introduced novel operations of PLTSs by a linguistic scale function. Later, using this function, Mao et al. (2019) introduced new operation laws and developed the GPLHWA operator and GPLHOWA operator based on Archimedean t-norms and s-norms. Zhang (2018) investigated the large group decision making and proposed a PL-WAA operator. Recently, diverse types of operators, such as power operators  and dependent weighted average operators , have been proposed one after another. Furthermore, Mi et al. (2020) conducted a survey of existing achievements on operations of PLTSs and aggregation operators.
(2) Distance measure of PLTSs. Pang et al. (2016) and Zhang et al. (2016) defined Euclidean distance and Hamming distance measures, respectively, of PLTSs based on the probabilities and subscripts of possible LTs. Afterwards,  argued that the results derived by these two distance measures (Pang et al., 2016;Zhang et al., 2016) occasionally are against human intuition. To overcome this defect, they proposed an improved distance measure. However, the calculation of this improved distance measure is very complex. Lately, Mao et al. (2019) proposed a simple Euclidean distance measure by a linguistic scale function (Gou and Xu, 2016). Unfortunately, some counterintuitive results still appear based on this Euclidean distance, please see Table 2 in Section 3.2.
(3) Entropy and cross entropy of PLTSs. Due to the uncertainty of information with PLTSs, how to measure such a certainty is important. Entropy of PLTSs, as an efficient tool to measure the uncertainty of information, is an interesting topic but has not gained wide attention. Only Lin et al. (2019) developed an information entropy of PLTSs by probabilities of possible LTs in a PLTS. Liu et al. (2018) proposed three kinds of entropy of PLTSs by extending the entropy of HFSs into the PLTS context, including the fuzzy entropy, hesitancy entropy and total entropy. On the other hand, the cross entropy is effective for measuring the differences between PLTSs. Up to now, only  presented two cross entropy measures of PLTSs based on sine and tangent trigonometric functions.
(4) Decision methods for solving MAGDM problems with PLTSs. It is important to choose proper decision methods for selecting best alternatives. At present, pools of methods have been proposed to solve MAGDM problems with PLTSs. For example, Pang et al. (2016) presented a distance-based extended TOPSIS method with PLTSs. Zhang X.F. et al. (2019) argued that this extended TOPSIS only considered the distance proximity of alternatives with respect to ideal and negative ideal solutions, while ignored the proximity in direction. In response to this problem, Zhang X.F. et al. (2019) developed a projection method to solve MAGDM with PLTSs. Later, by improving the classical MULTI-MOORA method in different angles,  and  developed a PL-MULTIMOORA method and an extended MULTIMOORA method, respectively. Unfortunately, these two methods did not take negative ideal solutions into account. Considering both positive and negative ideal solutions, Li et al. (2020) put forward a PP-MULTIMOORA method. In addition, diverse types of outranking methods (Lin et al., 2019;Liao et al., 2019;Xu et al., 2019;Peng et al., 2020;Liu and Li, 2018), such as ELEC-TRE and PROMETHEE, are also presented. To facilitate the study on decision methods of PLTSs, Liao et al. (2020) provided a survey which includes most decision methods and their applications.
Although many achievements have been achieved, they suffer from some limitations.
(1) Though the GPLHWA operator (Mao et al., 2019) improved the PLWA operator (Pang et al., 2016) and the PL-WAA operator  to some extent, it has some flaws that the aggregated result obtained by the GPLHWA operator is not a PLTS in a strict sense. The reason is that the sum of probabilities of all possible LTs in the aggregated result is more than 1. Therefore, proposing new aggregation operators with desirable properties is helpful for aggregating evaluation values of alternatives.
(2) By existing distance measures of PLTSs (Pang et al., 2016;Mao et al., 2019;Zhang et al., 2016), some counter-intuitive results appears. Although the results derived by the distance measure  almost agree with human intuition, the computation of this distance measure is too complex to use easily. Therefore, it is necessary to develop a new simple distance measure of PLTSs by which computation results coincide with human intuition.
(3) The research on entropy and cross entropy of PLTSs is very little. Up to date, only Lin et al. (2019) and Liu et al. (2018) addressed the entropy of PLTSs. However, the distinguishing power of the entropy proposed by Lin et al. (2019) is not high enough. Although the fuzzy entropy and hesitancy entropy proposed by Liu et al. (2018) can neatly measure the uncertainty of PLTSs, some properties of this fuzzy entropy are counterintuitive and the computation of this hesitancy entropy is not simple. In addition, the cross entropy proposed by  fails for the symmetric linguistic term sets in PLTSs. Thereby, it is valuable to study deeply on the entropy and cross entropy of PLTSs.
(4) Existing decision methods are fruitful to solve decision problems with PLTSs. Nevertheless, some methods (Lin et al., 2019;Xu et al., 2019;Li et al., 2020) can only solve MADM problems, but fail for MAGDM problems. Although these methods (Pang et al., 2016;Zhang X.F. et al., 2019; are capable to solve MAGDM problems with PLTSs, DMs' weights or attribute weights are not considered or assigned in advance, which may result in subjective random weights. Hence, it is interesting to seek a novel method which not only can solve MAGDM problems, but also determine DMs' weights and attribute weights objectively. To make up for above limitations, this paper proposes a novel method for solving MAGDM problems with PLTSs. First, two aggregation operators, including PLWAM (probabilistic linguistic weighted arithmetic mean) operator and PLWGM (probabilistic linguistic weighted geometric mean) operator, respectively, are proposed and some desirable properties are studied. To measure the hesitancy degree of a PLTS, a hesitancy index of the PLTS is introduced. Then a general distance measure of PLTSs is defined to measure the deviation between two PLTSs. Considering the fact that the uncertainty of a PLTS includes the fuzziness and hesitancy, a fuzzy entropy and hesitancy entropy of a PLTS are defined and then a total entropy of a PLTS is derived to measure such uncertainties. Meanwhile, a cross entropy of PLTSs is defined to measure the distinction between PLTSs. Afterwards, by minimizing the total entropy and the cross entropy of a DM, an objective program model is built to determine DMs' weights. Subsequently, individual decision matrices are aggregated into a collective one by the PLWAM operator. To derive attributes weights objectively, a bi-objective program is built by minimizing the total entropy of attribute values with respect to each attribute as well as maximizing the cross entropy between attribute values of alternatives. By defining a new preference function in the form of PLTSs, an improved PL-PROMETHEE method is developed to rank alternatives. Thereby, a novel method is proposed for solving MAGDM problems with PLTSs. A case of car sharing platform selection is applied to show the effectiveness and advantages of the proposed method at length. The primary features of the proposed method are outlined as follows: (1) Two new probabilistic linguistic average aggregation operators of PLTSs (i.e. PLWAM and PLWAG operators) are proposed. A prominent characteristic of them is that the aggregated result obtained by these two operators is not only a PLTS with the sum of probabilities of possible LTs being equal to 1, but also is consistent with human intuition.
(2) A new generalized distance measure of PLTSs is defined. It is worth mentioning that the hesitancy degree of a PLTS is considered in this distance measure. Thus, the new distance has a stronger distinguishing power. Moreover, a ranking approach is presented to rank PLTSs.
(3) A fuzzy entropy, a hesitancy entropy and a cross entropy of PLTSs are introduced. The fuzzy entropy has desirable properties and the computation of the hesitancy entropy is simple. Meanwhile, the cross entropy can distinguish the deviations between PLTSs with symmetric linguistic term sets.
(4) Based on entropy and cross entropy of PLTSs, distinct objective programs are established to determine DMs' weights and attribute weights objectively. Finally, an improved PL-PROMETHEE method is developed to rank alternatives.
The remainder of this paper is organized as follows: In Section 2, some basic concepts of PLTSs are reviewed. Moreover, PLWAM and PLWGM operators are proposed. Section 3 introduces a hesitancy index of a PLTS and then develops a generalized distance measure of PLTSs. Based on this distance measure, a ranking approach is presented to sort PLTSs. Section 4 defines several types of entropy of PLTSs, including the fuzzy entropy, hesitancy entropy, total entropy and cross entropy of PLTSs. In Section 5, an improved PL-PROMETHEE method is developed to solve MAGDM problems with PLTSs. Section 6 provides a case study of car sharing platform selection to illustrate the application of the proposed method. Furthermore, comparison analyses are conducted to show advantages of the proposed method. Some conclusions are made in Section 7.

Preliminaries
In this section, some definitions and notions related to the PLTS are reviewed. Furthermore, two aggregation operators of PLTSs are proposed and some desirable properties of them are investigated.
Definition 1 (See Xu, 2005). Let S = {s α |α = −τ, . . . , −1, 0, 1, . . . , τ } be a finite and totally ordered discrete LTS, where s α represents a possible value for a linguistic term, and τ is a positive integer. Especially, the mid-linguistic label s 0 represents an assessment of "indifference", and the rest of them are placed symmetrically around it. s −τ and s τ are lower and upper bounds of linguistic labels.
To preserve all given linguistic information, Xu (2004) extended the discrete LTS S into a continuous LTSS = {s α |α ∈ [−τ, τ ]. If s α ∈ S, then s α is called an original linguistic term; Otherwise, s α is called a virtual linguistic term.
Definition 2 (See Gou and Xu, 2016). Let S = {s α |α = −τ, . . . , −1, 0, 1, . . . , τ } be a LTS, the linguistic term s α that expresses the equivalent information to the membership degree γ is obtained by a linguistic scale function g: Additionally, the membership degree γ that expresses the equivalent information to the linguistic term s α is obtained by the following function Definition 3 (See Pang et al., 2016). Let S = {s α |α = −τ, . . . , −1, 0, 1, . . . , τ } be a linguistic term set, a PLTS is defined as where L (k) (p (k) ) represents the linguistic term L (k) associated with the probability p (k) , and #L(p) is the number of all different linguistic terms in L(p). Definition 4 (See Pang et al., 2016). Given a PLTS L(p) with #L(p) k=1 p (k) 1, the normalized PLTSL(p) is defined as: To ensure that the operational results among PLTSs can be straightforwardly determined, Mao et al. (2019) defined an ascending ordered PLTS below.
In real decision making problems, the numbers of linguistic terms in two different PLTSs are often different. This makes trouble to operate. In order to make them have the same number of linguistic terms, Pang et al. (2016) provided a method to add the number of linguistic terms for PLTSs in which the number of linguistic terms is relatively small as follows.

Aggregation Operators of PLTSs
This subsection reviews the operational laws of PLTSs presented by Gou and Xu (2016), and then proposes PLWAM and PLWGM aggregation operators, which are used to aggregate the decision information in GDM discussed in later sections.
Remark 1. As reviewed in Mi et al. (2020), up to now existing operations of PLTSs can be divided into five types: symbolic-based computation (Pang et al., 2016), semantic-based computation (Gou and Xu, 2016;Wu and Liao, 2019), triangular norms-based computation (Klement et al., 2000), evidence reasoning-based computation (Li and Wei, 2019) and double alpha-cut-based computation (Jiang and Liao, 2020). Among these operations, only the semantic-based computation (Gou and Xu, 2016;Wu and Liao, 2019) and the evidence reasoning-based computation (Li and Wei, 2019) satisfy the closure of operation of PLTSs. The semantic-based computations (Gou and Xu, 2016;Wu and Liao, 2019) defined common operations, such as additive, subtraction and numerical operations of PLTSs. However, the computation (Li and Wei, 2019) only defined the additive operation of PLTSs. To sum up, the semantic-based computation (Gou and Xu, 2016;Wu and Liao, 2019) is more reasonable compared with other computations. Therefore, it is natural to choose this kind of computation to complete the operations of PLTSs. The differences between operations in Gou and Xu (2016) and Wu and Liao (2019) are the selections of linguistic scale functions and the tools for handling two PLTSs with different numbers of LTs. Compared with the computation in Wu and Liao (2019) (see Definition 3 in Wu and Liao, 2019), the computation in Gou and Xu (2016) is easier to operate and more widely used. For example, many achievements (Bai et al., 2017;Liang et al., 2018;Feng et al., 2020) have been gained based on the computation in Gou and Xu (2016), while only one method  was developed by the computation in Wu and Liao (2019). At this point, we choose the computation in Gou and Xu (2016) to conduct the operations of PLTSs in this paper.
In virtue of Definition 9, a PLWAM operator of PLTS is presented.
Theorem 1. Given n PLTSs L 1 (p), L 2 (p), . . . , L n (p), the aggregated value by the PLWAM operator is also a PLTS as Furthermore, Analogously, a PLWGM operator of PLTSs is presented. where In a way similar to the proving process of Theorem 1, Theorem 2 can be proved based on the operational laws (ii) and (iv) in Definition 9.
Theorem 2. Given n PLTSs L 1 (p), L 2 (p), . . . , L n (p), the aggregated value by the PLWGM operator is also a PLTS as follows and Remark 2. Compared with existing aggregation operators (Pang et al., 2016;Mao et al., 2019;Zhang, 2018), the proposed operators have significant advantages. For instance, they satisfy the closure of operations and the computation results are closer to human intuition.  (1) The PLWAM operator satisfies the closure of operations, i.e. the aggregated result of PLTSs obtained by the PLWAM operator is still a PLTS. However, the aggregated result obtained by the PLWA operator (Pang et al., 2016) is a linguistic term rather than a PLTS. Thus, the probability information is lost. Although the aggregated result obtained by the GPLHWA operator (Mao et al., 2019) seems to be a PLTS, the sum of probabilities of all possible LTs in this result is more than 1. Therefore, this result is not a PLTS in a strict sense. Hence, there exists information distortion while using the GPLHWA operator.
(2) The aggregated result obtained by the proposed PLWAM operator is more consistent with human intuition. Observing L 1 (p), L 2 (p) and L 3 (p), it is deduced that the aggregated result should be larger than s 1 and smaller than s 3 . Furthermore, the probability of s 3 should be more than that of s 1 because three PLTSs include s 3 but only L 2 (p) contains s 1 with a small probability. By the PLWAM operator proposed in this paper, the aggregated result (see Table 1) is in accordance with this analysis indeed. On the other hand, the aggregated result obtained by the PL-WAA operator (Zhang, 2018) is counterintuitive because LTs s −3 , s −2 and s −1 are all included in this result, which means that the information is distorted in aggregating process.
(3) Observing Table 1, the aggregated results obtained by the proposed PLWAM operator and the PL-WAA (Zhang, 2018) operator are both PLTSs. Nevertheless, this paper proves that all the aggregated results by the PLWAM operator are all PLTSs for any PLTSs (see Theorem 1), but Zhang (2018) failed to do so.
Analogously, the PLWGM operator has above advantages.

A Hesitancy Index of a PLTS and an Approach to Ranking PLTSs
In this section, a hesitancy index of a PLTS is introduced to measure the hesitancy degree of a PLTS, and then a new generalized distance measure between two PLTSs is developed.
Finally, based on the proposed distance measure, a TOPSIS-based approach is presented to rank PLTSs.

A Hesitancy Index of a PLTS
This result is in accordance with human intuition.

A General Distance Measure Between PLTSs
First, this subsection reviews existing distance measures between PLTSs (Pang et al., 2016;Mao et al., 2019;Zhang et al., 2016), and then develops a general distance measure between PLTSs based on the proposed hesitancy index. In the end, comparative analyses of the proposed distance measure and existing ones are conducted.
where p(r Mao et al. (2019) argued some limitations of above distance measures and developed an Euclidean distance as This paper defines a general distance measure between PLTSs by the proposed hesitancy index.
Definition 13. Given two PLTSsL where q 1. Specially, when q = 1, a Manhattan distance is obtained as when q = 2, a novel Euclidean distance of PLTSs is developed as According to Eqs. (14) and (15), it is easily proved that d M and d E have desirable properties described in Theorem 3.
Theorem 3. For two PLTSsL 1 (p) andL 2 (p) with #L 1 (p) = #L 2 (p), d M and d E are the novel Manhattan and Euclidean distances of PLTSs, respectively. They satisfy: Zhang's distance (Zhang et al., 2016) The Manhattan distance proposed in this paper The Euclidean distance proposed in this paper Pang's distance (Pang et al., 2016) The Euclidean distance proposed in this paper The Mao's Euclidean distance (Mao et al., 2019) d(L 5 (p),L 6 (p)) = 0 The Manhattan distance proposed in this paper The Euclidean distance proposed in this paper Table 2 Remark 3. Compared with existing distance measures in Pang et al. (2016), Mao et al. (2019) and Zhang et al. (2016), the distance measures defined in Definition 13 have a stronger distinguishing power. Table 2 shows this advantage explicitly.
It is shown from Table 2 that PLTSs in each pair are not identical completely. Therefore, the distance between each pair PLTSs should be more than zero intuitively. The proposed distance measures do this exactly. However, the distances by existing distance measures (Pang et al., 2016;Mao et al., 2019;Zhang et al., 2016) are all zeros, so the difference of PLTSs in each pair cannot be distinguished. This means that the proposed distance measures have stronger distinguishing powers, which can be attributed to the following factors: (1) According to Zhang's distance (Zhang et al., 2016) (i.e. Eq. (11)), the distance between two PLTSs is zero if these two PLTSs have the same possible LTs regardless of their probabilities. Hence, such PLTSs, likeL 1 (p) andL 2 (p) in Table 2, cannot be distinguished even if they are not remarkably identical.
(4) The distance measures in this paper are developed by considering each possible element of a PLTSL(p) as a two-dimensional vector (i.e.L (k) (p (k) ) = (g(L (k) ),p (k) )) and then extending the corresponding classical distance measures into the PLTS context. Furthermore, it is worth noticing that the proposed distances take the hesitancy index of a PLTS into consideration. Thus, the distances proposed in this paper overcome above limitations and strengthen the distinguishing power.

A Ranking Approach for PLTSs
Based on the developed distance measures, this subsection introduces a TOPSIS based ranking approach to comparing PLTSs.
Let L(p) + = {s τ (1)} be a positive ideal PLTS and L(p) − = {s −τ (1)} be a negative ideal PLTS. According to the TOPSIS method, the closer a PLTS L(p) = {L (k) (p (k) )|k = 1, 2, . . . , #L(p)} is to L(p) + and at the same time the farther L(p) to L(p) − , the better the PLTS L(p). Thus, a closeness degree of the PLTS L(p) is defined as whereL(p),L(p) + andL(p) − are the normalized ordered PLTSs corresponding to L(p), In virtue of the closeness degree, a ranking approach for PLTSs is introduced.

Entropy and Cross Entropy of PLTSs
In order to judge the quality of decision information provided by DMs and the discriminations between DMs in the sequel, this section develops some entropy measures and a cross entropy measure.

Entropy measures of PLTSs
This subsection addresses the entropy of PLTSs and proposes some new entropy measures of PLTSs. The main motivations are outlined as follows: (1) Measure the uncertainty of a PLTS neatly. Although Lin et al. (2019) defined an entropy measure of a PLTS, the distinguishing power of this measure is not high and cannot neatly measure the uncertainty of a PLTS. Therefore, it is necessary to develop some new entropy measures of a PLTS.
(2) Overcome the shortcomings of entropy measures proposed by Liu et al. (2018). As Liu et al. (2018) stated that, the uncertainty of a PLTS includes the fuzziness and the hesitancy.
To measure such fuzziness and hesitancy, Liu et al. (2018) defined a fuzzy entropy measure and a hesitancy entropy measure of a PLTS respectively. However, these entropy measures have two shortcomings: one is that some properties of the fuzzy entropy measure are not consistent with intuition. The other is that the computation of the hesitancy entropy is time-consuming. To overcome these shortcomings, this subsection proposes a new fuzzy entropy measure with desirable properties and a hesitancy entropy measure with simple computation, respectively. Finally, the total entropy of a PLTS is defined by combining the proposed fuzzy entropy and the hesitancy entropy.
According to Definition 15, some properties of the fuzzy entropyĒ F are not in accordance with intuitions (please see Remark 4). Therefore, we define a new fuzzy entropy of a PLTS with desirable properties as below.
Proof. Please see Appendix B.
Remark 4. Compared with the fuzzy entropyĒ F , the proposed fuzzy entropy E F has the following merits: (1) The properties of the latter are more consistent with intuitions. Let's compare the property (i) ofĒ F with that of E F . From (i) in Definition 15, we haveĒ F (s 0 (p), s g (1 − p)) = 0. Specially, p = 0.5, the equalityĒ F (s 0 (0.5), s g (0.5)) = 0 means that the fuzzy uncertainty of the PLTS L * (P ) = {s 0 (0.5), s g (0.5)} is zero. However, in real decision making, the decision information represented by L * (P ) is very fuzzy because L * (P ) indicates that the probability of DM dissatisfying an alternative totally and that of DM satisfying an alternative completely are both 0.5. Therefore, the prop-ertyĒ F (s 0 (0.5), s g (0.5)) = 0 is counter-intuitive. However, according to the proposed fuzzy entropy E F defined in Eq. (17), one getsĒ F (s 0 (p), s g (1 − p)) = f (p). This indicates that the fuzzy entropy of L * (P ) varies with p. Employing Theorem 4, we have E F (s 0 (0.5), s g (0.5)) = f (0.5) = 1. At this point, the proposed entropy E F is more consistent with intuition.
(2) The latter can be regarded as an extension of the former in some sense. It is worth mentioning that a normalized PLTS can be considered as a type of discrete random variables. Therefore, the expectation value E(L(P )) represents the average level ofL(P ). Based on this idea, it is reasonable that the decision information described byL(P ) and that of {s g/2 (1)} in Definition 15 are equivalent if E(L(P )) = 0.5. Hence, the prop-ertyĒ F (s g/2 (1)) = 1 in Definition 15 can be extended to the one E F (L(p)) = 1 ⇔ e(L(p)) = 0.5 in Definition 16. Likewise, the property (iii) in Definition 15 is also extended into the one (iii) in Definition 16.
According to Theorem 4, let f (t) = ln(1+t−t 2 ) ln 5−ln 4 . A fuzzy entropy of a PLTS is defined as We call E H a hesitancy entropy of a PLTS if it satisfies the following properties:

Theorem 5. Given a normalized ordered PLTSL(p). The hesitancy index
is a hesitancy entropy, whereḡ = 1 Proof. Please see Appendix C.
where γ ij = |I (l i )−I (l j )| g and I (l i ) is the subscript of the linguistic term l i .
ComparedĒ h with E h , the difference between them is that the former measures the deviations between different linguistic terms in a PLTS, while the latter measures the deviations between all the transformation function values of possible linguistic terms and their mean value. Although the expressions ofĒ h and E h are different, they have similar properties (please see Definition 5 in Liu et al. (2018). However, the computation of E h needs #L(P ) + 1 times, whileĒ h needs #L(P )(#L(P ) − 1)/2 times. Thus, the computations of E h is more time-saving when #L(P ) 4.
Combining the fuzzy entropy and hesitancy entropy with an adjusted coefficient θ , a total entropy is determined as As 0 E F (L(p)), E H (L(p)) 1, one has 0 E T (L(p)) 1. Furthermore, the total entropy of a probabilistic linguistic matrix is derived.
Definition 18. LetŪ = (L ij (p)) m×n be a matrix with PLTSs. The total entropy ofŪ is defined as Remark 6. Lin et al. (2019) proposed an information entropy as where L(p) is a PLTS, #L(p) is the number of all different linguistic terms in L(p) and p (k) is the probability of kth possible linguistic value of L(p).
It is obvious from Eq. (23) that this information entropy only considered the probabilities of possible linguistic terms in a PLTS, but ignored all possible linguistic terms. However, the total entropy presented in this paper takes them into consideration together. Therefore, compared with μ(L(p)) in Eq. (23), the presented total entropy E T (L(p)) has a stronger distinguishing power while measuring the uncertainty of a PLTS, which can be verified by Example 2. Obviously, the uncertainty of L 1 (p) is larger than that of L 2 (p). In fact, the total entropy is calculated as E T (L 1 (p)) = 0.552 > E T (L 2 (p)) = 0.398, which is in accordance with our intuition. However, by Eq. (23), it obtains that μ(L 1 (p)) = μ(L 2 (p)) = 1.5395, by which the uncertainties of L 1 (p) and L 2 (p) are considered the same. In other words, μ(L(p)) in Eq. (23) is unable to distinguish the degrees of uncertainty of L 1 (p) and L 2 (p). Therefore, the total entropy proposed in this paper has a stronger distinguishing power.

The Cross Entropy of PLTSs
In this subsection, to measure the discrimination degree between PLTSs, a cross entropy is defined. Motivated by the cross entropy in the hesitant fuzzy linguistic environment (Gou et al., 2017), the axiomatic definition of cross-entropy measure for PLTSs is given below.
Proof. Please see Appendix D.
It is observed from Eq. (24) that the cross entropy, CE(L 1 (p),L 2 (p)), is not symmetric. A symmetric cross entropy between PLTSs can be obtained as Moreover, a symmetric cross entropy between two probabilistic linguistic matrices is defined below.

A Novel Method for MAGDM with PLTSs
In this section, we first describe the problems of MAGDM with PLTSs, and then propose a novel method for solving such problems.

Determine DMs' Weights Based on the Total Entropy and the Symmetric Cross Entropy
In this subsection, an approach is developed to determine DMs' weights objectively by using the proposed total entropy and symmetric cross entropy of PLTSs. As we know, the less uncertainty (i.e. a smaller total entropy) of an individual matrix provided by a DM, the better the quality of decision information reflected by this matrix is. Thus, the bigger the weight of this DM should be assigned. In virtue of this criterion, a programming model for deriving DMs' weights is built by minimizing the total entropy of decision matrices, i.e.
where E T (Ū k ) are the total entropy ofŪ k (k = 1, 2, . . . , t). Solving Eq. (27) with Lagrange multiplier approach, it is obtained as Normalizing λ * k 1 , one has On the other hand, a closer degree of a DM to other DMs means that the information supplied by this DM is much closer to that of the group. In this case, this DM should be assigned a larger weight. From this viewpoint, another optimal model for determining DMs' weights is constructed based on the symmetric cross entropy as Eq. (30) can also be solved with Lagrange multiplier approach, and the optimal solutions are derived as 1, 2, . . . , t).

Construct Bi-Objective Programs For Deriving Attribute Weights
In group decision making process, attribute weights play an important role because different attribute weights may result in diverse ranking of alternatives. In this subsection, considering the information of attribute weights is completely unknown or partially known, two bi-objective programs are constructed respectively for deriving attribute weights.
(i) Aggregating individual decision matrices into a collective one By employing DMs' weights determined in Section 5.2, individual decision matrices are integrated into a collective one with the proposed PLWAM operator. Denote the collective decision matrix by U = (L ij (p)) m×n , where Then, the normalized ordered collective decision matrixŪ = (L ij (p)) m×n is obtained from U = (L ij (p)) m×n by Definition 4, whereL ij (p) = {L (ii) Constructing bi-objective programs for deriving attribute weights It is known to us that an attribute plays a more important role if the performance values of alternatives on it have distinct differences. Thus, such an attribute should be given a bigger weight. Conversely, if the evaluation values of alternatives with respect to a certain attribute have little difference, this attribute should be given a smaller weight. In addition, the credibility of decision information on an attribute should also be taken into account while determining attribute weights. The more credibility of evaluation values on an attribute, the bigger weight of this attribute should be assessed. As mentioned before, the symmetric cross entropy and the total entropy can reflect the differences between alternatives and the credibility of evaluation values, respectively. Keeping this idea in mind, we can determine attribute weights by maximizing the symmetric cross entropy as well as minimizing the total entropy. Thus, a bi-objective program is established if the information of attribute weights is completely unknown, i.e.
As 0 E T (L ij (p)) 1, Eq. (35) can be converted into the following single objective program: Solving Eq. (36) with Lagrange multiplier approach, the weights of attributes are derived as The normalized weights w j are obtained by normalizing w * j , i.e.
Similarly, for the situations where the information of attribute weights is incomplete, another programming model is obtained by modifying Eq. (36), i.e.

Ranking Alternatives by an Improved PL-PROMETHEE Method
The PROMETHEE method is a popular outranking method in decision making. It has been extended to different fuzzy environments, such as intuitionistic fuzzy sets (Krishankumar et al., 2017) and linguistic variables Halouani et al. (2009). In the PLTS environment, although Xu et al. (2019) proposed a PL-PROMETHEE method, the distinguishing power of PL-PROMETHEE is not high, which can be verified by Example 3 in the sequel. To improve the distinguishing power, this subsection constructs new probabilistic linguistic preference relations, the integrated preference index and positive/negative outranking flows to adapt the main framework of PROMETHEE in the PLTS context. Thus, an improved PL-PROMETHEE method is proposed to solve MAGDM problems with PLTSs.
(i) The probabilistic linguistic preference function Suppose that A i and A r are two alternatives in the alternative set.L ij (p) andL rj (p) are respectively normalized ordered attribute values of A i and A r with respect to attribute u j . Preference function P j X (A i , A r ) ofL ij (p) with respect toL rj (p) is defined as where In Eq. (41), parameters p( 0) and q( 0) are the indifference threshold and the strict preference threshold, respectively. Functions g(x) and g −1 (x) are defined in Definition 2.
Obviously, the preference function P j X (A i , A r ) constructed by Eqs. (40) and (41) is a PLTS.
(ii) The integrated preference index According to the probabilistic linguistic preference index P j X (A i , A r ), an integrated preference index is defined as where T (P j (A i , A r )) is the closeness degrees of P j X (A i , A r ) and computed by Eq. (16). (iii) The positive and negative outranking flows Employing the integrated preference indices, the positive/negative outranking flows, φ + (A i ) and φ − (A i ), are defined as follows.
In virtue of Eqs. (43) and (44), the net flow of alternative A i is calculated as Finally, alternatives are ranked by the descending order of net flows.
Remark 7. In PL-PROMETHEE method , the preference function and the total preference index are defined respectively as and where P (A i,j > A r,j ) is the probability of A i preferred to A r with respect to u j . Obviously, the value of the preference function P j (A i , A r ) in Eq. (46) only takes 0 or 1, which may result in the loss of decision information. However, the preference function P j X (A i , A r ) in Eq. (40) is described by the PLTS which is more flexible than the crisp number for representing decision information. Therefore, compared with PL-PROMETHEE, the distinguishing power of the improved PL-PROMETHEE is stronger and the sensitivity to the ratings of alternatives on attributes is increased. To illustrate this advantage, Example 3 is given below.

A Novel Method for MAGDM with Probabilistic Linguistic Information
A novel method is generated for MAGDM problems with PLTSs. The main procedure of this method is outlined below.

A Case Study
In this section, an example of a car sharing platform selection is provided to illustrate the application of the proposed method. Furthermore, the comparative analyses are performed to show the merits of the proposed method.  Table 4 Probabilistic linguistic decision matrices U 1 , U 2 and U 3 .

Car Sharing Platform Selection
With the rapid development of internet technology and the deep advocation of green travel, the car sharing has sprung up over the last three years. Up to now, several car sharing platforms have emerged in China, such as Evcard, Gofun, Togo and so on. The popularization of car sharing greatly facilitates peoples' travel and relieves the traffic pressure. As a famous tourist city in China, Guilin cannot satisfy the travel of tourists due to the limited operational capacity of public transportation. So it is necessary for the government to introduce a car sharing platform to resolve the traffic problem. Now, the government invites three DMs to select the best car sharing platform from four candidate platforms (alternatives), including Evcard, Gofun, Togo and Greengo. Four attributes are considered, including safety (u 1 ), convenience (u 2 ), service (u 3 ) and car brand (u 4 ). Suppose a linguistic term set S = {s −4 , s −3 , s −2 , s −1 , s 0 , s 1 , s 2 , s 3 , s 4 } is provided in Table 3. Three invited DMs, d 1 , d 2 and d 3 , utilize PLTSs to evaluate four alternatives with respect to four attributes, and construct three probabilistic linguistic decision matrices U 1 , U 2 and U 3 .
Step 1. Each DM establishes decision matrices U 1 , U 2 and U 3 as shown in Table 4, and the corresponding normalized ordered matricesŪ 1 ,Ū 2 andŪ 3 are listed in Table 5.
Step 3. Determine the weight vector of DMs.
Step 4. Aggregate individual decision matrices into a collective one.
Steps 5. Calculate the total entropy and symmetric cross entropy of alternatives.
Take θ = 0.5, and the total entropy matrix E = E T (L ij (p)) is obtained by Eqs. (18)  Step 6. Determine the weight vector of attributes. The weight vector of attributes is determined by Eq. (38) as w = (0.2932, 0.2632, 0.2070, 0.2367) T .
Step 7. Compute the preference function values. Using Eqs. (40) and (41), preference function values are computed. In view of limited space, the concrete preference function values are not listed.
Step 8. Obtain the integrated preference index.
By Eq. (42), the total preference index matrix is derived as Step 9. Calculate positive and negative flows of alternatives. In virtue of Eqs. (43) and (44), positive and negative flows of alternatives are respectively calculated as Step 10. Rank alternatives.
Net flows of alternatives are obtained by Eq. (45) as , alternatives are ranked as A 2 A 1 A 4 A 3 . Therefore, the alternative A 2 (Gofun platform) is the best one.

Comparative Analyses
To show the advantages of the proposed method, comparative analyses with Mao's method (Mao et al., 2019), PL-PROMETHEE method  and other existing decision methods are performed in the sequel. Mao et al. (2019) presented a new method integrating ELECTRE and TOPSIS to solve MAGDM problems with PLTSs. Now we use Mao's method (Mao et al., 2019) to solve the above problem. Set γ = 1 and λ = 1. DMs' subjective weight vector is assigned as λ 1 = (0.04, 0.88, 0.08) T . Thus, the corresponding ordering results are derived as A 1 A 2 A 4 A 3 . Clearly, these results are not consistent with those obtained by the proposed method, which can be explained as follows:

Comparison with Mao's method
(1) The proposed method determines DMs' weights objectively by the total entropy and cross entropy of decision information. However, method (Mao et al., 2019) assigned DMs' weights subjectively, by which the random cannot be avoided. Although an adjusted coefficient obtained from the consistency degree between DMs is used to adjust the subjective weights of DMs, the adjusted DMs' weights may be irrational because the consistency degree ρ ir indicating the consistency between DM d i and d r (i.e. Eq. (40) in method Mao et al., 2019) does not satisfy the symmetry. For example, the consistency degree ρ 12 = 0.5125, while ρ 21 = 0.4799, So ρ 21 = ρ 12 . In fact, the consistency degree should be symmetrical theoretically. Hence, DMs' weights derived by the proposed method are more reasonable.

Compared with the PL-PROMETHEE Method
The PL-PROMETHEE method  is applied to solve the problem in Section 6.1. Suppose attribute weight vector as w = (0.3, 0.2, 0.3, 0.2) T . The Borda's scores of alternatives for each DM are shown as Table 8. It finds from Table 8 that the ranking result is A 1 ∼ A 3 A 2 ∼ A 4 , which is different from the one derived by the proposed method, i.e, A 2 A 1 A 4 A 3 . The primary reasons may result from the following two aspects.
(1) PL-PROMETHEE method ignored the determination of attribute weights and DMs' weights, but assigned them subjectively. Thus, the arbitrariness cannot be avoided. By contrast, the proposed method derives attribute weights and DMs' weights objectively based on the total entropy and symmetric cross entropy of decision information. Therefore, the subjectivity is effectively reduced and the decision results are more credible.
(2) In the preference function of PL-PROMETHEE method , the preference value is taken as 1 if the probability of one alternative preferred to the other with respect to an attribute is bigger than 0.5. Otherwise, the preference value is taken as 0. Therefore, the distinguishing power of this preference function is not STRONG enough, which may be the reason why alternatives A 1 and A 3 cannot be discriminated. However, in the proposed method, the preference function is neatly defined based on the deviations of attribute values and described by PLTSs, please see Eqs. (40) and (41). Hence, the proposed method has a stronger distinguishing power.

Compared with Other Existing Decision Making Methods
To further demonstrate the superiority of the proposed method, this subsection conducts a theoretical analysis and a practical analysis with other existing methods (Pang et al., 2016;Gou et al., 2017;Liu and Li, 2018;Mao et al., 2019;Peng et al., 2020).
(1) Theoretical analysis (i) The proposed method introduces two new probabilistic linguistic weighted average operators (i.e. PLWAM operator and PLWAG operator). Compared with existing operators mentioned in Pang et al. (2016) and Gou et al. (2017), the proposed operators have some advantages, such as the closure of operations and operation values being consistent with intuition (please see Remark 2). Therefore, the aggregated information obtained by the proposed operators is more reliable, and thus the decision results based on such aggregated information are more reasonable.
(ii) In the proposed method, two cases including the attribute weights with unknown or partially known values are both taken into account. To determine attribute weights, the proposed method builds two different bi-objective programming models by maximizing the cross entropy and minimizing the total entropy of the collective evaluation values. As analysed in Section 5.3, this model considers the quality of the collective evaluations and the deviations between evaluations. However, in the extended MULTIMOORA method , the attribute weights were given in advance, which may be not able to avoid the subjective arbitrariness in decision process. Although the TOPSIS method (Pang et al., 2016) derived attribute weights objectively by the maximum deviation approach, it ignored the quality of decision information. Thus, the creditability of attribute weights derived by the proposed method is higher.
(iii) Although method  and the proposed method can solve decision making problems with PLTSs, the former failed to handle group decision making problems while the latter can manage both single decision making problems and the group ones. Thereby, the latter has wider application fields.
(2) Practical analysis The distinguishing power of the proposed method is stronger than that of method , which can be verified by Section 6.2. This is largely due to the fact that the preference functions of the improved PL-PROMETHEE approach proposed in this paper are in the form of PLTSs and can distinguish alternatives neatly, whereas the preference functions defined in method  and the other PL-PROMETHEE approach (Liu and Li, 2018) are all crisp numbers.
(2) The stableness of the proposed method is better than method . The former determines DMs' weights and attribute weights by building objective programming model. Thereby, the decision result is unique under the given decision information. In contrast, the latter assigned DMs' weights or attribute weights in advance. Hence, the change of such weights may result in different decision results.
(3) Method Gou et al. (2017) is suitable for the environment of hesitant fuzzy linguistic variables. Nevertheless, it cannot handle MAGDM problems with PLTSs which can be solved by the proposed method.

Conclusions
In today's internet age, car sharing is more and more popular. The car sharing platform selection is important for tourists, which can be regarded as a MAGDM problem. The PLTS is a powerful tool to represent the evaluation information of DMs in complex MAGDM problems. This paper introduces PLWAM and PLWGM operators firstly and studies some desirable properties of them. Subsequently, a hesitancy index of a PLTS and a general distance measure between PLTSs are defined. Then, a new approach is proposed to rank PLTSs. To measure the fuzziness and hesitancy of a PLTS, a fuzzy entropy and a hesitancy entropy of PLTSs are presented. Afterwards, a total entropy of PLTSs is defined to measure the uncertainty of a PLTS. Meanwhile, a cross entropy between PLTSs is presented. Based on the total entropy and the cross entropy of PLTSs, DMs' weights and attribute weights are determined objectively. Finally, an improved PL-PROMETHEE method is developed by defining new preference functions and a total preference index. A car sharing platform selection is operated at length to illustrate the applications and advantages of the proposed method.
Apart from solving the selection of car sharing platform, the proposed method can be applied into many decision making fields, such as financial management (Kou 2019a, (i) The sufficiency obviously holds. Next, we prove the necessity. When h(L(p)) = 0, suppose #L(p) > 1. In virtue of Eq. (9), we have g L (k) −ḡ = 0 for any k = 1, 2, . . . , #L(p).
(iv) This is obviously true.
G. Xu is an associate professor at College of Science, Guilin University of Technology. She received her master's degree from school of Mathematics and Information Science, Guangxi University, China, 2007, and earned her PhD degree from College of Information Technology, Jiangxi University of Finance and Economics, China, 2017. Up until now, she has contributed 10 journal articles to professional journals. Her research interests include fuzzy mathematics and decision making.
S. Wan is a professor at College of Information Technology, Jiangxi University of Finance and Economics. He received his PhD degree from College of Information Technology, Nankai University, China, 2005. Up until now, he has contributed over 80 journal articles to professional journals. His research interests include decision making, information fusion and supply chain management.
J. Dong is a professor at School of Statistics, Jiangxi University of Finance and Economics. She received her PhD degree from School of Mathematics Sciences, Nankai University, China, 2013. Up until now, she has contributed over 30 journal articles to professional journals. Her research interests include decision making and graph theory.