An extended minimax absolute and relative disparity approach to obtain the OWA operator weights

Abstract

Determining the ordered weighted averaging (OWA) operator weights is important in decision making applications. Several approaches have been proposed in the literature to obtain the associated weights. This paper provides two new disparity models to obtain the associated weights, which is determined by considering the absolute deviation and relative deviation of any distinct pairs of weights. The proposed mathematical models improve the existing minimize disparity approach and chi-square method, which is suggested by Amin and Emrouzenjad (2006, 2010) and Wang (2007). A numerical example and an application in search engines prove the usefulness of the generated OWA operator weights.

Keywords

OWA operator absolute difference relative difference

1 Introduction

The ordered weighted averaging (OWA) operator defined by Yager [1] provides a general class of parametric aggregation operators that includes the min, max and average, and has shown to be useful for modeling many different kinds of aggregation problems. The OWA operator has been used in a wide range of applications such as neural networks [2, 3], fuzzy logic controllers [4, 5], decision making [6 –10] and data mining [11, 12], etc. To apply the OWA operator, a very crucial issue is the determination of the weights of the operator.

A number of techniques have been suggested for generating the weights. O’Hagan [13] first determined the OWA operator weights and suggested a maximum entropy method. Fullér and Majlender [14] showed that the maximum entropy model could be transformed into a polynomial equation that can be solved analytically. Fullér and Majlender [15] also suggested a minimum variance method to obtain the minimal variability OWA operator weights.

Liu [16, 17] presented the minimal variability OWA operator generating method with the equidifferent OWA operator. Wang and Parkan [18] proposed a linear programming model with a minimax disparity approach to obtain the OWA operator. Majlender [19] extended the maximum entropy method to Rényi entropy and proposed a maximum Rényi entropy OWA operator. Wang et al. [20] proposed minimize the relative disparity method to determine the OWA operator weights.

As Wang and Parkan [18] mentioned to determine the OWA operator weights when the α= orness(W) constraint is taken into consideration, this models above mentioned could be interpreted as making all the weights as close to each other as possible for the given degree of orness. So they imposed the disparities between two adjacent weights are made as small as possible. Minimizing the difference between two adjacent weights is a very strong constraint in the current disparity approach that renders the weights less usable in practice [21]. Thus, the minimax disparity approach was also extended by Wang et al. [20], Amin and Emrouznejad [21, 22]. Amin [23] explored some properties of the extended minimax disparity model. This paper aims to relax this condition and improves the relative disparity method by minimax absolute and relative disparity of any distinct pairs of weights instead of adjacent weights.

The remainder of this paper is organized as follows. In Section 2, we give a review of the OWA operator disparity models for determining its weights. In Section 3, we present two new models for determining OWA weights. In Section 4, we present one numerical example to examine the two methods. In Section 5, we use search engine selection problem to illustrate the proposed method. The conclusions are discussed in Section 6.

2 The OWA operator disparity models

An ordered weighted aggregation (OWA) operator to aggregate information was introduced by Yager [1]. This operator has been investigated in many documents and used in an astonishingly wide range of applications.

Definition 1. An OWA operator of dimension n is a mapping, OWA : Rⁿ → R, that has an associated n vector = (w₁, w₂, …, w_n) such that w_i ∈ [0, 1] and $\sum_{i = 1}^{n} w_{i} = 1$ . Furthermore, ${OWA}_{w} (a_{1}, a_{2}, \dots, a_{n}) = \sum_{j = 1}^{n} w_{j} b_{j}$ (1) where b_j is the jth largest element of the collection of the aggregated objects a₁, a₂, …, a_n. Yager further introduced a characterizing measure called orness measure associated with the weighting vector w of an OWA operator, where the orness measure of the aggregation is defined as $orness (w) = \frac{1}{n - 1} \sum_{i = 1}^{n} (n - i) w_{i}$ (2)

Which lies in the unit interval [0,1] and characterizes the degree to which the aggregation is like an or operation.

One important issue in the OWA operator is to determine its associated weights. To determine OWA operator weights, O’Hagan [13] suggested a maximum entropy method (MEM), which requires the solution of the following constrained nonlinear optimization model (3).

Fullér and Majlender [15] proposed a minimum variance method (MVM), which minimizes the variance of OWA operator weights under a given level of orness. Their method requires the solution of the following mathematical programming model (4).

$\begin{matrix} Max Disp (W) = - \sum_{i = 1}^{n} w_{i} ln w_{i} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, \\ 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (3) and

$\begin{matrix} {Min D}^{2} (W) = \frac{1}{n} \sum_{i = 1}^{n} {(w_{i} - \frac{1}{n})}^{2} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, \\ 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (4)

Wang and Parkan [18] proposed a minimax disparity approach (MDA) for generating OWA operator weights. They presented the following model for minimizing the maximum disparity:

$\begin{matrix} Min {\underset{i = 1, \dots, n - 1}{Max} | w_{i} - w_{i + 1} |} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, \\ 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (5)

There is a common characteristic for above mentioned models. That is the OWA operator weights should be made as equally important as possible. So, based on such a characteristic, Wang and Luo [20] propose the chi-square model (CSM) for obtaining OWA operator weights.

$\begin{matrix} Min J (W) = \sum_{i = 1}^{n - 1} (\frac{w_{i}}{w_{i + 1}} + \frac{w_{i + 1}}{w_{i}} - 2) \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (6)

Amin and Emrouznejad [21] propose an extended minimax disparity model (EMDM) by minimax disparity of any distinct pairs of weights instead of adjacent weights. That is:

$\begin{matrix} Min {\underset{\begin{matrix} i = 1, \dots, n - 1 \\ j = i + 1, \dots, n \end{matrix}}{Max} | w_{i} - w_{j} |} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (7)

Emrouznejad and Amin [22] propose a new disparity model imposing less restriction of disparity between w_j and w_i as follows:

$\begin{matrix} Min \sum_{i = 1}^{n - 1} \sum_{j = i + 1}^{n} δ_{ij} \\ s . t . - δ_{ij} \leq w_{i} - w_{j} \leq δ_{ij} \\ orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \\ δ_{ij} \geq 0, i = 1, \dots, n - 1, j = i + 1, \dots, n \end{matrix}$ (8)

The following section provides two new disparity OWA weight determination models as an extension to the above models.

3 The extended minimax absolute and relative disparity model for obtaining OWA operator weights

As Wang et al. [20] mentioned to determine the OWA operator weights, when the orness constraint is taken into consideration, models (5) and (8) could be understood as making all the weights as close to each other as possible under a given degree of orness. In this paper, two difference models between any distinct pairs of weights are made. We propose the follow model for minimizing the maximum disparity:

$\begin{matrix} Min J_{1} (W) = min {max_{\begin{matrix} i = 1, \dots, n - 1 \\ j = i + 1, \dots, n \end{matrix}} | \frac{w_{i}}{w_{j}} - 1 |} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (9) and

$\begin{matrix} Min J_{2} (W) = min {max_{\begin{matrix} i = 1, \dots, n - 1 \\ j = i + 1, \dots, n \end{matrix}} | \frac{w_{i}}{w_{j}} - 1 | \\ + max_{\begin{matrix} i = 1, \dots, n - 1 \\ j = i + 1, \dots, n \end{matrix}} | w_{i} - w_{j} |} \\ s . t . orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (10)

The two new models both produce as equally important OWA operator weights as possible for a given degree of orness. Note that model (9) is an extension of model (6), which minimizes the maximum disparity any two weights under a given level of orness. The model (10) is an extension of model (7), which takes into account the absolute deviation, relative deviation, and any distinct pairs of OWA operator weights. For convenience, we refer to model (9) as the extended relative difference model (ERDM) for determining the OWA operator weights and model (10) as the extended relative and absolute difference model (ERADM). The ERDM (9) and ERADM (10) can be reconfigured as the follow models:

$\begin{matrix} Min ρ \\ s . t . 1 - ρ \leq \frac{w_{i}}{w_{j}} \leq 1 + ρ, \\ i = 1, \dots, n - 1, j = i + 1, \dots, n \\ orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (11) and

$\begin{matrix} Min {ρ + ξ} \\ s . t . 1 - ρ \leq \frac{w_{i}}{w_{j}} \leq 1 + ρ, \\ i = 1, \dots, n - 1, j = i + 1, \dots, n \\ - ξ \leq w_{i} - w_{j} \leq ξ, \\ i = 1, \dots, n - 1, j = i + 1, \dots, n \\ orness (W) = \sum_{i = 1}^{n} \frac{n - i}{n - 1} \cdot w_{i} = α, 0 \leq α \leq 1 \\ \sum_{i = 1}^{n} w_{i} = 1, 0 \leq w_{i} \leq 1, i = 1, \dots, n \end{matrix}$ (12)

Models (11) and (12) are both nonlinear and can be solved by using LINGO or MATLAB software package. Note that model (11) and (12) are not applicable to α= 0 and α= 1, which are two extreme degrees of orness, but can be approached by using α= 0.0001 and α= 0.9999. We will use these two values to represent the two extreme cases in the next section.

4 Numerical example

In this section, we examine the two new models with a numerical example and verify their applicability in determining OWA operator weights. The weights determined by the two models are also compared with above mentioned models and the NOWA weights [24]. Suppose n = 5 and it is needed to determine the OWA operator weights satisfying different degrees of orness: α= 0.5, 0.6, … , 1, which are provided by the decision maker.

Tables 1 and 2 show the OWA operator weights determined by models (9) and (10), respectively, which are also depicted in Figs. 1 and 2. The models are solved by using MATLAB software package.

It is easy to see from Tables 1 –4 that the OWA operator weights determined by different methods are slightly different from one another, but there are no significant differences among them. This shows the use of the two new models to determine OWA operator weights is not only feasible but alsoeffective.

To find the slight differences among the above mentioned methods, we consider the distributions of the OWA operator weights in Table 5. Take α= 0.7 for example. Table 5 shows the distribution of the OWA operator weights determined by the eight methods under the given degree of orness (α) = 0.7, from which it can be seen very clearly that the weights determined by the two new methods vary neither in the form of exponential nor in the form of arithmetical progression, but in a very general way. This is what we are expecting because there is no evidence to support the OWA operator weights to follow a very regular distribution, and must be adjacent weight minimum requirements [20].

5 Application of the proposed OWA weights in internet search engine

Keyhanipour et al. [25] used optimistic OWA (O-OWA) operator to measure the performance of the search engines by factors such as average click rate, total relevancy of the returned results and the variance of the clicked results.

Relevance judgments were done according to users’ preferences, e.g., overall value of the returned results. Relevant items were scored as 2, irrelevant ones as 0 and undecided items as 1. One sample from the judgment of a user which is extracted from Keyhanipour is shown in the following matrix.

In their method they used O-OWA operator as introduced in Yager [1] to combine the decisions of underlying search engines which have different degrees of importance in a certain category. Emrouznejad [26] also used most preferred OWA (MP-OWA) operator to rank the search engines.

To compare the results with other models, Table 7 shows the score is given to each search engine using ERDM and ERADM as developed in this paper. According to the models (9) and (10), for an orness level of 0.75, the highest score is assigned to the WebFusion-OWA search engine, and the lowest score is given to the LookSmart or Lycos search engine, which is consistent with the method of O-OWA and MP-OWA. Some disagreements exist between the models-for example search engine AltaVista is given a very high rank (3rd place) by models (9) and (10), and a very low rank (9th or 6th place) using O-OWA and MP-OWA. In fact, from the matrix of judgment of a user for a sample query results (Table 6), we can find search engines AltaVista and ProFusion have no significant difference. Therefore, the rank of our method is more reasonable.

6 Conclusions

In this paper, we have proposed two new models for determining OWA operator weights. The two new models prove to be practical and effective and can produce the OWA operator weights that are very close to those obtained by existing methods. But the weights determined by considering the absolute deviation and relative deviation of any distinct pairs of weights, thus the two new models do not follow a regular distribution and therefore make more sense.

In future research, we will further investigate the related properties of the absolute deviation and relative deviation model through an analytic solution. A detailed application study of different models may also be a topic of interest.

Footnotes

Acknowledgments

The author is very grateful to the editor and the anonymous reviewers for their constructive comments and suggestions that have led to an improved version of this paper. The work was supported in part by the Natural Science Foundation of Jiangsu Province of China (No. BK20130242), the Fundamental Research Funds for the Central Universities (Nos. 2015B28014, 2015B23914).

References

Yager

R.R.

, On ordered weighted averaging aggregation operators in multi-criteria decision making, IEEE Transactions on Systems, Man and Cybernetics18 (1988), 183–190.

Amin

G.R.

, and Emrouznejad

, Parametric aggregation in ordered weighted averaging, International Journal of Approximate Reasoning52 (2011), 819–827.

Amin

G.R.

, and Emrouznejad

, Optimizing search engines results using linear programming, Expert System with Application38 (2011), 11534–11537.

Liu

, Some properties of the weighted OWA operator, IEEE Transactions on Systems, Man and Cybernetics, Part B36 (2006), 118–127.

Yager

R.R.

, Time series smoothing and OWA aggregation, IEEE Transactions on Fuzzy Systems16 (2008), 994–1007.

Zhou

L.G.

, Chen

H.Y.

, and Liu

J.P.

, Generalized multiple averaging operators and their applications to group decision making, Group Decision and Negotiation22 (2013), 331–358.

Zhou

L.G.

, Chen

H.Y.

, and Liu

J.P.

, Generalized logarithmic proportional averaging operators and their applications to group decision making, Knowledge-Based Systems36 (2012), 268–279.

Zhou

L.G.

, Tao

Z.F.

, Chen

H.Y.

, and Liu

J.P.

, Generalized ordered weighted logarithmic harmonic averaging operators and their applications to group decision making, Soft Computing19 (2015), 715–730.

Zeng

S.Z.

, and Xiao

, TOPSIS method for intuitionistic fuzzy multiple-criteria decision making and its application to investment selection, Kybernetes45 (2016), 282–296.

10.

Zeng

S.Z.

, Chen

J.P.

, and Li

X.S.

, A hybrid method for pythagorean fuzzy multiple-criteria decision making, International Journal of Information Technology & Decision Making15 (2016), 403–422.

11.

Sicilia

M.A.

, and García

, Barriocanal and S. Sánchez Alonso, Empirical assessment of a collaborative filtering algorithm based on OWA operators, International Journal of Intelligent Systems23 (2008), 1251–1263.

12.

Torra

, OWA operators in data modeling and reidentification, IEEE Transactions on Fuzzy Systems12 (2004), 652–660.

13.

O’Hagan

, Using maximum entropy-based weighted averaging to construct a fuzzy neuron, Proceedings of the 24th Annual IEEE Asilomar Conference on Signals Systems Computation, Pacific Grove1990, pp. 618–623.

14.

Fullér

, and Majlender

, An analytic approach for obtaining maximal entropy OWA operator weights, Fuzzy Sets and Systems124 (2001), 53–57.

15.

Fullér

, and Majlender

, On obtaining minimal variability OWA operator weights, Fuzzy Sets and Systems136 (2003), 203–215.

16.

Liu

, On the properties of equidifferent OWA operator, International Journal of Approximate Reasoning43 (2006), 90–107.

17.

Liu

, The solution equivalence of minimax disparity and minimum variance problems for OWA operators, International Journal of Approximate Reasoning45 (2007), 68–81.

18.

Wang

Y.M.

, and Parkan

, A minimax disparity approach for obtaining OWA operator weights, Information Sciences175 (2005), 20–29.

19.

Majlender

, OWA operators with maximal Rényi entropy, Fuzzy Sets and Systems155 (2005), 340–360.

20.

Wang

Y.M.

, Luo

, and Liu

, Two new models for determining OWA operator weights, Computers & Industrial Engineering52 (2007), 203–209.

21.

Amin

, and Emrouznejad

, An extended minimax disparity to determine the OWA operator weights, Computers & Industrial Engineering50 (2006), 312–316.

22.

Emrouznejad

, and Amin

, Improving minimax disparity model to determine the OWA operator weights, Information Sciences180 (2010), 1477–1485.

23.

Amin

, Notes on properties of the OWA weights determination model, Computers & Industrial Engineering52 (2007), 533–538.

24.

Z.S.

, An overview of methods for determining OWA weights, International Journal of Intelligent Systems20 (2005), 843–865.

25.

Keyhanipour

A.H.

, Moshiri

, Kazemian

, Piroozmand

, and Lucas

, Aggregation of web search engines based on users’ preferences in WebFusion, Knowledge-Based Systems20 (2007), 321–328.

26.

Emrouznejad

, MP-OWA: The most preferred OWA operator, Knowledge-Based Systems21 (2008), 847–851.