Next Article in Journal
Regional Inequality in Underdeveloped Areas: A Case Study of Guizhou Province in China
Next Article in Special Issue
Using Goal-Programming to Model the Effect of Stakeholder Determined Policy and Industry Changes on the Future Management of and Ecosystem Services Provision by Ireland’s Western Peatland Forests
Previous Article in Journal
De-[Constructing] Growth
Previous Article in Special Issue
Assessing Urban Forest Structure, Ecosystem Services, and Economic Benefits on Vacant Land
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying Data Envelopment Analysis and Grey Model for the Productivity Evaluation of Vietnamese Agroforestry Industry

1
Department of Industrial Engineering and Management, Fortune University, Kaohsiung 83158, Taiwan
2
Department of Industrial Engineering and Management, National Kaohsiung University of Applied Sciences, Kaohsiung 80778, Taiwan
3
Department of Supply Chain Management, National Kaohsiung Marine University, Kaohsiung 80778, Taiwan
4
Department of Management Information Systems, National Chengchi University, Taipei 11605, Taiwan
*
Authors to whom correspondence should be addressed.
Sustainability 2016, 8(11), 1139; https://doi.org/10.3390/su8111139
Submission received: 12 June 2016 / Revised: 10 October 2016 / Accepted: 1 November 2016 / Published: 5 November 2016
(This article belongs to the Special Issue Decision Support for Forest Ecosystem Management Sustainability)

Abstract

:
Agriculture and forestry play important roles in Vietnam, particularly as they contribute to the creation of food, conservation of forest resources, and improvement of soil fertility. Therefore, understanding the performances of relevant enterprises in this field contributes to the sustainable development of this country’s agroforestry industry. This research proposes a hybrid model, which includes a grey model (GM) and a Malmquist productivity index (MPI), to assess the performances of Vietnamese agroforestry enterprises over several time periods. After collecting the data of selected input and output variables for 10 Vietnam agroforestry enterprises in the period of 2011–2014, GM is used to forecast the future values of these input and output variables for the 10 agroforestry enterprises in 2015 and 2016. Following the results of GM, the MPI is used to measure the performance of these enterprises. The MPI scores showed some enterprises will become more efficient, while others will become less efficient. The proposed model gives past–present–future insights in order for decision-makers to sustain agroforestry development in Vietnam. This hybrid approach can be applied to performance analysis of other industries as well.

1. Introduction

Agriculture and forestry play essential roles in a country, contributing to the creation of food, conservation of forest resources, and improvement of soil fertility. Agroforestry is a land-use system and technology [1]. Agroforestry systems can be advantageous over conventional agricultural and forest production methods. Typically, it increases productivity, economic benefits, and diversity in ecological goods and services. An agroforestry system includes the following features: (1) two or more species (possibly both trees and animals), which must be perennial; (2) at least two or more products in the system; (3) a production cycle of more than one year; (4) ecological and economical diversity; and (5) the need for perennials and other compatible component relationships [2]. Every country has unique natural and ecological conditions that lead it to create different ways for developing an agroforestry industry to satisfy societal needs. For example, Thailand developed methods for holding water agroforestry, retaining moisture, and improving the ecological environment and people’s living standards [3]. India implemented a famous “green revolution” [4] to produce enough food not only to overcome hunger but also to support and increase extra rice exports.
Agroforestry studies in Vietnam began about three decades ago. To sustain the use of land and secure food sources, many programs have been initiated. One of them involved cooperation with Finland, which introduced various agroforestry systems into the areas around the Ban Kan province for the sustainable development of vast uplands [5]. All relevant programs have adapted to the local climate conditions, capacity of local peoples, and market availability. In 2003, Vietnam implemented a project called “Capacity Building of Agriculture and Forestry in Vietnam”, which aimed to establish projects in Vietnam and Southeast Asia to support the development of alternatives to slash and burn, and to strengthening research capacity and agricultural and forestry education. In this project, for example, Vietnam developed an amplification system to replace the slash and burn system.
In short, agroforestry in Vietnam is a traditional industry, and its importance has been steadily increasing. However, according to statistical data, the export value of agroforestry products in the first quarter in 2014 declined by approximately 13.2%. Coffee exports dropped sharply during the first quarter of this year. Data showed that around 350,000 tons of coffee were exported in the period, worth $734 million, decreasing by 41.4% in volume and 37.3% in value compared with last year. Rice exports slipped by 28.1% and 32% in export volume and value, respectively. In March, 517,000 tons of rice were exported, taking the total rice export volume in the first quarter to 1.01 million tons, with a total value of $440 million [6]. However, some products happened to increase in export value despite a drop in export volume, such as with tea, cashews, and pepper. The agricultural and forestry sectors needed to enhance the quality of products and build brands to become more competitive in the world market. Vietnam has the benefit of exporting agricultural and forestry products, which reached $30.8 billion last year. The Vietnamese annual growth rate is dedicated to accelerating the growth momentum in the coming years [7]. Therefore, new discoveries and a more efficient system need to be emphasized. Agroforestry has had a wide impact on the Vietnamese community in the past centuries, and has been adopted in Vietnamese mainstream forestry and agricultural development.
The sustainable development of the agroforestry industry in Vietnam, and other countries, requires an effective approach to change productivity performance. Moreover, from the views of decision-makers and stakeholders in agroforestry, good productivity is an important issue for them to deploy their operational strategy and to judge their investment. The performance index is not just counted for a single time period; they should be combined with several time periods for catching development tends. Moreover, people not only want to know the historical performance but want to find potential good candidates for the future [8]. Therefore, this research proposes a hybrid approach, which combines a grey model (GM) with a Malmquist productivity index (MPI), is used to assess and foresee “past-current-future” performances of the agroforestry enterprises in Vietnam. Specifically, based on the historical data collected from the stock market from 2011 to 2014, the GM was first used to forecast the values of selected input and output variables for Vietnam agroforestry enterprises in 2015 and 2016. Then, the MPI was used to measure efficiency changes and productivity changes for these selected enterprises from 2011 to 2016. The proposed method provides insight for decision-makers to adjust their strategies. For government, the proposed method also guides policy directions toward sustainable development of the agroforestry industry in Vietnam. For global investors or stakeholders, the proposed method provides a channel to obtain performance information about an enterprise, as well [9]. Moreover, the proposed method may be applicable to other industries for performance evaluation.
The remainder of this paper is organized as follows: Section 2 is a literature review; Section 3 presents the research procedure; the methodology introduction is depicted in Section 4; Section 5 conducts an empirical study and results analysis; and conclusions and suggestions for future directions are presented in Section 6.

2. Literature Review

Data envelopment analysis (DEA) is a useful method for estimating production frontiers [10]. It has been applied to various application domains, including operations research, management science, economics, etc. In essence, DEA is a non-parametric data analytical technique whose domain of inquiry is a set of decision-making units (DMUs), which can receive multiple inputs and express multiple outputs. Given a set of DMUs, the DEA can establish the relative efficiencies of each DMU within this set. When applying DEA, only limited data is required to measure performance. The selection of input and output variables is important for DEA as they can affect decision-making. One prerequisite for applying DEA is that the selected input and output variables should keep an isotonic relationship, which can be validated by correlation analysis [11]. If the correlation between input and output variables are positive, this means the variables maintain an isotonic relationship and can be used by the DEA model. Otherwise, it needs to re-examine these variables. In past studies, Liu [12] used the number of employees, assets, and purchased funds as input variables, and demand deposits, short-term loans, and long-term loans as output variables. Romano and Guirrini [13] used the cost of labor, cost of material, cost of service, and cost of lease as inputs, and the population served and water delivered as outputs. The input and output variables used in these studies are referred to in this present research. In addition to comparing the relative performance of DMUs at a specific period, DEA can also be used to calculate the productivity changes of a DMU over time to examine the relative progress among competitors. In the past two decades, increasing publications have adopted MPI as the evaluation techniques to assess the efficiency of an enterprise. For example, Briec et al. [14] employed MPI to analyze the productivity growth and technological change of a Portuguese hydroelectric plant in 2001–2008. Yang et al. [15] thought forestry enterprises were important elements of the forestry economy in China. The increasing investment of the government in science and technology accelerates the development of the forestry economy in China. Qazi and Yulin [16] had to use MPI to measure changes in China’s high-tech industry’s productivity in 2000 with fifteen high-tech industries. It concluded that the electronic component and office equipment industries are considered to be valid.
Grey system theory was introduced by Deng [17]. It has become popular due to the advantage of managing a system with unknown parameters. Compared with traditional statistical methods, grey system theory requires little data for forecasting. Under the missing information or partially unknown parameters or uncertainty problems, this theory has become a popular method. Superior to conventional statistical models, grey forecasting can only use a limited amount of data to evaluate the unknown systems [18]. The GM (1,1) is one of the most popular grey forecasting models. Ren proved that using GM (1,N) to predict the yield of bio-hydrogen under scanty data conditions could give a better predictability result than the use of an artificial neural network [19]. There are some studies to have employed grey system theory. Kuo et al. [20] applied grey relational analysis for multiple attributes decision making. These studies indicate that grey system theory as been employed in various application areas. Lin et al. [21] proposed a gray forecast model to deal with the factor analysis of multivariate series forecasting problems. The results show that the model is better than other existing models.
Chen and Chen [22] used DEA and MPI to explore Taiwanese chip manufacturing company operating performance. The results show that if the wafer manufacturing company in Taiwan wished to improve their operating performance, they should improve their efficiency of constant returns to scale (CRS) and variable returns to scale (VRS), as there is no economy of scale. Wang et al. evaluated the performance of the Indian energy industry under multiple different input and output criteria. The DEA and grey theory are used to conduct this study [23].
Referring to the operational characteristics of Vietnamese agroforestry companies and summarizing the DEA literature mentioned in Section 2 [11,12,13,14,15,16,22,23], this research selects total assets, liabilities, and equity as three input factors because they are key financial indicators contributing to the performance of DMUs in the agroforestry industry. On the other hand, net revenue and gross profit are selected as output factors, as they are important indices for measuring the performance of DMUs in the agroforestry industry.

3. Research Procedure

In this research, a seven-step procedure is used, and the details of each step are described below and in Figure 1.
Step 1: Choose DMUs and collect their data: This step focuses on choosing DMUs in the Vietnam agroforestry industry and collecting their relevant information. The study investigates the relevant enterprises to find all potential candidates’ DMU lists. A total of 10 Vietnamese agroforestry companies with financial reports on the stock market during the period 2011–2014 are selected.
Step 2: Choose input/output variables: This step focuses on selecting the input and output variables. The selected input and output variables are critical, as they can affect the results. In this research, we refer to previous studies in order to make an appropriate choice. The Pearson correlation test will be used in Step 6 to check the suitability of these selected variables. However, according to the rule of thumb from Golany and Roll [11], the total DMUs need to be more or equal than double the number of inputs plus outputs.
Step 3: Grey prediction: In this step, the GM is used to predict the future data in 2015 and 2016. Since the prediction always exists with errors, the mean absolute percent error (MAPE) is used in the next step to check the prediction accuracy of the GM. Even though there are several prediction methodologies, almost all of them need to prepare lots of historical data. The GM is the only one that uses a minimum of four historical data types for prediction. The historical data are difficult to be collected in the Vietnamese agroforestry industry. Therefore, this study uses GM as a prediction method for overcoming a shortage of data.
Step 4: Check forecasting accuracy: As stated in the previous step, the MAPE is employed to check the prediction accuracy of the GM (1,1). Lewis [24] gave MAPE an assessment standard: a MAPE value < 10% is considered “good”; a MAPE value between 10%–20% is considered “qualified”; a MAPE value between 20%–50% is considered “just”; and a MAPE value > 50% is considered “unqualified”. If the forecasting error is too high, we have to go back to the Step 2 to reselect the input and output variables.
Step 5: Pearson correlation test: DEA requires the selected input and output variables to have an isotonic relationship. Therefore, to test that the data matches this isotonic prerequisite, the variable correlation analysis is calculated to verify a positive correlation between the selected inputs and outputs. If there are negative coefficient variables, the input and output variables need to be changed. Thus, we have to revert to the Step 2, until this prerequisite is met. In this research, we will employ the Pearson correlation test for this purpose.
Step 6: Choose the DEA model: This step includes many DEA models. Most of them can only analyze a single time period. However, only MPI can combine several time periods for integrated analysis and can separately discuss the changes of efficiency and technology. These characteristics match our data of several time periods and match our requirement of integrated analysis. Therefore, this study chooses a MPI model to assess and rank the efficiencies of DMUs for our analysis.
Step 7: Performance analysis: In this step, the forecasting performances of DMUs will be thoroughly analyzed.

4. Methodology

4.1. Grey Forecasting Model

Due to its popularity and computational efficiency, the GM (1,1) is used in this research for forecasting. The GM (1,1) is a type of time-series forecasting model with multiple inputs and outputs. One advantage of the GM (1,1) is that it does not need an entire set of historical data, but at least a set of four series data is required [18]. These series data have to be taken at equal intervals and bypassing any data in a sequential order is not allowed. The procedure of applying the GM (1,1) to forecasting is detailed as follows [25].
From the historical data we can derive the original time series dataset X ( 0 ) , as shown in Equation (1):
X ( 0 ) = ( X ( 0 ) ( 1 ) , X ( 0 ) ( 2 ) , ... , X ( 0 ) ( n ) ) ,   n 4
where X ( 0 ) is a non-negative sequence, and n is the total number of data observations.
From X ( 0 ) , we can derive the accumulated time series data X ( 1 ) , as shown in Equation (2):
X ( 1 ) = ( X ( 1 ) ( 1 ) , X ( 1 ) ( 2 ) , ... , X ( 1 ) ( n ) ) ,   n 4
where:
X ( 1 ) ( 1 ) = X ( 0 ) ( 1 )   if   k = 1
X ( 1 ) ( k ) = i = 1 k X ( 0 ) ( i )   if   k = 2 , , n
Then, the mean sequence Z ( 1 ) of X ( 1 ) can be derived as shown in Equation (5):
Z ( 1 ) = ( Z ( 1 ) ( 2 ) , Z ( 1 ) ( 3 ) , , Z ( 1 ) ( n ) )
where Z(1)(k) is the mean value of two adjacent data, i.e.,:
Z ( 1 ) ( k ) = 1 2 ( X ( 1 ) ( k ) + X ( 1 ) ( k 1 ) ) ,   k = 2 , 3 , ... , n
From X ( 1 ) , a GM (1,1) model which corresponds to the first order differential equation X ( 1 ) ( k ) (Equation (7)) can be constructed as follows:
d X ( 1 ) ( k ) d k + a X ( 1 ) ( k ) = b
where parameters a and b are called the developing coefficient and grey input, respectively.
In practice, the parameters a and b are not directly calculated from Equation (7); they are actually derived from Equation (8):
[ a , b ] T = ( B T B ) 1 B T Y
where Y = [ X ( 0 ) ( 2 ) X ( 0 ) ( 3 ) ... ... X ( 0 ) ( n ) ] and B = [ Z ( 1 ) ( 2 ) 1 Z ( 1 ) ( 3 ) 1 ... ... Z ( 1 ) ( n ) 1 ] .
The vector Y is called the data series, vector B is called the data matrix, and [ a , b ] T is called the parameter series.
After substituting the derived values of a and b into Equation (9), this resulting equation can then be used to generate the accumulated time series values of prediction, which can be organized into a set of values as shown in Equation (10).
X ^ ( 1 ) ( k + 1 ) = [ X ( 0 ) ( 1 ) b a ] e a k + b a
where X ^ ( 1 ) ( k + 1 ) denotes the accumulated value of prediction of variable X at the time point k + 1.
X ^ ( 1 ) = ( X ^ ( 1 ) ( 1 ) , X ^ ( 1 ) ( 2 ) , ... , X ^ ( 1 ) ( n ) )
where X ^ ( 1 ) ( 1 ) = X ( 0 ) ( 1 ) .
Substituting the k value (k = 1, …, n) into Equation (11) we can then find the set of forecasting values as shown in Equation (12):
X ^ ( 0 ) ( k + 1 ) = [ X ( 0 ) ( 1 ) b a ] e a k ( 1 e a )
X ^ ( 0 ) = ( X ^ ( 0 ) ( 1 ) , X ^ ( 0 ) ( 2 ) , ... , X ^ ( 0 ) ( n ) )
Errors always exist in forecasting, which can be classified into bias or random types. A bias error is a consistent result of mistakes—predicting too high or low. Such errors are often demand patterns being inaccurately estimated, such as trend, seasonal, or cyclic pattern results. Random errors are unpredictable factors. To understand the accuracy of forecasting models one needs to check the forecast errors. For this purpose the measurement, MAPE, is used. The equation of MAPE is defined in Equation (13), where n is the total number of data observations. A smaller MAPE means a higher forecasting accuracy.
MAPE = | A c t u a l t F o r e c a s t t | A c t u a l t × 100 n
Lewis [24] gave an assessment standard for MAPE as follows: a MAPE value < 10% is considered “good”; a MAPE value between 10%–20% is considered “qualified”; a MAPE value between 20%–50% is considered “just”; a MAPE value > 50% is considered “unqualified”.

4.2. Malmquist Productivity Index

Caves et al. developed MPI in 1982 and the major DEA model for productivity measurement is MPI [26]. The MPI calculated the relative performance of a DMU in different periods of time using the technology of the base period [27]. The MPI is used to evaluate the change in efficiency of a DMU between two time periods [28]:
MPI = [ d o s ( x t , y t ) d o s ( x s , y s ) × d o t ( x t , y t ) d o t ( x s , y s ) ] 1 / 2
where d o s is a distance function measuring the efficiency of the conversion of inputs xs to outputs ys in the period s (note that DEA efficiency is considered a distance measure in the literature as it reflects the efficiency of converting inputs to outputs).
Importantly, if there is a technological change in period t, then, d o t ( x s , y s ) = Efficiency   of   conversion   of   input   in   period   s   to   output   in   period   s d o s ( x s , y s ) .
MPI considers two-period efficiency and technological change in the geometric mean. MPI was defined in Equation (15) and was being used to estimate the change of overall productivity of a pharmaceutical company over time. If MPI > 1 it means overall productivity increases; if MPI = 1 it indicates no change; a MPI < 1, it indicates a decrease of overall productivity:
MPI = d o t ( x t , y t ) d o s ( x s , y s ) [ d o s ( x s , y s ) d o t ( x s , y s ) × d o s ( x t , y t ) d o s ( x t , y t ) ] 1 / 2 = Efficiency   change × Technological   change
The “efficiency change” is also called the “catch-up” effect, indicating the extent that a DMU improves or deteriorates its efficiency. If the change in efficiency > 1 then it represents relative efficiency progress from time period s to t; if change in efficiency = 1 or < 1 it indicates there is no change or recession, respectively. The “technological change” is also called the “frontier shift effect” (or innovative effect), reflecting the frontier change between the two periods. If technological change > 1 it represents technological progression; technological change < 1 shows technological recession.

5. Empirical Results

Following the research procedure proposed in Section 3, Table 1 shows 10 Vietnamese agroforestry companies selected for this research. These companies are the top companies in the agroforestry industry in Vietnam, and their data are posted publicly and prestigiously. These companies are qualified with transparent financial data. Information about these DMUs was collected from the market observation posting system of Vietstock.vn, which is a premier site for business and financial market news in Vietnam.
Referring to the literature review, this research selects total assets, liabilities, and equity as the three inputs factors, and net revenue and gross profit as the output factors. The total has 10 DMUs with three inputs and two outputs. Matching the rule of thumb from Golany and Roll [11], which states that the total number of DMUs need to be more or equal than double the number of inputs plus outputs. The correlation of these selected input and output variables will be validated by the Pearson correlation test.
Having determined these input and output variables, we collect the historical data of these DMUs from 2011 to 2014 [6,29]. However, due to space limitations, only the data of the year 2014 are listed in Table 2.
Following the research procedure, the GM (1,1) is used in this section to predict the future data from 2015 and 2016 for the selected DMUs. Table 3 shows the DMU1’s historical data (2011–2014) used by the GM (1,1).
Below, we use the “total asset data” of DMU1 in Table 3 to illustrate the generation of the forecasting data of DMU1 step-by-step.
(1)
Create the primitive series
X ( 0 ) = ( 2 , 430 , 078 ,   2 , 821 , 172 ,   3 , 328 , 252 ,   3 , 312 , 062 )
(2)
Generate the accumulated series X ( 1 )
X ( 1 ) = ( 2 , 430 , 078 ,   5 , 251 , 250 ,   8 , 579 , 502 ,   11 , 891 , 564 )
Where
X ( 1 ) ( 1 ) = X ( 0 ) ( 1 ) = 2 , 430 , 078 X ( 1 ) ( 2 ) = X ( 0 ) ( 1 ) + X ( 0 ) ( 2 ) = 5 , 251 , 250 X ( 1 ) ( 3 ) = X ( 0 ) ( 1 ) + X ( 0 ) ( 2 ) + X ( 0 ) ( 3 ) = 8 , 579 , 502 X ( 1 ) ( 4 ) = X ( 0 ) ( 1 ) + X ( 0 ) ( 2 ) + X ( 0 ) ( 3 ) + X ( 0 ) ( 4 ) = 11 , 891 , 564
(3)
Create mean series dataset Z ( 1 )
To find the mean series dataset Z ( 1 ) , Equation (6) is used and we derive the following data.
X ( 1 ) ( 1 ) = X ( 0 ) ( 1 ) = 2 , 430 , 078 Z ( 1 ) ( 2 ) = 1 2 ( 2 , 430 , 078 + 5 , 251 , 250 ) = 3 , 840 , 664 Z ( 1 ) ( 3 ) = 1 2 ( 5 , 251 , 250 + 8 , 579 , 502 ) = 6 , 915 , 376 Z ( 1 ) ( 4 ) = 1 2 ( 8 , 579 , 502 + 11 , 891 , 564 ) = 10 , 235 , 533
Then, Z ( 1 ) = ( Z ( 1 ) ( 2 ) , Z ( 1 ) ( 3 ) , Z ( 1 ) ( 4 ) ) = ( 3 , 840 , 660 ,   6 , 915 , 376 ,   10 , 235 , 533 ) .
(4)
Find the values for coefficients a and b
Let B = [ 3 , 840 , 664 6 , 915 , 376 10 , 235 , 533 1 1 1 ] , θ ^ = [ a b ] , y N = [ 2 , 821 , 172 3 , 328 , 252 3 , 312 , 062 ] .
Then, using the Equation (8), we can derived the values of a and b as below.
[ a b ] = θ ^ = ( B T B ) 1 B T y N = [ 0 . 0756791 2 , 624 , 287 . 78 ]
(5)
Generate the accumulated data series
Substitute the two coefficients a and b, as well as the k values (k = 0, …, 7), into the Equation (9); then we can derive the accumulated data series in the third column in Table 4.
(6)
Generate the series values of prediction
Substitute the two coefficients a and b, as well as the k values (k = 0, …, 7) into the Equation (11); then we can the derive the fifth column in Table 4.
With the same computational process, the study could obtain the forecasting data of all DMUs in 2015 and 2016, as shown in Table 5.
In this step we check the prediction accuracy of the GM (1,1) based on the MAPE values. When a MAPE value is small it means the predicted value is close to the actual value. The MAPE values obtained are shown in Table 6.
As the MAPE values obtained are, mostly, smaller than 10%, especially as the average MAPE of the 10 DMUs reaches 5.111% (below 10%, as well), it confirms that the GM (1,1) provides a good prediction accuracy in this research.
However, before using this DEA model, we have to ensure that input and output variables have isotonic relationships. This means that if the input quantity increases, then the output quantity could not decrease under the same condition. Thus, a Pearson correlation test is first used to ensure this prerequisite. A higher correlation coefficient refers to a closer relationship between two variables, and a low correlation coefficient means a low relationship. The interpretation of the correlation coefficient is explained in more detail as follows: The correlation coefficient is always between −1 and +1. When the correlation coefficient is closer to ± 1, this means two groups are closer to a perfect linear relationship. When a correlation value is less than 0.2 then the degree of correlation is considered “very low”; a value between 0.4–0.6 is considered “average high”; a value between 0.6–0.8 is considered “efficient high”; a value between 0.8–1 is considered “upper high”. Our empirical results (Table 7) indicate that these variables have strong positive correlations, indicating that these selected input and output factors are qualified.
Though some DEA models are available for assessing DMUs, the efficiencies of DMUs of several time periods cannot be compared for most of the DEA models [30]. In this research, the Malmquist model is used to evaluate the performances of DMUs. Table 8 shows the results derived from the Malmquist O-V model. To facilitate the analysis, the values of efficiency change, technological change, and MPI are depicted in Figure 2, Figure 3 and Figure 4, respectively.
Figure 2 shows the “efficiency changes” (catch-up effects) of the 10 DMUs over the time period 2011–2016. It is found that the efficiency changes of these DMUs fluctuated over these years, especially DMU3, DMU4, and DMU8. DMU3 had experienced a dramatic efficiency rise during the time period 2011–2015, and DMU4 experienced a dramatic efficiency drop during the time period 2015–2016, while the DMU8 had experienced a dramatic drop during the time period 2014–2015. Except for DMU2 and DMU6, which had experienced slight efficiency changes, the other DMUs all had experienced some degree of efficiency fluctuations. In terms of “efficiency change”, DMU9, DMU4, and DMU6 are the top three best companies, while the top three worst companies are DMU3, DMU7, and DMU1. DMU2 is found with a stable efficiency, but it also had not improved its efficiency in the current and near future. From Figure 2, we can determine and predict the efficiency changes, or the “catch-up” effects, of each DMU.
Figure 3 show the “technological changes” (frontier-shift) of the 10 DMUs from 2011–2016. It is found that, except for DMU5, other DMUs had experienced an upward technological change, though most of them had experienced an efficiency drop during 2013–2014. This implies most of these DMUs had continued to improve their technological capabilities. DMU5 always keeps a lower “technological changes” score, indicating that it has not actively improved its technology. Thus, DMU5 requires further investigation into the causes that lead to its tardiness on technology development. In terms of “technological change”, DMU10, DMU4, and DMU9 are the top three best companies, while the top three worst companies are DMU5, DMU1, and DMU8.
Figure 4 shows the MPI, the total productivity of the 10 DMUs from 2011–2016. It is noted that DMU6 and DMU9 had a long-term upward trend during 2011–2016, though DMU6 had experienced a drop during 2013–2014. In terms of MPI, DMU9, DMU4, and DMU10 are the top three best companies, while the top three worst companies are DMU5, DMU1, and DMU7. DMU8 appears unstable due to having the largest fluctuation; thus, it requires intense care for its development. In summary, to sustain the development of the agroforestry industry, the Vietnamese government needs to especially focus on DMU5, DMU1, DMU7, and DMU8.

6. Conclusions

For sustainable development of the Vietnamese agroforestry industry, this research has proposed a hybrid approach combining GM and DEA to assess the current, and predict the future, performance of agroforestry companies. This study first selects total assets, liabilities, and equity as input variables, and net revenue and gross profit as output variables. After collecting the historical data (2011–2014) from the Vietnam stock market, the GM (1,1) model is used to forecast the future data in 2015–2016 for selected agroforestry enterprises. The MAPE value of 5.11% confirms the desired accuracy of the GM (1,1). Based on the historical and predicted data, the Malmquist model is used to assess the performance of these DMUs. The derived results provide insight views into the agroforestry companies in terms of “efficiency changes”, “technological changes”, and MPI (total productivity). These results can also provide insight into the “past-current-future” performances of these DMUs. In summary, to sustain the development of Vietnamese agroforestry industry, the government should take more care in the developments of these companies, which are DMU5, DMU1, DMU7, and DMU8.
This research also demonstrates that the proposed approach can help decision-makers. Policy-makers can use the approach to develop strategies to sustain the development of the Vietnamese agroforestry industry. Investors can use of this approach to discover the good companies for investments. This mathematical approach reduces the errors and risks in decision-making. In addition, the approach can be used in other industries to extend its contributions. This study provides many significant and noticeable results. This numerical study gives us better “past-present-future” insights through the integration method, and this work could be used as a better model for performance analysis among the decision-makers of varied industries.
The proposed approach has been applied to the Vietnamese agroforestry industry; however, it only includes 10 companies listed on the stock market. Including more companies can provide further focus. In addition, some other input and output variables (such as the number of staff, number of branches, revenue, research and development, etc.) can be taken to measure the performance of these companies. Moreover, different DEA models can be further investigated for comparisons. Furthermore, this approach can be applied to other industries.

Acknowledgments

This research was partially supported by the Ministry of Science and Technology of Taiwan under the grant MOST 105-2221-E-151-039 and MOST M105-2622-E-151-006-CC3. The authors also appreciate the support from National Kaohsiung University of Applied Sciences in Taiwan and Minister of Industry and Trade in Vietnam.

Author Contributions

Chia-Nan Wang guided the research direction and found the solutions; Han-Sung Lin conducted the literature review and designed the experiments; Hsien-Pin Hsu summarized the data and edited the paper; Van-Tinh Le contributed to the experiment and analyzed the results. Tsung-Fu Lin worked on the writing and formatting. All authors contributed to produce results.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
DEAdata envelopment analysis
GMgrey model
MPIMalmquist productivity index
DMUdecision making unit
CCRCharnes, Cooper and Rhodes
BCCBanker, Charnes and Cooper
MAPEmean absolute percent error
CRSconstant returns to scale
VRSvariable returns to scale

References

  1. Torquebiau, E.F. A renewed perspective on agroforestry concepts and classification. C. R. Acad. Sci. 2000, 323, 11009–11017. [Google Scholar] [CrossRef]
  2. DeWit, R.; Leibreich, J.; Vernier, F.; Delmas, F.; Beuffe, H.; Maison, P.; Capdeville, P. Relationship between land-use in the agro-forestry system of les Landes, nitrogen loading to and risk of macro-algal blooming in the Bassin d’Arcachon coastal lagoon (SW France). Estuar. Coast. Shelf Sci. 2005, 62, 353–465. [Google Scholar]
  3. Hurni, H.; Nuntapong, S. Agro-forestry improvements for shifting cultivation systems: Soil conservation research in northern Thailand. Mount. Res. Dev. 1983, 3, 338–345. [Google Scholar] [CrossRef]
  4. Sharma, V.P.; Kumar, A. Factors influencing adoption of agro-forestry programme: A case study from North-West India. Indian J. Agric. Econ. 2000, 55, 300–510. [Google Scholar]
  5. Craswell, E.T.; Sajjapongse, A.; Howlett, D.J.B.; Dowling, A.J. Agroforestry in the management of sloping lands in Asia and the Pacific. Agrofor. Syst. 1997, 38, 121–137. [Google Scholar] [CrossRef]
  6. Annual Report, Vietnamese Minister of Industry and Trade. Available online: http://www.moit.gov.vn/en/Pages/default.aspx (accessed on 15 January 2016).
  7. Fritzen, S. Growth, inequality and the future of poverty reduction in Vietnam. J. Asian Econ. 2002, 13, 635–657. [Google Scholar] [CrossRef]
  8. Wang, C.N.; Nguyen, N.T.; Tran, T.T. Integrated DEA models and grey system theory to evaluate past-to-future performance: A case of Indian electricity industry. Sci. World J. 2015, 2015, 1–17. [Google Scholar] [CrossRef] [PubMed]
  9. Wang, C.N.; Nguyen, X.T.; Wag, Y.H. Automobile industry strategic alliance partner selection: The application of a hybrid DEA and Grey theory model. Sustainability 2016, 8, 273. [Google Scholar] [CrossRef]
  10. Charnes, A.; Cooper, W.W.; Rhodes, E. Measuring the efficiency of decision making units. Eur. J. Oper. Res. 1978, 2, 429–444. [Google Scholar] [CrossRef]
  11. Golany, B.; Roll, Y. An application procedure for DEA. Omega 1989, 17, 237–250. [Google Scholar] [CrossRef]
  12. Liu, S.T. Measuring and categorizing technical efficiency and productivity change of commercial banks in Taiwan. Expert Syst. Appl. 2010, 37, 2783–2789. [Google Scholar] [CrossRef]
  13. Romano, G.; Guerrini, A. Measuring and comparing the efficiency of water utility companies: A data envelopment analysis approach. Util. Policy 2011, 19, 202–209. [Google Scholar] [CrossRef]
  14. Briec, W.; Peypoch, N.; Ratsimbanierana, H. Productivity growth and biased technological change in hydroelectric dams. Energy Econ. 2011, 33, 853–858. [Google Scholar] [CrossRef]
  15. Yang, H.; Yuan, T.; Zhang, X.; Li, S. A decade trend of total factor productivity of key state-owned forestry enterprises in China. Forests 2016, 7, 57. [Google Scholar] [CrossRef]
  16. Qazi, A.Q.; Yulin, Z. Productivity measurement of hi-tech industry of China malmquist productivity index–DEA approach. Procedia Econ. Financ. 2012, 1, 330–336. [Google Scholar] [CrossRef]
  17. Deng, J.L. Control problems of grey systems. Syst. Control Lett. 1982, 1, 288–294. [Google Scholar]
  18. Deng, J.L. Introduction to grey system theory. J. Grey Syst. 1989, 1, 1–24. [Google Scholar]
  19. Ren, J.; Gao, S.; Tan, S.; Dong, L. Prediction of the yield of biohydrogen under scanty data conditions based on GM (1,N). Int. J. Hydrog. Energy 2013, 38, 13198–13203. [Google Scholar] [CrossRef]
  20. Kuo, Y.; Yang, T.; Huang, G.W. The use of grey relational analysis in solving multiple attribute decision-making problems. Comput. Ind. Eng. 2008, 55, 80–93. [Google Scholar] [CrossRef]
  21. Lin, Y.H.; Wang, J.S.; Pai, P.F. A grey prediction model with factor analysis technique. J. Chin. Inst. Ind. Eng. 2004, 21, 535–542. [Google Scholar] [CrossRef]
  22. Chen, Y.S.; Chen, B.Y. Applying DEA, MPI, and grey model to explore the operation performance of the Taiwanese wafer fabrication industry. Technol. Forecast. Soc. Chang. 2011, 78, 536–546. [Google Scholar] [CrossRef]
  23. Wang, C.N.; Lin, L.C.; Murugesan, D. Analyzing PSU’s performance: A case from ministry of petroleum and natural gas of India. Math. Probl. Eng. 2013, 2013, 802690. [Google Scholar] [CrossRef]
  24. Lewis, C.D. International and Business Forecasting Methods; Butterworths: London, UK, 1982. [Google Scholar]
  25. Liu, S.; Forrest, J.Y.L. Grey Systems: Theory and Applications; Springer: Chennai, India, 2010. [Google Scholar]
  26. 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]
  27. Färe, R.; Grosskopf, S.; Norris, M.; Zhang, Z. Productivity growth, technical progress, and efficiency change in industrialized countries. Am. Econ. Rev. 1994, 84, 66–83. [Google Scholar]
  28. Chen, Y. A non-radial Malmquist productivity index with an illustrative application to Chinese major industries. Int. J. Prod. Econ. 2003, 83, 27–35. [Google Scholar] [CrossRef]
  29. Vietstock. Available online: http://vietstock.vn/ (accessed on 15 January 2016).
  30. Banker, R.D.; Charnes, A.; Cooper, W.W. Some models for estimating technological and scale inefficiencies in data envelopment analysis. Manag. Sci. 1984, 30, 1078–1092. [Google Scholar] [CrossRef]
Figure 1. Research development flow.
Figure 1. Research development flow.
Sustainability 08 01139 g001
Figure 2. Efficiency change or catch-up (2011–2016).
Figure 2. Efficiency change or catch-up (2011–2016).
Sustainability 08 01139 g002
Figure 3. Technological change or Frontier shift (2011–2016).
Figure 3. Technological change or Frontier shift (2011–2016).
Sustainability 08 01139 g003
Figure 4. MPI scores (2011–2016).
Figure 4. MPI scores (2011–2016).
Sustainability 08 01139 g004
Table 1. The list of 10 DMUs selected from Vietnam agroforestry industry.
Table 1. The list of 10 DMUs selected from Vietnam agroforestry industry.
Number OrderDMUCompany NameStock Market Code
1DMU1Dong Phu Rubber Joint Stock CompanyHOSE: DPR
2DMU2HAGL Joint Stock CompanyHOSE: HAG
3DMU3Hoa Binh Rubber Joint Stock CompanyHOSE: HRC
4DMU4National Seed Joint Stock CompanyHOSE: NSC
5DMU5Phuoc Hoa Rubber Joint Stock CompanyHOSE: PHR
6DMU6Southern Seed CorporationHOSE: SSC
7DMU7Thong Nhat Rubber Joint Stock CompanyHOSE: TNC
8DMU8Tay Ninh Rubber Joint Stock CompanyHOSE: TRC
9DMU9Quang Nam Rubber Investment JSCHOSE: VHG
10DMU10Vinacafe Son Thanh Joint Stock CompanyOTC: AUM
Table 2. Inputs and outputs data of DMUs in 2014.
Table 2. Inputs and outputs data of DMUs in 2014.
DMUsInputs ( Millions of VND)Outputs (Millions of VND)
(I) Total Asset(I) Liabilities(I) Equity(O) Net Revenue(O) Gross Profit
DMU13,312,062658,0392,195,471938,682285,480
DMU236,368,86420,978,62414,237,7283,054,3311,226,993
DMU3686,336181,918504,418161,3947059
DMU4953,008202,337744,946719,136287,323
DMU53,428,8301,129,5282,253,0821,604,981275,288
DMU6537,355162,978373,160608,625190,717
DMU7332,60530,901301,70463,873−355
DMU81,587,255140,3081,446,947500,638121,259
DMU91,071,970126,832900,392394,87346,077
DMU1013,283267510,60821,611158
Table 3. The historical data (2011–2014) of the input and output variables of DMU1.
Table 3. The historical data (2011–2014) of the input and output variables of DMU1.
DMU1Inputs (Millions of VND)Outputs (Millions of VND)
(I) Total Asset(I) Liabilities(I) Equity(O) Net Revenue(O) Gross Profit
20112,430,078491,0561,910,1131,837,202866,224
20122,821,172556,9602,179,9901,376,951553,307
20133,328,252606,6002,298,5991,100,122424,953
20143,312,062658,0392,195,471938,682285,480
Table 4. The derived prediction values X(k).
Table 4. The derived prediction values X(k).
kX(k)ValueX(k)Value
k = 0X(1)2,430,078X(1)2,430,078
k = 1X(2)5,347,264.63X(2)2,917,186.63
k = 2X(3)8,493,789.83X(3)3,146,525.2
k = 3X(4)11,887,683.4X(4)3,393,893.54
k = 4X(5)15,548,392.5X(5)3,660,709.08
k = 5X(6)19,496,893.1X(6)3,948,500.68
k = 6X(7)23,755,810.6X(7)4,258,917.41
k = 7X(8)28,349,548.5X(8)4,593,737.97
Table 5. Predicted input and output data of all DMUs in 2015 and 2016.
Table 5. Predicted input and output data of all DMUs in 2015 and 2016.
YearDMUsInputs (Currency Unit: Millions of VND)Outputs (Currency Unit: Millions of VND)
(I) Total Asset(I) Liabilities(I) Equity(O) Net Revenue(O) Gross Profit
2015DMU13,660,709.08715,051.352,239,955.23759,042.36216,764.72
DMU238,126,475.919,810,053.217,301,061.42,169,114.511,234,682.16
DMU3710,528.52192,011.93518,491.25140,856.832963.25
DMU41,427,295.43256,101.601,160,150.10816,038.05343,084.35
DMU53,646,880.391,210,168.952,382,125.121,369,690.16209,289.74
DMU6601,752.38169,376.87436,089.19715,599.93213,766.23
DMU7308,391.5821,024.41288,886.8049,637.401344.39
DMU81,765,695.56259,012.541,529,264.81370,704.2480,840.40
DMU91,534,558.85108,114.201,366,325.23584,979.7143,400.28
DMU1013,221.752551.0810,749.1524,795.80177.04
2016DMU13,948,500.68777,109.402,247,634.11624,477.72158,370.95
DMU241,364,890.2320,100,095.7620,670,001.671,748,357.561,248,329.47
DMU3725,873.82196,885.59528,955.1491,601.171261.76
DMU42,418,284.18349,367.532,116,736.15939,274.77412,839.94
DMU53,829,240.641,250,205.832,512,384.371,167,365.28142,525.83
DMU6678,367.07173,410.58517,540.57813,804.44242,919.54
DMU7288,308.5915,960.23276,386.2531,491.09398.52
DMU81,816,554.47248,283.061,606,736.41275,684.5352,543.82
DMU92,386,716.78101,837.682,318,982.81971,069.14138,457.14
DMU1012,975.192287.6310,853.2529,651.57109.37
Table 6. Average MAPE errors of DMUs.
Table 6. Average MAPE errors of DMUs.
DMUsMAPEDMUsMAPE
DMU11.626%DMU67.792%
DMU25.663%DMU72.069%
DMU32.256%DMU81.124%
DMU44.001%DMU914.513%
DMU55.021%DMU107.046%
Average MAPE: 5.111%
Table 7. The correlations between input and output variables in 2011–2016.
Table 7. The correlations between input and output variables in 2011–2016.
2011Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9979510.9933870.7551610.781175
Liabilities0.99795110.984010.7191370.744403
Equity0.9933870.9840110.813370.840045
Net Revenue0.7551610.7191370.8133710.99058
Gross profit0.7811750.7444030.8400450.990581
2012Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9975930.9877210.9169870.865799
Liabilities0.99759310.974520.8922730.832981
Equity0.9877210.9745210.9575690.924983
Net Revenue0.9169870.8922730.95756910.983508
Gross profit0.8657990.8329810.9249830.9835081
2013Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9972830.9970840.8472840.927948
Liabilities0.99728310.9889450.8210530.90642
Equity0.9970840.98894510.8741790.946157
Net Revenue0.8472840.8210530.87417910.967438
Gross profit0.9279480.906420.9461570.9674381
2014Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9981390.996820.9015230.964055
Liabilities0.99813910.9903010.883170.955216
Equity0.996820.99030110.9232120.970017
Net Revenue0.9015230.883170.92321210.956117
Gross profit0.9640550.9552160.9700170.9561171
2015Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9981820.9985550.8291650.960343
Liabilities0.99818210.9939880.8101350.955803
Equity0.9985550.99398810.8469130.961551
Net Revenue0.8291650.8101350.84691310.895381
Gross profit0.9603430.9558030.9615510.8953811
2016Total AssetLiabilitiesEquityNet RevenueGross Profit
Total asset10.9977540.9981230.7225420.95035
Liabilities0.99775410.9928050.6939720.94331
Equity0.9981230.99280510.7491760.955682
Net Revenue0.7225420.6939720.74917610.818206
Gross profit0.950350.943310.9556820.8182061
Table 8. The catch up, frontier, and MPI of DMUs from 2011 to 2016.
Table 8. The catch up, frontier, and MPI of DMUs from 2011 to 2016.
Catch-up
DMUs2011–20122012–20132013–20142014–20152015–2016Average
DMU10.8290641.0175340.4811220.4802250.5886110.679311
DMU2111111
DMU30.5074930.7501960.491170.7696740.5609380.615894
DMU41.500321.5117641.2192311.6392570.9246591.359046
DMU50.7343991.2958570.8783221.1304130.690980.945994
DMU60.6542291.362341.364081.3530881.361181.218983
DMU70.4648720.9381840.5228820.8275530.4314610.63699
DMU81.3111970.2517781.5193960.4531930.6940190.845917
DMU90.6424121.1250741.685042.3130332.0639881.565909
DMU101.1463580.8600331.0234630.8978240.9353410.972604
Frontier
DMUs2011–20122012–20132013–20142014–20152015–2016Average
DMU10.2942750.4921470.5728741.0425271.1664150.713647
DMU211.0468660.925847110.994543
DMU31.0693731.0169910.8256231.080881.0724511.013064
DMU40.8838120.9635770.7765331.2869131.6033981.102847
DMU50.4023610.3589680.5213820.3927110.4558770.42626
DMU61.1076211.0015660.740051.0239851.104370.995519
DMU71.0377611.0673160.6433541.1566751.3341361.047848
DMU80.7304530.9547310.7918981.0934951.059640.926043
DMU91.1157981.0247820.7963851.0723491.3437071.070604
DMU101.0637731.0492611.1782191.1818821.2233641.1393
MPI
DMUs2011–20122012–20132013–20142014–20152015–2016Average
DMU10.2439730.5007770.2756220.5006480.6865640.441517
DMU211.0468660.925847110.994543
DMU30.5426990.7629430.4055210.8319250.6015790.628934
DMU41.3260011.4567010.9467732.1095811.4825961.464331
DMU50.2954940.4651710.4579410.4439250.3150020.395506
DMU60.7246371.3644731.0094881.3855421.5032461.197477
DMU70.4824261.0013380.3363980.9572090.5756280.6706
DMU80.9577680.240381.2032070.4955640.735410.726466
DMU90.7168021.1529551.341942.4803792.7733951.693094
DMU101.2194640.9023991.2058641.0611221.1442621.106622

Share and Cite

MDPI and ACS Style

Wang, C.-N.; Lin, H.-S.; Hsu, H.-P.; Le, V.-T.; Lin, T.-F. Applying Data Envelopment Analysis and Grey Model for the Productivity Evaluation of Vietnamese Agroforestry Industry. Sustainability 2016, 8, 1139. https://doi.org/10.3390/su8111139

AMA Style

Wang C-N, Lin H-S, Hsu H-P, Le V-T, Lin T-F. Applying Data Envelopment Analysis and Grey Model for the Productivity Evaluation of Vietnamese Agroforestry Industry. Sustainability. 2016; 8(11):1139. https://doi.org/10.3390/su8111139

Chicago/Turabian Style

Wang, Chia-Nan, Han-Sung Lin, Hsien-Pin Hsu, Van-Tinh Le, and Tsung-Fu Lin. 2016. "Applying Data Envelopment Analysis and Grey Model for the Productivity Evaluation of Vietnamese Agroforestry Industry" Sustainability 8, no. 11: 1139. https://doi.org/10.3390/su8111139

APA Style

Wang, C. -N., Lin, H. -S., Hsu, H. -P., Le, V. -T., & Lin, T. -F. (2016). Applying Data Envelopment Analysis and Grey Model for the Productivity Evaluation of Vietnamese Agroforestry Industry. Sustainability, 8(11), 1139. https://doi.org/10.3390/su8111139

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop