Next Article in Journal
On Convergence Rates of Some Limits
Next Article in Special Issue
On the Advent of Fractional Calculus in Econophysics via Continuous-Time Random Walk
Previous Article in Journal
Some Identities Involving Two-Variable Partially Degenerate Hermite Polynomials Induced from Differential Equations and Structure of Their Roots
Previous Article in Special Issue
Fractional Dynamics and Pseudo-Phase Space of Country Economic Processes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Gross Domestic Product per Capita of Countries

Informatics Institute, Istanbul Technical University, 34467 Istanbul, Turkey
*
Author to whom correspondence should be addressed.
Mathematics 2020, 8(4), 633; https://doi.org/10.3390/math8040633
Submission received: 19 February 2020 / Revised: 29 March 2020 / Accepted: 9 April 2020 / Published: 20 April 2020
(This article belongs to the Special Issue Mathematical Economics: Application of Fractional Calculus)

Abstract

:
In this study, a new approach for time series modeling and prediction, “deep assessment methodology,” is proposed and the performance is reported on modeling and prediction for upcoming years of Gross Domestic Product (GDP) per capita. The proposed methodology expresses a function with the finite summation of its previous values and derivatives combining fractional calculus and the Least Square Method to find unknown coefficients. The dataset of GDP per capita used in this study includes nine countries (Brazil, China, India, Italy, Japan, the UK, the USA, Spain and Turkey) and the European Union. The modeling performance of the proposed model is compared with the Polynomial model and the Fractional model and prediction performance is compared to a special type of neural network, Long Short-Term Memory (LSTM), that used for time series. Results show that using Deep Assessment Methodology yields promising modeling and prediction results for GDP per capita. The proposed method is outperforming Polynomial model and Fractional model by 1.538% and by 1.899% average error rates, respectively. We also show that Deep Assessment Method (DAM) is superior to plain LSTM on prediction for upcoming GDP per capita values by 1.21% average error.

1. Introduction

In the last quarter of the century, the data exchange with not only person to person but also, machine to machine has increased tremendously. Developments in technology and informatics in parallel with the development of data science lead the companies, institutions, universities and especially, the countries to give priority to evaluating produced data and predicting what can be forthcoming. The modeling of all technical, economic, social events and data has been the interest of scientists for many years [1,2,3,4]. Many authors have been investigating the modeling and predicting events, options, choices and data. Especially, there is a huge research interest in finding any relation between telecommunication, economic growth and financial development [5,6,7,8,9,10,11,12]. One of the approaches to model a physical phenomenon or a mathematical study is to model the dependent variable satisfying differential equation with respect to the independent variable. However, the differential equations with an integer-order proposed for mathematical economics or data modeling cannot describe processes with memory and non-locality because the integer-order derivatives have the property of the locality. On the other hand, the fractional-order differential equation is a branch of mathematics that focuses on fractional-order differential and integral operators and can be used to address the limitations of integer order differential models. Using the fractional calculus or converting the integer-order differential equation into the non-integer order differential equations lead to a very essential advantage which is memory property of the fractional-order derivative. This is very crucial for models related to economics which in general, deal with the past and the effect of the past and now on future [12,13]. The memory capability of the fractional differential approach is the foundation of our motivation.
Fractional calculus (FC) as a question to Wilhelm Leibniz (1646–1716) first arose in 1695 from French Mathematician Marquis de L’Hopital (1661–1704) [11]. The main question of interest was what if the order of derivative were a real number instead of an integer. After that, the FC idea has been developed by many mathematicians and researchers throughout the eighteenth and nineteenth centuries. Now, there exist several definitions of the fractional-order derivative, including Grünwald-Letnikov, Riemann-Liouville, Weyl, Riesz and the Caputo representation. The fractional approach is used in many studies because the fractional derivative represents the intermediate states between two known states. For example, zero order-derivative of the function means the function itself while the first-order derivative represents the first derivative of the function. Between these known states, there are infinite intermediate states [11]. The use of semi-derivatives and integrals in the mass and heat transfer become an important instant in the field of fractional calculus due to employing the mathematical definitions into physical phenomena [12,13]. In the last decade, using fractional operators which explain the events, situations or modes between two different stages or the phenomena with memory provide more accurate models in many branches of science and engineering including chemistry, biology, biomedical devices, nanotechnology, diffusion, diffraction and economics [12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31]. In References [25,26,27,28,29,30,31], the modeling and comparison of the countries and trends in the sense of economics and its parameters are implemented. In References [25,26], economic processes with memory are discussed and modeling is obtained by using the fractional calculus. The studies with similar purposes as we aim such as modeling or prediction exist. In these studies, the fractional calculus is employed to model the given dataset and to predict for the forthcoming. In Reference [28], the orthogonal distance fitting method is used. The study is trying to minimize the sum of the orthogonal distance of data points in order to obtain an optimized continuous curve representing the data points. In Reference [32], the one-parameter fractional linear prediction is studied using the memory of two, three or four samples, without increasing the number of predictor coefficients defined in the study. In Reference [33], the generalized formulation of the optimal low-order linear prediction by using the fractional calculus approach is developed with restricted memory. All these studies focus on modeling or prediction for a phenomenon with fractional calculus. Also, in our previous studies, methods based on FC that works for modeling were introduced. In these studies, the children’s physical growth, subscriber’s numbers of operators, GDP per capita were modeled and compared with other modeling approaches such as Fractional Model-1 and Polynomial Models [34,35,36]. According to the results, proposed fractional models had better results compared to the results obtained from Linear and Polynomial Models [34,35,36]. Our previous works do not take into account the previous values of the dataset for any time instant. Their purpose is to model the dataset with minimum error and faster way compared with classical methods such as Polynomial and Linear Regression.
In this study, we extend our prior works by predicting the next incoming values as well as modeling the data itself. We introduce a new mathematical model, namely “Deep Assessment,” based on the fractional differential equation for modeling and prediction by using the properties of fractional calculus. Different to the literature and our early studies mentioned above, this model can be used for prediction as well as modeling. The proposed approach is built on the fractional-order differential equation and corresponding Laplace transform properties are utilized. Here, the modeling is implemented with mathematical tools similar to those developed in the previous study [4] with a different approach in which the finite numbers of previous values and the derivatives are taken into account. Then, the prediction is obtained by assuming a value in a specific time can be expressed as the summation of the previous values weighted by unknown coefficients and the function to be modeled is continuous and differentiable. In this way, the proposed method takes previous values and variation rates between different time samples (derivative) of the dataset into account while modeling the data itself and predicting upcoming values. Combining the previous values with the variations weighted by the unknown coefficients lead to calling the method “deep assessment.”
In this study, we assessed the proposed method by the modeling, testing and predicting GDP per capita of the following countries and the European Union: Brazil, China, European Union, India, Italy, Japan, the UK, the USA, Spain and Turkey. GDP per capita is a measure of a country’s total economic output divided by the number of the population of the country. In general, it is a reasonable and good measurement of a country’s living quality and standard [37]. Therefore, the modeling of GDP per capita is crucial and predicting GDP per capita is very essential not only for researchers but also for companies, investors, manufacturers and institutions. To assess the performance of Deep Assessment in modeling, we compare the proposed model with Polynomial Regression and Fractional Model-1 [34]. Besides, in the same way, for the prediction, we compared the model with Long-Short Term Memory (LSTM), a special type of neural networks used in time series problems.
The structure of the study is the following. Section 2 explains the formulation of the problem. After that, Section 3, namely Our Approach, is devoted to explaining how to obtain modeling, simulation, testing and prediction. Then, in Section 4, the results are presented. Lastly, Section 5 highlights the conclusion of the study.

2. Formulation of the Problem

In this section, the mathematical foundation of the proposed method is given. Before going into the mathematical manipulations, it is better to explain the approach and the main steps for the formulation. The study aims to model and then, to predict GDP per capita data at any time t by using the previous GDP per capita values of the countries. Here, we assume that countries’ historical data and the change of these data over time create an eco-genetics for the forthcoming. In other words, mathematically, GDP per capita at a time t is assumed to be the summation of both its previous values and the changes in time with unknown constant coefficients. In the second stage, we express a function for the GDP per capita as a series expansion by using Taylor expansion of a continuous and bounded function. Then, the differential equation obtained from this series expansion is defined. After that, the unknown constant coefficients are found by the least-squares method. The method aims to minimize the error between the proposed GDP per capita function and the dataset.
First, it is a reasonable idea to approximate a function g ( x ) as the finite summation of the previous values of the same function weighted with unknown coefficients   α k and the summation of the derivatives of the previous values of the same function weighted with unknown coefficients   β k because, intuitively, the recent value of data, in general, is related to and correlated with its previous values and the change rates. The purpose is to find the upcoming values of any dataset with a minimum error by employing the previously inherited features of the dataset. As a starting point, an arbitrary function is assumed to be approximately the finite summation of the previous values and the change rates weighted with some constant coefficients. To use the heritability of fractional calculus, this presupposition for modeling of the function itself and predicting future values is done [6,28,34].
g ( x ) k = 1 l α k g ( x k ) + k = 1 l β k g ( x k ) .
Here, g is the first derivative of g ( x k ) with respect to   x . After assuming Equation (1), the function g ( x ) can be expanded as the summation of polynomials with unknown constant coefficients, a n as given in Equation (2). Here, g ( x ) is assumed to be a continuous and differentiable function.
g ( x ) = n = 0 a n x n .
Then, g ( x k ) becomes as Equation (3)
g ( x k ) = n = 0 a n ( x k ) n
The final form of g ( x ) is given as Equation (4).
g ( x ) k = 1 l α k n = 0 a n ( x k ) n + k = 1 l β k n = 0 a n n ( x k ) n 1 .
After combining α k a n as a k n , β k a n as b k n and approximating Equation (4), Equation (5) is obtained. Here, truncation of to M is performed. After truncation, the first derivative of g ( x ) is taken and given in Equation (6).
g ( x ) k = 1 l n = 0 M a k n ( x k ) n + k = 1 l n = 0 M b k n n ( x k ) n 1
d g ( x ) d x k = 1 l n = 1 M a k n n ( x k ) n 1 + k = 1 l n = 1 M b k n n ( n 1 ) ( x k ) n 2 .
The expression given in Equation (7) is the definition of Caputo’s fractional derivative [11]. Throughout the study, Caputo’s description of the fractional derivative is employed.
D x γ g ( x ) = d γ g ( x ) d x γ = 1 Γ ( n γ ) 0 x g ( n ) ( k ) d k ( x k ) γ n + 1 , ( n 1 < γ < n ) .
In Equation (7), Γ ( 1 γ ) is the Gamma function, the fractional derivative is taken with respect to x in the order of γ and g ( n ) corresponds to the n t h derivative again, with respect to   x . In our study, n = 1 is assumed and the fractional-order spans between 0 and 1. Here, two expansions are done to express g ( x ) , approximately. The first one is to express the function as the finite summation of the previous values of the function. Second, expressing the function g ( x ) as the summation of polynomials known as Taylor Expansion assuming that   g ( x ) is a continuous and differentiable function.
Finally, the mathematical background is enough to go further in the proposed methodology. As a summary, above, we mentioned three important tools. First, a function is expressed as the summation of its previous samples. Second, Taylor expansion for a continuous and differentiable function is defined. After that, the Caputo definition of the fractional derivative is given. Now, it is time to express Deep Assessment Methodology by using fractional calculus for the modeling and prediction. Apart from above, there is an assumption that the fractional derivative f ( x ) in the order of γ is equal to Equation (8). After this assumption, it is required to find unknown f ( x ) which satisfies the fractional differential equation below and models the discrete dataset.
d γ f ( x ) d x γ k = 1 l n = 1 a k n n ( x k ) n 1 + k = 1 l n = 1 b k n n ( n 1 ) ( x k ) n 2 ,
where, f ( x ) stands for the GDP per capita of the countries and x corresponds to the time.
Note that, in (6), allowing the order of the derivation in the left-hand side of Equation (6) to be non-integer gives a more general model [28]. This generalization is employed in Deep Assessment Methodology for f ( x ) which stands for the GDP per capita.
Here, the motivation is to find a k n and b k n   given in Equation (8). To find the unknowns, the differential equation needs to be solved. The strategy is as follows—first, it is required to take the Laplace transform which leads to having an algebraic equation instead of a differential equation. In other words, the Laplace transform is taken for Equation (8) to reduce the differential equation to algebraic equation, then, by using inverse Laplace transform properties, the final form of   f ( x ) is obtained as Equation (9) [11].
f ( x , γ ) f ( 0 ) + k = 1 l n = 1 a k n C k n ( x , γ ) + k = 1 l n = 1 b k n D k n ( x , γ ) , where , C k n ( x , γ ) Γ ( n + 1 ) Γ ( n + γ ) ( x k ) n + γ 1 D k n ( x , γ ) Γ ( n + 1 ) Γ ( n + γ 1 ) ( x k ) n + γ 2 .
To obtain the numerical calculation, the infinite summation of polynomials is approximated as a finite summation given in Equation (10).
f ( x , γ ) f ( 0 ) + k = 1 l n = 1 M a k n C k n ( x , γ ) + k = 1 l n = 1 M b k n D k n ( x , γ ) .
Here,   f ( 0 ) ,   a k n and b k n are unknown coefficients that need to be determined. Note that, below, properties of the Laplace transform ( ) are given to find Equations (9) and (10) [11].
[ ( x k ) n 1 ] = Γ ( n ) s n e k s   and   [ d γ f ( x ) d x γ ] = s γ F ( s ) s γ 1 f ( 0 )   for   0 < γ < 1 .
where, stands for the Laplace transform and [ f ( x ) ] = F ( s ) .
For the numerical calculation, the infinite summation is converted into a finite summation, as given in Equation (10).

3. Our Approach

3.1. Modeling with Deep Assessment

In this part, the methodology for the modeling of the problem is given in detail. To predict the upcoming years, the problem has four regions as given in Figure 1. Dataset spans in Region 1, 2 and 3. Note that, there is no data for Region 4 where the prediction is aimed. Region 1 is called “before modeling region” which consists of historical data. Each of the coefficients   ( x k ) n + γ 1 and derivative coming from previous values of GDP per capita for different values of k and multiplication by different weights as given in Equation (10) will add the contribution to the recent data. For modeling, the historical data is employed directly for the modeling of the data located in Region 2. Region 2 and 3 are named as modeling and testing, respectively. In the modeling region, the GDP per capita is tried to be modeled and the unknown coefficients are found. Note that, the approach uses the previous l values ( P i 1 , P i 2 , …, P i l and corresponding f ( i 1 ) ,   f ( i 2 ) ,   f ( i l ) ) for arbitrary P i located in Region 2. The third region consists of the data used to test for upcoming predictions. Finally, Region 4 is called the “prediction region” where the aim is to find the GDP per capita values for the time that the actual values have not known yet and implement prediction. The region division is required because there are parameters given in the previous section (Equation (10)) such as M ,   l ,   γ which need to be found before the prediction. In Region 2, the modeling is done to find the optimum values of coefficients a k n ,   M ,   l ,   γ in Equation (10) for modeling. To model the data, Least Squares Method is employed, which is explained later in this section. After that, one of the purposes of the study is achieved. This is the modeling of the data using the fractional approach. Then, the second purpose comes which is to predict the values of GDP per capita for the upcoming unknown years. In order to find optimum M ,   l ,   γ values for the prediction, Region 3, namely testing is needed. In the region, there is an iterative solution where the real discrete data is again known. For instance, in Region 3, it is required to find f ( m 1 + 1 ) . Then, by using the proposed method employing the fractional calculus and Least Squares Method, f ( m 1 + 1 ) is obtained with a minimum error by optimizing M ,   l ,   γ values for f ( m 1 + 1 ) itself. Then, f ( m 1 + 1 ) is included the dataset for the next test which is done for f ( m 1 + 2 ) . This continues up to   f ( m ) . Then, with optimized   M ,   l ,   γ , the predicted f ( m x ) is found in Region 4.
To model the known data, f ( x ) representing the data optimally should be obtained. In other words, the unknowns   a k n , b k n and   f ( 0 ) in Equation (10) or Equation (11) should be determined. For this, the Least Squares Method is employed.
f ( i , γ ) = f ( 0 ) + k = 1 l n = 1 M a k n C k n ( i , γ ) + k = 1 l n = 1 M b k n D k n ( i , γ ) .
In Equation (12), the squares of total error ϵ T 2 is given. The main purpose of the modeling region is to minimize   ϵ T 2 by a gradient-based approach which requires minimization of the square of the total error as the following.
   ϵ T 2 = i = l m 1 ( P i f ( i , γ ) ) 2
ϵ T 2 f ( 0 ) = 0 ,
ϵ T 2 a r t = 0 ,
and
ϵ T 2 b r t = 0 .
where, r = 1 ,   2 ,   3 ,   l and t = 1 ,   2 ,   3 ,   M .
It is better to give an example of how to obtain ϵ T 2 a r t = 0 and ϵ T 2 b r t = 0 .
ϵ T 2 a r t = 0 a r t i = l m 1 ( P i f ( i , γ ) ) 2 = 0
Then,
2 i = l m 1 [ P i f ( i , γ ) ] C r t ( i , γ ) = 0
i = l m 1 C r t ( i , γ ) P i = f ( 0 ) i = l m 1 C r t ( i , γ ) + i = l m 1 { k = 1 l n = 1 M a k n C k n ( i , γ ) } C r t ( i , γ )
The same procedure is followed for ϵ T 2 b r t = 0 .
ϵ T 2 b r t = 0 b r t i = l m 1 ( P i f ( i , γ ) ) 2 = 0
Then,
i = l m 1 [ P i f ( i , γ ) ] D r t ( i , γ ) = 0
i = l m 1 D r t ( i , γ ) P i = f ( 0 ) i = l m 1 D r t ( i , γ ) + i = l m 1 { k = 1 l n = 1 M a k n D k n ( i , γ ) } D r t ( i , γ )
This leads to having a system of linear algebraic equations (SLAE) as given in (13).
[ A ] . [ B ] = [ C ]
where, [A], [B] and [C] is shown in Equations (14), (15) and (16), respectively.
To find f ( x ) the continuous curve modeling with a minimum error, the optimum fractional-order γ is inquired between (0, 1). Then, with optimum fractional-order γ , the unknown coefficients are determined. In the study, the GDP per capita of Brazil, China, the European Union, India, Italy, Japan, the UK, the USA, Spain and Turkey were used from 1960 until 2018 [38]. The dataset is shown in Table A1 and Table A2.
Among them, the year 2018 is in Region 3 as testing to predict for the next years.
Here,
t (years): [1960, 1961, ..., 2018]
i (points): [1, 2, ..., 59]
P i (value of   i ): [ P 1 , P 2 , , P 59 ]
P i : It shows the actual GDP per capita of each country in each i t h year. For example, P 2 is the GDP per capita of the country in 1961.
i : It stands for the number for each year. For example,   i = 1 for 1960, i = 3 for 1962 and i = 59 for 2018.
A = [ A 1 , 1 A 1 , 2 A 2 , 1 A 2 , 2 ]
A matrix consists of the matrix set below, where,
C k n ( x , γ ) = C k n   and   D k n ( x , γ ) = D k n A 1 , 1 = m 1 l + 1 i = l m 1 C 11 i = l m 1 C 1 M i = l m 1 C 21 i = l m 1 C 2 M i = l m 1 C l 1 i = l m 1 C l M C 11 i = l m 1 C 11 C 11 i = l m 1 C 1 M C 11 i = l m 1 C 21 C 11 i = l m 1 C 2 M C 11 i = l m 1 C l 1 C 11 i = l m 1 C l M C 11 C 12 i = l m 1 C 11 C 12 i = l m 1 C 1 M C 12 i = l m 1 C 21 C 12 i = l m 1 C 2 M C 12 i = l m 1 C 11 C 12 i = l m 1 C l M C 12 C l m i = l m 1 C 11 C l M i = l m 1 C 1 M C l M i = l m 1 C 21 C l M i = l m 1 C 2 M C l M i = l m 1 C l 1 C l M i = l m 1 C l M C l M A 2 , 1 = D 11 i = l m 1 C 11 D 11 i = l m 1 C 1 M D 11 i = l m 1 C 21 D 11 i = l m 1 C 2 M D 11 i = l m 1 C l 1 D 11 i = l m 1 C l M D 11 i = l m 1 D 12 i = l m 1 C 11 D 12 i = l m 1 C 1 M D 12 i = l m 1 C 21 D 12 i = l m 1 C 2 M D 12 i = l m 1 C l 1 D 12 i = l m 1 C l M D 12 i = l m 1 D l M i = l m 1 C 11 D l M i = l m 1 C 1 M D l M i = l m 1 C 21 D l M i = l m 1 C 2 M D 1 M i = l m 1 C l 1 D l M i = l m 1 C l M D l M A 1 , 2 = i = l m 1 D 11 i = l m 1 D 1 M i = l m 1 D l M i = l m 1 D 11 C 11 i = l m 1 D 1 M C 11 i = l m 1 D l M C 11 i = l m 1 D 11 C 12 i = l m 1 D 1 M C 12 i = l m 1 D l M C 12 i = l m 1 D 11 C l M i = l m 1 D 1 M C l M i = l m 1 D l M C l M A 2 , 2 = i = l m 1 D 11 D 11 i = l m 1 D 1 M D 11 i = l m 1 D l M D 11 i = l m 1 D 11 D 12 i = l m 1 D 1 M D 12 i = l m 1 D l M D 12 i = l m 1 D 11 D l M i = l m 1 D 1 M D l M i = l m 1 D l M D l M
[ B ] = [ f ( 0 ) a 11 a 12 a 1 M   a 21 a 22     a 2 M a l 1 a l M   b 11 b 12   b 1 M b 21   b 2 M b l 1 b l 2 b l M ] T
[ C ] = [ i = l m 1 P i i = l m 1 P i C 11 i = l m 1 P i C 12 i = l m 1 P i C l M     i = l m 1 P i D 11 i = l m 1 P i D 12 i = l m 1 P i D l M ] T

3.2. Prediction with Deep Assessment

To find the optimized values of the unknowns for the prediction, the testing region (3rd region) is required. The predictions obtained in the test region ( m 1 < i < m ) are also given in Table 1. For testing, the data up to m 1 = 58 have been taken into consideration in the operations. The f ( m 1 + 1 ) value was found from the obtained modeling. Then, the value is kept, and the next step was started again for ( f ( m 1 + 2 ) ) . These operations are done until the last value of the test zone. In our case, m = 59 .
The last region is called the “Prediction Region.” Here, using Region 1, 2 and 3, the prediction for the upcoming years is obtained. After having modeled and tested regions, the unknowns in Equation (11) have already found in an optimal manner. After testing, Region 4 is started. In the region, the first prediction f ( m + 1 ) is found by using the coefficients and unknowns found by the testing region. After that, the first predicted value ( f ( m + 1 ) ) is included in Region 3 (testing) for the consecutive prediction f ( m + 2 ) . This procedure is reiterated and recycled up to   f ( m x ) .
The prediction results for 2019 are given in Table 2. For example, as of the end of 2019 ( f ( m + 2 ) ), Brazil, China, European Union, India, Italy, Japan, the UK, the USA, Spain and Turkey’s GDP per capita values are expected as listed.
In Figure 2, the algorithm for prediction with DAM is illustrated. The first step of the algorithm is to initialize the parameters ( l ,   M ,   x 1 ,   x 2 ,   x m and   P 1 ,   P 2 ,     P m ). Then, the counter variable N is introduced, which counts the number of prediction steps. The total number of required predicted steps is denoted as n 0 . As an initial value, the fractional-order γ is assigned 0 and the increment is 0.01 for each loop to find the optimized value. For each value of γ between 0 and 1, matrix A given as Equation (14) is created and then, the unknown coefficients given in Equation (10) are calculated. After that, using the actual data in Region 1 and Region 2, the modeling of data between P l and P m is actualized for Region 2. Then, the error defined in Equation (12) is calculated. The value of the error is analyzed and compared to previously obtained values. If it is smaller than the previous one, the corresponding fractional-order value is memorized. At the end of Loop II, the optimal value of the fractional-order, which coincides with the optimal modeling is found and corresponding coefficients given in Equation (10) is determined. Then, the prediction for the next forthcoming value is made with Equation (10). After that, all the procedures starting from the increment of N is repeated so that the previously predicted value is added to the initial data for the next step prediction. This process is repeated up to the termination of Loop I. Finally, n 0 the number of predictions is obtained. Keep in mind that, for the parameters l and M , there exist two loops starting from 1 to L 0 and 1 to M 0 searching the optimum values of the parameters in order to get the outcomes with a minimum error for the testing region, respectively. Here, L 0 and M 0 are pre-defined some constant values.

3.3. Long Short-Term Memory

In our study, we compare the modeling with the polynomial curve fitting method and in the prediction, we compare Deep Assessment with the LSTM method. Conventional neural networks are insufficient for modeling the content of temporal data. Recursive neural networks (RNN) model the sequential structure of data by feeding itself with the output of the previous time step. LSTMs are special types of RNNs that operate over sequences and are used in time series analysis [39]. An LSTM cell has four gates: input, forget, output and gate. With these gates, LSTMs optionally inherit the information from the previous time steps. Forget gate   ( f ) , input gate ( i ) and output gate ( o ) are sigmoid functions ( σ ) and they take values between 0 and 1. Gate g has hyperbolic tangent ( t a n h ) activation and is between −1 and 1. The Gate and forward propagation equations are listed below as Equations (17)–(22). Here c t l and h t l refer to cell state and hidden state of layer l at time step t , respectively. Each gate takes input from the previous time step ( h t 1 l ) and previous layer ( h t l 1 ) and has its own set of learnable parameters W ’s and b ’s.
f t = σ ( W f [ h t 1 l , h t l 1 ] + b f )
i t = σ ( W i [ h t 1 l , h t l 1 ] + b i )
o t = σ ( W o [ h t 1 l , h t l 1 ] + b o )
g t = t a n h ( W g [ h t 1 l , h t l 1 ] + b g )
c t l = f   c t 1 l + i   g
h t l = o   tanh ( c t l )
Here, is the Hadamard product. Each LSTM neuron in a network may consist of one or more cells. In every time step, every cell updates its own cell state,     c t l . Equation (22) describes how these cells get updated with forget gate and input gate; f gate decides how much of previous cell state that cell should remember while i gate decides how much it should consider the new input from the previous layer. Then, LSTM neuron updates its internal hidden state by multiplying output and squashed version of c t l . An LSTM neuron gives outputs only in hidden state information to another LSTM neuron. Gate o and c t are used internally in the computation of forward time steps [40]. To forecast time series and compare our proposed approach to neural networks, we employed a stacked LSTM model with 2 layers of LSTMs (each having 50 hidden units) and a linear prediction layer. LSTM model is trained with the Adam optimizer [40].

4. Numerical Results

In this section, we report the modeling and prediction performance of the Deep Assessment methodology. Further, we compare the proposed method to other modeling and prediction approaches such as Polynomial Model, Fractional Model-1 [34,35] and LSTM. In this section, results are reported with the Mean Average Precision Error (MAPE) metric and calculated as follows:
M A P E = 1 k i = 1 k | v ( i ) v ~ ( i ) v ( i ) | × 100 ,
where k is the total number of samples, v ( i ) is the actual value and v ~ ( i )   is the predicted value for i th sample.
Before presenting the results, it is important to highlight that for modeling, M 0 and l 0 are taken 20 and 10, respectively whereas for prediction, M 0 and l 0 are taken 8 and 25, respectively. The number of prediction, n 0 is equal to 1.

4.1. Modeling Results

In this part, we compare the modeling performance with Polynomial, Fractional Model-1 and Deep Assessment models.
To achieve modeling, l value needs to be investigated. For the modeling of the GDP per capita of each country, the required previous data l of past years used in the algorithm differs after optimization. In order to make a fair evaluation, l value is fixed among all countries to 10. Modeling results for Deep Assessment, Polynomial Model and Fractional Model-1 are shown in Table 1. Optimized M values after processing can be seen in the last column. The Deep Assessment model has a %4.308 average MAPE and outperforms Polynomial and Fractional Model-1 by %1.538 and %1.899 average error rates. All three methods model the US best with %0.81, %1.01 and %1.06 error. Further, in the case of Italy, Fractional Model-1 uses the fractional-order value of 1 and produces %8.81 MAPE, equal to the Polynomial method as expected because for the fractional-order value of 1 is the same with the Polynomial Method. However, DAM yields fractional order of 0.39, decreasing the error to 4.70%, justifying the advantage of employing fractional calculus and previous values of the data itself.
M A P E M o d e l i n g = 1 m l + 1 i = l m | P ( i ) f ( i , γ ) P ( i ) | × 100 .
GDP per capita data, Deep Assessment, Polynomial and Fractional Model-1 modeling results are shown in Figure 1 for each country. One can conclude that when data points have high variance all models produce high error rates, as in Turkey and Italy. For Japan and Brazil, DAM (Deep Assessment Method) and Polynomial models produce similar results. Also, it can be seen from the Figure 3, both Deep Assessment and Fractional Model-1 have a low bias when compared to the Polynomial model and overfits to dataset less. This is possible because of the memory property of the proposed approach. Except for Brazil, India and the EU, the proposed method yields superior results compared to other models.

4.2. Prediction Results

In this section, we compare the accuracy rate of the prediction of Deep Assessment and Deep Learning models. As in modeling, the GDP per capita dataset is used to assess the performance of the proposed method. Table 2 illustrates optimized   γ , l , M values and the corresponding performance of DAM and LSTM. Here, column 6 reports the performance of DAM while column 7 represents LSTM. Column 5 shows that the Deep Assessment methodology predicts GDP per capita with an average 0.29% error with predicting all countries with 1.< (less than 1 percent) of error. The best-predicted country is Spain while UK’s prediction is the least accurate with 0.91% error. On the other hand, LSTM yields 1.51% error on average. For both DAM and LSTM, UK yields the highest error. Table 2 demonstrates that in the implemented setting, DAM outperforms LSTM by 1.21% average error and produces fair results.
M A P E P r e d i c t i o n = 1 m 1 l + 1 i = m 1 m | P ( i ) f ( i , γ ) P ( i ) | × 100 .
Table 3 reports the prediction of GDP per capita for the year 2019 is illustrated in Table 2 for both DAM and LSTM methods. For countries Brazil, China, India, Turkey, the UK and the US, predictions obtained by the two models are similar. On the other hand, Italy and Spain yield different results.

5. Conclusions

In this study, a model called “Deep Assessment” is introduced which employs Fractional Calculus to model discrete data as the summation of previous values and derivatives. Different to the literature and our previous work, the proposed approach also predicts the incoming values of the discrete data in addition to modeling. The method is evaluated on modeling and predicting GDP per capita, using a dataset including the period of 1960–2018 for nine countries (Brazil, China, European Union, India, Italy, Japan, UK, the USA, Spain and Turkey) and the European Union. Using the fractional differential equation and the summation of previous values for the modeling of GDP per capita at a specific time instant bring non-locality, memory and generalization of the problem for different fractional order. In experiments, first, GDP per capita is modeled. The Deep Assessment model has a 4.308% average MAPE and outperforms Polynomial and Fractional Model-1 by 1.538% and 1.899% average error rates for modeling. For prediction, LSTM, a special type of neural network is used to assess the performance of the model. In the selected test region, it is shown that Deep Assessment is superior to LSTM by 1.51% average error. Results illustrate that the proposed method yields promising results and demonstrates the benefits of combining fractional calculus and differential equations. Evaluation of multivariable and multifunctional problems, analyzing time windows, randomness, noise and error changes are left to future work.

Author Contributions

The contribution of each author is listed as follows. E.K. has contributed to supervision, conceptualization, investigation, methodology, and administration. V.T. plays an important role in resources, supervision, and validation. K.K. supported conceptualization, writing, and editing. N.Ö.Ö. was the key person about visualization, investigation, administration, validation, and writing. E.E. has contributed to validation, visualization, writing, and editing. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Istanbul Technical University (ITU) Vodafone Future Lab with the grant number ITUVF20180901P11.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses or interpretation of data; in the writing of the manuscript or in the decision to publish the results.

Appendix A

Table A1. GDP per capita (US dollars) values of the countries.
Table A1. GDP per capita (US dollars) values of the countries.
i YearsBrazilChinaEUIndiaItaly
11960210.109989.52054890.405682.1886804.4926
21961205.040875.80584959.7185.3543887.3367
31962260.425770.909411037.32689.88176990.2602
41963292.252174.313641135.194101.12641126.019
51964261.666685.498561245.499115.53751222.545
61965261.354498.486781346.058119.31891304.454
71966315.7972104.32461448.55189.997311402.442
81967347.493196.589531546.80496.339141533.693
91968374.786891.472721602.0699.875961651.939
101969403.8843100.12991762.472107.62231813.388
111970445.0231113.1631950.732112.43452106.864
121971504.7495118.65462195.145118.60322305.61
131972586.2144131.88362611.729122.98192671.137
141973775.2733157.09043296.935143.77873205.252
1519741004.105160.14013685.596163.47813621.146
1619751153.831178.34184274.046158.03624106.994
1719761390.625165.40554406.238161.09214033.099
1819771567.006185.42284968.988186.21354603.6
1919781744.257156.39646064.883205.69345610.498
2019791908.488183.98327377.165224.0016990.286
2119801947.276194.80478384.718266.57788456.919
2219812132.883197.07157391.077270.47067622.833
2319822226.767203.33497093.702274.11137556.523
2419831570.54225.43196859.966291.23817832.575
2519841578.926250.7146572.019276.6687739.715
2619851648.082294.45886775.647296.43527990.687
2719861941.491281.92819265.924310.465911,315.02
2819872087.308251.81211,432.23340.416814,234.73
2919882300.377283.537712,711.96354.149315,744.66
3019892908.496310.881912,936.46346.112916,386.66
3119903100.28317.884715,989.22367.556620,825.78
3219913975.39333.142116,496.51303.055621,956.53
3319922596.92366.460717,919.02316.953923,243.47
3419932791.209377.389816,256.42301.15918,738.76
3519943500.611473.492317,194.12346.10319,337.63
3619954748.216609.656719,898.44373.766520,664.55
3719965166.164709.413820,295.17399.950123,081.6
3819975282.009781.744219,121.21415.493821,829.35
3919985087.152828.580519,763.51413.298922,318.14
4019993478.373873.287119,698.89441.998821,997.62
4120003749.753959.372518,261.97443.314220,087.59
4220013156.7991053.10818,457.89451.57320,483.22
4320022829.2831148.50820,055.33470.986822,270.14
4420033070.911288.64324,310.25546.726627,465.68
4520043637.4621508.66827,960.05627.774231,259.72
4620054790.4371753.41829,115.63714.86132,043.14
4720065886.4642099.22930,960.56806.753333,501.66
4820077348.0312693.9735,630.941028.33537,822.67
4920088831.0233468.30438,185.62998.522340,778.34
5020098597.9153832.23634,019.281101.96137,079.76
51201011,286.244550.45433,740.651357.56436,000.52
52201113,245.615618.13236,506.641458.10438,599.06
53201212,370.026316.91934,328.821443.8835,053.53
54201312,300.327050.64635,683.861449.60635,549.97
55201412,112.597651.36636,787.231573.88135,518.42
5620158814.0018033.38832,319.451605.60530,230.23
5720168712.8878078.7932,425.131729.26830,936.13
5820179880.9478759.04233,9081981.26932,326.84
5920188920.7629770.84736,569.732009.97934,483.2
Table A2. GDP per capita (US dollars) values of the countries.
Table A2. GDP per capita (US dollars) values of the countries.
i YearsJapanSpainUKUSTurkey
11960478.9953396.39231397.5953007.123509.4239
21961563.5868450.05331472.3863066.563283.8283
31962633.6403520.20611525.7763243.843309.4467
41963717.8669609.48741613.4573374.515350.6629
51964835.6573675.24161748.2883573.941369.5834
61965919.7767774.76161873.5683827.527386.3581
719661058.504889.65991986.7474146.317444.5494
819671228.909968.30682058.7824336.427481.6937
919681450.62950.54571951.7594695.923526.2135
1019691669.0981077.6792100.6685032.145571.6178
1119702037.561212.2892347.5445234.297489.9303
1219712272.0781362.1662649.8025609.383455.1049
1319722967.0421708.8093030.4336094.018558.421
1419733997.8412247.5533426.2766726.359686.4899
1519744353.8242749.9253665.8637225.691927.7991
1619754659.123209.8374299.7467801.4571136.375
1719765197.8073279.3134138.1688592.2541275.956
1819776335.7883627.5914681.449452.5771427.372
1919788821.8434356.4395976.93810,564.951549.644
2019799105.1365770.2157804.76211,674.192079.22
2119809465.386208.57810,032.0612,574.791564.247
22198110,361.325371.1669599.30613,976.111579.074
2319829578.1145159.7099146.07714,433.791402.406
24198310,425.414478.58691.51915,543.891310.256
25198410,984.874489.9898179.19417,121.231246.825
26198511,584.654699.6568652.21718,236.831368.401
27198617,111.856513.50310,611.1119,071.231510.677
28198720,745.258239.61413,118.5920,038.941705.895
29198825,051.859703.12415,987.1721,417.011745.365
30198924,813.310,681.9716,239.2822,857.152021.859
31199025,359.3513,804.8819,095.4723,888.62794.35
32199128,925.0414,811.919,900.7324,342.262735.708
33199231,464.5516,112.1920,487.1725,418.992842.37
34199335,765.9113,339.9118,389.0226,387.293180.188
35199439,268.5713,415.2919,709.2427,694.852270.338
36199543,440.3715,471.9623,123.1828,690.882897.866
37199638,436.9316,109.0824,332.729,967.713053.947
38199735,021.7214,730.826,734.5631,459.143144.386
39199831,902.7715,394.3528,214.2732,853.684496.497
40199936,026.5615,715.3328,669.5434,513.564108.123
41200038,532.0414,713.0728,149.8736,334.914316.549
42200133,846.4715,355.727,744.5137,133.243119.566
43200232,289.3517,025.5330,056.5938,023.163659.94
44200334,808.3921,463.4434,419.1539,496.494718.2
45200437,688.7224,861.2840,290.3141,712.86040.608
46200537,217.6526,419.342,030.2944,114.757384.252
47200635,433.9928,365.3144,599.746,298.738035.377
48200735,275.2332,549.9750,566.8347,975.979711.874
49200839,339.335,366.2647,28748,382.5610,854.17
50200940,855.1832,042.4738,713.1447,099.989038.52
51201044,507.6830,502.7239,435.8448,466.8210,672.39
52201148,16831,636.4542,038.549,883.1111,335.51
53201248,603.4828,324.4342,462.7151,603.511,707.26
54201340,454.4529,059.5543,444.5653,106.9112,519.39
55201438,109.4129,461.5547,417.6455,032.9612,095.85
56201534,524.4725,732.0244,966.156,803.4710,948.72
57201638,794.3326,505.6241,074.1757,904.210,820.63
58201738,331.9828,100.8540,361.4259,927.9310,513.65
59201839,289.9630,370.8942,943.962,794.599370.176

References

  1. Beaudry, P.; Portier, F. Stock prices, news, and economic fluctuations. Am. Econ. Rev. 2006, 96, 1293–1307. [Google Scholar] [CrossRef] [Green Version]
  2. Simpson, D.M.; Rockaway, T.D.; Weigel, T.A.; Coomes, P.A.; Holloman, C.O. Framing a new approach to critical infrastructure modelling and extreme events. Int. J. Crit. Infrastruct. 2005, 1, 125–143. [Google Scholar] [CrossRef]
  3. Ho, T.Q.; Strong, N.; Walker, M. Modelling analysts’ target price revisions following good and bad news? Account. Bus. Res. 2018, 48, 37–61. [Google Scholar] [CrossRef]
  4. Karaçuha, E.; Tabatadze, V.; Önal, N.Ö.; Karaçuha, K.; Bodur, D. Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Population of Countries. Authorea 2020. [Google Scholar] [CrossRef]
  5. Forestier, E.; Grace, J.; Kenny, C. Can information and communication technologies be pro-poor? Telecommun. Policy 2002, 26, 623–646. [Google Scholar] [CrossRef] [Green Version]
  6. Kretschmer, T. Information and Communication Technologies and Productivity Growth: A Survey of the Literature; OECD Digital Economy Papers; OECD Publishing:: Paris, France, 2012; Volume 195. [Google Scholar] [CrossRef]
  7. Han, E.J.; Sohn, S.Y. Technological convergence in standards for information and communication technologies. Technol. Forecast. Soc. Chang. 2016, 106, 1–10. [Google Scholar] [CrossRef]
  8. Grover, V.; Purvis, R.L.; Segars, A.H. Exploring ambidextrous innovation tendencies in the adoption of telecommunications technologies. IEEE Trans. Eng. Manag. 2007, 54, 268–285. [Google Scholar] [CrossRef]
  9. Torero, M.; Von Braun, J. (Eds.) Information and Communication Technologies for Development and Poverty Reduction: The Potential of Telecommunications; International Food Policy Research Institute: Washington, DC, USA, 2006. [Google Scholar]
  10. Pradhan, R.P.; Arvin, M.B.; Norman, N.R. The dynamics of information and communications technologies infrastructure, economic growth, and financial development: Evidence from Asian countries. Technol. Soc. 2015, 42, 135–149. [Google Scholar] [CrossRef]
  11. Podlubny, I. Fractional differential equations. Math. Sci. Eng. 1999, 198, 79. [Google Scholar]
  12. Machado, J.T.; Kiryakova, V.; Mainardi, F. Recent history of fractional calculus. Commun. Nonlinear Sci. Numer. Simul. 2019, 16, 1140–1153. [Google Scholar] [CrossRef] [Green Version]
  13. Tarasov, V.E. On history of mathematical economics: Application of fractional calculus. Mathematics 2019, 7, 509. [Google Scholar] [CrossRef] [Green Version]
  14. Baleanu, D.; Güvenç, Z.B.; Machado, J.T. (Eds.) New Trends in Nanotechnology and Fractional Calculus Applications; Springer: New York, NY, USA, 2010. [Google Scholar]
  15. Gómez, F.; Bernal, J.; Rosales, J.; Cordova, T. Modeling and simulation of equivalent circuits in description of biological systems-a fractional calculus approach. J. Electr. Bioimpedance 2019, 3, 2–11. [Google Scholar] [CrossRef] [Green Version]
  16. Machado, J.T.; Mata, M.E. Pseudo phase plane and fractional calculus modeling of western global economic downturn. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 396–406. [Google Scholar] [CrossRef] [Green Version]
  17. Moreles, M.A.; Lainez, R. Mathematical modelling of fractional order circuit elements and bioimpedance applications. Commun. Nonlinear Sci. Numer. Simul. 2017, 46, 81–88. [Google Scholar] [CrossRef]
  18. Tabatadze, V.; Karaçuha, K.; Veliev, E.I. The Fractional Derivative Approach for the Diffraction Problems: Plane Wave Diffraction by Two Strips with the Fractional Boundary Conditions. Prog. Electromagn. Res. 2019, 95, 251–264. [Google Scholar] [CrossRef] [Green Version]
  19. Sierociuk, D.; Dzieliński, A.; Sarwas, G.; Petras, I.; Podlubny, I.; Skovranek, T. Modelling heat transfer in heterogeneous media using fractional calculus. Philosophical Trans. R. Soc. A 2013, 371, 20120146. [Google Scholar] [CrossRef] [Green Version]
  20. Bogdan, P.; Jain, S.; Goyal, K.; Marculescu, R. Implantable pacemakers control and optimization via fractional calculus approaches: A cyber-physical systems perspective. In Proceedings of the Third International Conference on Cyber-Physical Systems, Beijing, China, 17–19 April 2012. [Google Scholar]
  21. Magin, R.L. Fractional calculus models of complex dynamics in biological tissues. Comput. Math. Appl. 2010, 59, 1586–1593. [Google Scholar] [CrossRef] [Green Version]
  22. Chen, Y.; Xue, D.; Dou, H. Fractional calculus and biomimetic control. In Proceedings of the 2004 IEEE International Conference on Robotics and Biomimetics, Shenyang, China, 22–26 August 2004; pp. 901–906. [Google Scholar]
  23. Dimeas, I.; Petras, I.; Psychalinos, C. New analog implementation technique for fractional-order controller: A DC motor control. Aeu-Int. J. Electron. Commun. 2004, 78, 192–200. [Google Scholar] [CrossRef]
  24. Sierociuk, D.; Skovranek, T.; Macias, M.; Podlubny, I.; Petras, I.; Dzielinski, A.; Ziubinski, P. Diffusion process modeling by using fractional-order models. Appl. Math. Comput. 2015, 257, 2–11. [Google Scholar] [CrossRef] [Green Version]
  25. Tarasova, V.; Tarasov, V. Elasticity for economic processes with memory: Fractional differential calculus approach. Fract. Differ. Calc. 2016, 6, 219–232. [Google Scholar] [CrossRef]
  26. Tarasov, V.; Tarasova, V. Macroeconomic models with long dynamic memory: Fractional calculus approach. Appl. Math. Comput. 2018, 338, 466–486. [Google Scholar] [CrossRef]
  27. Tejado, I.; Pérez, E.; Valério, D. Fractional calculus in economic growth modelling of the group of seven. Fract. Calc. Appl. Anal. 2019, 22, 139–157. [Google Scholar] [CrossRef]
  28. Škovránek, T.; Podlubny, I.; Petráš, I. Modeling of the national economies in state-space: A fractional calculus approach. Econ. Model. 2012, 29, 1322–1327. [Google Scholar] [CrossRef]
  29. Inés, T.; Valério, D.; Pérez, E.; Valério, N. Fractional calculus in economic growth modelling: The Spanish and Portuguese cases. Int. J. Dyn. Control 2017, 5, 208–222. [Google Scholar]
  30. Ming, H.; Wang, J.; Fečkan, M. The application of fractional calculus in Chinese economic growth models. Mathematics 2019, 7, 665. [Google Scholar] [CrossRef] [Green Version]
  31. Blackledge, J.; Kearney, D.; Lamphiere, M.; Rani, R.; Walsh, P. Econophysics and Fractional Calculus: Einstein’s Evolution Equation, the Fractal Market Hypothesis, Trend Analysis and Future Price Prediction. Mathematics 2019, 7, 1057. [Google Scholar] [CrossRef] [Green Version]
  32. Despotovic, V.; Skovranek, T.; Peric, Z. One-parameter fractional linear prediction. Comput. Electr. Eng. 2018, 69, 158–170. [Google Scholar] [CrossRef]
  33. Skovranek, T.; Despotovic, V.; Peric, Z. Optimal fractional linear prediction with restricted memory. Ieee Signal Process. Lett. 2019, 26, 760–764. [Google Scholar] [CrossRef]
  34. Önal, N.Ö.; Karaçuha, K.; Erdinç, G.H.; Karaçuha, B.B.; Karaçuha, E. A Mathematical Approach with Fractional Calculus for the Modelling of Children’s Physical Development. Comput. Math. Methods Med. 2019, 2019, 3081264. [Google Scholar] [CrossRef] [Green Version]
  35. Önal, N.Ö.; Karaçuha, K.; Erdinç, G.H.; Karaçuha, B.B.; Karaçuha, E. A mathematical model approach regarding the children’s height development with fractional calculus. Int. J. Biomed. Biol. Eng. 2019, 13, 252–260. [Google Scholar]
  36. Önal, N.Ö.; Karaçuha, K.; Karaçuha, E. A Comparison of Fractional and Polynomial Models: Modelling on Number of Subscribers in the Turkish Mobile Telecommunications Market. Int. J. Appl. Phys. Math. 2019, 10, 41. [Google Scholar] [CrossRef] [Green Version]
  37. Maddison, A. A comparison of levels of GDP per capita in developed and developing countries, 1700–1980. J. Econ. Hist. 1983, 43, 27–41. [Google Scholar] [CrossRef]
  38. World Bank Databank World Development Indicators. Available online: https://databank.worldbank.org/source/world-development-indicators# (accessed on 10 December 2019).
  39. Hochreiter, S.; Schmidhuber, J. Long short-term memory. Neural Comput. 1997, 9, 1735–1780. [Google Scholar] [CrossRef] [PubMed]
  40. Kingma, D.P.; Ba, J. Adam. A method for stochastic optimization. arXiv 2014, arXiv:1412.6980. [Google Scholar]
Figure 1. The regions of the dataset.
Figure 1. The regions of the dataset.
Mathematics 08 00633 g001
Figure 2. The algorithm for the prediction.
Figure 2. The algorithm for the prediction.
Mathematics 08 00633 g002
Figure 3. Modelling results of the countries (Brazil, China, India, Italy, Japan, the UK, the USA, Spain and Turkey) and the European Union or Deep Assessment (Blue), Fractional Model 1 (Yellow), Polynomial Model (Purple).
Figure 3. Modelling results of the countries (Brazil, China, India, Italy, Japan, the UK, the USA, Spain and Turkey) and the European Union or Deep Assessment (Blue), Fractional Model 1 (Yellow), Polynomial Model (Purple).
Mathematics 08 00633 g003aMathematics 08 00633 g003b
Table 1. Comparison of modeling results ( γ , M and MAPE values) of countries for l = 10 .
Table 1. Comparison of modeling results ( γ , M and MAPE values) of countries for l = 10 .
Country γ
Deep Assessment
γ
Fractional Model-1
Deep Assessment *
( l < i < m )
Polynomial Model *
( l < i < m )
Fractional Model-1 *
( l < i < m )
M
US0.440.540.81%1.01%1.06%15
UK0.140.855.38%7.03%6.61%15
Brazil0.060.587.26%7.13%9.00%17
China0.030.952.84%5.62%5.67%11
India0.150.023.09%2.51%4.10%16
Japan0.260.694.45%4.64%5.82%20
EU0.060.894.02% 3.41%5.71%20
Italy0.3914.70%8.81%8.81%9
Spain0.220.584.44%6.49%6.36%13
Turkey0.710.016.09%11.81%8.93%10
* M A P E M o d e l i n g values.
Table 2. Test ( m 1 < i < m ) results ( γ , l , M and MAPE) of GDP per capita for corresponding countries.
Table 2. Test ( m 1 < i < m ) results ( γ , l , M and MAPE) of GDP per capita for corresponding countries.
Country γ l M γ
Interpolation
Deep Assessment *Deep Learning *
Brazil0.182430.320.1303%0.4728%
China0.971130.50.7147%1.6365%
India0.96320.990.3379%
5
0.7203%
Italy0.432040.430.1048%3.0796%
Japan0.574310.3499%1.1091%
Spain0.99230.990.0560%1.5683%
Turkey0.391740.390.1167%2.3691%
EU0.322050.220.1044%0.2522%
US0.392520.180.1081%0.8424%
UK0.181870.050.9129%3.0508%
* M A P E P r e d i c t i o n values.
Table 3. GDP per Capita Prediction of Countries for 2019 (US dollars).
Table 3. GDP per Capita Prediction of Countries for 2019 (US dollars).
CountryDeep AssessmentDeep Learning
Brazil79328013
China10,31210,273
India21541967
Italy39,02835,141
Japan34,42137,994
Spain30,38535,372
Turkey82608920
US65,76763,844
UK44,89744,702
EU40,48736,487

Share and Cite

MDPI and ACS Style

Karaçuha, E.; Tabatadze, V.; Karaçuha, K.; Önal, N.Ö.; Ergün, E. Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Gross Domestic Product per Capita of Countries. Mathematics 2020, 8, 633. https://doi.org/10.3390/math8040633

AMA Style

Karaçuha E, Tabatadze V, Karaçuha K, Önal NÖ, Ergün E. Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Gross Domestic Product per Capita of Countries. Mathematics. 2020; 8(4):633. https://doi.org/10.3390/math8040633

Chicago/Turabian Style

Karaçuha, Ertuğrul, Vasil Tabatadze, Kamil Karaçuha, Nisa Özge Önal, and Esra Ergün. 2020. "Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Gross Domestic Product per Capita of Countries" Mathematics 8, no. 4: 633. https://doi.org/10.3390/math8040633

APA Style

Karaçuha, E., Tabatadze, V., Karaçuha, K., Önal, N. Ö., & Ergün, E. (2020). Deep Assessment Methodology Using Fractional Calculus on Mathematical Modeling and Prediction of Gross Domestic Product per Capita of Countries. Mathematics, 8(4), 633. https://doi.org/10.3390/math8040633

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