# An Application of Monte-Carlo-Based Sensitivity Analysis on the Overlap in Discriminant Analysis

Authors: S. Razmyan; F. Hosseinzadeh Lotfi
Journal: Journal of Applied Mathematics
 (2012)

---

## Abstract

Discriminant analysis (DA) is used for the measurement of estimates of a discriminant function by minimizing their group misclassifications to predict group membership of newly sampled data. A major source of misclassification in DA is due to the overlapping of groups. The uncertainty in the input variables and model parameters needs to be properly characterized in decision making. This study combines DEA-DA with a sensitivity analysis approach to an assessment of the influence of banks’ variables on the overall variance in overlap in a DA in order to determine which variables are most significant. A Monte-Carlo-based sensitivity analysis is considered for computing the set of first-order sensitivity indices of the variables to estimate the contribution of each uncertain variable. The results show that the uncertainties in the loans granted and different deposit variables are more significant than uncertainties in other banks’ variables in decision making.


---

## Body

## 1. Introduction

The classification problem of assigning observations to one of different groups plays an important role in decision making. When observations are restricted to one of two groups, the Binary classification has wide applicability in business environments.Discriminant analysis (DA) is a classification method that can distinguish the group membership of a new observation. A group of observations for which the memberships have already been identified is used for the estimation of a discriminant function by some criteria, such as the minimization of misclassification. A new sample is classified into one of the groups based on the gained results [1].Mangasarian [2] identified that linear programming (LP) could be used to determine separating hyperplanes, namely, when two set of observations are linearly separable use linear discriminant function. Freed and Glover [3] and Hand [4] using objectives such as minimization of the sum of deviations (MSD) or maximization of the minimum deviation (MMD) of misclassified observations from the separating hyperplane, when sets of observations that are not necessarily linearly separable, proposed LP methods for generating linear discriminant function.Then a model is based on the goal programming (GP) extension of LP by choosing different criteria, such as minimizing the maximum deviation, maximizing the minimum deviation, minimizing the sum of interior deviation, minimizing the sum of deviations, minimizing misclassified observations, minimizing external deviations, maximizing internal deviations, maximizing the ratio of internal to external, and hybrid models for which there are both advantages and deficiencies [5–9].In DA, LP and other mathematical programming (MP) based approaches are non-parametric and more flexible than statistical methods [6, 10]. Retzlaff-Roberts [11, 12] and Tofallis [13] proposed the use of DEA-ratio model for DA. Sueyoshi [14] using a data envelopment analysis (DEA) additive model, described a goal programming formulation of DA in which the proposed model is more directly linked to minimizing the sum of deviations from the separating hyperplane; this method was named DEA-DA to distinguish it from other DA and DEA approaches. The original GP version of DEA-DA could not deal with negative data. Therefore, Sueyoshi [15] extended DEA-DA to overcome this deficiency. This approach was designed to minimize the total distance of misclassified observations and formulated by two-stage GP formulations. The number of misclassifications can, however, be considered as measures of misclassification, in which binary variables indicate whether observations are correctly or in correctly classified. Bajgier and Hill [16] proposed a Mixed Integer Programming (MIP) model that included the number of misclassifications in the objective function for the two-group discriminant problem. Gehrlein [17] and Wilson [18] introduced a MIP approach for minimizing the number of misclassified observations in multigroup problems. Chang and Kuo [19] proposed a procedure based on benchmarking model of DEA to solve the two-group problems. Sueyoshi [20] reformulated DEA-DA by MIP to minimize the total number of misclassified observations. When an overlap between two groups is not a serious problem, dropping the first stage of the two-stage MIP approach simplifies the estimation process [21].Sensitivity analysis provides an understanding of how the model outputs is affected by changes in the inputs, therefore it can assist to increase the confidence in the model and its predictions. sensitivity analysis can use in deciding whether inputs estimates are sufficiently precise to give reliable predictions or we can find the model parameters that can be eliminated.Two classes in sensitivity analysis have been distinguished [22].Local Sensitivity Analysis
Studies how some small variations of inputs around a given value, change the value of the output. This approach is practical when the variation around a baseline of the input-output relationship to be assumed linear.Global Sensitivity Analysis
Takes into account all the variation range of the inputs, and has for aim to apportion output uncertainty to inputs’ ones. It quantified the output uncertainty due to the uncertainty in the input parameters. Global sensitivity analysis apportions the output uncertainty to the uncertainty in the input factors, described typically by probability distribution functions that cover the factors’ ranges of existence.
Local methods are less helpful when sensitivity analysis is used to compare the effect of various factors on the output, as in this case the relative uncertainty of each input should be weighted. A global sensitivity analysis technique thus incorporates the influence of the whole range of variation and the form of the probability density function of the input. The variance-based methods can be considered as a quantitative method for global sensitivity analysis. In this Study, the Sobol’ decomposition in the framework of Monte Carlo simulations (MCS), [22], which is from the family of quantitative methods for global sensitivity analysis, is applied to study of the effect of the variability in DA due to the uncertainty in the variables. The results of the sensitivity analysis can determine which of the variables have a more dominant influence on the uncertainty in the model output.This paper is organized as follows: Section2 briefly introduces the DEA-DA model; Section 3 describes the sensitivity analysis based on a Monte-Carlo simulation. Section 4 contains an example and the conclusion is provided in Section 5.


## 2. Data Envelopment Analysis-Discriminant Analysis (DEA-DA)

The two-stage MIP approach [20] is used in this study to describe DEA-DA. We considered two groups (G1 and G2) for which the sum of the two groups has n observations (j=1,…,n). Each observation has m independent factors (i=1,…,m), denoted by rij. It is necessary to identify the group membership of each observation before its computation. In the two stage approach, the computation process consists of classification and overlap identification, and handling overlap. The first stage is formulated as follows [20].Stage 1.
(2.1)
Minimize
s
,
Subject  to
∑
i
=
1
m
(
λ
i
+
-
λ
i
-
)
r
i
j
-
d
+
s
≥
0
,
j
∈
G
1
,
subject  to
∑
i
=
1
m
(
λ
i
+
-
λ
i
-
)
r
i
j
-
d
-
s
≤
-
ε
,
j
∈
G
2
,
subject  to
∑
i
=
1
m
(
λ
i
+
+
λ
i
-
)
=
1
,
subject  to
ξ
i
+
≥
λ
i
+
≥
ε
ξ
i
+
,
ξ
i
-
≥
λ
i
-
≥
ε
ξ
i
-
,
i
=
1
,
…
,
m
,
subject  to
ξ
i
+
+
ξ
i
-
≤
1
,
i
=
1
,
…
,
l
,
subject  to
∑
i
=
1
m
(
ξ
i
+
+
ξ
i
-
)
=
l
,
subject  to
d
,
s
:
unrestricted
,
ξ
i
+
,
ξ
i
-
∈
{
0,1
}
and
all
other
variables
≥
0
.
Here “d” indicates a discriminant score for group classification and “s” indicates the size of an overlap between two groups.
Letλi* (=λi+*-λi-*) and d* and s* be an optimal solution of the model (2.1). Then, the original data set (G) is classified into the following subsets G=G1∪G2=C1∪N1∪C2∪N2, where
(2.2)C1={j∈G1|∑i=1mλi*rij>d*+s*},C2={j∈G2|∑i=1mλi*rij<d*-s*},N1=G1-C1,N2=G2-C2.
Then, we determine that observations inC1 belong to G1 and the observations of C2 belong to G2 because their location is identified from model (2.1). The two subsets N1 and N2 consist of the observations have not yet been classified in the first stage.Stage 2.
Ifs*≥0, then the existence of an overlap is identified in the fist stage. In this stage, we reclassify all of the observations belonging to the overlap (N1∪N2) because the group membership of these observations is still undetermined. The second stage is reformulated as follows [20]:
(2.3)Minimize∑j∈D1lj+w∑j∈D2lj,Subject  to∑i=1m(λi+-λi-)rij-c+Mlj≥0,j∈N1,subject  to∑i=1m(λi+-λi-)rij-c-Mlj≤-ε,j∈N2,subject  to∑i=1m(λi++λi-)=1,subject  toξi+≥λi+≥εξi+,ξi-≥λi-≥εξi-,i=1,…,m,subject  toξi++ξi-≤1,i=1,…,m,subject  toλi++λi-≥ε,for  somei,subject  toc:unrestricted,ξi+,ξi-,yj∈{0,1},and  all  other  variables≥0.
Here, the binary variable (lj) counts the number of observations classified incorrectly. The objective function minimizes the number of such misclassifications. The weight (w) identifies the importance between G1 and G2 in terms of the number of observations. In the presented model (2.3), it is necessary to prescribe a large number (M) and a small number (ε). The equation λi++λi-≥ε(forsomei) indicates that some pairs avoid the occurrence of λi+=0 and λi-=0.
After gaining an optimal solution onλi* and c*, the second stage classifies observations in the overlap as follows: if ∑i=1mλi*rij≥c*, then the jth observation belongs to G1, or if ∑i=1mλi*rij≤c*-ε, then it belongs to G2. Thus, all of the observations in G are classified into G1 or G2 at the end of the second stage.


## 3. Sensitivity Analysis Based on Monte-Carlo Simulation (MCS)

Sensitivity analysis was created to deal simply with uncertainties in the input variables and model parameters [22]. The results of an sensitivity analysis can determine which of the input parameters have a more dominant influence on the uncertainty in the model output [23]. A variance-based sensitivity analysis, which addresses the inverse problem of attributing the output variance to uncertainty in the input, quantifies the contribution that each input factor makes to the variance in the output quantity of interest. A global sensitivity analysis of complex numerical models can be performed by calculating variance-based importance measures of the input variables, such as the Sobol’ indices. These indices are calculated by evaluating a multidimensional integral using a Monte-Carlo technique. This approach allows analyzing the influence of different variables and their subsets, the structure of f(x), and so forth.It is assumed that a mathematical model havingn input parameters gathered in an input vector with a joint probability density function (pdf) can be presented as a model function:
(3.1)s=f(r),r=(r1,r2,…,rm)∈Un≡[0;1]n,thatrij∈Un≡[0;1]ni=1,…,m;j=1,…,n,
where ri=(r1i,r2i,…,rni)T. Because of the variables are affected by several kinds of heterogeneous uncertainties that reflect the imperfect knowledge of the system, it is assumed that input variables are independent and that the probability density function p(r) is known, even if the rij are not actually random variables.The Sobol’ sensitivity method explores the multidimensional space of the unknown input parametersx with a certain number of MC samples. The sensitivity indices are generated by a decomposition of the model function f in an n-dimensional factor space into summands of increasing dimensionality [22]:
(3.2)f(r)=f(r1,…,rm)=f0+∑i=1m∑q1<⋯<qifq1,…,qi(rq1,rq2,…,rqi).
where the constant f0 is the mean value of the function, and the integral of each summand over any of its independent variables is zero. Due to this property, the summands are orthogonal to each other in the following form:
(3.3)∫Unfi1,…,itfj1,…,jqdr=0,(i1,…,it)≠(j1,…,jq).The sensitivity index,SIi, represents the fractional contribution of a given factor ri to the variance in a given output variable, s. To calculate the sensitivity indices, the total variance, Vs, in the model output, s, is apportioned to all of the input factors, r1,…, rm, as follows:
(3.4)Vs=∫Unf2(r)dr-f02.By integrating the square of (2.2) and with (2.3), it is possible to decompose the total variance (3.1) as follows [24]:
(3.5)Vs=∑iVi+∑i<hVih+∑i<h<kVk+⋯+V12⋯o,
where Vi=V[E(s|ri=ri*)],
(3.6)Vih=V[E(s|ri=ri*,rh=rh*)]-Vi-Vh,
and so on. V[E(s∣ri=ri*)] is referred to as the variance of the conditional expectation and is the variance over all of the values of ri* in the expectation of s given that ri has a fixed value  ri*. This is an intuitive measure of the sensitivity of s to a factor ri, as it measures the amount by which E(s∣ri=ri*) varies with the value of ri* whilst averaging over the rh’s, h≠i. Following the above definition for the partial variances, the sensitivity indices are defined as
(3.7)SIi=ViVsHigher order indices can be calculated with a similar approach. With regard to (3.2), the decomposition of the sensitivity indices can be written in the following form:
(3.8)∑i=1mSIi+∑1≤i≤h≤mSIih+⋯+SI1,2,…,m=1.The Sobol’ indices are usually computed with a MC simulation. The mean value and total and partial variance can be derived withNsim samples in the following [22]:
(3.9)f^0=1Nsim∑p=1Nsimf(rp),V^s=1Nsim∑p=1Nsimf2(rp)-f^02,V^i=1Nsim∑p=1Nsimf2(r(~i)p(1),rip(1))f(r(~i)p(2),rip(1))-f^02.
In the later equations, rm is a sampled variable in Un, and
(3.10)r(~i)p=(r1p,r2p,…,r(-1)p,r(i+1)p,…,rmp).
The superscripts (1) and (2) indicate that two different samples are generated and mixed.


## 4. Illustrative Examples

Classification methods are widely used in economic and finance. They are useful for classifying the sectors based upon their performance in different groups and predict the group memberships of new firms. Most of researchers used classification methods to classifying the firms based upon their performance assessment. DA is a classification method that is used in this study. The purpose of the first stage in DEA-DA is to determine whether there is an overlap between the two groups. The existence of an overlap is the main source of misclassification in DA. By identification of the overlap between two groups, it is possible to increase the number of observations classified correctly. If there is no overlap, any DA method may produce an almost perfect classification. However, if there is an overlap, an additional computation process is needed to deal with such an overlap [20]. So, there is a tradeoff between computational effort/time and a high level of classification capability.Misclassification can result as a consequence of an intersection between two groups. Many researchers have proposed approaches that try to reveal the advantage of identifying the minimized overlap of two groups for risk management on the classification problem [19, 20, 25, 26]. Given the importance of the banking sector for, in general, the whole economy and, in particular, for the financial system, in this section, we present an application of the sensitivity analysis to overlaps, s, on data from a commercial bank of Iran. This assertion is illustrated numerically for bank branches that have more than 20 and 30 personnel, in two different examples.If we wish to take into account the inherent randomness with respect to what the criteria might experience, we have to bring stochastic characterization into play. The stochastic efficiency assessment of banking branches normally requires performing a set of analyses on DMUs with a suit of variables as criteria.At first, we use the additive model to discriminate banking branches. Most models need to examine both a DEA efficiency score and slacks or an efficiency score measured by them, which depends upon input-based or output-based measurement. The additive model [27] aggregates input-oriented and output-oriented measures to produce the efficiency score. Consequently, the efficiency status is more easily determined by the additive model than by the radial model.In two different examples, the real data set consists of 78 and 18 banking branches. This study selects 31 and 8 branches as inefficient branches, and 47 and 10 branches as efficient branches, respectively, in example 1 and 2, documented in Tables1 and 2. For determining the classifications based on the additive model, three variables of personnel, payable interest, and non-performing loans are considered as inputs. Nine variables of loans granted, long-term deposit, current deposit, non-benefit deposit, short-term deposit, received interest, and received fee are assumed as outputs.Table 1
Personnel
Loans granted
Long-term deposits
Current deposits
Non-benefit deposit
Short-term deposits
Payable interest
Received interest
Received fee
Non-performing loans
Efficient branches
20
311
E
+
08
90
E
+
08
62
E
+
08
85
E
+
08
71
E
+
08
8
E
+
08
48
E
+
07
117
E
+
06
386
E
+
07
39
737
E
+
08
135
E
+
08
228
E
+
08
74
E
+
08
225
E
+
08
14
E
+
08
234
E
+
07
152
E
+
06
9
E
+
07
26
377
E
+
08
100
E
+
08
126
E
+
08
28
E
+
08
154
E
+
08
11
E
+
08
158
E
+
07
52
E
+
06
44
E
+
07
20
507
E
+
08
172
E
+
08
153
E
+
08
42
E
+
08
135
E
+
08
16
E
+
08
221
E
+
07
59
E
+
06
110
E
+
07
34
538
E
+
08
105
E
+
08
173
E
+
08
434
E
+
08
158
E
+
08
10
E
+
08
22
E
+
07
161
E
+
06
1226
E
+
07
27
1383
E
+
08
246
E
+
08
508
E
+
08
92
E
+
08
190
E
+
08
23
E
+
08
59
E
+
07
145
E
+
06
221
E
+
07
28
1740
E
+
08
590
E
+
08
253
E
+
08
86
E
+
08
242
E
+
08
47
E
+
08
445
E
+
07
294
E
+
06
1075
E
+
07
28
960
E
+
08
100
E
+
08
987
E
+
08
120
E
+
08
109
E
+
08
95
E
+
08
122
E
+
07
515
E
+
06
460
E
+
07
29
1036
E
+
08
225
E
+
08
148
E
+
08
118
E
+
08
199
E
+
08
18
E
+
08
49
E
+
07
186
E
+
06
895
E
+
07
27
1068
E
+
08
97
E
+
08
1066
E
+
08
68
E
+
08
125
E
+
08
10
E
+
08
123
E
+
07
245
E
+
06
3121
E
+
07
27
1593
E
+
08
109
E
+
08
1977
E
+
08
92
E
+
08
103
E
+
08
11
E
+
08
42
E
+
07
301
E
+
06
1135
E
+
07
20
734
E
+
08
366
E
+
08
142
E
+
08
203
E
+
08
350
E
+
08
33
E
+
08
40
E
+
07
176
E
+
06
537
E
+
07
24
572
E
+
08
189
E
+
08
239
E
+
08
87
E
+
08
268
E
+
08
23
E
+
08
197
E
+
07
628
E
+
06
1530
E
+
07
22
454
E
+
08
168
E
+
08
34
E
+
08
64
E
+
08
143
E
+
08
17
E
+
08
565
E
+
07
102
E
+
06
5130
E
+
07
27
350
E
+
08
144
E
+
08
609
E
+
08
95
E
+
08
236
E
+
08
17
E
+
08
30
E
+
07
423
E
+
06
140
E
+
07
20
268
E
+
08
158
E
+
08
111
E
+
08
49
E
+
08
163
E
+
08
15
E
+
08
16
E
+
07
120
E
+
06
100
E
+
07
21
624
E
+
08
78
E
+
08
125
E
+
08
98
E
+
08
164
E
+
08
97
E
+
08
21
E
+
07
193
E
+
06
362
E
+
07
21
288
E
+
08
19
E
+
08
1367
E
+
08
29
E
+
08
499
E
+
08
27
E
+
08
13
E
+
07
265
E
+
06
31
E
+
07
20
86
E
+
08
274
E
+
08
314
E
+
08
32
E
+
08
210
E
+
08
24
E
+
08
31
E
+
07
36
E
+
06
26
E
+
07
25
126
E
+
08
235
E
+
08
163
E
+
08
104
E
+
08
477
E
+
08
27
E
+
08
14
E
+
07
75
E
+
06
37
E
+
07
20
137
E
+
08
162
E
+
08
169
E
+
08
27
E
+
08
208
E
+
08
16
E
+
08
14
E
+
07
116
E
+
06
409
E
+
07
25
224
E
+
08
356
E
+
08
491
E
+
08
107
E
+
08
309
E
+
08
32
E
+
08
1
E
+
07
160
E
+
06
59
E
+
07
20
369
E
+
08
110
E
+
08
131
E
+
08
49
E
+
08
190
E
+
08
13
E
+
08
24
E
+
07
81
E
+
06
306
E
+
07
21
358
E
+
08
241
E
+
08
109
E
+
08
129
E
+
08
279
E
+
08
22
E
+
08
23
E
+
07
94
E
+
06
45
E
+
07
24
108
E
+
08
696
E
+
08
46
E
+
08
38
E
+
08
389
E
+
08
60
E
+
08
12
E
+
07
63
E
+
06
5
E
+
07
42
339
E
+
08
645
E
+
08
157
E
+
08
53
E
+
08
238
E
+
08
42
E
+
08
28
E
+
07
48
E
+
06
667
E
+
07
23
199
E
+
08
307
E
+
08
205
E
+
08
31
E
+
08
155
E
+
08
26
E
+
08
17
E
+
07
162
E
+
06
10
E
+
07
37
914
E
+
08
397
E
+
08
241
E
+
08
61
E
+
08
298
E
+
08
34
E
+
08
116
E
+
07
242
E
+
06
4950
E
+
07
33
1685
E
+
08
153
E
+
08
117
E
+
08
166
E
+
08
255
E
+
08
17
E
+
08
17
E
+
07
40
E
+
06
200
E
+
07
21
634
E
+
08
173
E
+
08
125
E
+
08
120
E
+
08
287
E
+
08
18
E
+
08
44
E
+
07
83
E
+
06
270
E
+
07
34
288
E
+
08
490
E
+
08
86
E
+
08
121
E
+
08
549
E
+
08
47
E
+
08
3
E
+
07
151
E
+
06
368
E
+
07
20
235
E
+
08
112
E
+
08
159
E
+
08
99
E
+
08
489
E
+
08
9
E
+
08
12
E
+
07
42
E
+
06
622
E
+
07
24
486
E
+
08
101
E
+
08
233
E
+
08
208
E
+
08
988
E
+
08
9
E
+
08
15
E
+
07
107
E
+
06
363
E
+
07
20
1361
E
+
08
1194
E
+
08
354
E
+
08
78
E
+
08
346
E
+
08
94
E
+
08
82
E
+
07
242
E
+
06
69
E
+
07
23
712
E
+
08
369
E
+
08
102
E
+
08
43
E
+
08
106
E
+
08
29
E
+
08
43
E
+
07
81
E
+
06
190
E
+
07
28
1574
E
+
08
249
E
+
08
204
E
+
08
86
E
+
08
132
E
+
08
21
E
+
08
64
E
+
07
192
E
+
06
1963
E
+
07
36
2163
E
+
08
323
E
+
08
781
E
+
08
132
E
+
08
184
E
+
08
27
E
+
08
394
E
+
07
323
E
+
06
995
E
+
07
21
352
E
+
08
231
E
+
08
108
E
+
08
93
E
+
08
107
E
+
08
17
E
+
08
20
E
+
07
64
E
+
06
261
E
+
07
22
255
E
+
08
558
E
+
08
188
E
+
08
46
E
+
08
178
E
+
08
41
E
+
08
48
E
+
07
68
E
+
06
551
E
+
07
20
191
E
+
08
171
E
+
08
90
E
+
08
88
E
+
08
159
E
+
08
13
E
+
08
27
E
+
07
104
E
+
06
250
E
+
07
27
1373
E
+
08
1002
E
+
08
393
E
+
08
45
E
+
08
597
E
+
08
76
E
+
08
269
E
+
07
482
E
+
06
5347
E
+
07
23
542
E
+
08
538
E
+
08
144
E
+
08
15
E
+
08
238
E
+
08
39
E
+
08
2
E
+
07
113
E
+
06
918
E
+
07
21
1454
E
+
08
431
E
+
08
336
E
+
08
34
E
+
08
563
E
+
08
41
E
+
08
312
E
+
07
124
E
+
06
20
E
+
07
24
1141
E
+
08
723
E
+
08
367
E
+
08
28
E
+
08
558
E
+
08
56
E
+
08
32
E
+
07
193
E
+
06
81
E
+
07
21
1796
E
+
08
1188
E
+
08
416
E
+
08
38
E
+
08
801
E
+
08
93
E
+
08
442
E
+
07
726
E
+
06
103
E
+
07
26
361
E
+
08
223
E
+
08
236
E
+
08
75
E
+
08
441
E
+
08
21
E
+
08
8
E
+
07
58
E
+
06
201
E
+
07
24
1367
E
+
08
730
E
+
08
129
E
+
08
28
E
+
08
105
E
+
08
55
E
+
08
213
E
+
07
538
E
+
06
496
E
+
07
Inefficient branches
33
1088
E
+
08
509
E
+
08
283
E
+
08
50
E
+
08
183
E
+
08
49
E
+
08
172
E
+
07
95
E
+
06
2198
E
+
07
39
744
E
+
08
282
E
+
08
248
E
+
08
341
E
+
08
282
E
+
08
25
E
+
08
69
E
+
07
188
E
+
06
1373
E
+
07
31
1723
E
+
08
565
E
+
08
313
E
+
08
55
E
+
08
481
E
+
08
61
E
+
08
11
E
+
07
578
E
+
06
523
E
+
07
29
1020
E
+
08
196
E
+
08
255
E
+
08
120
E
+
08
171
E
+
08
18
E
+
08
15
E
+
07
164
E
+
06
713
E
+
07
21
417
E
+
08
194
E
+
08
90
E
+
08
82
E
+
08
180
E
+
08
17
E
+
08
18
E
+
07
129
E
+
06
773
E
+
07
29
757
E
+
08
151
E
+
08
218
E
+
08
288
E
+
08
204
E
+
08
14
E
+
08
34
E
+
07
101
E
+
06
1599
E
+
07
34
461
E
+
08
319
E
+
08
183
E
+
08
64
E
+
08
309
E
+
08
29
E
+
08
65
E
+
07
199
E
+
06
144
E
+
07
25
461
E
+
08
120
E
+
08
184
E
+
08
70
E
+
08
168
E
+
08
12
E
+
08
4
E
+
07
114
E
+
06
699
E
+
07
29
1519
E
+
08
521
E
+
08
382
E
+
08
80
E
+
08
215
E
+
08
44
E
+
08
61
E
+
07
236
E
+
06
1960
E
+
07
22
226
E
+
08
245
E
+
08
821
E
+
08
90
E
+
08
188
E
+
08
21
E
+
08
5
E
+
07
73
E
+
06
172
E
+
07
24
769
E
+
08
258
E
+
08
72
E
+
08
28
E
+
08
182
E
+
08
22
E
+
08
6
E
+
07
56
E
+
06
446
E
+
07
24
629
E
+
08
310
E
+
08
289
E
+
08
28
E
+
08
439
E
+
08
33
E
+
08
85
E
+
07
291
E
+
06
217
E
+
07
30
67
E
+
08
505
E
+
08
201
E
+
08
51
E
+
08
464
E
+
08
49
E
+
08
53
E
+
07
184
E
+
06
94
E
+
07
21
1182
E
+
08
428
E
+
08
298
E
+
08
63
E
+
08
452
E
+
08
48
E
+
08
132
E
+
07
66
E
+
06
197
E
+
07
32
331
E
+
08
331
E
+
08
122
E
+
08
67
E
+
08
378
E
+
08
33
E
+
08
30
E
+
07
105
E
+
06
424
E
+
07
36
376
E
+
08
306
E
+
08
102
E
+
08
126
E
+
08
202
E
+
08
26
E
+
08
18
E
+
07
138
E
+
06
377
E
+
07
25
275
E
+
08
355
E
+
08
333
E
+
08
115
E
+
08
475
E
+
08
40
E
+
08
62
E
+
07
150
E
+
06
3097
E
+
07
24
367
E
+
08
218
E
+
08
66
E
+
08
117
E
+
08
227
E
+
08
21
E
+
08
59
E
+
07
49
E
+
06
59
E
+
07
31
257
E
+
08
147
E
+
08
380
E
+
08
74
E
+
08
237
E
+
08
16
E
+
08
22
E
+
07
548
E
+
06
182
E
+
07
20
170
E
+
08
220
E
+
08
74
E
+
08
30
E
+
08
198
E
+
08
21
E
+
08
7
E
+
07
63
E
+
06
329
E
+
07
27
452
E
+
08
156
E
+
08
323
E
+
08
62
E
+
08
314
E
+
08
17
E
+
08
23
E
+
07
90
E
+
06
768
E
+
07
27
406
E
+
08
227
E
+
08
176
E
+
08
110
E
+
08
251
E
+
08
22
E
+
08
22
E
+
07
119
E
+
06
237
E
+
07
26
334
E
+
08
166
E
+
08
78
E
+
08
80
E
+
08
146
E
+
08
14
E
+
08
10
E
+
07
49
E
+
06
246
E
+
07
34
381
E
+
08
218
E
+
08
469
E
+
08
129
E
+
08
228
E
+
08
21
E
+
08
10
E
+
07
250
E
+
06
1167
E
+
07
28
565
E
+
08
133
E
+
08
186
E
+
08
225
E
+
08
156
E
+
08
12
E
+
08
11
E
+
07
168
E
+
06
530
E
+
07
40
675
E
+
08
382
E
+
08
121
E
+
08
52
E
+
08
109
E
+
08
28
E
+
08
208
E
+
07
106
E
+
06
1694
E
+
07
20
313
E
+
08
467
E
+
08
98
E
+
08
31
E
+
08
199
E
+
08
40
E
+
08
73
E
+
07
30
E
+
06
133
E
+
07
24
121
E
+
08
849
E
+
08
154
E
+
08
23
E
+
08
455
E
+
08
73
E
+
08
31
E
+
07
24
E
+
06
682
E
+
07
25
327
E
+
08
109
E
+
08
237
E
+
08
46
E
+
08
196
E
+
08
13
E
+
08
44
E
+
07
251
E
+
06
547
E
+
07
24
1001
E
+
08
672
E
+
08
444
E
+
08
81
E
+
08
419
E
+
08
56
E
+
08
133
E
+
07
585
E
+
06
1007
E
+
07
34
482
E
+
08
192
E
+
08
173
E
+
08
53
E
+
08
199
E
+
08
17
E
+
08
94
E
+
07
87
E
+
06
252
E
+
07Table 2
Personnel
Loans granted
Long-term deposits
Current deposits
Non-benefit deposit
Short-term deposits
Payable interest
Received interest
Received fee
Non-performing loans
Efficient branches
34
538
E
+
08
105
E
+
08
173
E
+
08
434
E
+
08
158
E
+
08
10
E
+
08
22
E
+
07
161
E
+
06
1226
E
+
07
42
339
E
+
08
645
E
+
08
157
E
+
08
53
E
+
08
238
E
+
08
42
E
+
08
28
E
+
07
48
E
+
06
667
E
+
07
37
914
E
+
08
397
E
+
08
241
E
+
08
61
E
+
08
298
E
+
08
34
E
+
08
116
E
+
07
242
E
+
06
4950
E
+
07
34
288
E
+
08
490
E
+
08
86
E
+
08
121
E
+
08
549
E
+
08
47
E
+
08
3
E
+
07
151
E
+
06
368
E
+
07
36
2163
E
+
08
323
E
+
08
781
E
+
08
132
E
+
08
184
E
+
08
27
E
+
08
394
E
+
07
323
E
+
06
995
E
+
07
39
744
E
+
08
282
E
+
08
248
E
+
08
341
E
+
08
282
E
+
08
25
E
+
08
69
E
+
07
188
E
+
06
1373
E
+
07
31
1723
E
+
08
565
E
+
08
313
E
+
08
55
E
+
08
481
E
+
08
61
E
+
08
11
E
+
07
578
E
+
06
523
E
+
07
30
67
E
+
08
505
E
+
08
201
E
+
08
51
E
+
08
464
E
+
08
49
E
+
08
53
E
+
07
184
E
+
06
94
E
+
07
31
257
E
+
08
147
E
+
08
380
E
+
08
74
E
+
08
237
E
+
08
16
E
+
08
22
E
+
07
548
E
+
06
182
E
+
07
34
381
E
+
08
218
E
+
08
469
E
+
08
129
E
+
08
228
E
+
08
21
E
+
08
10
E
+
07
250
E
+
06
1167
E
+
07
Inefficient branches
33
1088
E
+
08
509
E
+
08
283
E
+
08
50
E
+
08
183
E
+
08
49
E
+
08
172
E
+
07
95
E
+
06
2198
E
+
07
34
461
E
+
08
319
E
+
08
183
E
+
08
64
E
+
08
309
E
+
08
29
E
+
08
65
E
+
07
199
E
+
06
144
E
+
07
32
331
E
+
08
331
E
+
08
122
E
+
08
67
E
+
08
378
E
+
08
33
E
+
08
30
E
+
07
105
E
+
06
424
E
+
07
36
376
E
+
08
306
E
+
08
102
E
+
08
126
E
+
08
202
E
+
08
26
E
+
08
18
E
+
07
138
E
+
06
377
E
+
07
40
675
E
+
08
382
E
+
08
121
E
+
08
52
E
+
08
109
E
+
08
28
E
+
08
208
E
+
07
106
E
+
06
1694
E
+
07
34
482
E
+
08
192
E
+
08
173
E
+
08
53
E
+
08
199
E
+
08
17
E
+
08
94
E
+
07
87
E
+
06
252
E
+
07
39
737
E
+
08
135
E
+
08
228
E
+
08
74
E
+
08
225
E
+
08
14
E
+
08
234
E
+
07
152
E
+
06
9
E
+
07
33
1685
E
+
08
153
E
+
08
117
E
+
08
166
E
+
08
255
E
+
08
17
E
+
08
17
E
+
07
40
E
+
06
200
E
+
07Then, for sensitivity analysis in DEA-DA, each observation is modeled as random parameters as follows [28]:
(4.1)rij=μrij(1+Vrijαl),
where μrij and Vrij are the mean value and the coefficients of variation (COV) of the random parameters, respectively, and αl is generated random parameter with a zero mean that is used in MC simulation. The determination of the bank branches’ parameters carries a high degree of uncertainty, and the specification of these parameters can involve a significant degree of expert judgment. Additionally, the COV of these variables plays an important role in the variation of the efficiency. Here, the COV of the all of the parameters are assumed to equal 0.05.To compute the sensitivity indices, the Sobol’ sampling scheme has been used. Sequences of Sobol’ sampling vectors are essentially quasi-random sequences that are defined as sequences of points that have no intrinsic random properties. In this study, a sensitivity analysis is applied to assess the influence of the banks’ variables on the overlap.After one hundred estimates with a sample size of 5000, convergence was seen in the first-order Sobol’ indices derived by a Sobol’ sampling of the uniform criteria spaces for different banks’ variables. The sensitivity indices,Si, are depicted in Figures 1 and 2. These figures present the comparison of the first-order indices of the banks’ variables. The total fraction of the variance captured by the first order functions is approximately 99%. This indicates that, for this problem, higher order contributions to the Sobol’ series are relatively small. The overall variance in banks’ efficiency is affected by the variances in each of the random variables. Figure 1 indicates that 54% and 33% of the overall variance in overlap in DEA-DA is attributable to the variance in loans granted and different deposits, respectively, while the personnel, received interest, fee, and non-performing loans variables have little effect. Also, Figure 2 indicates that uncertainties in loans granted and different deposits are the dominant variables in overlap in DEA-DA.Figure 1
The relative contributions of different variables on the overlap in DEA-DA in example 1 that the real data set consists of 78 banking branches.Figure 2
The relative contributions of different variables on the overlap in DEA-DA in example 1 that the real data set consists of 18 banking branches.


## 5. Conclusions

Due to the inherent complexity and randomness of the data in DEA and problems involving unpredictable or stochastic variables, a probabilistic analysis may be the most rational method of analysis. Therefore, in a probabilistic-based approach, the results open the door to understanding the appropriate estimation of the deciding variables in DEA. In the overlap in DA, the analytical results show that the loans granted and different deposit variables are the main sources of uncertainty, while other variables have a relatively small effect. The main advantage of the used sensitivity analysis approach is that it provides quantified evaluation of the influence of individual variables in DA, and results may be used for decision making.



---

*Source: 315868-2012-11-10.xml*