Next Article in Journal
Classes of N-Dimensional Nonlinear Fokker-Planck Equations Associated to Tsallis Entropy
Previous Article in Journal
Generalized Hamilton’s Principle for Inelastic Bodies Within Non-Equilibrium Thermodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Generalized Maximum Entropy Stochastic Frontier Measuring Productivity Accounting for Spatial Dependency

1
LEI, Part of Wageningen UR, PO Box 29703, The Hague 2585 DB, The Netherlands
2
International Rice Research Institute, DAPO Box 7777, Metro Manila 1301, The Philippines
*
Author to whom correspondence should be addressed.
Entropy 2011, 13(11), 1916-1927; https://doi.org/10.3390/e13111916
Submission received: 29 September 2011 / Revised: 17 October 2011 / Accepted: 24 October 2011 / Published: 28 October 2011

Abstract

:
In this paper, a stochastic frontier model accounting for spatial dependency is developed using generalized maximum entropy estimation. An application is made for measuring total factor productivity in European agriculture. The empirical results show that agricultural productivity growth in Europe is driven by upward movements of technology over time through technological developments. Results are then compared for a situation in which spatial dependency in the technical inefficiency effects is not accounted.
PACS Codes:
89.65.Gh; 89.70.Cf

1. Introduction

This paper develops a stochastic frontier model accounting for spatial dependency through spatial econometrics techniques and using generalized maximum entropy (GME) estimation. An application is provided computing Total Factor Productivity (TFP) growth for the agricultural sector of the EU-27 countries as well as three European Union (EU) candidate countries (CC): Croatia. Former Yugoslav Republic of Macedonia (FYROM) and Turkey.
Stochastic frontiers, simultaneously introduced by Aigner [1], Meeusen [2] and Battese [3], usually comprise dual error components, one accounting for the random error and the other for the technical inefficiency. In order to compute total factor productivity growth a time-varying inefficiency model is usually required [4]. A farm-specific temporal variation in technical inefficiency as developed by Cuesta [5] is applied to allow countries to have different temporal technical inefficiency patterns.
Spatial dependence is accounted for in the technical inefficiency component of the stochastic frontier. In the literature, [6] detected that the employment rate exhibits a correlation with economic growth in neighbouring counties in the American Midwest, while [7] found that technical efficiency scores can display spatial dependency given that their determinants are often spatially correlated. Therefore, it important to account for the spatial dependency between countries since economic conditions in neighbouring countries are likely to be similar as a result of spillovers of economic activities across borders and even more for the countries that are part of the EU. These knowledge externalities can create dependency in the technical efficiency. Moreover, several other unobserved or immeasurable characteristics may be similar across country boundaries, which may lead to dependency in the residuals of empirical models of technical efficiency. The spatial dependency is introduced through a spatial weight matrix exogenously defined to capture the spillover effects between countries included in the study. In so doing, the technical inefficiency component is allowed to follow a spatially autoregressive error process [8].
Standard maximum likelihood stochastic frontier estimations require ad-hoc assumptions on the distribution of the inefficiency component. The GME estimator avoids the need for making unnecessary assumptions on the inefficiency component. Maximum entropy estimators were introduced by Golan [9], and proved to be robust under ill-posed and ill-conditioned problems. A GME is particularly suited in this context given the very large number of parameters to be estimated. The data are from the Food and Agriculture Organization (FAO) for the period 1993–2006 comprising former EU-15 countries, EU-12 as well as three CC.
The paper is structured as follows: Section 2 contains the time-varying technical inefficiency stochastic frontier model accounting for spatial dependency. Section 3 describes the data and the empirical model. Section 4 presents the empirical results. Finally, in Section 5 we provide our conclusions.

2. The Empirical Model

A stochastic production frontier is estimated in order to compute TFP growth which measures the change in total output relative to changes in the use of all inputs. Approaches using stochastic-frontier production functions have been independently proposed by Aigner [1], Meeusen [2] and Battese [3]. Schmidt [10] extended to panel data the production frontier model. For a stochastic production frontier estimated by GME see Macedo [11]. The agricultural production of country i (i = 1, ..., N) at time t (t = 1, ..., T), Y i t , is assumed to be produced using the input array, X i t constituted by land ( X 1 i t ) , machinery ( X 2 i t ) , labor ( X 3 i t ) , fertilizer ( X 4 i t ) and livestock ( X 5 i t ) respectively.
A trascendetal logarithmic (i.e. translog) production frontier is selected for its flexibility in measuring TFP growth. It provides a second-order differential approximation which is linear in parameters. A translog production frontier allows production elasticities to vary at each data point and non-neutral Hicksian technical change. The empirical model is specified following the ‘true’ fixed effects model of Greene [12,13]. The ‘true’ fixed effects model allows to capture potential country-specific time invariant heterogeneity. However when the number of time periods available in each panel is not large enough the ‘true’fixed effects model may be subject to the incidental parameter problem. Countries can be characterized by different stage of development, liberalization and different conditions in terms of soil quality, drainage, altitude and climate. The empirical model is specified as follows:
y i t = α i + n =1 N β n x i t n + β t t + 1 2 n = 1 N m = 1 N β n m x i t n x i t m + 1 2 β t t t 2 + n = 1 N β n t x i t n t Z i t + ε i t
where lower case letters (y,x) indicate natural logs of upper case letters (Y,X), α i = ( α 0 , α 1 , , α I ) is the vector of country-specific intercepts capturing time invariant heterogeneity, β n = ( β 0 , β 1 , , β N ) is the vector of unknown parameters to be estimated, t a time trend used in order to account for technological change. The inefficiency Z i t = l n ( τ i t ) is assigned a non-negative random variable and the random error ε i t = l n ( ξ i t ) has a symmetric distribution with mean zero.
Regularity conditions are a priori imposed during the estimation to ensure that input elasticities are nonnegative at all observed output levels [14] with the following inequality restrictions:
E x i t , n y i t x i t n = ζ i t , n 0 ,  for  n = 1 , ... , 5.
A generalization of the time-varying inefficiency model of Battese [15] proposed by Cuesta [5] is implemented and extended to account for spatial dependency in the technical inefficiency effects. Cuesta’s [5] function is specified as:
Z i t = exp [ η i ( t T ) ] Z i
where the temporal pattern of inefficiency effects (i.e., ηi) is now a country-specific parameters responsive to different temporal variations among countries.
Spatial dependency is introduced in the inefficiency terms equation (3) by the following expression:
Z i = [ ( I ρ W ) 1 Ω i ]
where I is an identity matrix, ρ the parameter accounting for spatial dependency in the inefficiency effects, W the spatial weighting matrix and Ω i a stochastic term ensuring that the technical efficiency exp ( Z i t ) lies between zero and one. Substituting Equation (4) into (3):
Z i t = exp [ η i ( t T ) ] [ ( I ρ W ) 1 Ω i ]
First of all, different weight matrix forms are possible, such as contiguity weights and distance weights. The choice of the weight matrix depends on the study at hand and the subject being investigated. In this study, a distance-based weight matrix is more appropriate because not all countries involved are contiguous (sharing borders). The weight matrix is defined a priori and exogenously based on arc distances between the geographical midpoints of the countries. The neighbourhood is defined using a cut-off distance that enforces all countries to have at least one neighbour. Two countries are considered neighbors when the geographical distance (arc distance) between them is less than a threshold distance. The weight matrix therefore takes value one when countries are neighbors and zero when they are not. The threshold distance in this case is chosen in such a way to guarantee all countries have at least one neighbor. No country should be a neighbor to itself. This forces the weight matrix to have zero on the diagonal. The weight matrix is further standardized. Mathematical and statistical reasons for standardization can be found in Bell and Bockstael [16]. The technical efficiency of each country in each year is obtained through the conditional expectation of exp ( Z i t ) , given the value of ( ε i t Z i ) .
TFP growth is computed and decomposed into the efficiency, technical change and scale components through the identity with the so-called Malmquist index (MI) [17]. The MI can be defined by using either an input or an output orientation. For country-level analysis, an output orientation is the proper choice [4]. An output orientation looks for the maximal proportional expansion of an output vector, given an input or resource vector. The efficiency change component is given by:
E C i s = exp ( Z i t Z i s ) = Z i s / Z i t
The technical change component requires one to evaluate the partial derivatives of the production frontier with respect to time using the data for the i-th country in period s and t . Then the technical change between the adjacent periods s and t can be derived through the geometric mean of the aforementioned partial derivatives. In the case of a translog specification, this is equivalent to the exponential of the arithmetic mean of the log derivatives as given by:
T C i s = exp { 1 2 [ y i s s + y i t t ] }
To detect potential scale change effects, scale change is introduced in computing TFP following Orea [18], who uses Diewert’s quadratic identity to derive a MI. The scale change component is given by:
S C i s = exp { 1 2 n = 1 N ( ζ i s n S F i s + ζ i t n S F i t n ) ( x i s n / x i t n ) }
where S F i s = ( ζ i s 1 ) / ζ i s , ζ i s = n = 1 N ζ i s n (return to scale) and ζ i s n = y i s x i s n (input elasticities). Each single component is then summed up to recover the MI of TFP growth.

3. Data and the Estimation

Data are drawn from the FAOSTAT database of the Statistics Division of the Food and Agriculture Organization (FAO) of the United Nations. The data include the former EU-15 countries as well as the new member states EU-12, notably: Austria, Belgium-Luxembourg, Denmark, Finland, France, Germany, Greece, Ireland, Italy, The Netherlands, Portugal, Spain, Sweden, and United Kingdom, Bulgaria, Czech Republic, Cyprus, Estonia, Hungary, Latvia, Lithuania, Malta, Poland, Romania, Slovakia, and Slovenia. In addition three CC to the EU accession are considered: Croatia, the Former Yougoslav Republic of Macedonia (FYROM) and Turkey. The data cover the period 1993–2006. The sample includes 29 countries for a period of 14 years yielding 406 observations.
One agricultural aggregate output variable, defined as agricultural production net of quantities of various commodities used as feed and seed, is considered. The output series is based on international commodity average prices (1999–2001 1000 I$) and expressed in a single currency unit. Five input variables are taken into account. Land comprises arable land land and permanent crops and the area for permanent meadows and pastures. It is measured in 1000 Ha. Agricultural tractors refers to wheel and crawler or track-laying tractors (excluding garden tractors) used in agriculture without reference to the horsepower of the tractors. It is measured in number (1000 No) of tractors in use. Agricultural population defined as all persons depending for their livelihood on agriculture, hunting, fishing and forestry. It comprises all persons economically active in agriculture as well as their non-working dependents. It is measured in thousands of people (Pe). Fertiliser is an aggregate of Nitrogen (N), Potassium (P) and Phosphate (K) consumed in agriculture and expressed in tonnes of nutrients. Livestock is constructed by aggregating five categories of animals into sheep equivalent. The categories of animals considered are: buffaloes, cattle, pigs, sheep and goats. Numbers of these animals are converted into sheep equivalents (SE) using the same conversion factors as in Coelli [4]. Table 1 presents descriptive statistics for the data.
Table 1. Descriptive statistics, 1993–2006.
Table 1. Descriptive statistics, 1993–2006.
VariablesUnitsAvgMinMaxStd
Agr. Output×103 I$8281701.253350.038125640.09950171.2
Land×103 Ha8335.29.041223.010174.4
Machinery×103 No349.10.51900.0458.2
Labour×103 Pe875.92.010454.01854.0
Fertilizer×103 Mt740.60.95064.01019.5
Livestock×103 SE42007.8236.8191956.050067.3
A GME estimator is used. The advantages and properties of the GME estimator are extensively discussed in Golan [9], Mittelhammer [19] and Perloff [20]. Here we only list the main advantages for the problem at hand. First, the GME estimator efficiently exploits all the information contained in each data point, instead of using moment conditions as done by the more traditional estimators. Second, it is less influenced by outlying observations because of the weighting between signal and noise in the objective function. Third, it is a robust estimator, even when noise is not normally distributed and/or the design matrix exhibits high condition indexes. Multicollinearity is often problematic in stochastic frontier analysis affecting the sign of the estimated coefficients as well as the precision of the estimated coefficients. Fourth, the GME estimator enables the researcher to easily impose prior information on parameters [21]. Finally, the GME estimator does not require strong behavioural assumptions on the underlying data generating process.
The parameters of the stochastic production frontier are reparameterized according to the GME formalism as follows:
β k = m = 1 5 p k m v k m  and  η i = u = 1 3 p i u v i u
with p k m and p i m being the probability vector associated with the stochastic production frontier parameters to be estimated and the temporal pattern of inefficiency effects, v k m and v i m representing their support values. The supports for v k m are initially set to ( 1000 , 500 , 0 , 500 , 1000 ) and the supports for v i u to ( 0.5 , 0 , 0.5 ) . The inefficiency terms and random error are defined by:
Ω i = m = 1 5 w i n g i n and ε i t = l = 1 3 w i t l g i t l
with w i n and w i n l being the probability vectors to be estimated and g i n and g i t l representing the support values. With respect to the number of support points, Golan [9] show through simulation experiments that the greatest improvements in the estimate precision is obtained by selecting five or more support values. The number of support values is set to five for the stochastic frontier parameters and the inefficiency terms whereas it is set to three for the random error. For the signal space definition, when the researcher does not know the possible realizations, it is suggested that v be set to be symmetrical around zero with large negative and positive boundaries [22]. The supports for the stochastic terms Ω i are set following Campbell [23] to ( 0 , 0.005 , 0.01 , 0.0015 , 1 ) . The noise support bounds are specified following the widely accepted three sigma rule of Pukelsheim [24], see also Golan [9], defining a uniform zero-centred symmetric distribution with a support space specified as follows: ( 3 τ ˜ , 1.5 τ ˜ , 0 , 1.5 τ ˜ , 3 τ ˜ ) , where τ ˜ is the empirical standard deviation of the dependent variable. The parameter ρ accounting for spatial dependency is reparameterized as follows:
ρ = m = 1 5 r m s m
where r m being the probability vector to be estimated and s m representing the support values defining a uniform symmetric distribution with a support space specified as follows: ( 0 , 0.25 , 0.5 , 0.75 , 1 ) .
The GME objective function H to maximize is defined as:
H ( p , w , w , r ) = k = 1 29 m = 1 5 p k m ln p k m i = 1 29 u = 1 3 p i u l n p i u i = 1 29 n = 1 5 w i n ln w i n i = 1 29 t = 1 14 l = 1 3 w i t l ln w i t l m = 1 5 r m ln r m
subject to Equation (1), Equation (2), Equation (5), and Equations (9–11) and the following adding-up constraints
m = 1 5 p k m = 1 , m = 1 3 p i u = 1 , n = 1 5 w i n = 1 , l = 1 3 w i t l = 1 ,  and  m = 1 5 r m = 1
for all values of k, i and t.

4. Empirical Results and Discussion

In this section, we discuss the final estimates. The estimates for the time-varying technical inefficiency stochastic frontier model accounting for spatial dependency estimated by GME, are displayed in Table 2.
Table 2. Estimated coefficients.
Table 2. Estimated coefficients.
CoefGME1Support1GME2Support2GME3Support3
β10.2005m = (0, 0.1, 0.2, 0.3, 0.4)0.2005m = (0, 0.1, 0.2, 0.3, 0.4)0.2005m = (0, 0.1, 0.2, 0.3, 0.4)
β20.1996m = (0, 0.1, 0.2, 0.3, 0.4)0.1996m = (0, 0.1, 0.2, 0.3, 0.4)0.1996m = (0, 0.1, 0.2, 0.3, 0.4)
β30.1999m = (0, 0.1, 0.2, 0.3, 0.4)0.1999m = (0, 0.1, 0.2, 0.3, 0.4)0.1999m = (0, 0.1, 0.2, 0.3, 0.4)
β40.1991m = (0, 0.1, 0.2, 0.3, 0.4)0.1991m = (0, 0.1, 0.2, 0.3, 0.4)0.1991m = (0, 0.1, 0.2, 0.3, 0.4)
β50.2008m = (0, 0.1, 0.2, 0.3, 0.4)0.2008m = (0, 0.1, 0.2, 0.3, 0.4)0.2008m = (0, 0.1, 0.2, 0.3, 0.4)
βt0.0126m = (−1000, −500, 0, 500, 1000)0.0126m = (−500, −250, 0, 250, 500)0.0126m = (−250, −125, 0, 125, 250)
β110.2159m = (−1000, −500, 0, 500, 1000)0.2159m = (−500, −250, 0, 250, 500)0.2158m = (−250, −125, 0, 125, 250)
β12−0.1479m = (−1000, −500, 0, 500, 1000)−0.1478m = (−500, −250, 0, 250, 500)−0.1477m = (−250, −125, 0, 125, 250)
β13−0.0098m = (−1000, −500, 0, 500, 1000)−0.0098m = (−500, −250, 0, 250, 500)−0.0097m = (−250, −125, 0, 125, 250)
β14−0.0272m = (−1000, −500, 0, 500, 1000)−0.0272m = (−500, −250, 0, 250, 500)−0.0271m = (−250, −125, 0, 125, 250)
β150.1382m = (−1000, −500, 0, 500, 1000)0.1382m = (−500, −250, 0, 250, 500)0.1381m = (−250, −125, 0, 125, 250)
β1t0.0108m = (−1000, −500, 0, 500, 1000)0.0108m = (−500, −250, 0, 250, 500)0.0108m = (−250, −125, 0, 125, 250)
β220.1058m = (−1000, −500, 0, 500, 1000)0.1058m = (−500, −250, 0, 250, 500)0.1058m = (−250, −125, 0, 125, 250)
β230.1171m = (−1000, −500, 0, 500, 1000)0.1171m = (−500, −250, 0, 250, 500)0.1170m = (−250, −125, 0, 125, 250)
β24−0.0419m = (−1000, −500, 0, 500, 1000)−0.0419m = (−500, −250, 0, 250, 500)−0.0419m = (−250, −125, 0, 125, 250)
β250.0661m = (−1000, −500, 0, 500, 1000)0.0661m = (−500, −250, 0, 250, 500)0.0660m = (−250, −125, 0, 125, 250)
β2t0.0066m = (−1000, −500, 0, 500, 1000)0.0066m = (−500, −250, 0, 250, 500)0.0066m = (−250, −125, 0, 125, 250)
β330.0185m = (−1000, −500, 0, 500, 1000)0.0184m = (−500, −250, 0, 250, 500)0.0184m = (−250, −125, 0, 125, 250)
β34−0.0539m = (−1000, −500, 0, 500, 1000)−0.0539m = (−500, −250, 0, 250, 500)−0.0539m = (−250, −125, 0, 125, 250)
β35−0.0667m = (−1000, −500, 0, 500, 1000)−0.0667m = (−500, −250, 0, 250, 500)−0.0667m = (−250, −125, 0, 125, 250)
β3t−0.0100m = (−1000, −500, 0, 500, 1000)−0.0100m = (−500, −250, 0, 250, 500)−0.0100m = (−250, −125, 0, 125, 250)
β440.0228m = (−1000, −500, 0, 500, 1000)0.0228m = (−500, −250, 0, 250, 500)0.0229m = (−250, −125, 0, 125, 250)
β450.1690m = (−1000, −500, 0, 500, 1000)0.1689m = (−500, −250, 0, 250, 500)0.1688m = (−250, −125, 0, 125, 250)
β4t−0.0064m = (−1000, −500, 0, 500, 1000)−0.0064m = (−500, −250, 0, 250, 500)−0.0064m = (−250, −125, 0, 125, 250)
β55−0.4092m = (−1000, −500, 0, 500, 1000)−0.4091m = (−500, −250, 0, 250, 500)−0.4087m = (−250, −125, 0, 125, 250)
β5t−0.0012m = (−1000, −500, 0, 500, 1000)−0.0012m = (−500, −250, 0, 250, 500)−0.0012m = (−250, −125, 0, 125, 250)
βtt−0.0017m = (−1000, −500, 0, 500, 1000)−0.0017m = (−500, −250, 0, 250, 500)−0.0017m = (−250, −125, 0, 125, 250)
ρ0.5001m = (0, 0.25, 0.5, 0.75, 1)0.5001m = (0, 0.25, 0.5, 0.75, 1)0.5001m = (0, 0.25, 0.5, 0.75, 1)
All variables, before estimation, are rescaled in order to have unit means so that the estimated first-order coefficients represent the elasticities of output with respect to inputs when evaluated at the variable means. One of the advantages of maximum entropy is that it allows to include prior information on parameters either introducing additional consistency constraints characterizing the system at hand or modifying the support values of the parameters to be estimated. Three different support spaces are specified for the estimated models (GME1, GME2, and GME3). The support space for the first order coefficients is kept fixed at m = ( 0 , 0.1 , 0.2 , 0.3 , 0.4 ) since we attribute a prior expectation of 0.2 to all five first order coefficients in order to fulfill globally monotonicity ensuring a well behaved production frontier function. For the remaining coefficients the supports bounds are progressively shrinked although keeping a wide support definition introducing minimal prior information on the coefficient to be estimated. The estimated coefficients appear to be stable to different choices of support values. The first order coefficient of the time trend variable provides information on the average annual rate in neutral technical change. The annual percentage change in output due to technical change is estimated to be 1.3 percent. The coefficient of time squared is negative indicating that the rate of technical change is increasing at a decreasing rate through time. The coefficient of time interacted with the land, machinery, labour, fertilizer and livestock input are positive for the first two input variables and negative for the remaining three input variables. This suggests that technical change has been land-, machinery-saving but labour-, fertilizer- and livestock-using over the time considered. The country-specific intercepts and the coefficients capturing the temporal variation of the technical inefficiency for the estimated model are not reported to save space but are available upon request.
Table 3. Estimated technical efficiency scores for 2006.
Table 3. Estimated technical efficiency scores for 2006.
CountriesGME1ASupport1CountriesGME1BSupport2
Austria0.80499w = (0. 0.005. 0.010. 0.015. 1)Austria0.86293w = (0. 0.005. 0.010. 0.015. 0.67)
Denmark0.80447w = (0. 0.005. 0.010. 0.015. 1)Denmark0.86261w = (0. 0.005. 0.010. 0.015. 0.67)
Greece0.80180w = (0. 0.005. 0.010. 0.015. 1)Greece0.86062w = (0. 0.005. 0.010. 0.015. 0.67)
Germany0.80153w = (0. 0.005. 0.010. 0.015. 1)Germany0.86052w = (0. 0.005. 0.010. 0.015. 0.67)
Finland0.79926w = (0. 0.005. 0.010. 0.015. 1)Finland0.85875w = (0. 0.005. 0.010. 0.015. 0.67)
Sweden0.79923w = (0. 0.005. 0.010. 0.015. 1)Sweden0.85874w = (0. 0.005. 0.010. 0.015. 0.67)
France0.79763w = (0. 0.005. 0.010. 0.015. 1)France0.85758w = (0. 0.005. 0.010. 0.015. 0.67)
Italy0.79735w = (0. 0.005. 0.010. 0.015. 1)Italy0.85736w = (0. 0.005. 0.010. 0.015. 0.67)
Netherlands0.79659w = (0. 0.005. 0.010. 0.015. 1)Netherlands0.85682w = (0. 0.005. 0.010. 0.015. 0.67)
Belgium0.79392w = (0. 0.005. 0.010. 0.015. 1)Belgium0.85504w = (0. 0.005. 0.010. 0.015. 0.67)
Ireland0.75978w = (0. 0.005. 0.010. 0.015. 1)Ireland0.82970w = (0. 0.005. 0.010. 0.015. 0.67)
Portugal0.75977w = (0. 0.005. 0.010. 0.015. 1)Portugal0.82970w = (0. 0.005. 0.010. 0.015. 0.67)
Spain0.75973w = (0. 0.005. 0.010. 0.015. 1)Spain0.82968w = (0. 0.005. 0.010. 0.015. 0.67)
United Kingdom0.75968w = (0. 0.005. 0.010. 0.015. 1)United Kingdom0.82966w = (0. 0.005. 0.010. 0.015. 0.67)
Slovakia0.80520w = (0. 0.005. 0.010. 0.015. 1)Slovakia0.86310w = (0. 0.005. 0.010. 0.015. 0.67)
Poland0.80483w = (0. 0.005. 0.010. 0.015. 1)Poland0.86283w = (0. 0.005. 0.010. 0.015. 0.67)
Hungary0.80473w = (0. 0.005. 0.010. 0.015. 1)Hungary0.86275w = (0. 0.005. 0.010. 0.015. 0.67)
Czech Republic0.80462w = (0. 0.005. 0.010. 0.015. 1)Czech Republic0.86268w = (0. 0.005. 0.010. 0.015. 0.67)
Slovenia0.80454w = (0. 0.005. 0.010. 0.015. 1)Romania0.86265w = (0. 0.005. 0.010. 0.015. 0.67)
Romania0.80454w = (0. 0.005. 0.010. 0.015. 1)Slovenia0.86261w = (0. 0.005. 0.010. 0.015. 0.67)
Malta0.80235w = (0. 0.005. 0.010. 0.015. 1)Malta0.86102w = (0. 0.005. 0.010. 0.015. 0.67)
Bulgaria0.80212w = (0. 0.005. 0.010. 0.015. 1)Bulgaria0.86084w = (0. 0.005. 0.010. 0.015. 0.67)
Lithuania0.80070w = (0. 0.005. 0.010. 0.015. 1)Lithuania0.85981w = (0. 0.005. 0.010. 0.015. 0.67)
Latvia0.80002w = (0. 0.005. 0.010. 0.015. 1)Latvia0.85932w = (0. 0.005. 0.010. 0.015. 0.67)
Estonia0.79925w = (0. 0.005. 0.010. 0.015. 1)Estonia0.85875w = (0. 0.005. 0.010. 0.015. 0.67)
Cyprus0.79803w = (0. 0.005. 0.010. 0.015. 1)Cyprus0.85785w = (0. 0.005. 0.010. 0.015. 0.67)
Croatia0.80490w = (0. 0.005. 0.010. 0.015. 1)Croatia0.86287w = (0. 0.005. 0.010. 0.015. 0.67)
FYROM0.80440w = (0. 0.005. 0.010. 0.015. 1)FYROM0.86256w = (0. 0.005. 0.010. 0.015. 0.67)
Turkey0.79041w = (0. 0.005. 0.010. 0.015. 1)Turkey0.85228w = (0. 0.005. 0.010. 0.015. 0.67)
Standard Deviation0.01496Standard Deviation0.01100
Table 3 presents the estimated technical efficiency scores assuming two different upper bounds in the support of the inefficiency component. Following Campbell [23] we estimated two models (GME1A and GME1B) where the upper limit for the support points is changed. The model GME1A assumes a prior technical efficiency mean of 81.4 percent with a maximum of efficiency of 100 percent and a minimum of 36.8 percent. The model GME1B assumes a prior technical efficiency mean of 86.9 percent with a maximum of efficiency of 100 percent and a minimum of 51.2 percent. The ranking of the technical efficiency scores is preserved under the two different distributional assumptions with the exception made for Romania and Slovenia. GME1A displays higher standard deviation in the estimates and a larger entropy value.
The decomposition of TFP growth for each country is presented in Table 4. On average the EU-12 is characterized by higher agricultural productivity growth rates as compared to EU-15 and CC. TFP growth ranges from 0.48 for Turkey to 3.48 for Slovenia. All the countries under analysis showed productivity growth. Efficiency change ranges from −0.20 for Belgium to 0.17 for Germany. Technical change ranges from −0.34 for Croatia to 2.49 for Slovenia. Scale change ranges from -0.68 for France to 2.87 for Malta. The estimates suggest that productivity growth was driven by technical change. The estimates of technical efficiency change were then compared with a model not allowing for spatial dependency ( ρ = 0 ) . In Figure 1 cumulated technical efficiency change is reported for the model allowing for spatial dependency and for the more standard model not including spatial dependency. From this application, it appears that not accounting for spatial dependency has an impact on the final technical inefficiency estimates. The cumulated technical efficiency changes appear underestimated in the standard model not accounting for spatial dependency.
Figure 1. Cumulated technical efficiency changes, (percent).
Figure 1. Cumulated technical efficiency changes, (percent).
Entropy 13 01916 g001
Table 4. Decomposition of TFP growth for 2006, (percent).
Table 4. Decomposition of TFP growth for 2006, (percent).
CountriesEfficiency ChangeTechnical ChangeScale ChangeTFP Change
Austria0.01502.1045−0.10272.0168
Denmark0.11081.50950.13931.7595
Sweden−0.04791.8101−0.03621.7259
Spain0.07731.60650.00861.6924
Finland−0.00971.6997−0.01041.6796
Ireland0.00731.38440.00691.3986
United Kingdom−0.09911.5839−0.22591.2589
Belgium−0.19521.17070.17651.1519
France0.05041.7017−0.67951.0725
Italy−0.00861.5461−0.48851.0490
Greece−0.05971.1068−0.05470.9365
Germany0.16861.2156−0.61970.7645
Portugal0.03380.7013−0.01420.7209
Netherlands0.06040.32250.23510.6180
EU-15 average0.00741.3902−0.11901.2747
Slovenia0.03452.49050.96453.4894
Estonia−0.01281.75391.28143.0225
Malta0.0565−0.00442.86952.9216
Latvia−0.04861.73210.72462.4081
Cyprus−0.01030.49351.20801.6912
Lithuania−0.04581.56180.02131.5373
Bulgaria−0.02891.31120.03681.3192
Slovakia−0.06660.44430.42150.7992
Hungary0.04650.8790−0.13430.7912
Romania0.10220.8323−0.36320.5714
Poland−0.05640.58740.00090.5319
Czech Republic−0.06200.5649−0.01840.4846
EU-12 average−0.00761.05390.58441.6306
Fyrom0.11461.84250.32552.2826
Croatia0.0089−0.34941.38511.0446
Turkey−0.01870.13200.37070.4841
CC average0.03490.54170.69371.2704

5. Conclusions

A time-varying technical inefficiency stochastic frontier model was developed accounting for spatial dependency through spatial econometrics techniques and using a GME estimation procedure. The estimated model also controlled for country-specific time invariant heterogenity. An application was provided computing total factor productivity growth for the agricultural sector of the EU-27 countries as well as three CC: Croatia, FYROM and Turkey.
The GME estimator allowed us to estimate a model with a relatively large number of parameters and include prior information on parameters by defining support bounds and adding specific consistency constraints. The inclusion of prior information helped in decreasing multicollinearity issues among the exogenous variables The economic analysis showed that agricultural productivity growth in Europe is driven by upward movements of technology over time through technological developments. EU-12 countries are characterized by larger agricultural productivity growth rates as compared to the EU-15 countries and the CC. Most of the countries that were well below the frontier at the beginning of the sample period are displaying high TFP growth rates evidencing catch-up. The empirical results show that not allowing for spatial dependency underestimate the cumulated technical inefficiency changes.
Topics for future research are the programming of standard errors with bootstrapping techniques and to examine the consistency properties of the GME estimator. Care should be taken when interpreting the results since the input data used did not allow to control for differences in input qualities across countries. In additon, the potential incidental parameter problem and related bias should be limited in our application given the relatively large number of time periods available in our panel data set.

References

  1. Aigner, D.J.; Lovell, C.A.K.; Schmidt, P.J. Formulation and estimation of stochastic frontier production function models. J. Econom. 1977, 6, 21–37. [Google Scholar] [CrossRef]
  2. Meeusen, W.; van den Broeck, J. Efficiency estimation from Cobb-Douglas production functions with composed error. Int. Econ. Rev. 1977, 18, 435–444. [Google Scholar] [CrossRef]
  3. Battese, G.E.; Corra, G. Estimation of a production frontier model with application to the pastoral zone of Easter Australia. Aust. J. Agric. Econ. 1977, 21, 167–179. [Google Scholar]
  4. Coelli, T.J.; Prasada Rao, D.S.; O’Donnell, C.G.; Battese, G.E. An Introduction to Efficiency and Productivity Analysis, 2nd ed.; Kluwer Academic: Boston, MA, USA, 2005; p. 367. [Google Scholar]
  5. Cuesta, R. A production model with firm specific temporal variation in technical efficiency: With application to Spanish dairy farms. J. Prod. Anal. 2000, 13, 139–158. [Google Scholar] [CrossRef]
  6. Monchuk, D.C.; Miranovski, J.A. Spatial Labor Market and Technology Spillovers—Analysis from US Midwest, No. 22250. In Proceedings of the 2003 Annual Meeting, Montreal, Canada, 27–30 July 2003.
  7. Druska, V.; Horrace, W.C. Generalized moment estimation for spatial panel data: Indonesian rice farming. Am. J. Agric. Econ. 2004, 86, 185–198. [Google Scholar] [CrossRef]
  8. Anselin, L. Spatial Econometrics: Methods and Models, 1st ed.; Kluver: London, UK, 1988; p. 304. [Google Scholar]
  9. Golan, A.; Judge, J.; Miller, D. Maximum Entropy Econometrics: Robust Estimation with Limited Data, 1st ed.; John Wiley and Sons: Chichester, UK, 1996; p. 324. [Google Scholar]
  10. Schmidt, P.; Sickles, R.C. Production frontier and panel data. J. Bus. Econ. Stat. 1984, 2, 299–326. [Google Scholar]
  11. Macedo, P.; Silva, E. A stochastic production frontier model with a translog specification using the generalized maximum entropy estimator. Econ. Bull. 2010, 1, 587–596. [Google Scholar]
  12. Greene, W.H. Distinguishing between heterogeneity and inefficiency: Stochastic frontier analysis of the World Health Organization’s panel data on national health care systems. Health Econ. 2004, 13, 959–980. [Google Scholar] [CrossRef] [PubMed]
  13. Greene, W.H. Fixed and random effects in stochastic frontier models. J. Prod. Anal. 2005, 23, 7–32. [Google Scholar] [CrossRef]
  14. Terrell, D. Incorporating monotonicity and concavity conditions in flexible functional forms? J. Appl. Econ. 1996, 11, 179–194. [Google Scholar] [CrossRef]
  15. Battese, G.E.; Coelli, T. Frontier production functions, technical efficiency and panel data: With application to paddy farmers in India. J. Prod. Anal. 1992, 3, 153–169. [Google Scholar] [CrossRef]
  16. Bell, K.P.; Bockstael, N.E. Applying the generalized-moments estimation approach to spatial problems involving micro-level data. Rev. Econ. Statictics 2000, 82, 72–82. [Google Scholar] [CrossRef]
  17. Caves, D.W.; Christensen, L.R.; Diewert, W.E. The economic theory of index numbers and the measurement of input, output and productivity. Econometrica 1982, 50, 1393–1414. [Google Scholar] [CrossRef]
  18. Orea, L. Parametric decomposition of a generalized Malmquist index. J. Prod. Anal. 2002, 18, 5–22. [Google Scholar] [CrossRef]
  19. Mittelhammer, R.C.; Judge, G.G.; Douglas, J.M. Econometric Foundations, 1st ed.; Cambridge University Press: Cambridge, UK, 2000; p. 784. [Google Scholar]
  20. Perloff, J.M.; Karp, L.S.; Golan, A. Estimating Market Power and Strategies, 1st ed.; Cambridge University Press: Cambridge, UK, 2007; p. 352. [Google Scholar]
  21. Golan, A.; Perloff, J.M.; Shen, E.Z. Estimating a demand system with non-negativity constraints: Mexican meat demand. Rev. Econ. Statistics 2001, 83, 541–550. [Google Scholar] [CrossRef]
  22. Golan, A. Information theoretic approach for estimating nonlinear dynamic models. Stud. Nonlinear Dyn. E 2003, 7, 1–24. [Google Scholar] [CrossRef]
  23. Campbell, R.; Rogers, K.; Rezek, J. Efficient frontier estimation: A maximum entropy approach. J. Prod. Anal. 2008, 30, 213–221. [Google Scholar] [CrossRef]
  24. Pukelsheim, F. The three sigma rule. Am. Stat. 1994, 48, 88–91. [Google Scholar]

Share and Cite

MDPI and ACS Style

Tonini, A.; Pede, V. A Generalized Maximum Entropy Stochastic Frontier Measuring Productivity Accounting for Spatial Dependency. Entropy 2011, 13, 1916-1927. https://doi.org/10.3390/e13111916

AMA Style

Tonini A, Pede V. A Generalized Maximum Entropy Stochastic Frontier Measuring Productivity Accounting for Spatial Dependency. Entropy. 2011; 13(11):1916-1927. https://doi.org/10.3390/e13111916

Chicago/Turabian Style

Tonini, Axel, and Valerien Pede. 2011. "A Generalized Maximum Entropy Stochastic Frontier Measuring Productivity Accounting for Spatial Dependency" Entropy 13, no. 11: 1916-1927. https://doi.org/10.3390/e13111916

Article Metrics

Back to TopTop