Next Article in Journal
Effects of Exogenous Trans-Zeatin and Lovastatin on Abortion of Small Seeds in ‘Dawuxing’ Loquat (Eriobotrya Japonica Lindl.)
Next Article in Special Issue
Estimation of Frost Hazard for Tea Tree in Zhejiang Province Based on Machine Learning
Previous Article in Journal
Sugar Beet Root Yield and Quality with Leaf Seasonal Dynamics in Relation to Planting Densities and Nitrogen Fertilization
Previous Article in Special Issue
German Farmers’ Attitudes on Adopting Autonomous Field Robots: An Empirical Survey

Prediction of Food Production Using Machine Learning Algorithms of Multilayer Perceptron and ANFIS

Doctoral School of Economic and Regional Sciences, Hungarian University of Agriculture and Life Sciences, 2100 Godollo, Hungary
Department of Biosystem Engineering, University of Mohaghegh Ardabili, Ardabil 5619911367, Iran
Institute of Economic Sciences, Hungarian University of Agriculture and Life Sciences, 2100 Godollo, Hungary
Institute of Information Society, University of Public Service, 1083 Budapest, Hungary
Faculty of Civil Engineering, Technische Universitat Dresden, 01069 Dresden, Germany
John von Neumann Faculty of Informatics, Obuda University, 1034 Budapest, Hungary
Information Systems, University of Siegen, 57072 Siegen, Germany
Authors to whom correspondence should be addressed.
Academic Editors: Selma Boumerdassi, Eric Renault and Christopher Robin Bryant
Received: 2 February 2021 / Revised: 26 April 2021 / Accepted: 29 April 2021 / Published: 2 May 2021
(This article belongs to the Special Issue Latest Advances for Smart and Sustainable Agriculture)


Advancing models for accurate estimation of food production is essential for policymaking and managing national plans of action for food security. This research proposes two machine learning models for the prediction of food production. The adaptive network-based fuzzy inference system (ANFIS) and multilayer perceptron (MLP) methods are used to advance the prediction models. In the present study, two variables of livestock production and agricultural production were considered as the source of food production. Three variables were used to evaluate livestock production, namely livestock yield, live animals, and animal slaughtered, and two variables were used to assess agricultural production, namely agricultural production yields and losses. Iran was selected as the case study of the current study. Therefore, time-series data related to livestock and agricultural productions in Iran from 1961 to 2017 have been collected from the FAOSTAT database. First, 70% of this data was used to train ANFIS and MLP, and the remaining 30% of the data was used to test the models. The results disclosed that the ANFIS model with generalized bell-shaped (Gbell) built-in membership functions has the lowest error level in predicting food production. The findings of this study provide a suitable tool for policymakers who can use this model and predict the future of food production to provide a proper plan for the future of food security and food supply for the next generations.
Keywords: food production; machine learning; agricultural production; prediction; big data; data science; deep learning; forecasting; data-driven decision making; food demand; artificial intelligence food production; machine learning; agricultural production; prediction; big data; data science; deep learning; forecasting; data-driven decision making; food demand; artificial intelligence

1. Introduction

Climate change, natural hazards, drought, uncertainty in recourses, and population growth are increasingly threatening the food security of the global nations [1]. It is estimated that the world’s population will exceed 9.7 billion by 2050, which will encourage worldwide hunger and food insecurity [2]. In general, there are two means of the food supply, i.e., domestic production and imports [3]. Awareness of a region’s potential for producing food provides the foundation for developing informed policies for food security. Thus, advancing accurate prediction models is considered essential for food governance and business models [4]. Reliable food prediction models can be used by policymakers to reconsider the annual food import volumes and prices [5]. Furthermore, insight into the food production value to better manage the poverty and support vulnerable groups exposed to food insecurity [6]. Conventional time series and mathematical models had been often used to project food production [7]. Advanced data-driven methods based on artificial intelligence and machine learning have recently shown promising results in providing accurate prediction models. The research for the advancement of reliable artificial intelligence and machine learning methods to be used in a higher level of policymaking is still in the early stage [8,9,10].
A review of the literature for studies that predicted agricultural and livestock production, as the essential representatives of food production, shows that the available studies at the microlevel often focus on a specific crop or individual livestock. For instance, Nosratabadi et al. [7], Pantazi et al. [8], and Sengupta and Lee [9], used machine learning techniques to develop models for crop yield prediction. Nosratabadi et al. [7] developed a gray wolf optimizer of neural networks (GWO-ANN), a hybrid machine learning model, to predict the yield of wheat crops in Iran and they also state that this model has a lower error rate and higher predictive accuracy (with R = 0.48 and root mean square error (RMSE) = 3.19) compared to other models. Pantazi et al. [8] designed a supervised Kohonen networks (SNK) model to predict wheat yield. They report that the accuracy of their model in the prediction of wheat yield was 81.65%. Sengupta and Lee [9] using a support vector machine (SVM) tried to identify the number of immature green citrus and they report that the accuracy of their model was 80.4%. In addition, Morales et al. [10], Alonso, Villa, and Bahamonde [11], and Alonso, Castañón, and Bahamonde [12], for example, have employed machine learning techniques to design models for livestock production. Morales et al. [10] develop an SVM model for the early detection of problems in the production curves of hens’ eggs. They claim that the accuracy of their mode has been equal to 98%. Alonso et al. [11] developed an SVM model to forecast cattle weight trajectories with only one or a few weights. Additionally, they report that the level of error metrics of mean absolute percentage error (MAPE) for their model were between 3.9 and 9.3 for different datasets. Alonso et al. [12] develop an SVM/ support vector regression (SVR) to estimate the beef cattle’ carcass weight 150 days before slaughter. They used MAPE to test the accuracy of their model and they report that the average MAPE of their model was 4.27%. Although research has used advanced machine learning tools to predict agricultural and livestock production, the focus of the research has been on a specific product or livestock, and developed models are not designed to forecast different production at the macrolevel of a country. To address this gap in the literature, the present study intends to develop a model for predicting food production at the macro level of a country using machine learning models.
Since there is ample evidence that agriculture in Iran is facing many problems due to a lack of water resources (e.g., Karandish et al. [13] and Qasemipour and Abbasi [14]), with successive droughts (e.g., Paymard et al. [15]) and poor water management (e.g., Raeisi et al. [16] and Akhoundi and Nazif [17]) cited as reasons for Iran’s lack of water. Such problems have hampered food security at the macrolevel in Iran. On the other hand, Iran, with 79 million in 2015 [18], is one of the most populous countries in the world and is expected to have positive population growth in Iran in the future [18]. There are plenty of studies that explain that some Iranian households are exposed to food insecurity for reasons such as low levels of education and low levels of income (e.g., Ekhlaspour et al. [19], Esfarjani et al. [20], Fathi Beyranvand et al. [21], Najafi Alamdarlo et al. [22]). Therefore, in the present study, Iran was selected as a case study, and the time-series data of agricultural and livestock products related to Iran were used to develop and test the research model.
In the literature, there are advanced and accurate methods for predicting future trends using past data. Artificial intelligence models have the ability to learn from data and can predict non-linear phenomena with very high accuracy based on existing data. There is ample evidence that neural networks, as one of the tools of artificial intelligence, have a very high performance in predicting time series data. For example, Tealeb [23] conducts a review study detailing the articles that used artificial neural networks (ANNs) to predict time series data and shows that the results of ANN are promising in predicting time series data. On the other hand, Tealab, Hefny, and Badr [24] debate that it is better to use advanced and hybrid ANN models in predicting non-linear time series data. Adaptive network-based fuzzy inference system (ANFIS) is a hybrid ANN that is combined with fuzzy systems that can be applied for the time-series data. Hence, the main objective of the current study is to compare the predictive performance of multilayer perceptron (MLP), a type of ANN, and ANFIS in the prediction of the future of agricultural and livestock production in Iran to select the most accurate model. The output of the present study provides policymakers with a comprehensive picture of the future food supply in Iran. Information on predicting indigenous food production provides knowledge to macro-decision makers to design appropriate policies for food security and provide adequate food for future generations. The research has been designed based on a comparative analysis of MLP and ANFIS. Our study investigates the model performance of neural networks and neuro-fuzzy. The structure of the manuscript is represented as follows. First, the data, data source, and the data collection process are elaborated. The machine learning methods used in this paper are then described in detail. After that, the results of comparing MLP and ANFIS are presented. In the next stage, the most accurate model for predicting food production based on the results of accuracy metrics is presented.

1.1. Food Security in Iran

Iran is one of the countries exposed to drought [15] as climate change and inadequate agricultural irrigation systems are among the main reasons mentioned in the literature for the problem of drought in Iran [25]. Drought is a serious threat to food security and has created many challenges for food supply in Iran. Iran is a vast country with diverse climatic conditions that have led to the cultivation of various agricultural products in different parts of the country. Drought and rising population growth, nonetheless, have jeopardized food supply and food security in the country. Qasemipour and Abbasi [14] believe that intensive agricultural practices in Iran led to water scarcity of 206%. Of course, research solutions have been proposed to address water management in order to increase food security and improve food production in Iran. Raeisi et al. [16], for example, consider greenhouses as an alternative to traditional farming because of better water management and higher crop yields. On the other hand, Akhoundi and Nazif [15] propose a model by which wastewater is used to irrigate agricultural fields instead of using natural water. Besides, Esfahani et al. [26] introduce a more creative model to deal with water scarcity in Iran. They consider overseas cultivation as a solution to contribute to food security in Iran.

1.2. Application of Data Science in Food and Agriculture

Many researchers have used data science to solve research problems related to food and agriculture. Since machine learning and deep learning models have the ability to analyze big data, find trends, and make accurate predictions, they have become highly useful tools for researchers [27]. Sengupta and Lee [9] and Su, Xu, and Yan [28], for instance, have used the SVM model and Ali et al. [29] has used the ANFIS model to predict crop yield. The use of learning machine to detect diseases is one of the other applications of machine learning in agriculture. For example, Chung et al. [30] and Ebrahimi et al. [31] used the SVM model to detect diseases in rice and strawberry crops, respectively. The use of ANN models to detect wheat diseases has been very common. So, that Moshou et al. [32] has used the ANN/MLP model, Moshou et al. [33] employ the ANN/SOM model to detect wheat diseases. There are also studies that have used machine learning models to detect weeds. For example, Pantazi et al. [34] and Pantazi, Moshou, and Bravo [35] use an ANN model to detect weeds. Water management and soil management are other applications that have used machine learning models to improve agricultural production. For example, Feng et al. [36] and Patil and Deka [37] use the ANN model to estimate evapotranspiration. Estimation of soil temperature and humidity are also among the applications of machine learning models for soil management. In addition, the use of machine learning models to solve problems related to livestock management has become trendy. Craninx et al. [38], for example, has used the ANN model to forecast rumen fermentation pattern from milk fatty acids in cattle. Alonso, Villa, and Bahamonde [11] uses the SVM model to estimate the weight of cattle at different stages of growth with the least number of weights. Alonso, Castañón, and Bahamonde [12] also used the SVM model to predict carcass weight for beef cattle 150 days before slaughter.
Researchers have also used machine learning models in the food industry. The main applications of machine learning and deep learning in food are to estimate the quality of food. For example, Liu et al. [39] combined stacked sparse autoencoder (SSAE) with CNN to develop a model that detect the quality of vegetables. In addition, Rodriguez et al. [40] and Azizah et al. [41] use CNN to study the quality of fruits. There are studies that evaluate the quality of meat and aquatic products using deep learning models [42,43]. Using machine learning models to study food contaminations is another example of using machine learning in the food industry [44,45].

2. Materials and Methods

2.1. Data

The aim of this study is to develop a model for predicting food production for the next decade in Iran. In the present study, two subvariables of agricultural production and livestock production have been considered to evaluate food production. Three variables, livestock yield, live animals, and animal slaughtered, are used to measure livestock production. This study has also considered two variables, agricultural production yields, and losses, to evaluate the agricultural production. Figure 1 represents the model of the study. In this study, the production of barley, beans, dates, maize, millet, potatoes, rice, soybeans, wheat, rye, and olives is considered as agricultural production in Iran. According to this model, agricultural production yields and losses of the aforementioned products are evaluated as two input variables of agricultural production. Since the losses refer to the loss of productions the respective arrow is drawn outward. For the livestock production, the data related to the live animals such as beehives, buffalo, camel, cattle, chicken, duck, geese, goat, pig, sheep, and turkey, the data related to indigenous meat of buffalo, camel, cattle, chicken, duck, geese, goat, pig, sheep, and turkey, and the data related to milk of buffalo milk, cow, goat, and sheep are collected. These data are collected from the FAO database, i.e., FAOSTAT, that can be accessed on (accessed on 20 September 2020). The collected data covers the period of 1961–2017.
Figure 1 shows that the indigenous livestock production quantity and the indigenous agricultural production in Iran are considered as the country’s potential food production for this country. Two variables of yield and losses were used to evaluate and measure agricultural production, and three variables of live animals, livestock yield, and slaughtered animals were used to measure livestock production quantity.

2.2. Methods

For predicting the future trends of food production in Iran, two models of MLP and ANFIS are applied in the collected data, and the predictive performance of the models are compared based on the accuracy metrics. We trained the proposed models by minimizing a regularized loss function on the training set and evaluated the models by comparing the accuracy metrics on the test set.

2.2.1. MLP

Multilayer perceptron (MLP) is a type of neural network that has a supervised learning technique using the back-propagation method. Figure 2 shows that MLP benefits from a three-layer structure, including the input layer, hidden layer/s, and output layer/s, in which each neuron is connected to all the neurons in the next layer. It is frequently reported that MLP has a great function in non-linear problems [46,47].
Equation (1) shows how the output of input variables, bias values, and input values are calculated:
S i = i = 1 n w i j I i + β i
where I represent the input layer, Ii is the input variable i, n shows the total number of inputs, βj is a bias value, and ωij is the weight of connections in j level. The sigmoid function is mostly used as the activation functions in MLP and it can be calculated through Equation (2):
f j = 1 1 + e S j
where, S is the activation function. Therefore, the ultimate output neuron j can be measured in Equation (3):
y i = f i ( i = 1 n w i j I i + β i )
where, y presents the output value of the MLP method, which needs to be compared by the target values for calculating the model performance. MLP was trained by 70% of total data as a training dataset, which has been sorted randomly by the model. The training was performed by different sets of the neuron numbers in the hidden layer for finding the best architecture for the predictor model from 10 to 18 by interval 4. The activation function was selected to be the Tanh(x) due to its higher performance compared with other activation functions.

2.2.2. ANFIS

The adaptive network-based fuzzy inference system is a hybrid neural network in which a fuzzy logic (FL) is embedded to the artificial neural network (ANN) architecture to identify the optimal distribution of membership functions [48]. The inference system of ANFIS consists of five layers in which the input of each layer is the output of the previous layer. This method applies fuzzy if-then rules of Sugeno, and if an ANFIS model has two inputs (x, y) and one output (fi), for example, the two rules for a first-order two-rule are:
  • Rule 1: if x is A1 and y is B1 then z is f1(x, y);
  • Rule 2: if x is A2 and y is B2 then z is f2(x, y).
where x and y are the ANFIS inputs, A and B are the fuzzy sets, and fi(x, y) is the outputs of the first-order Sugeno fuzzy. The architecture of an ANFIS model constitutes adaptive nodes and fixed nodes (see Figure 3). The first layer of the model includes adaptive nodes that can be calculated through Equations (4)–(6).
O 1 ,   i μ A i ( x )   f o r   i = 1 , 2
O 1 ,   i μ B i ( y )   f o r   i = 1 , 2
μ ( x ) = 1 1 + ( x c i a i ) 2 b i
where x and y are the inputs, A and B are the linguistic labels, µ(x) and µ(y) are membership functions that take values between 0 and 1, and ai, bi, and ci are the parameter sets.
The second layer, which is shown in red circles in Figure 3, is a fixed node and can be calculated through Equation (7). It is worth mentioning that ωi is the firing strength of a rule.
O 2 , i = w i = μ A i ( x ) · μ B i ( y ) ,           for   i = 1 , 2
O2,i as the output of the second layer enters to the third layer. The third layer, which is presented in yellow circles in Figure 3, is also a fixed node. Its main goal is to normalize the firing strength by using Equation (8).
O 3 , i = w ¯ i = w i   w i = w i w 1 + w 2 ,       for   i = 1 , 2
The fourth layer is an adaptive node as well and depicted as green squares. Equation (9) is used to measure the fourth layer.
O 4 , i = w ¯ i · f i ,         for   i = 1 , 2
  • Rule 1: if x is A1 and y is B1 then f1 = p1x + q1y + r1;
  • Rule 2: if x is A2 and y is B2 then f2 = p2x + q2y + r2.
where pi, qi, and ri are the parameters sets.
The fifth layer is also a fixed node presented in the form of a blue circle in Figure 3 and can be calculated through Equation (10).
O 5 , i = f o u t =   w ¯ i · f i = O v e r a l o u t p u t ,       for   i = 1 , 2
The final output of an ANFIS structure, which is shown as fout in Figure 3, can be calculated through Equation (11).
c f o u t = w ¯ 1 f 1 + w ¯ 2 f 2 = w 1 w 1 + w 2 f 1 + w 2 w 1 + w 2 f 2 = ( w ¯ 1 x ) p 1 + ( w ¯ 1 y ) q 1 + ( w ¯ 1 ) r 1 + ( w ¯ 2 x ) p 2 + ( w ¯ 2 y ) q 2 + ( w ¯ 2 ) r 2
ANFIS was trained using 70% of the total data (randomly selected). Input variables were time-series data. The training parameter was the type of the membership function (MF). Since it has the maximum effect on the accuracy and performance of the ANFIS model. Triangular, trapezoidal, and G-bell types were selected as the frequently used and popular MF types for comparison purposes in the presence of linear output MF type (for its highest accuracy in comparison with constant type MF). Other parameters like the number of MF types and hybrids method were considered to be constant because they did not have any significant effect on the modeling procedure in the present study. One of the main reasons can be the dimension of the dataset in the present study. The rest of the data set (30% of the total dataset) was employed for the testing step.

2.3. Accuracy Metrics

To compare the predictive power and accuracy performance of MLP and ANFIS two evaluation criteria namely RMSE and determination coefficient (R) are measured for both models. Equations (12) and (13) respectively show how to calculate RMSE and R2.
R M S E = 1 N i = 1 N ( A P ) 2
R 2 = 1 ( i = 1 n ( A P ) 2 i = 1 n ( A ) 2 )
where A is the target values, P refers to the predicted values (output of models), and N is the number of data. Using these performance parameters, the accuracy of models can be calculated for comparison purposes.

3. Results

In this study, the process of selecting the appropriate model with better predictive power was designed in such a way that the models were first trained by 70% of the data. After the training phase, the predictive performance of the models was tested on the remaining 30% of the data, and then the accuracy of the models was measured and compared by accuracy metrics RMSE and R2. Table 1 shows that the variables of xt−1, xt−2, and xt−3, which are respectively the representation of live animals, animals slaughtered, and livestock yield, are the inputs variables of livestock production quantity and xt−4 and xt−5, which are respectively the representation of yield and losses of agricultural productions, are the input variables of agricultural productions. In other words, the current model constitutes two outputs: (1) livestock production and (2) agricultural production.

3.1. Training Results

As it is mentioned above, 70% of the data are used to train the models. The training phase was repeated three times, with each model being tested with a different number of neurons.
By changing the number of neurons, the accuracy of the MLP model can be controlled and it reveals the most accurate model. Table 2 shows that in the training phase of the MLP model with the number of neurons ten, fourteen, and eighteen were tested. At this stage, the model with ten neurons for predicting livestock production and the model with 18 neurons for predicting agricultural production had the best performance because the corresponding RMSEs were lower compared to other models.
On the other hand, to control the accuracy of the ANFIS model in the training phase, the predictive accuracy of different membership functions (MF) was tested. In this study, triangular-shaped (Tri.), trapezoidal-shaped (Trap.), and generalized bell-shaped (Gbell) built-in membership functions were evaluated. The results of the evaluation of the accuracy of MFs are presented in Table 3.
The results show that the model with the Trap. built-in membership function had the highest accuracy for predicting both livestock and agricultural production because the RMSE of this model was 4,080,579.79 for livestock production and 987,950.19 for agricultural production, which were lower than other membership functions. The comparison of Table 2 and Table 3 illustrates that the performance of the ANFIS model compared to the MLP model in predicting both agricultural and livestock production was higher. Since the values of RMSE of this model in all cases was lower than the MLP model in the training phase.

3.2. Testing Results

After training the models, the models were tested by 30 percent of the data to examine the predictivity power of models. The results of the testing phase of the MLP model were in accordance with the results of the training phase as the MLP model with ten neurons had the highest accuracy for predicting livestock production because the RMSE of this model was equal to 265,590,099.2, which was lower than other models with different neurons. In addition, the RMSE of the MLP model 18 neurons for testing agricultural production was 33,575,595.74 that was lower than the other models indicating the higher accuracy of this model compare to the other models (see Table 4).
However, Table 5 shows that in the testing phase, the ANFIS model with the Gbell membership function had more accurate results with less error levels in both livestock production prediction (with RMSE = 6,052,851.43) and agricultural production prediction (with RMSE = 1,724,426) while in the training phase the Trap. membership function model had the highest accuracy rate. Comparing the results of the testing phase with the training phase was the same, and in both phases, the ANFIS model provided higher performance than the MLP model due to the low level of error in predicting both livestock and agricultural production. Therefore, the present study proposes the ANFIS model for predicting food production.
The coefficient of determination of the ANFIS model was also tested. Figure 4 discloses that the coefficient of determination (R2) of the ANFIS model was very high for both livestock (Figure 4a) and agricultural production (Figure 4b) forecast so that R2 was equal to 0.99 for livestock production and 0.94 for agricultural production.

3.3. Prediction Results

The results showed that the ANFIS model with the Gbell membership function, due to the lower RMSE, not only had a better predictive performance in both agricultural and livestock production forecasting compared to the ANFIS model with other membership functions, but also it had a higher predictability power on the current data compared to the MLP model. Consequently, this model was selected to predict food production in Iran. The results of the prediction of Iranian agricultural and livestock products for 2018–2030 using the ANFIS model with the Gbell membership function are presented in Table 6.
In order to better represent the predicted trend for agricultural and livestock products in Iran, Figure 5 was designed based on the predicted data. Figure 5 shows that agricultural and livestock products in Iran are expected to have an upward trend with almost the same slope. This is because the predictive model of this study, using time series data, predicts that food production in Iran will increase in the upcoming decade.

4. Conclusions

As the world’s population grows, so does the demand for food, and in recent years the number of people exposed to hunger, and even severe hunger, is increasing daily. Governments and organizations active in the food industry are planning and preparing to prevent potential problems that may arise in the way of food security for future generations. To achieve food security goals, food is mainly supplied through domestic production and import. Therefore, studying a country’s potential for food supply is the first step in planning for food security. Food production prediction gives a realistic view to policy makers and activists in the agricultural and food industries for long-term and short-term planning. Therefore, the present study tried to provide a suitable model with high predictive performance for predicting food production. The present study predicted Iran’s agricultural and livestock production for the next ten years. According to the results, it is predicted that in the next ten years, the volume of both agricultural and livestock production in Iran will increase. The findings of this study provide a basis for planning the production volume required for the coming years, planning for budgeting and agricultural subsidies, planning for the active workforce in the agricultural and livestock sectors. In addition, according to forecasts, decision-makers can plan to import needed food production and export surplus domestic production. Using machine learning, researchers have come up with creative and precise solutions to a variety of food and agricultural problems, such as crop yields prediction. However, there is no research to predict food production. The present study used machine learning models to predict agricultural and livestock products in Iran. For this purpose, the performance of two models, MLP and ANFIS, was tested using time series data of agricultural and livestock production in Iran. The results of accuracy metrics revealed that the ANFIS model has higher predictive power than the MLP model due to its higher predictive accuracy. The current study contributes to food security research by providing a repayable tool to predict the future of agricultural and livestock production. Researchers and decision-makers can use this model to predict the future of food security in a region. Therefore, for future research, it is suggested that using the proposed model of the present study to predict food production in different countries and provide appropriate solutions to combat food insecurity. One of the limitations of this study is that forecasts for agricultural and livestock production are based only on time series data while other factors such as climate, government policies, and technological advances are considered constant. Another limitation of this article is the generalization of the finding that the ANFIS model outperforms the MLP model because this finding is limited to the time series data of Iran and the result may differ in data related to another country.

Author Contributions

Conceptualization, S.N., Z.L.; methodology, A.M., S.A.; investigation, S.N.; software, S.A.; formal analysis, A.M., S.A.; writing—original draft preparation, S.N.; writing—review and editing, C.M., Z.L.; visualization, A.M.; supervision, A.M., Z.L. All authors have read and agreed to the published version of the manuscript.


This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.


ANFISAdaptive network-based fuzzy inference system
ANNArtificial neural network
CNNConvolutional neural network
DSData science
MLMachine learning
MLPMultilayer perceptron
SOMSelf-organizing Map
SSAEStacked sparse autoencoder
SVMSupport vector machine


  1. World Health Organization. The State of Food Security and Nutrition in the World: Safeguarding against Economic Slowdowns and Downturns; FAO: Rome, Italy, 2019. [Google Scholar]
  2. Wunderlich, S.M.; Martinez, N.M. Conserving natural resources through food loss reduction: Production and consumption stages of the food supply chain. Int. Soil Water Conserv. Res. 2018, 6, 331–339. [Google Scholar] [CrossRef]
  3. Nosratabadi, S.; Mosavi, A.; Lakner, Z. Food Supply Chain and Business Model Innovation. Foods 2020, 9, 132. [Google Scholar] [CrossRef] [PubMed]
  4. Kang, Y.; Khan, S.; Ma, X. Climate change impacts on crop yield, crop water productivity and food security—A review. Prog. Nat. Sci. 2009, 19, 1665–1674. [Google Scholar] [CrossRef]
  5. Iizumi, T.; Sakuma, H.; Yokozawa, M.; Luo, J.-J.; Challinor, A.J.; Brown, M.E.; Sakurai, G.; Yamagata, T. Prediction of seasonal climate-induced variations in global food production. Nat. Clim. Chang. 2013, 3, 904–908. [Google Scholar] [CrossRef]
  6. Hertel, T.W.; Burke, M.B.; Lobell, D.B. The poverty implications of climate-induced crop yield changes by 2030. Glob. Environ. Chang. 2010, 20, 577–585. [Google Scholar] [CrossRef]
  7. Nosratabadi, S.; Karoly, S.; Beszedes, B.; Felde, I.; Ardabili, S.; Mosavi, A. Comparative Analysis of ANN-ICA and ANN-GWO for Crop Yield Prediction. In Proceedings of the 2020 RIVF International Conference on Computing and Communication Technologies (RIVF), Ho Chi Minh, Vietnam, 14–15 October 2020. [Google Scholar]
  8. Pantazi, X.E.; Moshou, D.; Alexandridis, T.; Whetton, R.L.; Mouazen, A.M. Wheat yield prediction using machine learning and advanced sensing techniques. Comput. Electron. Agric. 2016, 121, 57–65. [Google Scholar] [CrossRef]
  9. Sengupta, S.; Lee, W.S. Identification and determination of the number of immature green citrus fruit in a canopy under different ambient light conditions. Biosyst. Eng. 2014, 117, 51–61. [Google Scholar] [CrossRef]
  10. Morales, I.R.; Cebrián, D.R.; Blanco, E.F.; Sierra, A.P. Early warning in egg production curves from commercial hens: A SVM approach. Comput. Electron. Agric. 2016, 121, 169–179. [Google Scholar] [CrossRef]
  11. Alonso, J.; Villa, A.; Bahamonde, A. Improved estimation of bovine weight trajectories using Support Vector Machine Classification. Comput. Electron. Agric. 2015, 110, 36–41. [Google Scholar] [CrossRef]
  12. Alonso, J.; Castañón, Á.R.; Bahamonde, A. Support Vector Regression to predict carcass weight in beef cattle in advance of the slaughter. Comput. Electron. Agric. 2013, 91, 116–120. [Google Scholar] [CrossRef]
  13. Karandish, F.; Salari, S.; Darzi-Naftchali, A. Application of virtual water trade to evaluate cropping pattern in arid regions. Water Resour. Manag. 2015, 29, 4061–4074. [Google Scholar] [CrossRef]
  14. Qasemipour, E.; Abbasi, A. Virtual water flow and water footprint assessment of an arid region: A case study of South Khorasan province, Iran. Water 2019, 11, 1755. [Google Scholar] [CrossRef]
  15. Paymard, P.; Yaghoubi, F.; Nouri, M.; Bannayan, M. Projecting climate change impacts on rainfed wheat yield, water demand, and water use efficiency in northeast Iran. Theor. Appl. Climatol. 2019, 138, 1361–1373. [Google Scholar] [CrossRef]
  16. Raeisi, L.G.; Morid, S.; Delavar, M.; Srinivasan, R. Effect and side-effect assessment of different agricultural water saving measures in an integrated framework. Agric. Water Manag. 2019, 223, 105685. [Google Scholar] [CrossRef]
  17. Akhoundi, A.; Nazif, S. Sustainability assessment of wastewater reuse alternatives using the evidential reasoning approach. J. Clean. Prod. 2018, 195, 1350–1376. [Google Scholar] [CrossRef]
  18. Mehri, N.; Messkoub, M.; Kunkel, S. Trends, determinants and the implications of population aging in Iran. Ageing Int. 2020, 45, 1–17. [Google Scholar] [CrossRef]
  19. Ekhlaspour, P.; Foroumandi, E.; Ebrahimi-Mameghani, M.; Jafari-Koshki, T.; Arefhosseini, S.R. Household food security status and its associated factors in Baft-Kerman, IRAN: A cross-sectional study. Ecol. Food Nutr. 2019, 58, 608–619. [Google Scholar] [CrossRef]
  20. Esfarjani, F.; Hosseini, H.; Khaksar, R.; Roustaee, R.; Alikhanian, H.; Khalafi, M.; Khaneghah, A.M.; Mohammadi-Nasrabadi, F. Home Food Safety Practice and Household Food Insecurity: A Structural Equation Modeling Approach. Iran. J. Public Health 2019, 48, 1870. [Google Scholar] [PubMed]
  21. Fathi Beyranvand, H.; Eghtesadi, S.; Atai Jafari, A.; Movahedi, A. Prevalence of Food Insecurity in Pregnant Women in Khorramabad City and its Association with General Health and other Factors. Iran. J. Nutr. Sci. Food Technol. 2019, 14, 21–30. [Google Scholar]
  22. Alamdarlo, H.N.; Riyahi, F.; Vakilpoor, M.H. Wheat self-sufficiency, water restriction and virtual water trade in Iran. Netw. Spat. Econ. 2019, 19, 503–520. [Google Scholar] [CrossRef]
  23. Tealab, A. Time series forecasting using artificial neural networks methodologies: A systematic review. Future Comput. Inform. J. 2018, 3, 334–340. [Google Scholar] [CrossRef]
  24. Tealab, A.; Hefny, H.; Badr, A. Forecasting of non-linear time series using ANN. Future Comput. Inform. J. 2017, 2, 39–47. [Google Scholar] [CrossRef]
  25. Karandish, F.; Hoekstra, A. Informing national food and water security policy through water footprint assessment: The case of Iran. Water 2017, 9, 831. [Google Scholar] [CrossRef]
  26. Esfahani, A.K.; Mirdamadi, S.M.; Hosseini, S.J.; Lashgarara, F. Overseas cultivation: The complimentary approach for developing food security. Bulg. J. Agric. Sci. 2019, 25, 26–35. [Google Scholar]
  27. Vargas, R.; Mosavi, A.; Ruiz, R. Advances in Intelligent Systems and Computing; Springer: Berlin, Germany, 2017; Volume 7, pp. 122–148. [Google Scholar]
  28. Su, Y.; Xu, H.; Yan, L. Support vector machine-based open crop model (SBOCM): Case of rice production in China. Saudi J. Biol. Sci. 2017, 24, 537–547. [Google Scholar] [CrossRef] [PubMed]
  29. Ali, I.; Cawkwell, F.; Dwyer, E.; Green, S. Modeling Managed Grassland Biomass Estimation by Using Multitemporal Remote Sensing Data—A Machine Learning Approach. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2016, 10, 3254–3264. [Google Scholar] [CrossRef]
  30. Chung, C.L.; Huang, K.J.; Chen, S.Y.; Lai, M.H.; Chen, Y.C.; Kuo, Y.F. Detecting Bakanae disease in rice seedlings by machine vision. Comput. Electron. Agric. 2016, 121, 404–411. [Google Scholar] [CrossRef]
  31. Ebrahimi, M.A.; Khoshtaghaza, M.H.; Minaei, S.; Jamshidi, B. Vision-based pest detection based on SVM classification method. Comput. Comput. Electron. Agric. 2017, 137, 52–58. [Google Scholar] [CrossRef]
  32. Moshou, D.; Bravo, C.; West, J.; Wahlen, S.; McCartney, A.; Ramon, H. Automatic detection of “yellow rust” in wheat using reflectance measurements and neural networks. Comput. Electron. Agric. 2004, 44, 173–188. [Google Scholar] [CrossRef]
  33. Moshou, D.; Bravo, C.; Oberti, R.; West, J.; Bodria, L.; McCartney, A.; Ramon, H. Plant disease detection based on data fusion of hyper-spectral and multi-spectral fluorescence imaging using Kohonen maps. Real-Time Imaging 2005, 11, 75–83. [Google Scholar] [CrossRef]
  34. Pantazi, X.E.; Tamouridou, A.A.; Alexandridis, T.K.; Lagopodi, A.L.; Kashefi, J.; Moshou, D. Evaluation of hierarchical self-organising maps for weed mapping using UAS multispectral imagery. Comput. Electron. Agric. 2017, 139, 224–230. [Google Scholar] [CrossRef]
  35. Pantazi, X.-E.; Moshou, D.; Bravo, C. Active learning system for weed species recognition based on hyperspectral sensing. Biosyst. Eng. 2016, 146, 193–202. [Google Scholar] [CrossRef]
  36. Feng, Y.; Peng, Y.; Cui, N.; Gong, D.; Zhang, K. Modeling reference evapotranspiration using extreme learning machine and generalized regression neural network only with temperature data. Comput. Electron. Agric. 2017, 136, 71–78. [Google Scholar] [CrossRef]
  37. Patil, A.P.; Deka, P.C. An extreme learning machine approach for modeling evapotranspiration using extrinsic inputs. Comput. Electron. Agric. 2016, 121, 385–392. [Google Scholar] [CrossRef]
  38. Craninx, M.; Fievez, V.; Vlaeminck, B.; De Baets, B. Artificial neural network models of the rumen fermentation pattern in dairy cattle. Comput. Electron. Agric. 2008, 60, 226–238. [Google Scholar] [CrossRef]
  39. Liu, Z.; He, Y.; Cen, H.; Lu, R. Deep feature representation with stacked sparse auto-encoder and convolutional neural network for hyperspectral imaging-based detection of cucumber defects. Trans. ASABE 2018, 61, 425–436. [Google Scholar] [CrossRef]
  40. Rodriguez, F.J.; Garcia, A.; Pardo, P.J.; Chavez, F.; Luque-Baena, R.M. Study and classification of plum varieties using image analysis and deep learning techniques. Prog. Artif. Intell. 2018, 7, 119–127. [Google Scholar] [CrossRef]
  41. Azizah, L.M.; Umayah, S.F.; Riyadi, S.; Damarjati, C.; Utama, N.A. Deep learning implementation using convolutional neural network in mangosteen surface defect detection. In 7th IEEE International Conference on Control System. Comput. Eng. 2017, 242–246. [Google Scholar] [CrossRef]
  42. Cheng, J.H.; Sun, D.W. Partial least squares regression (PLSR) applied to NIR and HSI spectral data modeling to predict chemical properties of fish muscle. Food Eng. Rev. 2017, 9, 36–49. [Google Scholar] [CrossRef]
  43. Yu, X.J.; Wang, J.P.; Wen, S.T.; Yang, J.Q.; Zhang, F.F. A deep learning based feature extraction method on hyperspectral images for nondestructive prediction of TVB-N content in Pacific white shrimp (Litopenaeus vannamei). Biosyst. Eng. 2019, 178, 244–255. [Google Scholar] [CrossRef]
  44. Bisgin, H.; Bera, T.; Ding, H.J.; Semey, H.G.; Wu, L.H.; Liu, Z.; Barnes, A.E.; Langley, D.A.; Pava-Ripoll, M.; Vyas, H.J.; et al. Comparing SVM and ANN based machine learning methods for species identification of food contaminating beetles. Sci. Rep. 2018, 8, 12. [Google Scholar] [CrossRef] [PubMed]
  45. Song, Q.; Zheng, Y.J.; Xue, Y.; Sheng, W.G.; Zhao, M.R. An evolutionary deep neural network for predicting morbidity of gastrointestinal infections by food contamination. Neurocomputing 2017, 226, 16–22. [Google Scholar] [CrossRef]
  46. Boniecki, P.; Zaborowicz, M.; Pilarska, A.; Piekarska-Boniecka, H. Identification Process of Selected Graphic Features Apple Tree Pests by Neural Models Type MLP, RBF and DNN. Agriculture 2020, 10, 218. [Google Scholar] [CrossRef]
  47. Olier, I.; Sansom, A.; Lisboa, P.; Ortega-Martorell, S. Using MLP partial responses to explain in-hospital mortality in ICU. In IEEE International Conference on Data Analytics for Business and Industry: Way towards a Sustainable Economy (ICDABI); IEEE: Piscataway, NJ, USA, 2020; pp. 1–5. [Google Scholar]
  48. Ying, L.-C.; Pan, M.-C. Using adaptive network based fuzzy inference system to forecast regional electricity loads. Energy Convers. Manag. 2008, 49, 205–211. [Google Scholar] [CrossRef]
Figure 1. The proposed model of the study for indigenous food production in Iran.
Figure 1. The proposed model of the study for indigenous food production in Iran.
Agriculture 11 00408 g001
Figure 2. The architecture of the multilayer perceptron neural networks.
Figure 2. The architecture of the multilayer perceptron neural networks.
Agriculture 11 00408 g002
Figure 3. The architecture of adaptive network-based fuzzy inference system.
Figure 3. The architecture of adaptive network-based fuzzy inference system.
Agriculture 11 00408 g003
Figure 4. The coefficient of determination of the ANFIS model for livestock and agricultural production prediction. (a) Livestock Production; (b) Agricultural Production.
Figure 4. The coefficient of determination of the ANFIS model for livestock and agricultural production prediction. (a) Livestock Production; (b) Agricultural Production.
Agriculture 11 00408 g004
Figure 5. The result of predicting agricultural and livestock production for the next ten years in Iran.
Figure 5. The result of predicting agricultural and livestock production for the next ten years in Iran.
Agriculture 11 00408 g005
Table 1. The prepared dataset for time-series prediction.
Table 1. The prepared dataset for time-series prediction.
xt−1, xt−2, xt−3, xt−4, xt−5O1 = Livestock production
O2 = Agricultural production
Table 2. RMSE results for MLP models with different numbers of neurons in the training phase.
Table 2. RMSE results for MLP models with different numbers of neurons in the training phase.
VariableNeuron NumberRMSE
Livestock Production10275,284,878.3
Livestock Production14462,563,347.1
Livestock Production18320,412,824.4
Agri. Production1036,325,828
Agri. Production1477,746,693.65
Agri. Production1835,410,107.42
Table 3. RMSE results for ANFIS models with different MF types in the training phase.
Table 3. RMSE results for ANFIS models with different MF types in the training phase.
VariableMF TypeRMSE
Livestock ProductionTri.17,225,511.04
Livestock ProductionTrap.4,080,579.79
Livestock ProductionGbell6,750,734
Agri. ProductionTri.2,144,876.04
Agri. ProductionTrap.987,950.19
Agri. ProductionGbell9,751,562
Table 4. RMSE results for MLP models with different numbers of neurons in the testing phase.
Table 4. RMSE results for MLP models with different numbers of neurons in the testing phase.
VariableNeuron NumberRMSE
Livestock Production10265,590,099.2
Livestock Production14457,160,675.6
Livestock Production18311,543,277.9
Agri. Production1040,310,186.93
Agri. Production1482,380,698.29
Agri. Production1833,575,595.74
Table 5. RMSE results for ANFIS models with different MF types in the testing phase.
Table 5. RMSE results for ANFIS models with different MF types in the testing phase.
VariableMF TypeRMSE
Livestock ProductionTri.11,124,369
Livestock ProductionTrap.17,894,505.8
Livestock ProductionGbell6,052,851.43
Agri. ProductionTri.2,264,668
Agri. ProductionTrap.2,415,988
Agri. ProductionGbell1,724,426
Table 6. Prediction of agricultural and livestock production in Iran.
Table 6. Prediction of agricultural and livestock production in Iran.
Livestock ProductsAgricultural Products
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Back to TopTop