Accepted Manuscript
A non-radial directional distance method on classifying inputs and outputs in DEA: Application to banking industry Mehdi Toloo , Maryam Allahyar , Jana Hanˇclova PII: DOI: Reference:
S0957-4174(17)30638-3 10.1016/j.eswa.2017.09.034 ESWA 11555
To appear in:
Expert Systems With Applications
Received date: Revised date: Accepted date:
23 December 2016 14 May 2017 12 September 2017
Please cite this article as: Mehdi Toloo , Maryam Allahyar , Jana Hanˇclova , A non-radial directional distance method on classifying inputs and outputs in DEA: Application to banking industry, Expert Systems With Applications (2017), doi: 10.1016/j.eswa.2017.09.034
This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.
ACCEPTED MANUSCRIPT
Highlights A non-radial nor-oriented method is developed to deal with flexible measures.
Two optimistic and pessimistic approaches are proposed.
Each approach contains two individual and integrated models.
A case study of 61 banks in the Visegrad Four region validates the new models.
AC
CE
PT
ED
M
AN US
CR IP T
I
ACCEPTED MANUSCRIPT
A non-radial directional distance method on classifying inputs and outputs in DEA: Application to banking industry
Mehdi Toloo1
CR IP T
Department of Systems Engineering, Technical University of Ostrava, Sokolska tř. 33, 702 00 Ostrava 1, Ostrava, Czech Republic E-mail:
[email protected]
URL: http://homel.vsb.cz/~tol0013/
AN US
Tel: (+420) 792 272 272
Maryam Allahyar
Young Researchers and Elite Club, Yadegar-e-Imam Khomeini (RAH)Shahre Rey Branch, Islamic Azad University, Tehran, Iran
M
Email:
[email protected]
Jana Hančlova
ED
Department of Systems Engineering, Technical University of Ostrava, Sokolska tř. 33, 702 00 Ostrava 1, Ostrava, Czech Republic
AC
CE
PT
Email:
[email protected]
1
Corresponding Author
II
ACCEPTED MANUSCRIPT
Abstract The original Data Envelopment Analysis (DEA) models have required an assumption that the
CR IP T
status of all inputs and outputs be known exactly, whilst we may face a case with some flexible performance measures whose status is unknown. Some classifier approaches have been proposed in order to deal with flexible measures. This contribution develops a new classifier non-radial directional distance method with the aim of taking into account input contraction and output expansion, simultaneously, in the presence of flexible measures. To make the most appropriate decision for flexible measures, we suggest two pessimistic and optimistic
AN US
approaches from both individual and summative points of view. Finally, a numerical real example in the banking system in the countries of the Visegrad Four (i.e. Czech Republic, Hungary, Poland, and Slovakia) is presented to elaborate applicability of the proposed method. Keywords: Data Envelopment Analysis; Directional distance function; Non-radial non-oriented models;
1
M
Mixed integer linear programming; Flexible measure.
Introduction
ED
Data envelopment analysis (DEA) is a non-parametric approach based on mathematical programming for evaluating the relative performances of many different kinds of entities, commonly referred to as decision making units (DMUs), where the presence of multiple inputs
PT
and multiple outputs makes comparisons difficult. DEA was originally introduced by Charnes et al. (1978) (CCR model), built on the ideas of Farrell (1957), and extended by Banker et al. (1984) (BCC model). Nowadays, DEA is becoming a very important analysis tool and research
CE
method in various sciences such as management science, operational research, system engineering, decision analysis, etc. Input- and output-oriented models are two main approaches
AC
in DEA, which respectively measure the largest radial contraction of inputs and the largest radial expansion of outputs. Actually, the conventional DEA models developed with the assumption that the status of each measure is clearly stated as an input or an output variable. However, in the real world, certain variables, referred to as flexible measures, can play input role for some DMUs and output role for others. For example, “high-value customers” and “deposit” measures in a bank branch, “research income” measure in higher education application (Beasley, 1990), or “uptime” measure in evaluating robotics installations (Cook, Johnston and Mccutcheon, 1992) are considered as a flexible measure which could serve as either an input or an output. Bala and Cook (2003) presented a two-step procedure measurement tool for 1
ACCEPTED MANUSCRIPT
evaluating the performance of branches in the banking industry with flexible measures. In the first stage, all variables are assumed to be flexible and a discriminant model is used to designate the input/output status; the second stage performs the DEA analysis based on the variable designations chosen. Cook and Zhu (2007) introduced a method based on a fractional programming problem to determine whether a measure is an input or an output. However, Toloo (2009) claimed that their model may produce incorrect efficiency scores as result of
CR IP T
introducing a large positive number to the model and introduced a revised mixed integer linear programming (MILP) model with the aim of excluding the large positive number from the model. Amirteimoori and Emrouznejad (2011), Toloo (2012, 2014) declared that one drawback in this method is the requirement to enter extra information to decide about the role of each variable. On the other hand, in the presence of alternative optimal solutions in classifier models, the results of selecting a flexible measure as an input or an output are the same for some DMUs and it is logical that not be taken into account for classifying inputs and outputs. Accordingly,
AN US
Toloo (2012) referred to these cases as share cases and introduced a new classifier model that identifies share cases. In another study, Amirteimoori et al. (2013) developed a flexible slacksbased model to calculate the relative efficiency of DMUs in the presence of flexible measures. In order to allow the contraction in inputs and expansion in outputs, simultaneously, Chambers et al. (1996, 1998) introduced the directional distance function. In this approach, there is no need
M
to distinguish between input-oriented and output-oriented approaches. Indeed, the traditional input- and output-oriented models are special cases of this concept. However, unlike the distance function is
ED
traditional DEA models, the relative efficiency score computed using the directional technology for efficient DMUs and ranges between
and
for inefficient DMUs.
Among the classic DEA models, the non-radial models have a higher discriminating power in
PT
evaluating the efficiencies of DMUs and seem to be more effective in performance assessment. Considering this feature, the current paper, inspired by the general non-radial directional distance model (DD model), proposes a new non-radial non-oriented classifier method in both
CE
envelopment and multiplier form in the presence of flexible measures. The suggested classifier directional distance model treats each flexible factor as input, output or both to evaluate individual DMUs from two pessimistic and optimistic viewpoints. Practically, in some cases,
AC
determining the status of a flexible measure may tie and in order to break the tie, the adopted models are developed from the summative perspective as well. Moreover, the proposed approach will be applied to evaluate banks in the countries of the Visegrad Four (V4) in 2014. V4 is an alliance of four central European states – Czech Republic (CZ), Hungary (HU), Poland (PL), and Slovakia (SK) – for the purposes of furthering their European integration, as well as for advancing military, economic and energy cooperation with one another. The efficiency of the financial system is playing a significant role from the 2
ACCEPTED MANUSCRIPT
viewpoints of macroeconomic and microeconomic. The performance of the banking sector from a macroeconomic perspective was influenced by the cost of financial intermediation and the stability of the entire financial system. From a microeconomic point of view, the evaluation of bank operations is important because foreign bank entry to a domestic financial market increases the competitiveness and is necessary to improve an institutional regulation. Evaluating the bank efficiency in the literature distinguishes between two basic approaches, the
CR IP T
production and intermediation approaches that differ in defining the role of deposits as input or output in a production system of banks. This paper takes the deposits measure into consideration as a flexible measure.
The outline of the paper is organized as follows. In Section 2, we review two classifier methods based on classic DEA models. A new DEA-based classifier approach is introduced in both envelopment and multiplier forms in the presence of flexible measures in Section 3. The penultimate section illustrates the discrimination power of suggested models by a real case future research directions in Section 5.
1
Preliminaries
AN US
involving 61 bank branches in V4 countries. Finally, we have come to conclusion and some
This section contains two subsections. Firstly, we briefly review two classifier methods based on classic DEA models deal with flexible measures. Secondly, we provide some preliminary
M
concept about DD model that inspires us in developing the suggested method in the rest of the
1.1
Flexible measures
Assume that there are (
ED
paper.
DMUs
) to produce
each using
semi-positive outputs
semi-positive inputs . The following pair of
PT
multiplier and envelopment models under constant returns to scale (CRS) assumption can be , respectively:
CE
applied to assess the radial input-efficiency of ∑
AC
∑ ∑
(1)
∑
∑ ∑
(2)
3
ACCEPTED MANUSCRIPT
where decision variables
and
are the weight of
is the intensity vector corresponds to
input and
output, respectively, and
. These two models are mutual dual (for more
details see Cooper et al. 2007). Moreover, assume that there are
semi-positive flexible measures
(
) whose
input/output statuses are unknown. To handle such flexible measures, Cook & Zhu (2007) and Toloo (2012) established the following MILP models by modifying models (1) and (2),
∑ ∑ ∑
CR IP T
respectively: ∑ ∑
∑ ∑
∑
∑
AN US
(3)
(4)
ED
M
∑ ∑ ∑ ∑
is a large positive number and the indicator variable
measure
PT
where
is an input (
) or an output (
designates that flexible
). The output-oriented version of above
models can easily be formulated (see Toloo (2012)). It should be noted that models (3) and (4)
CE
are not mutual dual. Toloo (2012) also showed that the existence of the alternative optimal solutions of classifier models for some DMUs could cause incorrect results and referred to these
AC
units as share cases, which must not be taken into account in overall classifying inputs and outputs.
1.2 Directional distance function Chambers et al. (1996, 1998) with the aim of generalizing of the existing distance functions and introducing a non-oriented approach suggested the following generic DD model, which assess under CRS technology:
4
ACCEPTED MANUSCRIPT
∑ ∑
(5)
where the nonzero directional vector (
)
enables the model to
CR IP T
contract inputs and expand outputs simultaneously. Since data are semi-positive, a usual choice for the directional vector is the observed input and output levels. As a result, taking the direction vector (
)
into account, a well-known special case of model (5) can be
formulated as bellow:
AN US
∑ ∑
(6)
From now on, “*” indicates the optimal values. It can be readily demonstrated that however unlike the traditional DEA models, any positive number between
. Furthermore,
can be
for inefficient DMUs. In other words, the value
is the
. Note that model (6) is an extension of traditional DEA models: The
M
inefficiency score of
and
is efficient if
;
input-oriented CCR model (2) can be obtained if ( )
ED
output-oriented CCR model achieves when (
)
and analogously the .
Model (6) assumes that inputs are contracted and outputs are expanded at the same rate (proportionally) which it indeed signifies a radial measure of efficiency. However, in this study,
PT
we utilize the following non-radial version of model (6) which enables us to address nonproportionate improvement in inputs and outputs:
CE
∑
∑
∑ ∑
AC
(7)
In the formulation above the variables
and
are introduced to allow the non-proportionate
contraction of inputs and the non-proportionate expansion of outputs in the given directional vector
, respectively.
As inspection makes clear, we have score of
input and
and
. Indeed
and
are the inefficiency
output, respectively. Hence, the optimal solution of model (7) for 5
ACCEPTED MANUSCRIPT
tells us the amount of employed excess inputs shortfall outputs
and the amount of produced
. This model actually maximizes the sum of all input/output
inefficiencies. Here,
is efficient if
; otherwise it is inefficient in at least one of the components of
input or output vector.
∑ ∑
∑ ∑
CR IP T
Consider the following dual model of model (7):
(8)
Since this model is non-oriented naturally, the weighted sum of both inputs and outputs are ∑
AN US
included in the objective function which minimizes ∑
.
Models (7) and (8) are in fact envelopment and multiplier forms of non-radial DD models, respectively.
2
Proposed non-radial classifier directional distance approach
M
The main issue of this study is to facilitate the derivation of the input/output status of flexible measures. To pursue this objective, this section proposes a new classifier method based on non-
ED
radial DD model which combines both the input and output orientations. As a matter of fact, we adopt such model for reasons like the following: (i) Traditional classifier models are oriented and hence the role of a flexible measure in an
PT
input-oriented model may differ from an output-oriented model. However, DD models are non-oriented and hence extending classifier DD models accommodate the flexible measures by combining both input- and output-orientations.
CE
(ii) The special structure of formulations of envelopment and multiplier forms of DD models helps us to drive two classifier models and to decide the status of flexible measures from two different perspectives, i.e. pessimistic and optimistic.
AC
(iii) The main advantage of non-radial models is that they measure individual input reductions along with individual output expansions in order to provide information on the inefficiency score of each input and output.
In particular, this section develops a pair of non- radial classifier directional distance (CDD hereafter) models in both envelopment and multiplier forms, which introduce two pessimistic and optimistic approaches considering the individual and summative viewpoints. Individual models designate flexible measures in favor of each individual DMUs, meanwhile, summative 6
ACCEPTED MANUSCRIPT
models aim at accommodating flexible measures in favor of all DMUs, simultaneously. In other words, individual models are solved
times, one for each unit, while summative models are
applied just once to determine the overall role of flexible measures. We also provide two pessimistic and optimistic approaches for each individual and summative model. The former and latter approaches adjust the flexible measures status in line with maximizing and minimizing the possible distance from frontier in the given direction, respectively.
CR IP T
2.1 Pessimistic CDD approach
In order to develop a non-radial envelopment DEA model for determining the input/output status of flexible measures, we consider model (7) and the variable measure . There would be two possible cases. If ∑
corresponding to flexible
treats as an input role, then the constraint
must be active; otherwise the constraint ∑
must
hold. A method to model the either-or constraints is introducing the indicator variable designates that
is an input, and
designates it as an output. As a result, to
AN US
where
,
handle the either-or constraints, it is sufficient to impose the following constraints to model (7): ∑ ∑
If
is a large positive number. , then the constraint ∑
is active and the constraint ∑
M
where
is weakened to ∑ binding. In this case,
treats as an input. We can analogously conclude that the situation is . Hence, we formulate the following non-radial non-oriented MILP model,
ED
reversed when
, which is always non-
which classifies flexible measures in the envelopment form: ∑
∑
)
(9)
AC
CE
∑ ∑ ∑ ∑
(∑
PT
̅
It should be mentioned here that the direction vector in this model is where
and operator
element. In general, there are
multiples vectors
and
element by
various combinations for the direction vectors and the
following theorem proves that model (9) finds one direction with the maximum objective function value. Let
be the optimal objective value obtained by model (7) for
each combination. 7
ACCEPTED MANUSCRIPT
Theorem 1. ̅
{
}
Proof. By contradiction, suppose that ̅ Consider the
combination and let
{
where
}.
be the index of flexible measures, which treat as input.
Under these assumptions, model (7) can be rewritten as follows: ∑
)
∑ ∑ ∑ ∑
CR IP T
∑
(∑
(10)
is at hand. Clearly,
AN US
Suppose that model (10) is solved and the optimal solution ( (
is a feasible solution of model (9) where {
Corollary 1. If ̅
M
However, the objective function value of this feasible solution, i.e. optimal objective function value ̅ , which is a contradiction. ■
, is larger than the
then an identical inefficiency score is achieved when the flexible measure is efficient by model (9) then it
ED
is considered either as input or output. In other words, if
is also efficient by models (7) and (8), no matter flexible measure is designated as an input or an output.
PT
In fact, Theorem 1 proves that model (9) obtains the optimal direction with the maximum sum of input/output inefficiencies among all
possible combinations.
CE
Now is time to decide on overall input/output status of flexible measures for each individual
AC
DMU. As the simplest method, one criterion to make a decision would be based on the majority choice among the DMUs. For , let ̅ and be the optimal objective value and the optimal values of the indicator variable corresponding to flexible measure l in model (9), respectively. The following criterion is taken into account, mathematically, to handle the flexible measure.
For flexible measure l, assign whole DMUs to the following two groups: { |̅ { |̅
} }
{ |̅
} 8
ACCEPTED MANUSCRIPT
Set
includes the index of DMUs that are efficient in model (9). Referencing to Corollary 1,
these DMUs are efficient either the flexible measure plays the role of an input or an output. So, logically, these DMUs must not be considered for classifying inputs and outputs. Flexible measure l should be designated as an output if | | Moreover, ties for determining the status of flexible measure
| | and as an input if | | happen when | |
| |.
| |. In this
case, an alternative criterion is needed to examine the issue from the point of view of the may be
CR IP T
collection of all DMUs. Hence, ties for determining the status of flexible measure
broken by applying the following summative version of model (9) which carries out the simultaneous optimization the performance measure of all units and evaluates the pessimistic summative inefficiency corresponding each unit: ∑
∑
(∑
∑
)
AN US
∑ ∑ (
)
∑
(
)
where the intensity variables
(
.
Note
),
(
that
and
the
PT
.
),
(
) for
will end up either as an input, if
CE
Theorem 2. Let corresponding to
solution (
),
, and in a single stage. As a , or an output, if
be the part of optimal solution ∑ ∑ ) ̅ .
(∑
in model (11) then
̅
Proof: By contradiction, we suppose that
AC
are corresponding to optimal
for summative model (11) accommodates flexible measure matter of fact, the flexible measure
(11)
where
ED
for
and variables
M
∑
for
is a feasible solution to model (9) when
. It is easy to show that is evaluated. However, the
objective function value of this feasible solution is larger than the optimal objective function value ̅ , which is impossible. ■ 2.2 Optimistic CDD approach Referencing to Theorem 1, the formulated non-radial envelopment CDD model obtains the maximum possible inefficiency among all possible combinations for each DMU. Alternatively, 9
ACCEPTED MANUSCRIPT
we formulate non-radial multiplier CDD models, which give the minimum possible inefficiency for each DMU with the aim of increasing the discrimination power, which is needed for discriminating among efficient DMUs. To do this, we establish the following mixed integer nonlinear programming model extending model (8): ∑
(∑ ∑
)
∑
CR IP T
∑
∑
(12)
Where the decision variable
, the following method can be applied to eliminate the product of a
AN US
linear due to term
is the weight of flexible measure . Although model (12) is non‐
binary and a continuous variable: Let
be a indicator variable, and
be a continuous variable which
continuous variable is introduced to replace the product added to the model for forcing ̅ to take the value of
. Now a non-zero
. The following constraints must be
:
ED
M
̅ ̅ ̅ ̅ The validity of these constraints can be checked by examining all following two possible situations: If ̅
If
, then from the constraint ̅ and ̅
PT
we have ̅
are redundant.
, then the constraint ̅
or equivalently ̅
; subsequently the constraint ̅
CE
lead to ̅
We, therefore, utilize this method and let ̅
10
, and also ̅
is redundant. for
AC
(12):
. In this case, other constraints
in order to linearize model
ACCEPTED MANUSCRIPT
∑
∑
(∑
∑
∑
∑
∑
̅
)
̅
(13)
In this model, if
, then
CR IP T
̅ ̅ ̅
treats as an input and otherwise it treats as an output. The
following theorem proves that model (13) finds a direction with the minimum objective value possible combinations. Let
obtained by model (8) for the
} {
Proof. Suppose the
combination.
{
Theorem 3.
combination and let
be the optimal objective value
AN US
among all
} and contrary to our claim
. Consider
be the index set of flexible measures which treat as input. Under
∑
(∑
∑
∑
Let (
)
∑
ED
∑
∑
(14)
PT
∑
M
these assumptions, model (8) can be written as bellow:
be the optimal solution of model (14). Clearly, (
̅ is a feasible solution
AC
CE
of model (13) where
̅
{
However, the objective function value of the feasible solution objective function value
is smaller than the optimal
, which is a contradiction. ■
The following theorem proves that the optimal objective value for the proposed envelopment CDD model (9) is greater than or equal to the optimal objective value for the proposed multiplier CDD model (13): Theorem 4.
̅ . 11
ACCEPTED MANUSCRIPT
Proof: According to Theorems 1 and 2, we have ̅ combination
and
for
. Since, for each combination, models (7) and (8) are mutually dual, it ̅ . which implies
then follows that
Theorem 4 lays emphasis on the fact that the proposed non-radial CDD models (9) and (13) are formulated with pessimistic and optimistic points of view, respectively, and hence if a unit is
Corollary 2. If ̅
, then
CR IP T
efficient in the envelopment model (9), then it is also efficient in the multiplier model (13). .
It should be noted that the criterion based on the majority choice can be carried out here in order to designate the role (input or output) for each flexible measure by the optimal solution obtained in model (13). In order to break the possible ties, likewise to the aforementioned
AN US
summative envelopment model (11), the following envelopment form is suggested with the aim of evaluating the optimistic summative inefficiency corresponding each unit: ∑ ∑
∑
̅
)
̅
(15)
PT
ED
̅ ̅ ̅
,
and
represent the
CE
where
∑
M
∑
∑
(∑
input,
output, and
flexible measure weight for
. In fact, model (15) is an extended version of model (13) where the
AC
performance measure of all the DMUs is simultaneously considered. The optimal solution flexible measure (
for summative model (15) accommodates
in a single stage where
) for
, and
(
)
. The flexible measure
(
)
is an input, if
and otherwise it is an output. Theorem 4. Let ∑
be the optimal solution and ∑
̅
corresponding to
12
in model (15) then
∑ .
ACCEPTED MANUSCRIPT
Proof: By contradiction, we suppose that
for
is a feasible solution to model (10) when
. It is easy to show that is evaluated. However, the
objective function value of this feasible solution is smaller than the optimal objective function value
3
, which is impossible. ■
Application
CR IP T
To demonstrate the applications of proposed CDD models in both optimistic and pessimistic approaches (i.e. envelopment and multiplier forms as well as individual and summative models), we analyze the performance of banking industry in V4 countries. The economies of these countries have in common that have seen some changes after communism's collapse; the transition to a market economy, joining the European Union in May 2004 and especially the transformation of banking system. The transition from centrally planned economy to market economy had been accompanied with restructuring and liberalization of the banking system. It
AN US
had been associated with the privatization of some banks, the entry of foreign-owned banks, deregulation of interest rates and changes in legislation.
Exploring Banking Efficiency is important from both macroeconomic and microeconomic standpoints. From the macroeconomic point of view, the efficiency of the banking sector affects the cost of financial intermediation and the stability of the entire financial system. From the
M
microeconomic point of view, banking efficiency is particularly important for improving institutional regulation and supervision, and in particular for improving the competitiveness of banks. Increasing the efficiency of banks leads to a better distribution of financial resources, and
ED
thus better investment support and economic growth. Analysis of the banking performance considering a collection of inputs and outputs, which are
PT
clearly specified, discovers some problems in the system. However, the most controversial issue is the role of deposits, whether they should be classified as inputs or inputs. Boda and Zimkova (2015) investigate the bank efficiency in the Slovak banking industry using three the best-known
CE
approaches, the production, the intermediation, and the value-added ones. The production approach was founded by Benston (1965) and is based on assumption that the aim of banks is to
AC
produce deposits (liabilities) as well as loans (assets) and other services to customers. This service-oriented approach considers deposits as an output together with loans and the interest income. Hancock (1991) modified the production approach into the user-cost approach, where deposits are specified as both inputs and outputs of the cost/profit function of a bank. The second approach is the intermediation approach and was proposed by Sealey and Lindley (1977). This approach assumes that the main aim of a bank is to produce the intermediation services through the collection of deposits or other liabilities. The main banks are seen as production units, which transmute deposits into loans. They are also interested in the use of interest-earning assets and loans, securities and similar investments. The profit-oriented 13
ACCEPTED MANUSCRIPT
approach is a modification of the intermediation approach as the profit-maximizing tendency of the banks. The last value-added approach differs from previous ones in considering all liability and asset categories to have some output characteristics rather than distinguishing inputs from outputs in a mutually exclusive way. The major categories of produced deposits and loans are viewed as important outputs because they form a significant proportion of value added. The current study adopts the intermediation approach for evaluating bank efficiency in the
CR IP T
commercial banks of V4 countries recorded in 2014. The economic environment has not been dramatically changed since 2014. The most V4 economies are slowly recovering from the global crisis and optimism among investors and consumers is growing. A positive effect on the banking sector was limited for a number of reasons. The effects of regulation are not only direct costs of implementation of regulation but also the changing business environment requesting of actual new information using data mining from disposal data at banks with help of
AN US
technologies innovations.
There are 61 banks including 15 banks from CZ, 14 banks from HU, 23 banks from PL, and 9 banks from SK whose complete name and the country where they are located are provided in Table A.1 (see Appendix) in detail. We employed three groups of variables: three inputs (labor, physical capital, and credit risk ratio), three outputs (loans and advances, investments, and non-
M
interest income), and one flexible measure (deposit), which are summarized in Table 1. =====[Please insert Table 1 around here]======
ED
The annual balance sheet and income statement data have been collected from the world banking information source Bankscope2. The data for bank branches are listed in Table 2. Mean values of inputs, outputs, and deposits along with the number of banks corresponding each
PT
country are summarized in Table A.2. in Appendix. =====[Please insert Table 2 around here]======
CE
In order to make a comparison between input- and output-oriented classifier models of Cook & Zhu (2007) and Toloo (2012), we apply these models to the data set. Table 3 reports the results.
AC
As can be seen, different assignations for the flexible measure are obtained through input- and output-oriented models of Cook & Zhu (2007) for banks
. To cope with this issue,
we now apply the proposed non-radial non-oriented CDD models by which the flexible measures are accommodated without the need for distinguishing between input- and outputorientation. 2
Bankscope combines comprehensive financial statements with a wide range of other banking intelligence including ratings, an analysis model, bank structures, news, Anti-Money-Laundering (AML) documentation and banking research. The website has information on 32000 banks and is the definitive tool for bank research and analysis. We refer the readers for more details visit https://bankscope.bvdinfo.com/.
14
ACCEPTED MANUSCRIPT
=====[Please insert Table 3 around here]====== Table 4 provides the obtained result from both suggested individual models (9) and (13). To be more specific, the calculated input/output inefficiencies through model (9) are also reported in Table 5. =====[Please insert Tables 4 and 5 around here]====== of Table 4 show the value of
as the sum of all
CR IP T
The second and third columns
input/output inefficiencies for each DMU when deposit , i.e. , is considered as an input and an output, respectively (the columns labeled “
(input)” and “
(output)”). The optimal objective
function value of proposed individual models (9) and (13) along with the related optimal shown in the last four columns. We note that equalities ̅ =max and =min
AN US
assured by Theorem 1 and Theorem 3, respectively.
are are
As can be extracted from the obtained results, models (9) and (13) identify 12 and 17 efficient banks, respectively. Moreover, 12 efficient DMUs are identified by both models which follows from Corollary 2. Essentially, more efficient units introduced by model (13) compared to model (9) which reveals the optimistic/pessimistic viewpoint of these models. For these 12 banks (six financial institutions from the Czech Republic and six banks from Poland), i.e.
M
, we have
that means these
banks are efficient no matter the deposit measure is designated as an input or an output. Therefore, these banks are not taken into account in the overall decision on the role of deposit as
ED
input or output. Large financial institutions such as ČSOB, Česká spořitelna, private PPF, Komerční banka, Raiffeisenbank and the Czech-Moravian Guarantee and Development Bank are in a group of efficient banks in the Czech Republic. Moreover, in Poland, the efficient group
PT
consists large banks (i.e. Bank Zachodni WBK, ING Bank Slaski and Bank Handlowy w Warszawie) and also some small special banks (SGB Bank, RBS Bank, and Bank of Tokyo).
CE
Boxplot is a convenient way to graphically depict groups of data through their quartiles. Figure 1 shows boxplot for ̅ =max and =min . In terms of inefficiency comparison by model (9), the best banks and financial institutions for year 2014 are in the Czech Republic
AC
(median 3.57). They are also followed by Hungarian banks (27.39), Polish banks (41.91) and last are banks in Slovakia (60.02). A similar analysis for model (13) testifies that the best situation is in the Czech Banks (median is 0.000), followed by Hungarian banks (5.91), Polish banks (13.78) and least effective banks are in Slovakia (24.77). The obtained results of the comparison of banking inefficiency in the EU countries using the proposed pessimistic and optimistic models are consistent with empirical studies and developments in financial systems in these countries, which have been affected by the global financial crisis since 2008. The Czech banking system is one of the most resilient systems in the 15
ACCEPTED MANUSCRIPT
EU due to the high liquidity of the sector, which is dominated by deposits from clients and is less dependent on interbank market financing. The Hungarian banking system was one of the smallest banking system in the EU with a large bank concentration, but with the development of this system, new types of credit and financial institutions were entering the market. The macroeconomic function of the transfer of deposits into loans and the microeconomic function of providing a complex of banking services were the most successful in the Czech and
CR IP T
Hungarian banks. The second group consists of Polish and Slovak banks. Poland has the largest banking industry out of the Visegrad countries and this banking system is focused on domestic business and plays an important role in financing private households, SMEs, big infrastructure projects. Polish banking sector is also owned by foreign-owned institution. The main reason for the low-efficiency of Polish banks is the loan quality and not a well-developed payment system. The high level of current earnings means that companies had no need to take out new loans, being able to fund investment and going operations internally. The size of Slovak banking
AN US
system is rather small and influenced by largest banks. The most of the banks have the
AC
CE
PT
ED
M
universal banking license.
Figure 1. Boxplot for (Ineff_M9) = ̅ =max
and (Ineff_M13) =
=min
As can be seen from the results in Table 5, the variations of inefficiency are significant. Accordingly, for the sake of illustration, we only consider boxplot with the inefficiency axis in the range 0-250 in Figure 1 to be more readable. On the one hand, through model (9), ignoring 12 share cases and based on the majority choice among 49 banks, all these banks treat the flexible measure as an output and hence deposit 16
ACCEPTED MANUSCRIPT
should be considered as an output with a pessimistic standpoint. On the other hand, applying model (13), in all 49 banks the flexible measure plays an input role and hence deposit should be considered as an input with an optimistic point of view. All in all, a profit-oriented banking system is desirable from the pessimistic standpoint; meanwhile, the optimistic viewpoint follows a service-oriented banking system. Returning to Table 4, we observe that there are some banks possess very high inefficiency score
CR IP T
from the pessimistic standpoint whereas these banks have better performance from the optimistic view. This significant variance is logical since, following from Theorems 1 and 3, the pessimistic model (9) finds a specific combination (among
various combinations) which
maximizes the sum of inefficiency scores while the optimistic model (13) introduces the combination with minimum sum of inefficiency scores. Moreover, some banks have recived high inefficiency score from both optimistic and pessimistic views. In order to find the main
AN US
sources of inefficiency, one can refer to Table 5 which provides useful information as a management guide to improve efficiency in performance of each inefficient bank. For example, Bank29 (Post Bank, SK) is a bank with high inefficiency score from both optimistic and pessimistic perspectives and referencing to Table 5 reveals that its first output (loans and advances) is the main source of inefficiency. In fact, the value of the first output of this bank is
M
400 which is significantly low compared to the values of other its input/ouput measures. Although there is no tie for determining the status of deposit in this example, we apply the summative models as well in order to see the results from the summative inefficiency
ED
perspective. After solving models (11) and (15), ( , the optimal objective function value and
)
and (
)
, are obtained, respectively. By comparing
the obtained results through both individual and summative envelopment (or multiplier) forms
PT
for this case study, deposit ends up the same role. However, for 61 banks, the final results are obtained by running 122 individual models (envelopment ant multiplier form), whereas,
CE
summative inefficiency needs to solve only two models for any number of DMU. It should be noted that there is no guarantee that the same results are always achieved through individual and summative version, as mentioned in Cook & Zhu (2007), because the summative
AC
approach might be overly sensitive to extreme DMUs, or possibly to the larger DMUs.
4
Conclusion
One of the controversial issues in many applications in DEA literature is the existence of flexible measure which can serve as either an input or output. Some studies have been done to classify inputs and outputs in DEA. This paper deployed a new classifier method based on directional distance function. As the main feature of this study, we developed a pair of non-radial CDD 17
ACCEPTED MANUSCRIPT
models in both envelopment and multiplier forms which introduce two pessimistic and optimistic approaches from both individual and summative viewpoints. Individual approach accommodates the flexible measure in favor of each DMU; nevertheless, in order to break the possible ties for classifying inputs and outputs, we have extended summative approach which designates the flexible measures in favor of all DMUs, simultaneously. Finally, to validate the discrimination power of the new procedure, an empirical study on the banking industry in V4
CR IP T
countries has been presented in which the deposit factor is taken into consideration as the flexible measure. It is found that deposit is treated as an output and as an input from the pessimistic and optimistic viewpoints, respectively, by both individual and summative models. Finally, extending the classifier models in the presence of imprecise, negative, or stochastic data can be considered as some interesting future research directions. Moreover, dealing with Malmquist index in the presence of selective measures is an alternative research direction.
AN US
Estimating the state of returns to scale in the company of flexible measure under variable returns to scale technology also could be a motivating topic for further study. Acknowledgements The
research
was
supported
by
the
European
Social
Fund
within
the
project
M
CZ.1.07/2.3.00/20.0296, the Czech Science Foundation through project No. 16-17810S as well as the SGS Project No. SP2017/141 VŠB-Technical University of Ostrava. All support is greatly
ED
acknowledged. References
PT
Amirteimoori, A., & Emrouznejad, A. (2011). Flexible measures in production process: A DEAbased approach. RAIRO - Operations Research, 45(1), 63–74.
CE
Amirteimoori, A., Emrouznejad, A., & Khoshandam, L. (2013). Classifying flexible measures in data envelopment analysis: A slack-based measure. Measurement, 46(10), 4100–4107.
AC
Bala, K., & Cook, W. D. (2003). Performance measurement with classification information: An enhanced additive DEA model. Omega, 31, 439–450. Banker, R. D., Charnes, A., & and Cooper, W. W. (1984). Some Models for Estimating Technical and Scale Inefficiencies in Data Envelopment Analysis. Management Science, 30(9), 1078– 1092. Beasley, J. E. (1990). Comparing university departments, 18(2), 171–183. Benston, G. J. (1965). Branch banking and economies of scale. The Journal of Finance, 20(2), 312– 331.
18
ACCEPTED MANUSCRIPT
Boda, M., & Zimková, E. (2015). Efficiency in the Slovak Banking Industry: A Comparison of Three Approaches. Prague Economic Papers, 2015(4), 434–451. Chambers, R. G., Chung, Y., & Färe, R. (1996). Benefit and Distance Functions. Journal of Economic Theory, 70(2), 407–419.
CR IP T
Chambers, R. G., Chung, Y., & Färe, R. (1998). Profit, Directional Distance Functions, and Nerlovian Efficiency. Journal of Optimization Theory and Applications, 98(2), 351–364. Charnes, A., Cooper, W. W., & Rhodes, E. (1978). Measuring the efficiency of decision making units. European Journal of Operational Research, 2(6), 429–444. Cook, W. D., Johnston, D. A., & Mccutcheon, D. (1992). Implementations of Robotics Identifying Efficient Implementers. Omega-International Journal of Management Science, 20(2), 227–239.
AN US
Cook, W. D., & Zhu, J. (2007). Classifying inputs and outputs in data envelopment analysis. European Journal of Operational Research, 180(2), 692–699. Cooper, W. W., Seiford, L. M., & Tone, K. (2007). Data envelopment analysis: A comprehensive text with models, applications, references and DEA-solver software (2nd edi.). Springer US.
M
Farrell, M. J. (1957). The Measurement of Productive Efficiency. Journal of the Royal Statistical Society. Series A (General), 120(3), 253.
ED
Hancock, D. (1991). A Theory of Production for the Financial Firm. Dordrecht: Kluwer Academic Publishers. Sealey, C. W., & Lindley, J. T. (1977). Inputs, outputs, and A theory of prouduction and cost at depository financial institutions. The Journal of Finance, 32(4), 1251–1266.
PT
Toloo, M. (2009). On classifying inputs and outputs in DEA: A revised model. European Journal of Operational Research (Vol. 198).
CE
Toloo, M. (2012). Alternative solutions for classifying inputs and outputs in data envelopment analysis. Computers & Mathematics with Applications, 63(6), 1104–1110.
AC
Toloo, M. (2014). Notes on classifying inputs and outputs in data envelopment analysis : a comment. European Journal of Operational Research, 198(2009), 1625–1628.
19
ACCEPTED MANUSCRIPT
APPENDIX Table A.1. Commercial banks of V4 banking sector subject to the analysis Country
Bank name
DMUs
Ceskoslovenska Obchodni Banka-
32
CZ
Ceska Sporitelna
33
PL
OTP Bank Plc
34
1
CZ
2
CZ
3
HU
4
PL
5
CZ
Komercni Banka
36
6
PL
Bank Zachodni WBK
37
7
PL
mBank
38
8
PL
ING Bank Slaski-Capital Group
9
CZ
10
PL
Bank Polska Kasa Opieki SA-Bank Pekao
35
PL
HU CZ
Unicredit Bank Czech Republic and Slovakia Getin Noble Bank
Czech Moravian Guarantee and Develpoment Bank Santander Consumer Bank SGB Bank
Budapest Bank Nyrt-Budapest Hitel-és Fejleszési Bank Nyrt Modra pyramida stavebni sporitelna as Raiffeisen
PL
Euro Bank
39
SK
Sberbank Slovensko
40
SK
Prima banka Slovensko
41
CZ
Air Bank
42
Bank Millennium
Bank name
CZ
AN US
PL
CSOB
M
11
Country
CR IP T
DMUs
HU
Bank of Hungarian Savings Cooperatives LimitedTAKAREKBANK
SK
Vseobecna Uverova Banka
43
PL
Idea Bank
13
PL
Bank Handlowy w Warszawie
44
HU
FHB Kereskedelmi Bank Zrt
14
SK
Tatra Banka
45
PL
Bank Pocztowy
ED
12
HU
K&H Bank Zrt
46
SK
OTP Banka Slovensko
16
PL
Bank BPH
47
PL
MBank Hipoteczny
17
PL
Nordea Bank Polska
48
CZ
Expobank CZ
18
HU
Erste Bank Hungary Nyrt
49
PL
HSBC Bank Polska
19
HU
MKB Bank Zrt
50
PL
RBS Bank (Polska)
51
HU
KDB Bank Europe
21
Ceskoslovenska obchodna banka-
CE
20
PT
15
SK
CSOB
HU
Raiffeisen Bank Zrt
52
PL
FM Bank PBP
HU
CIB Bank Ltd-CIB Bank Zrt
53
CZ
Equa Bank
23 24
PL
Alior Bank Spólka Akcyjna
54
SK
Privatbanka
HU
UniCredit Bank Hungary Zrt
55
PL
Pekao Bank Hipoteczny
25
PL
BNP Paribas Bank Polska
56
HU
Bank of China (Hungary)
26
CZ
GE Money Bank
57
HU
27
PL
Bank Ochrony Srodowiska
58
PL
28
CZ
J&T Banka
59
HU
AC
22
20
MagNet Hungarian Community Bank Bank of Tokyo - Mitsubishi UFJ (Polska) Magyar Cetelem Bank Rt
ACCEPTED MANUSCRIPT
29
SK
Post Bank JSC-Postova Banka
60
SK
CSOB Stavebna Sporitelna
30
CZ
PPF banka
61
CZ
Evropsko-ruska banka
31
CZ
Stavební Sporitelna Ceské Sporitelny
Table A.2. Mean value of data for each country in 2014 Country
banks
Inputs PE
FA
Flexible CRR
DEP
15
75519
84461
49.43
8340805
HU
14
92061
112030
53.64
4773693
PL
23
94523
58398
65.12
6806943
SK
9
43178
41233
61.91
3311167
V4 group
61
81709
74583
58.16
6201704
LA
M ED PT CE AC
21
SEC
NEA
1542066
2430671
597703
318278
1542982
496632
308651
1928066
573831
149533
1039011
216800
590683
1832105
509307
AN US
CZ
Output
CR IP T
No. of
ACCEPTED MANUSCRIPT
TABLES Table 1. Description of three groups of variables Measures
Name
Description
Input variables =PE
Total personnel expenses
Physical capital
=FA
Fixed assets
Credit risk
=CR
Net loans/Total assets
=DEP
Deposit and short-term funding
Loans and Advances
=LA
Loans and advances to banks
Investments
=SEC
Other securities
Non-interest income
=NEA
Non-earning assets
CR IP T
Labor
Flexible Variable Deposit
AN US
Output variables
Table 2. Data for the 61 commercial banks in V4 countries in 2014 256231
283822
49.21
2
328507
516325
50.5
3
686796
879261
60.16
4
447771
382676
5
262864
286920
6
331134
7
194588
8
223067
1038336
9760384
3328338
26480895
919459
8683189
3083983
63.94
30924048
2438522
8512511
1878279
54.74
25475452
4587331
6096395
2016568
152290
64.18
20419727
542477
5651880
2473068
170827
65.41
19595973
835646
6284417
827197
147307
55.97
17622988
444088
5920254
2047469
94073
69871
62.4
12979049
2405106
3159656
448750
91023
77806
75.38
13336015
332169
2222990
1043667
131683
39288
73.25
11518691
424106
2088996
966800
CE
11
10033454
2746297
PT
10
6977311
29304187
M
1
9
31888107
ED
DMU
12
107800
122100
65.54
8831600
771700
2806200
252900
13
152122
92587
33.55
7931612
852185
4780990
683453
14
93100
67.41
7579000
151600
1937000
785500
97719
141144
45.8
6598928
263473
3423997
682015
16
140494
77325
68.77
3292296
118130
1617585
614698
17
55418
24363
82.45
7040428
192445
562362
584028
18
88419
36216
60.03
6622321
444410
2152468
377061
19
70382
205181
64.23
4884371
330247
785638
977851
20
66600
71800
67.58
4870600
42200
1682800
173900
21
86637
37672
66.12
5075290
93241
1238500
608247
22
78387
75933
67.73
4962492
429687
944721
510011
AC
110300
15
22
110811
51783
76.94
5140868
61195
816315
428223
24
54822
89401
51.45
5004387
572462
1949398
87976
25
67093
30261
78.53
3340130
19066
627818
391222
26
80332
26972
72.4
3418511
40603
832440
454181
27
42803
28094
65.45
2947753
36688
1131536
297671
28
25441
7691
45.92
3313358
471603
1297045
360764
29
43300
21700
48.8
3331200
400
1536900
406200
30
8456
1531
30.12
3174855
925601
936682
755753
31
6816
14069
36.8
3450439
1204393
975134
92214
32
7946
5030
14.5
3258977
2328455
698347
20338
33
45475
16226
80.45
2500825
131370
358193
155612
34
15070
2504
38.49
2806127
145116
1502488
231758
35
74221
54960
62.58
2408961
197956
733082
152612
36
10657
13865
51.95
2663425
1030343
324648
69838
37
8310
2078
47.8
2760803
55037
1407155
94182
38
41792
17261
76.26
2373643
224945
264787
132887
39
19900
19400
70.28
1589000
225100
167200
185600
40
15900
20200
62.1
1639000
104600
516500
78400
41
12720
8675
44.45
1661198
70345
782980
146121
42
16346
21817
24.23
1694606
385220
791736
175860
43
18513
16463
65.62
7051
475974
103701
44
11898
14269
36.23
1576177
476847
590476
56752
45
22244
10761
68.48
1590237
8739
443438
96921
46
16500
21400
74.33
1220300
34100
254100
55100
14681
49
10087
50
AN US
M
2191
84.58
394712
5368
145839
20992
23228
75.36
973367
185857
55460
12965
1180
46.51
758177
34883
340402
88471
PT
5489
48
1472071
5754
361
39.2
614265
66203
317845
74219
6697
1634
52.71
659852
139002
36815
186582
16502
CE
47
CR IP T
23
ED
ACCEPTED MANUSCRIPT
51
1305
57.2
576290
7697
226439
47018
53
12746
3470
69.92
596810
57705
70709
67101
54
5400
1300
33.13
547500
11100
387300
10900
55
2961
120
92.93
114783
10761
16346
1926
56
3221
8846
35.32
311966
170943
63159
13390
57
3729
978
36.65
354212
13818
206997
28652
58
2143
169
44.34
248031
160114
36111
8835
59
9573
1115
87.75
197242
19134
1567
11860
60
2900
100
67.98
192300
5000
63100
2700
61
3007
3364
35.37
193540
45006
29585
60096
AC
52
23
ACCEPTED MANUSCRIPT
Table 3. Input-oriented v.s. output-oriented classifier models (Input-oriented)
Cook & Zhu
Toloo
0 0 0 1 1 0 0 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 1 1 0 0 1 0 1 0
0 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 1 0 1 1 0
0 0 0 1 1 0 1 0 1 1 1 0 0 0 0 0 1 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 0 0
0 0 1 0 0 0 0 0 0 0 0 1 0 1 1 1 0 0 1 1 1 1 1 1 1 1 1 1 1 0 0 0 1 0 1 1 0
3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19
ED
20 21 22
PT
23
AC
CE
25 26 27 28 29 30 31 32 33 34 35 36 37
AN US
2
CR IP T
Toloo
1
24
(Output-oriented)
Cook & Zhu
M
DMU
24
ACCEPTED MANUSCRIPT
41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60
AC
CE
PT
ED
61
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0
1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 0 1 1 0 1 0 1
CR IP T
40
1 1 1 1 1 1 1 1 1 1 1 1 0 1 1 1 1 0 1 1 0 1 0 1
AN US
39
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 1 0 1 0
M
38
25
ACCEPTED MANUSCRIPT
Table 4. Results and input/output status of deposit for each bank Model (9)
Model (7) (input)
(output)
Model (13)
̅ =max
=min
1
0
0
0
0
2
0
0
0
0
3
3.79
6.03
6.03
1
4
2.78
2.78
2.78
0
5
0.59
0
0.59
0
6
0
0
0
0
7
8.85
10.83
10.83
1
0
0
0
0
9
3.57
3.57
10
22.79
25.67
11
16.35
18.36
12
9.17
15.34
0
0
14
27.33
60.02
15
4.53
25.03
16
0
17 18 19
7.23
20 21
3.79
0
2.78
0
0
1
0
1
8.85
0
0
1
0
3.52
0
25.67
1
22.79
0
18.36
1
16.35
0
15.34
1
9.17
0
0
0
0
0
60.02
1
27.33
0
25.03
1
4.53
0
86.24
86.24
1
0
0
23.85
61.74
61.74
1
23.85
0
7.69
20
20
1
7.69
0
26.91
26.91
1
7.23
0
59.16
258.68
258.68
1
59.16
0
28.83
91.58
91.58
1
28.83
0
ED
PT
22
0
3.57
M
13
1
0
AN US
8
0
CR IP T
DMU
7.81
27.87
27.87
1
7.81
0
50.09
190.01
190.01
1
50.09
0
12.25
20.01
20.01
1
12.25
0
25
91.74
596.76
596.76
1
91.74
0
26
42.19
248.19
248.19
1
42.19
0
27
30.08
265.43
265.43
1
30.08
0
28
2.41
9.73
9.73
1
2.41
0
29
1975.06
15484.94
15484.94
1
1975.06
0
30
0
0
0
0
0
0
31
0
0
0
0
0
1
32
0
0
0
0
0
0
33
13.78
76.62
76.62
1
13.78
0
34
0
0
0
0
0
0
35
7.8
60.03
60.03
1
7.8
0
23
AC
CE
24
26
ACCEPTED MANUSCRIPT
9.14
17.04
17.04
1
9.14
0
37
0
0
0
0
0
0
38
9.46
56.58
56.58
1
9.46
0
39
6.58
38.25
38.25
1
6.58
0
40
9.44
48.51
48.51
1
9.44
0
41
5.54
51.73
51.73
1
5.54
0
42
0.83
10.47
10.47
1
0.83
0
43
102.59
758.9
758.9
1
102.59
0
44
4.58
19.73
19.73
1
4.58
0
45
102.42
605.33
605.33
1
102.42
0
46
24.77
150.37
150.37
1
24.77
0
47
33.55
220.96
220.96
1
33.55
0
48
24.26
142.03
142.03
1
24.26
0
49
5.51
29.24
29.24
1
5.51
0
50
0
0
0
0
0
0
51
0
30.74
30.74
1
0
0
52
28.51
125.13
125.13
1
28.51
0
53
8.15
69.73
69.73
1
8.15
0
54
0
105.35
105.35
1
0
0
55
18.79
41.91
41.91
1
18.79
0
56
4.15
30.7
30.7
1
4.15
0
57
2.9
50.65
50.65
1
2.9
0
58
0
0
0
0
0
0
59
78.81
565.52
565.52
1
78.81
0
60
25.62
28.82
28.82
1
25.62
0
61
0
31.13
31.13
1
0
0
AN US
M
AC
CE
PT
ED
CR IP T
36
27
ACCEPTED MANUSCRIPT
Table 5. The calculated input/output inefficiencies DMU 0
0
0
0
0
0
0
2
0
0
0
0
0
0
0
3
0.47
0.4
0
0.37
4.43
4
0.47
0.15
0
0.09
2.07
5
0.31
0.05
0
0.21
0
6
0
0
0
0
0
7
0.25
0
0
0.19
10.36
8
0
0
0
0
0
0.44
0.38
0
0.23
0
0
11
0.65
0
0
12
0
0.13
0
13
0
0
0
14
0.29
0
0
15
0
0.12
0
16
0.51
0
0
17
0.26
0
18
0.53
0
19
0
0.59
20
0.16
21
0.5
22
0.13
0.2
0.12
0.36
0
0
0
0.02
0
0
0
0.03
0
0
0
0
2.43
AN US
9 10
CR IP T
1
0.22
24.38
0.84
0
0.19
16.91
0.6
0
1.41
12.26
0.98
0.57
0
0
0
57.39
1.14
0
1.34
23.34
0.23
0
4.06
80.19
1.48
0
1.31
54.65
5.51
0
0
1
18.07
0.4
0
0
1.98
20.56
3.77
0
0.37
0
2.47
254.37
1.31
0
0
0
1.68
87.83
1.58
0
M
0
1.2
ED
0
0
2.39
22.12
3.23
0
0
0
2.41
183.22
3.9
0
0
0.23
0
0.7
2.03
0.36
16.68
0.34
0
0
3.98
587.54
4.89
0
26
0.51
0
0
3.36
241.3
3.02
0
27
0.1
0
0
3.81
259.74
1.78
0
28
0.3
0
0
1.11
6.97
0.34
1.01
29
0.34
0
0
1.95
15482.17
0.48
0
30
0
0
0
0
0
0
0
31
0
0
0
0
0
0
0
32
0
0
0
0
0
0
0
33
0.26
0
0
4.35
58.34
7.94
5.73
34
0
0
0
0
0
0
0
35
0.15
0
0
5.92
49.18
4.41
0.37
36
0
0.86
0.27
0.5
0.13
2.64
12.64
24
AC
CE
25
PT
0
0.49
23
28
ACCEPTED MANUSCRIPT
0
0
0
0
0
0
0
38
0.2
0
0
4.66
35.52
10.87
5.34
39
0
0.39
0.44
4.11
23.64
9.67
0
40
0
0.51
0.52
2.97
42.78
1.72
0
41
0
0.14
0.43
2.12
48.57
0.47
0
42
0
0.35
0
2.07
7.53
0.51
0
43
0
0.31
0.47
4.18
751.5
2.45
0
44
0
0.67
0
1.51
1.35
0.98
15.22
45
0
0
0.22
4.58
590.58
3.79
6.17
46
0
0.52
0.59
4.55
139.99
4.72
0
47
0
0
0.82
4.45
201.51
2.75
11.42
48
0
0.89
0.31
4.66
7.65
28.32
100.2
49
0.35
0
0.5
2.23
19.45
1.12
5.58
50
0
0
0
0
0
0
0
51
0
0.04
0.33
2.46
0.37
27.55
0
52
0.51
0
0
3.56
107.29
2.3
11.48
53
0
0
0.41
7.17
31.36
18.14
12.65
54
0
0
0.46
2.75
64.04
0.5
37.6
55
0.78
0
0.97
1.17
5.74
3.49
29.76
56
0
0.93
0.68
2.88
1.06
4.65
20.5
57
0
0
0.67
3.02
58
0
0
59
0
60
0.78
61
0
AN US
CR IP T
37
0.92
8.6
0
0
0
0
0
0.42
10.1
27.01
490.54
37.45
0
0.96
0.06
10.49
0
16.52
5.28
16.12
8.4
0
ED
M
37.44
0
0.82
AC
CE
PT
0.52
29