# Malmquist Productivity Index by Extended VIKOR Method Using Interval Numbers

Authors: Mohammad Fallah; Amir Mohajeri; Esmaeil Najafi
Journal: Abstract and Applied Analysis
 (2013)

---

## Abstract

The VIKOR method was developed for multicriteria optimization of complex systems. It determines the compromise ranking list and the compromise solution obtained with the given weights. This method focuses on ranking and selecting from a set of alternatives in the presence of conflicting criteria. Here, the VIKOR method is used for two timest and t+1. In order to calculate the progress or regression via Malmquist productivity index, the positive and negative ideals at times t and t+1 are calculated first. Then we introduce the multi-criteria ranking index based on the particular measure of “closeness” to the ideal solution and calculate the separation of each alternative from the ideal solution at times t and t+1. Then we use the Malmquist productivity index to calculate the progress or regression of all alternatives. In this paper, productivity of alternatives available in decision matrix with interval numbers and their improvement or deterioration is researched. To achieve this practical goal, use of extended VIKOR is made to calculate Malmquist productivity index for multicriteria decision-making (MCDM) problem with interval numbers, and by applying Malmquist productivity index, productivity rate of growth for alternatives is calculated. Finally, a numerical example illustrates and clarifies the main results developed in this paper.


---

## Body

## 1. Introduction

Multi criteria optimization is the process of determining the best feasible solution according to the established criteria. Practical problems are often characterized by several conflicting criteria, and there may be no solution satisfying all criteria simultaneously. Thus, the solution is a set of noninferior solutions or a compromise solution according to the decision maker’s preferences. The compromise solution was established by Yu [1] and Zeleny [2] for a problem with conflicting criteria, and it can be helping the decision makers to reach a final solution. The compromise solution is a feasible solution, which is the closest to the ideal, and compromise means an agreement established by mutual concessions.A multiattribute decision-making (MADM) problem can be defined as shown in Table1 where A1,A2,…,Am are possible alternatives among which decision makers have to choose, C1,C2,…,Cn are criteria with which alternative performance is measured, fij is the rating of alternative Ai with respect to criterion Cj, and wj is the weight of criterion Cj [3–5].Table 1
Matrix format of a MADM problem.
C
1
C
2
…
C
n
A
1
f
11
f
12
…
f
1
n
A
2
f
21
f
22
…
f
2
n
⋮
⋮
⋮
⋮
⋮
A
m
f
m
1
f
m
2
…
f
m
n
W
=
[
w
1
,
w
2
,
…
,
w
n
].In classical MCDM methods, the ratings and the weights of the criteria are known precisely, whereas in the real world, in an imprecise and uncertain environment, it is an unrealistic assumption that the knowledge and representation of a decision maker or expert are so precise. To describe and treat imprecise and uncertain elements present in a decision problem, fuzzy and stochastic approaches are frequently used. In the literature, in the works of fuzzy decision making, fuzzy parameters are assumed to be with known membership functions [6–12], and in stochastic decision making, parameters are assumed to have known probability distributions [13–16]. However, in reality to a decision maker (DM), it is not always easy to specify the membership function or probability distribution in an inexact environment. At least in some of the cases, the use of interval numbers may serve the purpose better. An interval number can be thought as an extension of the concept of a real number [17].The TOPSIS method to solve decision-making problems with interval data has been extended by Jahanshahloo et al. [18]. Sayadi et al. [19] have extended the concept of VIKOR method to develop a methodology for solving MADM problems with interval numbers.In this paper, we extend the concept of VIKOR method to develop a methodology for the calculation of progress or regression via Malmquist productivity index for MADM problems with interval numbers.


## 2. Extended VIKOR Method for Decision-Making Problem with Interval Numbers

As it was said in the introduction, the interval numbers are more suitable to deal with the decision-making problems in the imprecise and uncertain environment, because they are the simplest form of representing uncertainty in the decision matrix. The interval numbers require the minimum amount of information about the values of attributes. Specifying an interval for a parameter in decision matrix indicates that the parameter can take any value within the interval. Note that, the interval numbers does not indicate how probable it is to the value to be in the interval, nor does it indicate which of the many values in the interval is the most likely to occur [20]. An interval number signifies the extent of the tolerance or a region that the parameter can possibly take. An extensive research on interval mathematics and its applications can be found in [17, 21, 22]. More information about the interval numbers and its differences with other methods of representing uncertainty such as probability and fuzzy theory can be found in [23–25].According to these facts, when determining the exact values of the attributes is difficult or impossible, it is more appropriate to consider them as interval numbers. Therefore, in the present paper, we extend the VIKOR method to calculate Malmquist productivity index for MADM problem with interval numbers.  Then, applying proposed Malmquist productivity index, we calculate the progress or regression of all alternatives. To do this, we explain the extended VIKOR method for MADM problem with interval numbers proposed by Sayadi et al. [19] in the first place. Then, we introduce a new algorithm for the calculation of progress or regression via extended VIKOR and Malmquist productivity index for MADM problem with interval numbers.Now, we suppose that a decision matrix with interval numbers has the form shown in Table2 where fij is not known exactly, and only we know thatfij∈[fijL,fijU], fijL, and fijU are lowest and highest value of fij, respectively. The extended VIKOR method consists of the following steps:(a)
Determine the positive ideal solution (PIS) and negative ideal solution (NIS):(1)A*={f1*,…,fn*}={(maxifijU∣j∈I)or(minifijL∣j∈J)},j=1,2,…,n,A-={f1-,…,fn-}={(minifijL∣j∈I)or(maxifijU∣j∈J)},j=1,2,…,n,
where A*and A- represent the positive ideal and the negative ideal, respectively. I is associated with the benefit criteria, and J is associated with the cost criteria.(b)
In this step, compute[SiL,SiU] and [RiL,RiU] intervals as follows:
(2)SiL=∑j∈Iwj(fj*-fijU)(fj*-fj-)+∑j∈Jwj(fijL-fj*)(fj--fj*),i=1,2,…,m,SiU=∑j∈Iwj(fj*-fijL)(fj*-fj-)+∑j∈Jwj(fijU-fj*)(fj--fj*),i=1,2,…,m,RiL=max{wj(fj*-fijU)(fj*-fj-)∣j∈I,wj(fijL-fj*)(fj--fj*)∣j∈J},i=1,2,…,m,RiU=max{wj(fj*-fijL)(fj*-fj-)∣j∈I,wj(fijU-fj*)(fj--fj*)∣j∈J},i=1,2,…,m,
where wj are the weights of criteria, expressing their relative importance.(c)
Compute the intervalQi=[QiL,QiU]; i=1,2,…,m, by the following relations:
(3)QiL=v(SiL-S*)(S--S*)+(1-v)(RiL-R*)(R--R*),QiU=v(SiU-S*)(S--S*)+(1-v)(RiU-R*)(R--R*),
where
(4)S*=miniSiL,S-=maxiSiU,R*=miniRiL,R-=maxiRiU,
where v is introduced as weight of the strategy of “the majority of criteria” (or “the maximum group utility”); here suppose that v=0.5.Table 2
Matrix format of a MADM problem with interval data.
C
1
C
2
…
C
n
A
1
[
f
11
L
,
f
11
U
]
[
f
12
L
,
f
12
U
]
…
[
f
1
n
L
,
f
1
n
U
]
A
2
[
f
21
L
,
f
21
U
]
[
f
22
L
,
f
22
U
]
…
[
f
2
n
L
,
f
2
n
U
]
⋮
⋮
⋮
⋮
⋮
A
m
[
f
m
1
L
,
f
m
1
U
]
[
f
m
2
L
,
f
m
2
U
]
…
[
f
m
n
L
,
f
m
n
U
]
W
=
[
w
1
,
w
2
,
…
,
w
n
].


## 3. Calculate Progress or Regression via Extended VIKOR and Malmquist Index for MADM Problem with Interval Numbers

Productivity management is one of the major sources of sustainable organizational effectiveness and a systematic understanding of the factors that affecting productivity is very important. The measurement and analysis of productivity change are always a controversial topic and have enjoyed a great deal of interest among organizational researchers and practitioners.Originally, Sten Malmquist, a Swedish economist and statistician, proposed a quantity index for measuring the standard of living for the purpose of consumption analysis [26], but later on the Malmquist index and its variations have mainly been used in the field of production analysis. However, most of these studies have concentrated on total factor productivity (TFP) measurement, although in the spirit of the original proposition, Malmquist index could be applied in other areas equally well.In early work in this field, productivity change was explained in terms of technical change, but recently it has become widely accepted that efficiency change can also contribute to it. In this framework, the Malmquist index was first introduced in productivity literature by Caves et al. [27]. The Malmquist productivity index was introduced as a theoretical index based on Shephard’s distance function (1970) and is widely developed and used in many fields. Nishimizu and Page [28] used a parametric programming approach to compute the index for the first time in the empirical context which then is further developed and popularized as an empirical index by Fare et al. [29, 30]. They decomposed productivity change into a part attributable to change of technical efficiency and technical change and used nonparametric mathematical programming models for its computation.So, Malmquist productivity index is defined by patching together efficiency variations of each alternative and technical variation. In addition to comparing the efficiency of each alternative at present and in the past to calculate the corresponding advancement or regression, bases itself on the society transcendence movement in these two-time intervals. Meaning that in the competition, the advancement or regression of the alternatives under consideration is defined on the advancement or regression of the society transcendence. The variations of the society transcendence are called technology variations. Therefore it is possible to examine the causes of productivity change and to see whether the productivity has improved, through a more efficient use of existing technology or through technological progress. Here, the Malmquist index is defined using the distance function and used to identify productivity differences between one firm over two-time periods. To estimate technical efficiency changes and technical changes over the period in question, the decomposed Malmquist productivity index was used. According to Fare et al. [31, 32], the Malmquist productivity index between time periods t and t+1 can be decomposed into two components and defined as
(5)Mt,t+1(yt,yt+1,xt,xt+1)︸M.P.I=Dt+1(yt+1,xt+1)Dt(yt,xt)︸Efficiencyvariations×[Dt(yt,xt)Dt+1(yt,xt)×Dt(yt+1,xt+1)Dt+1(yt+1,xt+1)]1/2︸Technologyvariations,
where the notation D represents the distance function and the value of M is the Malmquist productivity index. The term outside the brackets (efficiency variations) is a ratio of two distance functions, which measures the change in the measure of the Farell technical efficiency between period t and t+1. The square root term (technology variations) is a measure of the technical change in the production technology. It is an indicator of the distance covered by the efficient frontier from one period to another and thus a measure of technological improvements between the periods. The term (efficiency variations) is greater than, equal to, or less than one if the producer is moving closer to, unchanging, or diverging from the production frontier, respectively. The square root term (technology variations) is greater than, equal to, or less than one when the technological best practice is improving, unchanged, or deteriorating, respectively.In Figure1 the calculation of Malmquist productivity index is illustrated, where a single output y is produced using a single input x. It is assumed that the constant returns to scale (CRS) technology. The firm A produces at the point A1 in the first period and at the point A2 in the second period. The firm A is technically inefficientin the first period as the point A1 is below the frontier for that period. In the second period the point A2 is on the frontier and thereby firm A is technically efficient. The technical change includes a time component and involves advances in technology, which is represented by an upward shift in the production frontier from first period to the second period.Figure 1
Malmquist productivity index [33].Using (5) it is
(6)M1,2(yA1,yA2,xA1,xA2)︸M.P.I=yA2/yA2yA1/y1︸Efficiencyvariations[yA1/y1yA1/y2×yA2/y2yA2/yA2]1/2︸Technologyvariations.
It is easy to calculate that the value of M is greater than one (i.e., M>1), which implies productivity growth. In (6) the terms efficiency variations and technology variations are also greater than one (i.e., efficiency variations > 1 and technology variations > 1). Therefore the firm A experienced the positive technical efficiency change and technological change from one period to another.Efficiency variation is the individual observations shifting relative to the best frontier (i.e., how far the average firm is from the best-practice frontier), and technology variation is the shift in the best-practice frontier. Efficiency variation above unity means that the firm has moved closer to the efficient (best-practice) frontier and thus measures “catching up” (or “falling behind” if it is less than unity). Technology variation above unity indicates technological progress, meaning that the efficient frontier has shifted out compared to the previous period, and a value below unity suggests technology regress.As we mentioned previously, Malmquist productivity index measures the progress or regression in the period of(t,t+1). If Malmquist index on the basis of minimization of production factors was less than one, it indicates productivity decrease; on the other hand, if on the basis of maximization of production factors, the Malmquist index or any of its elements were less than one, and it signifies productivity getting worse, while if the index is bigger than one, it indicates productivity improvement.The importance of the above technique is that, in an industry occasionally, companies that have faced similar productivity decrease in an specific period of time; by evaluating productivity elements, it can be observed that productivity decrease of one company was mainly due to lack of technological advancements and nonexistence of necessary investments and for the other corporation was because of the decline of the size of activities and the limitation of productivity scale. Therefore equal decline in productivity does not signify a common reason, and it might have a specific reason for each company.Malmquist productivity index on interval data was first introduced by Hosseinzadeh Lotfi and Ghasemi [34]. They used Data Envelopment Analysis (DEA) to calculate Malmquist productivity index and by applying Malmquist productivity index, productivity rate of growth for telecommunication companies was calculated. If the ratings of each alternative with respect to each criterion are interval, we obtained interval models. Furthermore, the value of Malmquist productivity index for each alternative will be related to one interval that the lowest and highest values of this interval have the following form:
(7)M.P.Ii=[M.P.IiL,M.P.IiU].
In the present paper, we propose the extended VIKOR method to calculate Malmquist productivity index for MADM problem with interval numbers. Then, applying Malmquist productivity index, we calculate progress or regression of all alternatives. To do this, suppose that a decision matrix with interval data in two times of t and t+1 has the form shown in Table 3, where fijt=[fijt,L,fijt,U] and fijt+1=[fijt+1,L,fijt+1,U] are the lowest and highest values of fij in times t and t+1, respectively (Table 3).Table 3
Matrix format of a MADM problem with interval data in timest and t+1.
t
t
+
1
C
1
C
2
…
C
n
C
1
C
2
…
C
n
A
1
[
f
11
t
,
L
,
f
11
t
,
U
]
[
f
12
t
,
L
,
f
12
t
,
U
]
…
[
f
1
n
t
,
L
,
f
1
n
t
,
U
]
[
f
11
t
+
1
,
L
,
f
11
t
+
1
,
U
]
[
f
12
t
+
1
,
L
,
f
12
t
+
1
,
U
]
…
[
f
1
n
t
+
1
,
L
,
f
1
n
t
+
1
,
U
]
A
2
[
f
21
t
,
L
,
f
21
t
,
U
]
[
f
22
t
,
L
,
f
22
t
,
U
]
…
[
f
2
n
t
,
L
,
f
2
n
t
,
U
]
[
f
21
t
+
1
,
L
,
f
21
t
+
1
,
U
]
[
f
22
t
+
1
,
L
,
f
22
t
+
1
,
U
]
…
[
f
2
n
t
+
1
,
L
,
f
2
n
t
+
1
,
U
]
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
⋮
A
m
[
f
m
1
t
,
L
,
f
m
1
t
,
U
]
[
f
m
2
t
,
L
,
f
m
2
t
,
U
]
…
[
f
m
n
t
,
L
,
f
m
n
t
,
U
]
[
f
m
1
t
+
1
,
L
,
f
m
1
t
+
1
,
U
]
[
f
m
2
t
+
1
,
L
,
f
m
2
t
+
1
,
U
]
…
[
f
m
n
t
+
1
,
L
,
f
m
n
t
+
1
,
U
]
W
=
[
w
1
,
w
2
,
…
,
w
n
].Now, use of extended VIKOR is made to calculate Malmquist productivity index for MADM problem with interval numbers. The following steps are defined as follows.Step 1.
Determine the PIS and NIS as follows:(8)A*t={f1*t,…,fn*t}={(maxifijt,U∣j∈I)or(minifijt,L∣j∈J)},j=1,2,…,n,A-t={f1-t,…,fn-t}={(minifijt,L∣j∈I)or(maxifijt,U∣j∈J)},j=1,2,…,n,A*t+1={f1*t+1,…,fn*t+1}={(maxifijt+1,U∣j∈I)or(minifijt+1,L∣j∈J)},j=1,2,…,n,A-t+1={f1-t+1,…,fn-t+1}={(minifijt+1,L∣j∈I)or(maxifijt+1,U∣j∈J)},j=1,2,…,n,
where (A*t,A-t) and (A*t+1,A-t+1) represent the positive ideal and the negative ideal in time t and t+1,  respectively.Step 2.
Compute the separation measuring. The separation of each alternative from the ideal solution is given as(9)Si(t)(t),L=∑j∈Iwj(fj*t-fijt,U)(fj*t-fj-t)+∑j∈Jwj(fijt,L-fj*t)(fj-t-fj*t),i=1,2,…,m,Si(t)(t),U=∑j∈Iwj(fj*t-fijt,L)(fj*t-fj-t)+∑j∈Jwj(fijt,U-fj*t)(fj-t-fj*t),i=1,2,…,m,Si(t+1)(t+1),L=∑j∈Iwj(fj*t+1-fijt+1,U)(fj*t+1-fj-t+1)+∑j∈Jwj(fijt+1,L-fj*t+1)(fj-t+1-fj*t+1),i=1,2,…,m,Si(t+1)(t+1),U=∑j∈Iwj(fj*t+1-fijt+1,L)(fj*t+1-fj-t+1)+∑j∈Jwj(fijt+1,U-fj*t+1)(fj-t+1-fj*t+1),i=1,2,…,m,Ri(t)(t),L=max{wj(fj*t-fijt,U)(fj*t-fj-t)∣j∈I,wj(fijt,L-fj*t)(fj-t-fj*t)∣j∈J},i=1,2,…,m,Ri(t)(t),U=max{wj(fj*t-fijt,L)(fj*t-fj-t)∣j∈I,wj(fijt,U-fj*t)(fj-t-fj*t)∣j∈J},i=1,2,…,m,Ri(t+1)(t+1),L=max{wj(fj*t+1-fijt+1,U)(fj*t+1-fj-t+1)∣j∈I,wj(fijt+1,L-fj*t+1)(fj-t+1-fj*t+1)∣j∈J(fijt+1,L-fj*t+1)(fj-t+1-fj*t+1)},i=1,2,…,m,Ri(t+1)(t+1),U=max{wj(fj*t+1-fijt+1,L)(fj*t+1-fj-t+1)∣j∈I,wj(fijt+1,U-fj*t+1)(fj-t+1-fj*t+1)∣j∈J(fijt+1,L-fj*t+1)(fj-t+1-fj*t+1)},i=1,2,…,m,
where Si(t)(t)=[Si(t)(t),L,Si(t)(t),U] and Si(t+1)(t+1)=[Si(t+1)(t+1),L,Si(t+1)(t+1),U] represent the lowest and highest values of distance to the positive ideal and Ri(t)(t)=[Ri(t)(t),L,Ri(t)(t),U] and Ri(t+1)(t+1)=[Ri(t+1)(t+1),L,Ri(t+1)(t+1),U] represent the lowest and highest values of most inferior distance to the positive ideal, respectively, in time, t and t+1.
Now, we have(10)Si(t)(t+1),L=∑j∈Iwj(fj*t+1-fijt,U)(fj*t+1-fj-t+1)+∑j∈Jwj(fijt,L-fj*t+1)(fj-t+1-fj*t+1),i=1,2,…,m,Si(t)(t+1),U=∑j∈Iwj(fj*t+1-fijt,L)(fj*t+1-fj-t+1)+∑j∈Jwj(fijt,U-fj*t+1)(fj-t+1-fj*t+1),i=1,2,…,m,Si(t+1)(t),L=∑j∈Iwj(fj*t-fijt+1,U)(fj*t-fj-t)+∑j∈Jwj(fijt+1,L-fj*t)(fj-t-fj*t),i=1,2,…,m,Si(t+1)(t),U=∑j∈Iwj(fj*t-fijt+1,L)(fj*t-fj-t)+∑j∈Jwj(fijt+1,U-fj*t)(fj-t-fj*t),i=1,2,…,m,Ri(t)(t+1),L=max{wj(fj*t+1-fijt,U)(fj*t+1-fj-t+1)∣j∈I,wj(fijt,L-fj*t+1)(fj-t+1-fj*t+1)∣j∈J},i=1,2,…,m,Ri(t)(t+1),U=max{wj(fj*t+1-fijt,L)(fj*t+1-fj-t+1)∣j∈I,wj(fijt,U-fj*t+1)(fj-t+1-fj*t+1)j∈J},i=1,2,…,m,Ri(t+1)(t),L=max{wj(fj*t-fijt+1,U)(fj*t-fj-t)∣j∈I,wj(fijt+1,L-fj*t)(fj-t-fj*t)∣j∈J},i=1,2,…,m,Ri(t+1)(t),U=max{wj(fj*t-fijt+1,L)(fj*t-fj-t)∣j∈I,wj(fijt+1,U-fj*t)(fj-t-fj*t)∣j∈J},i=1,2,…,m,
where Si(t)(t+1)=[Si(t)(t+1),L,Si(t)(t+1),U] represent the lowest and highest values of distance Ai in time t from the positive ideal in time t+1 and Si(t+1)(t)=[Si(t+1)(t),L,Si(t+1)(t),U] represent the lowest and highest values of distance Ai in time t+1 from the positive ideal in time t. Also, Ri(t)(t+1)=[Ri(t)(t+1),L,Ri(t)(t+1),U] represent the lowest and highest value of most inferior distance Ai in time t from the positive ideal in time t+1, and Ri(t+1)(t)=[Ri(t+1)(t),L,Ri(t+1)(t),U] represent the lowest and highest values of most inferior distance Ai in time t+1, from the positive ideal in time t.Step 3.
Compute the valuesQi=[QiL,QiU],i=1,2,…,m, in forms of simple and mixed model by the following relations:
(11)Qi(t)(t),L=v(Si(t)(t),L-S*(t)(t))(S-(t)(t)-S*(t)(t))+(1-v)(Ri(t)(t),L-R*(t)(t))(R-(t)(t)-R*(t)(t)),i=1,2,…,m,Qi(t)(t),U=v(Si(t)(t),U-S*(t)(t))(S-(t)(t)-S*(t)(t))+(1-v)(Ri(t)(t),U-R*(t)(t))(R-(t)(t)-R*(t)(t)),i=1,2,…,m,Qi(t+1)(t+1),L=v(Si(t+1)(t+1),L-S*(t+1)(t+1))(S-(t+1)(t+1)-S*(t+1)(t+1))+(1-v)(Ri(t+1)(t+1),L-R*(t+1)(t+1))(R-(t+1)(t+1)-R*(t+1)(t+1)),i=1,2,…,m,Qi(t+1)(t+1),U=v(Si(t+1)(t+1),U-S*(t+1)(t+1))(S-(t+1)(t+1)-S*(t+1)(t+1))+(1-v)(Ri(t+1)(t+1),U-R*(t+1)(t+1))(R-(t+1)(t+1)-R*(t+1)(t+1)),i=1,2,…,m,Qi(t)(t+1),L=v(Si(t)(t+1),L-S*(t)(t+1))(S-(t)(t+1)-S*(t)(t+1))+(1-v)(Ri(t)(t+1),L-R*(t)(t+1))(R-(t)(t+1)-R*(t)(t+1)),i=1,2,…,m,Qi(t)(t+1),U=v(Si(t)(t+1),U-S*(t)(t+1))(S-(t)(t+1)-S*(t)(t+1))+(1-v)(Ri(t)(t+1),U-R*(t)(t+1))(R-(t)(t+1)-R*(t)(t+1)),i=1,2,…,m,Qi(t+1)(t),L=v(Si(t+1)(t),L-S*(t+1)(t))(S-(t+1)(t)-S*(t+1)(t))+(1-v)(Ri(t+1)(t),L-R*(t+1)(t))(R-(t+1)(t)-R*(t+1)(t)),i=1,2,…,m,Qi(t+1)(t),U=v(Si(t+1)(t),U-S*(t+1)(t))(S-(t+1)(t)-S*(t+1)(t))+(1-v)(Ri(t+1)(t),U-R*(t+1)(t))(R-(t+1)(t)-R*(t+1)(t)),i=1,2,…,m,
where
(12)S*(t)(t)=miniSi(t)(t),L,S-(t)(t)=maxiSi(t)(t),U,R*(t)(t)=miniRi(t)(t),L,R-(t)(t)=maxiRi(t)(t),U,S*(t+1)(t+1)=miniSi(t+1)(t+1),L,S-(t+1)(t+1)=maxiSi(t+1)(t+1),U,R*(t+1)(t+1)=miniRi(t+1)(t+1),L,R-(t+1)(t+1)=maxiRi(t+1)(t+1),U,S*(t)(t+1)=miniSi(t)(t+1),L,S-(t)(t+1)=maxiSi(t)(t+1),U,R*(t)(t+1)=miniRi(t)(t+1),L,R-(t)(t+1)=maxiRi(t)(t+1),U,S*(t+1)(t)=miniSi(t+1)(t),L,S-(t+1)(t)=maxiSi(t+1)(t),U,R-(t+1)(t)=miniRi(t+1)(t),L,R-(t+1)(t)=maxiRi(t+1)(t),U.Step 4 (Compute the Malmquist Productivity Index).
There are many different methods that could be used to measure the distance function, which makes up the Malmquist productivity index. These required distance functions can be calculated using either mathematical programming or econometric techniques. The extended VIKOR method to construct Malmquist indices defining distance function based on separation measures is used in the empirical part of this study. So, the Malmquist index is defined using the proposed distance function. An alternative definition of Malmquist productivity index can be presented as the Malmquist productivity index based on separation measure of each alternative from the ideal solution in which distance and the most inferior distance to the positive ideal solution are combined (Qi). The distance function makes it possible to describe a multiinput, multioutput production technology and does not require the profit maximization or cost minimization assumption. So, the distance function (Qi) is presented as a combination of distance to the positive ideal solution (Si) and the most inferior distance to the positive ideal solution (Ri). In the present paper the Malmquist productivity index will be defined using this distance function.To define a proposed distance function, there is a sample ofM firms using Xt∈ℜ+N inputs and Yt∈ℜ+K outputs in the time period t=1,…,T. Using these inputs and outputs data, we calculate the distance and most inferior distance of each firm to the positive ideal solution, (Si) and (Ri), respectively. Afterwards, Qi is defined as a combination of (Si) and (Ri). Multiple inputs and multiple outputs production technology may be defined using the set, P, which represents the set of all Q vectors, Qt=(Q1t,…,Qmt), which can be produced using the (x) and (y) vectors, xt=(x1t,…,xnt) and yt=(y1t,…,ykt) in the time period t=1,…,T. That is,
(13)Pt(xt,yt)={Qt:xt,ytcanproduceQtattimet}t=1,…,T.
The proposed distance function is less than or equal to one (i.e., Q≤1), if and only if Q belongs to the production possibility set of x, y (i.e., Q∈P(x,y)). Note that the distance function is equal to the unit (i.e., Q=1), if Q belongs to the “frontier” of the production possibility set. A firm is considered technically efficient if the distance function equals one.The proposed distance function and production possibility are shown in Figure2, where Q1 and Q2 are produced using (x,y).Figure 2
Separation-based distance function and production possibility set.To construct the VIKOR-based Malmquist productivity index for adjacent periods and MADM problem with interval numbers, it is needed to calculate eight different proposed distance functions,Qi(t)(t),L, Qi(t)(t),U, Qi(t+1)(t+1),L, Qi(t+1)(t+1),U, Qi(t)(t+1),L, Qi(t)(t+1),U, Qi(t+1)(t),L, and Qi(t+1)(t),U.In this section using the extended VIKOR method, we present an interval Malmquist productivity index. At first time we should evaluate the anti-ideal firm distance for each of time periodst and t+1. Taking time period t as the reference period, Qi(t)(t)and Qi(t+1)(t) are measured as proposed distance function.Similarly, it can be seen that we can obtain the anti-ideal firm distance for time periodst and t+1 taking time period t+1 as the reference period (Qi(t+1)(t+1) and Qi(t)(t+1)).Suppose that,Qi(t)(t),Land Qi(t)(t),U are the lower and upper bounds of the interval Qi(t)(t). Then distance interval is denoted by [Qi(t)(t),L,Qi(t)(t),U].Similarly, usingt+1 instead of t, we get the interval distance Qi(t+1)(t+1) as [Qi(t+1)(t+1),L,Qi(t+1)(t+1),U].The interval distance for mixed period measures is given as follows.The interval distance forQi(t+1)(t) is given as [Qi(t+1)(t),L,Qi(t+1)(t),U], and similarly, using t+1 instead of t, we get the interval distance for Qi(t)(t+1)  as  [Qi(t)(t+1),L,Qi(t)(t+1),U]. We know that the Malmquist productivity index can be decomposed into two components and is defined as follows:
(14)(5)⟹M.P.I=Q(t+1)(t+1)Q(t)(t)︸Efficiencyvariations×[Q(t)(t)Q(t)(t+1)×Q(t+1)(t)Q(t+1)(t+1)]1/2︸Technologyvariations.
Note that we obtain an interval for each of factors in above formula. Measuring the Malmquist index as an interval for any firm, the lower and upper bounds of Malmquist productivity index are given as follows:
(15)M.P.IiL=EfficiencyvariationsiL·TechnologyvariationsiL,i=1,2,…,m,M.P.IiU=EfficiencyvariationsiU·TechnologyvariationsiU,i=1,2,…,m,
where the first component measures the change in technical efficiency. The value of efficiency variations for each alternative will be related to one interval that the lowest and highest values of this interval [EfficiencyvariationsiL,EfficiencyvariationsiU] computed as follows:
(16)EfficiencyvariationsiL=Qi(t+1)(t+1),LQi(t)(t),U,i=1,2,…,m,EfficiencyvariationsiU=Qi(t+1)(t+1),UQi(t)(t),L,i=1,2,…,m.
And the second component measures the technology frontier shift between time periods t and t+1. The value of technology variations for each alternative will be related to one interval that the lowest and highest values of this interval [TechnologyvariationsiL,TechnologyvariationsiU] computed as follows:
(17)TechnologyvariationsiL=[Qi(t)(t),UQi(t)(t+1),U·Qi(t+1)(t),LQi(t+1)(t+1),L]1/2,i=1,2,…,m,TechnologyvariationsiU=[Qi(t)(t),LQi(t)(t+1),L·Qi(t+1)(t),UQi(t+1)(t+1),U]1/2,i=1,2,…,m.
Now we have
(18)M.P.IiL=Qi(t+1)(t+1),LQi(t)(t),U·[Qi(t)(t),UQi(t)(t+1),U·Qi(t+1)(t),LQi(t+1)(t+1),L]1/2=[Qi(t+1)(t+1),LQi(t)(t+1),U·Qi(t+1)(t),LQi(t)(t),U]1/2,i=1,2,…,m,M.P.IiU=Qi(t+1)(t+1),UQi(t)(t),L·[Qi(t)(t),LQi(t)(t+1),L·Qi(t+1)(t),UQi(t+1)(t+1),U]1/2=[Qi(t+1)(t+1),UQi(t)(t+1),L·Qi(t+1)(t),UQi(t)(t),L]1/2,i=1,2,…,m.Describing Interval Malmquist Productivity Index.For describing the interval M.P.I, the following relations are considered:(1)
IfM.P.IiL=M.P.IiU=1, then Alternativei did not have progress and regression.(2)
IfM.P.IiL=1,M.P.IiU>1, then Alternativei did have progress.(3)
IfM.P.IiL<1,M.P.IiU=1, then Alternativei did have regression.(4)
IfM.P.IiL<1,M.P.IiU>1, then ρ is applied to determine the progress or regression for Alternativei as shown in
(19)ρ=M.P.IiU-11-M.P.IiL,0<ρ<∞.(i)
Ifρ>1, then the percent of progress is more than the regression.
(ii)
Ifρ<1, then the percent of regression is more than the progress.


## 4. Numerical Example

In this section, we intend to apply the above proposed model to validate and verify our methodology in one illustration example. In this study, we consider ten bank branches of decision maker with respect to their performance in two-time periods (2004-2005). Also, we gather and analyze the efficient indices that can be influence on their efficiency. The exact definition of input and output variables in banking is a disputable issue [35]. There are two main approaches to the choice of how to measure the flow of services provided by banks. The majority of banking studies can be categorized as users of the intermediation model or of the production model. The intermediation approach characterizes banks as financial intermediaries whose function is to collect funds in the form of deposits and other lendable funds and to offer them as loans or other assets that earn income (Figure 3). With this approach, the data is typically assumed to be in the numbers of dollars of loans, deposits, or insurance in force [36].Figure 3
Intermediation model.Under the alternative production approach, banks are the institutions providing fee-based products and services to customers. Products and services such as loans and deposits are outputs in this model, and the resource consumed such as labour, capital, and operating expenses are inputs (Figure4). Under this approach, output is best measured by the number and type of transactions or documents processed over a given time period [35, 37].Figure 4
Production model.These two approaches, production and intermediation approachs have both advantage and disadvantages and cannot fully capture the role of banks. The selection of model, input, and output data is essential because they have a direct influence on the results of empirical analysis. The final choice of model depends upon the concept of what banks do, the stated problem, and the availability of data. Here, the second model, production approach, is chosen to capture the role of banks.


### 4.1. Data

The purpose of this section is to present the empirical variable of banks in Iran and to describe the choice procedure of data for modeling. The data used in this study covers two-time periods (2004-2005). This data is from annual balance sheets and income statements of the banks involved.The Central Bank of Iran has regularly collected the public reports by individual banks. This data was used as a database for the present paper. All variables, with the exception of number of customer and number of personnel, are reported in billions of Rials (Iranian currency unit) and corrected to the 2004 price level using the consumer price index. The example includes ten domestic branches of a commercial bank operating in Iran.To evaluate the productivity performance of bank branches in Iran, the essential element is the selection of variables. The first selection of model variables was on the basis of the research aim, which focuses on productivity of bank branches performance and on the availability of data. Here, we select eight possible groups of data for the productivity model. They are as follows:(i)
loan to customers;(ii)
deposits from customers;(iii)
commissions received;(iv)
incomes;(v)
liabilities/deposits to credit institutions;(vi)
total costs;(vii)
number of personnel;(viii)
number of customers.These data groups, with the exception of personnel and customers, are from the bank branches balance sheets and income statements. For each essay/model six variables were selected. The selection was made on the basis of correlation matrix. The correlation matrix displays the correlation coefficient between the variables.Various goals of bank performance can be presented. For instance, the bank is devoted to increasing employee motivation, while only dedicated employees can create value to bank customers and thus to bank shareholders. For selection of variables it was presumed that the goal of bank performance is to receive income. Therefore, the first step in correlation analysis was to determine the relationship between income and other rations. The data is represented in Table4.Table 4
Correlation matrix for the variables.The strongest correlation appeared between income and the following variables: loan to customers (0.75), deposits from customers (0.81), total costs (0.65), number of personnel (0.72), and number of customers (0.69). Table4 shows remarkably weak correlation between income and the following data set: commission received (0.32) and liabilities/deposits to credit institutions (0.53). Based on weak correlation, these are excluded from the data of the productivity model.Based on correlation analysis results, we have selected the variables for productivity analysis for the present paper. In the first essay the production approach was used, and the variables were defined as follows.(i)
As a benefit criteria,C1 are incomes, C2 are the loans,  C3 are the deposits, and C4 are the numbers of customers.(ii)
As a cost criteria,C5 are the numbers of personnel, and C6 are the total costs.Our decision matrix is founded over data of ten branches of a commercial bank in Iran. Hereby, there were two-time periods between December 31, 2003, and December, 31 2005. Tables8 and 9 contain some information on the variables used. Note that, the values of decision matrix are not precise, and interval numbers are used to describe and treat the uncertainty of the decision problem.In this example, six known variables have relative importance weightswi>0,∑i=16wi=1, that can be calculated by the following steps.Step 1 (Design an Expert Questionnaire).
The fixed weights related to the variables are calculated based on expert opinion. In this way, we design a questionnaire with aim to provide an information profile of the expert which allows his/her diagnosis in determining of relative importance of variables with respect to productivity performance of bank branches. The data which is captured by the questionnaire represents all the information we have about an expert in order to calculate importance weights. The questionnaire consists of 6 questions to which the expert responds on a five-point scale (very low, low, medium, high, and very high). A sample of the filled up proposed questionnaire to determine the importance weights of the variables is given in Table5. The bank branches employees (office employee with diploma, experts, supervisors, middle and top managers, and faculty members) are the members of statistical population of this research. For this purpose 65 questionnaires were distributed and 50 questionnaires were returned with responses.Table 5
A sample of filled up questionnaire.
Question
Very low
Low
Medium
High
Very high
(1) What is the importance of theincomes criterion in a productivity performance of bank?
✓
(2) What is the importance of theloans criterion in a productivity performance of bank?
✓
(3) What is the importance of thedeposits criterion in a productivity performance of bank?
✓
(4) What is the importance of thenumber of customers criterion in a productivity performance of bank?
✓
(5) What is the importance of thenumber of personnel criterion in a productivity performance of bank?
✓
(6) What is the importance of thetotal costs criterion in a productivity performance of bank?
✓Step 2 (Group Integration).
After the polling process, we can obtain the weights of each variable. By using the semantics (Table6), we can transform these into corresponding fuzzy numbers which were introduced by Zadeh [38] to deal with vague, imprecise and uncertain problems and have been used as a modeling tool for complex systems that can be controlled by humans but hard to define precisely.Table 6
Semantically statements conversion table.
Semantic statementsA
Extremely unimportant
Unimportant
Ordinary
Important
Extremely important
Very low
Low
Medium
High
Very high
Fuzzy numbers(α,m,β)
(0.5, 0, 0.2)
(0.15, 0.25, 0.15)
(0.2, 0.5, 0.2)
(0.15, 0.75, 0.15)
(0.15, 1, 0)
Defuzzified fuzzy numbersμT(A)
0.083
0.282
0.5
0.717
0.934In this paper, we use triangular fuzzy numbers for the semantic statements based on Chen and Hwangs’ scaled conversion measurement [39]; see Figure 6 and Table 6. Triangular fuzzy number is a special kind of fuzzy sets. A triangular fuzzy number can be denoted as N=(a,b,c). This representation is interpreted as membership functions and holds the following conditions:(i)
a to b is increasing function;(ii)
b to c is decreasing function;(iii)
a≤b≤c.Figure5 is an illustration of the membership function of a triangular fuzzy number.Figure 5
Triangular fuzzy number.Figure 6
The relationship between semantic statements and fuzzy numbers.The membership function of triangular fuzzy number is(20)μ(x)={x-ab-aifa≤x≤b,c-xc-bifb≤x≤c,0else.Step 3 (Calculate the Variables’ Weights).
First, we defuzzify the fuzzy numbers related to the semantic statements according to (21) and report them as a membership function (μT(A)) in Table 6. After the completed questionnaires have been returned, the frequency analysis of the results for each variable with respect to the five-point scale was done and is reported in Table 7. Afterwards, based on these findings, we calculate the score of each variable according to (22) and report it in Table 7. Finally, by definition of (23), we calculate the variables’ weights:
(21)μT(A)=m+β2(1+β)+m2(1+α),(22)Scorei=∑j=1Snj·μT∑j=1Snj,(23)wi=scorei∑i=1nscorei.
In the third essay, the Malmquist indices of productivity change and their components for the production model are researched. To solve this example using the extended VIKOR method and the proposed Malmquist productivity index, we go through the following steps.(a)
The PIS and NIS in two-time periods are computed by (8) and shown in Tables 10 and 11.(b)
In this step, we compute[SiL,SiU] and [RiL,RiU] for two time periods using (9) and (10). The result is presented in Tables 12, 13, 14, and 15.(c)
We compute theQi=[QiL,QiU] for two-time periods by (11). Now, let us suppose that v=0.5.  The results are shown in Tables 16 and 17.(d)
Finally, with applying models (18), the Malmquist productivity index for ten bank branches is computed. For each branch, two indices of M.P.I under lower and upper bound are obtained and reported in Table 18. For example, the branch B1 is in the worst and best condition having 0.88 and 1.06 Malmquist productivity index, respectively. Therefore, this branch has regression and progress when it is under these conditions and ρ scale in (19) model is computed as follows:
(24)ρ=1.06-11-0.88=0.5.
Also, ρ scale shows that the percent of regression B1 is more than the percent of progress. Hence ρ scale for other bank branches is shown in Table 18 (column 3). Five alternatives have been the percent of progress more than the percent of regression, corresponding to achieved results, the branches which have more in comparison with others are in following order. 												(25)ρB8>ρB10>ρB7>ρB9>ρB3.
Four alternatives (B1, B2, B5, and B6) have been the percent of regression more than the percent of progress, and B4 has been only regression.Table 7
Importance weights related to the six known criteria.
Very lowμT=0.083
LowμT=0.282
MediumμT=0.5
HighμT=0.717
Very highμT=0.934
Score
w
C
1
13
10
10
5
12
0.47
0.1715
≅
0.17
C
2
6
5
10
12
17
0.63
0.229
≅
0.23
C
3
9
10
6
13
12
0.54
0.197
≅
0.2
C
4
18
10
10
7
5
0.38
0.1386
≅
0.14
C
5
22
12
11
3
2
0.29
0.105
≅
0.1
C
6
13
13
10
5
9
0.43
0.1569
≅
0.16Table 8
Interval decision matrix in the first-time period (year of 2004).
C
1
C
2
C
3
C
4
C
5
C
6
B
1
[7041, 7841]
[40003.8, 50003.8]
[3935, 5335]
[740, 840]
[7, 11]
[3352, 4352]
B
2
[1464, 2264]
[31346, 41346]
[2998, 4398]
[321, 421]
[9, 13]
[1425, 2425]
B
3
[4538, 5338]
[55029.56, 62029.56]
[2022, 3422]
[377, 477]
[13, 17]
[713, 1713]
B
4
[116535, 117335]
[515090.2, 52509.02]
[19317, 20717]
[1069, 1169]
[8, 12]
[2167, 3167]
B
5
[8122, 8922]
[36655.05, 46655.05]
[7973, 9373]
[608, 708]
[9, 13]
[1350, 2350]
B
6
[1659, 2459]
[15664.4, 25664.4]
[2189, 3589]
[440, 540]
[7, 11]
[1429, 2429]
B
7
[3562, 4362]
[27000.5, 37000.5]
[1903, 3303]
[306, 406]
[8, 12]
[693, 1693]
B
8
[51887, 52687]
[887196.2, 897196]
[30284, 31684]
[470, 570]
[15, 19]
[4523, 5523]
B
9
[6139, 6939]
[48326.2, 58326.2]
[4234, 5634]
[361, 461]
[3, 7]
[1265, 2265]
B
10
[885, 1685]
[20969.35, 30969.35]
[2965, 4365]
[172, 272]
[5, 9]
[811, 1811]Table 9
Interval decision matrix in the second-time period (year of 2005).
C
1
C
2
C
3
C
4
C
5
C
6
B
1
[7251, 8051]
[44563.8, 54563.8]
[4935, 6335]
[752, 852]
[8, 12]
[2352, 3352]
B
2
[1484, 2284]
[33516, 43516]
[2499, 3899]
[421, 521]
[11, 15]
[1225, 2225]
B
3
[3658, 4458]
[53125.99, 63125.99]
[1302, 2702]
[417, 517]
[11, 15]
[1113, 2113]
B
4
[130793, 131593]
[480171.2, 580171.2]
[19520, 20920]
[1269, 1369]
[9, 13]
[3664, 2664]
B
5
[7622, 8422]
[37934.05, 47934.05]
[8248, 9648]
[628, 738]
[9, 13]
[1223, 2223]
B
6
[1759, 2559]
[16936.4, 26936.4]
[1629, 3029]
[435, 535]
[5, 9]
[823, 1823]
B
7
[3108, 3908]
[31132.3, 41132.3]
[2213, 3613]
[254, 354]
[9, 13]
[1073, 2073]
B
8
[54207, 55007]
[890215.2, 900215]
[29595, 30995]
[532, 632]
[11, 15]
[5677, 4677]
B
9
[6219, 7019]
[58297.2, 68297.2]
[4277, 5677]
[387, 487]
[5, 9]
[1167, 2167]
B
10
[655, 1455]
[23010.35, 33010.35]
[2634, 4034]
[221, 321]
[7, 11]
[1151, 2151]Table 10
PIS and NIS (year of 2004).
C
1
C
2
C
3
C
4
C
5
C
6
A
*
t
117335
897196.2
31684
1169
3
693
A
-
t
885
15664
1903
172
17
5523Table 11
PIS and NIS (year of 2005).
C
1
C
2
C
3
C
4
C
5
C
6
A
*
t
+
1
131993
900215
30995
1369
5
823
A
-
t
+
1
655
16936.4
1302
221
15
5677Table 12
S and R interval numbers (year of 2004-year of 2004).
S
i
(
t
)
(
t
)
,
L
S
i
(
t
)
(
t
)
,
U
R
i
(
t
)
(
t
)
,
L
R
i
(
t
)
(
t
)
,
U
B
1
0.72
0.81
0.22
0.22
B
2
0.75
0.83
0.22
0.23
B
3
0.74
0.83
0.22
0.22
B
4
0.26
0.34
0.1
0.1
B
5
0.66
0.75
0.22
0.23
B
6
0.73
0.81
0.23
0.23
B
7
0.73
0.81
0.22
0.23
B
8
0.32
0.41
0.13
0.16
B
9
0.67
0.76
0.22
0.22
B
10
0.55
0.78
0.23
0.23Table 13
S and R interval numbers (year of 2005-year of 2005).
S
i
(
t
+
1
)
(
t
+
1
)
,
L
S
i
(
t
+
1
)
(
t
+
1
)
,
U
R
i
(
t
+
1
)
(
t
+
1
)
,
L
R
i
(
t
+
1
)
(
t
+
1
)
,
U
B
1
0.69
0.79
0.22
0.22
B
2
0.75
0.85
0.22
0.22
B
3
0.75
0.84
0.22
0.22
B
4
0.25
0.37
0.08
0.11
B
5
0.66
0.75
0.22
0.22
B
6
0.68
0.78
0.22
0.23
B
7
0.75
0.84
0.22
0.23
B
8
0.38
0.47
0.13
0.16
B
9
0.67
0.77
0.22
0.22
B
10
0.73
0.83
0.23
0.23Table 14
S and R interval numbers (year of 2004-year of 2005).
S
i
(
t
)
(
t
+
1
)
,
L
S
i
(
t
)
(
t
+
1
)
,
U
R
i
(
t
)
(
t
+
1
)
,
L
R
i
(
t
)
(
t
+
1
)
,
U
B
1
0.72
0.82
0.22
0.22
B
2
0.75
0.84
0.22
0.23
B
3
0.76
0.85
0.22
0.22
B
4
0.28
0.38
0.1
0.1
B
5
0.67
0.76
0.22
0.22
B
6
0.72
0.82
0.23
0.23
B
7
0.73
0.82
0.22
0.23
B
8
0.33
0.42
0.12
0.15
B
9
0.7
0.75
0.22
0.22
B
10
0.71
0.76
0.23
0.23Table 15
S and R interval numbers (year of 2005-year of 2004).
S
i
(
t
+
1
)
(
t
)
,
L
S
i
(
t
+
1
)
(
t
)
,
U
R
i
(
t
+
1
)
(
t
)
,
L
R
i
(
t
+
1
)
(
t
)
,
U
B
1
0.68
0.77
0.22
0.22
B
2
0.743
0.83
0.22
0.23
B
3
0.74
0.83
0.22
0.22
B
4
0.31
0.39
0.08
0.1
B
5
0.65
0.74
0.22
0.22
B
6
0.69
0.78
0.23
0.23
B
7
0.75
0.83
0.22
0.23
B
8
0.36
0.45
0.13
0.17
B
9
0.68
0.77
0.22
0.22
B
10
0.74
0.83
0.22
0.23Table 16
Q interval numbers.
Q
i
(
t
)
(
t
)
,
L
Q
i
(
t
)
(
t
)
,
U
Q
i
(
t
+
1
)
(
t
+
1
)
,
L
Q
i
(
t
+
1
)
(
t
+
1
)
,
U
B
1
0.87
0.94
0.83
0.92
B
2
0.89
1
0.88
0.92
B
3
0.88
0.96
0.88
0.96
B
4
0
0.07
0
0.2
B
5
0.81
0.93
0.81
0.88
B
6
0.91
0.98
0.83
0.94
B
7
0.87
0.98
0.88
0.99
B
8
0.17
0.36
0.28
0.45
B
9
0.82
0.90
0.82
0.9
B
10
0.75
0.96
0.9
0.98Table 17
Q interval numbers.
Q
i
(
t
)
(
t
+
1
)
,
L
Q
i
(
t
)
(
t
+
1
)
,
U
Q
i
(
t
+
1
)
(
t
)
,
L
Q
i
(
t
+
1
)
(
t
)
,
U
B
1
0.85
0.94
0.82
0.91
B
2
0.87
0.99
0.88
1
B
3
0.88
0.96
0.88
0.97
B
4
0
0.09
0
0.14
B
5
0.80
0.88
0.79
0.88
B
6
0.89
0.97
0.87
0.95
B
7
0.86
0.97
0.89
1
B
8
0.12
0.31
0.21
0.43
B
9
0.83
0.87
0.82
0.91
B
10
0.88
0.92
0.88
1Table 18
Result of M.P.I andρ scale.
M.P.I
L
M.P.I
U
ρ
B
1
0.88
1.06
0.5
B
2
0.88
1.09
0.75
B
3
0.92
1.1
1.25
B
4
0
0
—
B
5
0.88
1.09
0.75
B
6
0.87
1.05
0.38
B
7
0.91
1.15
1.67
B
8
0.73
3.08
7.70
B
9
0.93
1.1
1.43
B
10
0.95
1.22
4.4


## 4.1. Data

The purpose of this section is to present the empirical variable of banks in Iran and to describe the choice procedure of data for modeling. The data used in this study covers two-time periods (2004-2005). This data is from annual balance sheets and income statements of the banks involved.The Central Bank of Iran has regularly collected the public reports by individual banks. This data was used as a database for the present paper. All variables, with the exception of number of customer and number of personnel, are reported in billions of Rials (Iranian currency unit) and corrected to the 2004 price level using the consumer price index. The example includes ten domestic branches of a commercial bank operating in Iran.To evaluate the productivity performance of bank branches in Iran, the essential element is the selection of variables. The first selection of model variables was on the basis of the research aim, which focuses on productivity of bank branches performance and on the availability of data. Here, we select eight possible groups of data for the productivity model. They are as follows:(i)
loan to customers;(ii)
deposits from customers;(iii)
commissions received;(iv)
incomes;(v)
liabilities/deposits to credit institutions;(vi)
total costs;(vii)
number of personnel;(viii)
number of customers.These data groups, with the exception of personnel and customers, are from the bank branches balance sheets and income statements. For each essay/model six variables were selected. The selection was made on the basis of correlation matrix. The correlation matrix displays the correlation coefficient between the variables.Various goals of bank performance can be presented. For instance, the bank is devoted to increasing employee motivation, while only dedicated employees can create value to bank customers and thus to bank shareholders. For selection of variables it was presumed that the goal of bank performance is to receive income. Therefore, the first step in correlation analysis was to determine the relationship between income and other rations. The data is represented in Table4.Table 4
Correlation matrix for the variables.The strongest correlation appeared between income and the following variables: loan to customers (0.75), deposits from customers (0.81), total costs (0.65), number of personnel (0.72), and number of customers (0.69). Table4 shows remarkably weak correlation between income and the following data set: commission received (0.32) and liabilities/deposits to credit institutions (0.53). Based on weak correlation, these are excluded from the data of the productivity model.Based on correlation analysis results, we have selected the variables for productivity analysis for the present paper. In the first essay the production approach was used, and the variables were defined as follows.(i)
As a benefit criteria,C1 are incomes, C2 are the loans,  C3 are the deposits, and C4 are the numbers of customers.(ii)
As a cost criteria,C5 are the numbers of personnel, and C6 are the total costs.Our decision matrix is founded over data of ten branches of a commercial bank in Iran. Hereby, there were two-time periods between December 31, 2003, and December, 31 2005. Tables8 and 9 contain some information on the variables used. Note that, the values of decision matrix are not precise, and interval numbers are used to describe and treat the uncertainty of the decision problem.In this example, six known variables have relative importance weightswi>0,∑i=16wi=1, that can be calculated by the following steps.Step 1 (Design an Expert Questionnaire).
The fixed weights related to the variables are calculated based on expert opinion. In this way, we design a questionnaire with aim to provide an information profile of the expert which allows his/her diagnosis in determining of relative importance of variables with respect to productivity performance of bank branches. The data which is captured by the questionnaire represents all the information we have about an expert in order to calculate importance weights. The questionnaire consists of 6 questions to which the expert responds on a five-point scale (very low, low, medium, high, and very high). A sample of the filled up proposed questionnaire to determine the importance weights of the variables is given in Table5. The bank branches employees (office employee with diploma, experts, supervisors, middle and top managers, and faculty members) are the members of statistical population of this research. For this purpose 65 questionnaires were distributed and 50 questionnaires were returned with responses.Table 5
A sample of filled up questionnaire.
Question
Very low
Low
Medium
High
Very high
(1) What is the importance of theincomes criterion in a productivity performance of bank?
✓
(2) What is the importance of theloans criterion in a productivity performance of bank?
✓
(3) What is the importance of thedeposits criterion in a productivity performance of bank?
✓
(4) What is the importance of thenumber of customers criterion in a productivity performance of bank?
✓
(5) What is the importance of thenumber of personnel criterion in a productivity performance of bank?
✓
(6) What is the importance of thetotal costs criterion in a productivity performance of bank?
✓Step 2 (Group Integration).
After the polling process, we can obtain the weights of each variable. By using the semantics (Table6), we can transform these into corresponding fuzzy numbers which were introduced by Zadeh [38] to deal with vague, imprecise and uncertain problems and have been used as a modeling tool for complex systems that can be controlled by humans but hard to define precisely.Table 6
Semantically statements conversion table.
Semantic statementsA
Extremely unimportant
Unimportant
Ordinary
Important
Extremely important
Very low
Low
Medium
High
Very high
Fuzzy numbers(α,m,β)
(0.5, 0, 0.2)
(0.15, 0.25, 0.15)
(0.2, 0.5, 0.2)
(0.15, 0.75, 0.15)
(0.15, 1, 0)
Defuzzified fuzzy numbersμT(A)
0.083
0.282
0.5
0.717
0.934In this paper, we use triangular fuzzy numbers for the semantic statements based on Chen and Hwangs’ scaled conversion measurement [39]; see Figure 6 and Table 6. Triangular fuzzy number is a special kind of fuzzy sets. A triangular fuzzy number can be denoted as N=(a,b,c). This representation is interpreted as membership functions and holds the following conditions:(i)
a to b is increasing function;(ii)
b to c is decreasing function;(iii)
a≤b≤c.Figure5 is an illustration of the membership function of a triangular fuzzy number.Figure 5
Triangular fuzzy number.Figure 6
The relationship between semantic statements and fuzzy numbers.The membership function of triangular fuzzy number is(20)μ(x)={x-ab-aifa≤x≤b,c-xc-bifb≤x≤c,0else.Step 3 (Calculate the Variables’ Weights).
First, we defuzzify the fuzzy numbers related to the semantic statements according to (21) and report them as a membership function (μT(A)) in Table 6. After the completed questionnaires have been returned, the frequency analysis of the results for each variable with respect to the five-point scale was done and is reported in Table 7. Afterwards, based on these findings, we calculate the score of each variable according to (22) and report it in Table 7. Finally, by definition of (23), we calculate the variables’ weights:
(21)μT(A)=m+β2(1+β)+m2(1+α),(22)Scorei=∑j=1Snj·μT∑j=1Snj,(23)wi=scorei∑i=1nscorei.
In the third essay, the Malmquist indices of productivity change and their components for the production model are researched. To solve this example using the extended VIKOR method and the proposed Malmquist productivity index, we go through the following steps.(a)
The PIS and NIS in two-time periods are computed by (8) and shown in Tables 10 and 11.(b)
In this step, we compute[SiL,SiU] and [RiL,RiU] for two time periods using (9) and (10). The result is presented in Tables 12, 13, 14, and 15.(c)
We compute theQi=[QiL,QiU] for two-time periods by (11). Now, let us suppose that v=0.5.  The results are shown in Tables 16 and 17.(d)
Finally, with applying models (18), the Malmquist productivity index for ten bank branches is computed. For each branch, two indices of M.P.I under lower and upper bound are obtained and reported in Table 18. For example, the branch B1 is in the worst and best condition having 0.88 and 1.06 Malmquist productivity index, respectively. Therefore, this branch has regression and progress when it is under these conditions and ρ scale in (19) model is computed as follows:
(24)ρ=1.06-11-0.88=0.5.
Also, ρ scale shows that the percent of regression B1 is more than the percent of progress. Hence ρ scale for other bank branches is shown in Table 18 (column 3). Five alternatives have been the percent of progress more than the percent of regression, corresponding to achieved results, the branches which have more in comparison with others are in following order. 												(25)ρB8>ρB10>ρB7>ρB9>ρB3.
Four alternatives (B1, B2, B5, and B6) have been the percent of regression more than the percent of progress, and B4 has been only regression.Table 7
Importance weights related to the six known criteria.
Very lowμT=0.083
LowμT=0.282
MediumμT=0.5
HighμT=0.717
Very highμT=0.934
Score
w
C
1
13
10
10
5
12
0.47
0.1715
≅
0.17
C
2
6
5
10
12
17
0.63
0.229
≅
0.23
C
3
9
10
6
13
12
0.54
0.197
≅
0.2
C
4
18
10
10
7
5
0.38
0.1386
≅
0.14
C
5
22
12
11
3
2
0.29
0.105
≅
0.1
C
6
13
13
10
5
9
0.43
0.1569
≅
0.16Table 8
Interval decision matrix in the first-time period (year of 2004).
C
1
C
2
C
3
C
4
C
5
C
6
B
1
[7041, 7841]
[40003.8, 50003.8]
[3935, 5335]
[740, 840]
[7, 11]
[3352, 4352]
B
2
[1464, 2264]
[31346, 41346]
[2998, 4398]
[321, 421]
[9, 13]
[1425, 2425]
B
3
[4538, 5338]
[55029.56, 62029.56]
[2022, 3422]
[377, 477]
[13, 17]
[713, 1713]
B
4
[116535, 117335]
[515090.2, 52509.02]
[19317, 20717]
[1069, 1169]
[8, 12]
[2167, 3167]
B
5
[8122, 8922]
[36655.05, 46655.05]
[7973, 9373]
[608, 708]
[9, 13]
[1350, 2350]
B
6
[1659, 2459]
[15664.4, 25664.4]
[2189, 3589]
[440, 540]
[7, 11]
[1429, 2429]
B
7
[3562, 4362]
[27000.5, 37000.5]
[1903, 3303]
[306, 406]
[8, 12]
[693, 1693]
B
8
[51887, 52687]
[887196.2, 897196]
[30284, 31684]
[470, 570]
[15, 19]
[4523, 5523]
B
9
[6139, 6939]
[48326.2, 58326.2]
[4234, 5634]
[361, 461]
[3, 7]
[1265, 2265]
B
10
[885, 1685]
[20969.35, 30969.35]
[2965, 4365]
[172, 272]
[5, 9]
[811, 1811]Table 9
Interval decision matrix in the second-time period (year of 2005).
C
1
C
2
C
3
C
4
C
5
C
6
B
1
[7251, 8051]
[44563.8, 54563.8]
[4935, 6335]
[752, 852]
[8, 12]
[2352, 3352]
B
2
[1484, 2284]
[33516, 43516]
[2499, 3899]
[421, 521]
[11, 15]
[1225, 2225]
B
3
[3658, 4458]
[53125.99, 63125.99]
[1302, 2702]
[417, 517]
[11, 15]
[1113, 2113]
B
4
[130793, 131593]
[480171.2, 580171.2]
[19520, 20920]
[1269, 1369]
[9, 13]
[3664, 2664]
B
5
[7622, 8422]
[37934.05, 47934.05]
[8248, 9648]
[628, 738]
[9, 13]
[1223, 2223]
B
6
[1759, 2559]
[16936.4, 26936.4]
[1629, 3029]
[435, 535]
[5, 9]
[823, 1823]
B
7
[3108, 3908]
[31132.3, 41132.3]
[2213, 3613]
[254, 354]
[9, 13]
[1073, 2073]
B
8
[54207, 55007]
[890215.2, 900215]
[29595, 30995]
[532, 632]
[11, 15]
[5677, 4677]
B
9
[6219, 7019]
[58297.2, 68297.2]
[4277, 5677]
[387, 487]
[5, 9]
[1167, 2167]
B
10
[655, 1455]
[23010.35, 33010.35]
[2634, 4034]
[221, 321]
[7, 11]
[1151, 2151]Table 10
PIS and NIS (year of 2004).
C
1
C
2
C
3
C
4
C
5
C
6
A
*
t
117335
897196.2
31684
1169
3
693
A
-
t
885
15664
1903
172
17
5523Table 11
PIS and NIS (year of 2005).
C
1
C
2
C
3
C
4
C
5
C
6
A
*
t
+
1
131993
900215
30995
1369
5
823
A
-
t
+
1
655
16936.4
1302
221
15
5677Table 12
S and R interval numbers (year of 2004-year of 2004).
S
i
(
t
)
(
t
)
,
L
S
i
(
t
)
(
t
)
,
U
R
i
(
t
)
(
t
)
,
L
R
i
(
t
)
(
t
)
,
U
B
1
0.72
0.81
0.22
0.22
B
2
0.75
0.83
0.22
0.23
B
3
0.74
0.83
0.22
0.22
B
4
0.26
0.34
0.1
0.1
B
5
0.66
0.75
0.22
0.23
B
6
0.73
0.81
0.23
0.23
B
7
0.73
0.81
0.22
0.23
B
8
0.32
0.41
0.13
0.16
B
9
0.67
0.76
0.22
0.22
B
10
0.55
0.78
0.23
0.23Table 13
S and R interval numbers (year of 2005-year of 2005).
S
i
(
t
+
1
)
(
t
+
1
)
,
L
S
i
(
t
+
1
)
(
t
+
1
)
,
U
R
i
(
t
+
1
)
(
t
+
1
)
,
L
R
i
(
t
+
1
)
(
t
+
1
)
,
U
B
1
0.69
0.79
0.22
0.22
B
2
0.75
0.85
0.22
0.22
B
3
0.75
0.84
0.22
0.22
B
4
0.25
0.37
0.08
0.11
B
5
0.66
0.75
0.22
0.22
B
6
0.68
0.78
0.22
0.23
B
7
0.75
0.84
0.22
0.23
B
8
0.38
0.47
0.13
0.16
B
9
0.67
0.77
0.22
0.22
B
10
0.73
0.83
0.23
0.23Table 14
S and R interval numbers (year of 2004-year of 2005).
S
i
(
t
)
(
t
+
1
)
,
L
S
i
(
t
)
(
t
+
1
)
,
U
R
i
(
t
)
(
t
+
1
)
,
L
R
i
(
t
)
(
t
+
1
)
,
U
B
1
0.72
0.82
0.22
0.22
B
2
0.75
0.84
0.22
0.23
B
3
0.76
0.85
0.22
0.22
B
4
0.28
0.38
0.1
0.1
B
5
0.67
0.76
0.22
0.22
B
6
0.72
0.82
0.23
0.23
B
7
0.73
0.82
0.22
0.23
B
8
0.33
0.42
0.12
0.15
B
9
0.7
0.75
0.22
0.22
B
10
0.71
0.76
0.23
0.23Table 15
S and R interval numbers (year of 2005-year of 2004).
S
i
(
t
+
1
)
(
t
)
,
L
S
i
(
t
+
1
)
(
t
)
,
U
R
i
(
t
+
1
)
(
t
)
,
L
R
i
(
t
+
1
)
(
t
)
,
U
B
1
0.68
0.77
0.22
0.22
B
2
0.743
0.83
0.22
0.23
B
3
0.74
0.83
0.22
0.22
B
4
0.31
0.39
0.08
0.1
B
5
0.65
0.74
0.22
0.22
B
6
0.69
0.78
0.23
0.23
B
7
0.75
0.83
0.22
0.23
B
8
0.36
0.45
0.13
0.17
B
9
0.68
0.77
0.22
0.22
B
10
0.74
0.83
0.22
0.23Table 16
Q interval numbers.
Q
i
(
t
)
(
t
)
,
L
Q
i
(
t
)
(
t
)
,
U
Q
i
(
t
+
1
)
(
t
+
1
)
,
L
Q
i
(
t
+
1
)
(
t
+
1
)
,
U
B
1
0.87
0.94
0.83
0.92
B
2
0.89
1
0.88
0.92
B
3
0.88
0.96
0.88
0.96
B
4
0
0.07
0
0.2
B
5
0.81
0.93
0.81
0.88
B
6
0.91
0.98
0.83
0.94
B
7
0.87
0.98
0.88
0.99
B
8
0.17
0.36
0.28
0.45
B
9
0.82
0.90
0.82
0.9
B
10
0.75
0.96
0.9
0.98Table 17
Q interval numbers.
Q
i
(
t
)
(
t
+
1
)
,
L
Q
i
(
t
)
(
t
+
1
)
,
U
Q
i
(
t
+
1
)
(
t
)
,
L
Q
i
(
t
+
1
)
(
t
)
,
U
B
1
0.85
0.94
0.82
0.91
B
2
0.87
0.99
0.88
1
B
3
0.88
0.96
0.88
0.97
B
4
0
0.09
0
0.14
B
5
0.80
0.88
0.79
0.88
B
6
0.89
0.97
0.87
0.95
B
7
0.86
0.97
0.89
1
B
8
0.12
0.31
0.21
0.43
B
9
0.83
0.87
0.82
0.91
B
10
0.88
0.92
0.88
1Table 18
Result of M.P.I andρ scale.
M.P.I
L
M.P.I
U
ρ
B
1
0.88
1.06
0.5
B
2
0.88
1.09
0.75
B
3
0.92
1.1
1.25
B
4
0
0
—
B
5
0.88
1.09
0.75
B
6
0.87
1.05
0.38
B
7
0.91
1.15
1.67
B
8
0.73
3.08
7.70
B
9
0.93
1.1
1.43
B
10
0.95
1.22
4.4


## 5. Conclusion

Because of the fact that determining the exact values of the attributes is difficult or impossible, it was more appropriate to consider them as interval numbers. In this paper, we extended VIKOR method to calculate Malmquist productivity index for MADM problem with interval numbers. This method introduced the multicriteria ranking index based on the particular measure of “closeness” to the ideal solution and calculated the separation of each alternative from the ideal solution at two-time periods.Here, productivity of alternatives available in decision matrix with interval numbers and their improvement or deterioration was researched.  Productivity rate of growth for alternatives was calculated by applying proposed Malmquist productivity index. Finally, effectiveness of our methodology was illustrated by a numerical example in this paper.



---

*Source: 824316-2013-10-01.xml*