Modeling and Application of a New Nonlinear Fractional Financial Model

The paper proposes a new nonlinear dynamic econometric model with fractional derivative. The fractional derivative is defined in the Jumarie type.The corresponding discrete financial system is considered by removing the limit operation in Jumarie derivative’s. We estimate the coefficients and parameters of themodel by using the least squared principle.The new approach to financial system modeling is illustrated by an application to model the behavior of Japanese national financial system which consists of interest rate, investment, and inflation. The empirical results with different time step sizes of discretization are shown, and a comparison of the actual data against the data estimated by empirical model is illustrated. We find that our discrete financial model can describe the actual data that include interest rate, investment, and inflation accurately.


Introduction
The financial and economical systems exhibit ubiquitous complex dynamics evidenced by large-amplitude and periodic fluctuations, which have attracted much attention recently [1][2][3][4]. However, studying the dynamics of finance behavior and economics life is very challenging mainly because of the fact that it is very difficult to identify the inner relationship between different financial and economic variables quantitatively. In the recent decade, several nonlinear models have been proposed to study periodic or chaotic behaviors in financial and economic systems since it is believed that finance and economics always illustrate nonlinearity. In particular, the complex dynamics of economic cycles using the van der Pol model is studied in [5][6][7]. The advantage of employing forced van der Pol equation to depict economic variable is that the introduction of a forcing function enables us to model the complex interdependence between an individual national economy and the international economy in an era of globalization and the impact of climate variations such as the annual solar cycle on seasonable fluctuations of various markets. In [8,9], the problem of the existence of a finite lag between the accrual and the payment of taxes in a framework, where this type of lag has never been considered by the well-known IS-LM model, is proposed. The qualitative study of the system of functional differential equations shows that the finite lag may give rise to a wide variety of dynamic behaviors. These works, together with many recent papers, lead the studying of the nonlinear dynamics of financial and economic systems to be more attractive. Most of the studies about nonlinear dynamics of financial and economic systems concentrate on the periodic or chaotic behaviors in it. References [10,11] propose a simplified macrofinancial model which consists of interest rate, investment demand, and inflation based on the economic theory and analyzes the balance, stable periodic, chaotic motion, and so forth. In [12], the complex motion in nonlinear dynamical systems by reconsidering the Goodwin's nonlinear accelerator model with periodic investment outlays is studied, and it is believed that transient chaotic dynamics are prevalent in nonlinear economic model. In general, since the chaotic dynamics has adverse influences on predicting economy reasonably and effectively, more and more related papers pay attention to analyze and control the chaotic phenomenon of nonlinear dynamics in financial and economic systems. However, in these above works, only integer order differential equations are investigated. In recent years, modeling the real-world problems by using fractional differential equations becomes more and more popular. Fractional differential equation means that the order of derivative in differential equation can be noninteger and even complex number. As an excellent methodology of modeling, fractional calculus gained considerable development in the recent forty years. It has been applied to many scientific and engineering fields and verified to be a powerful tool in modeling most physical processes with memory effect, which cannot be described well by integer-order integral and differential equations. For a comprehensive review of theory and applications of fractional calculus, we refer to [13][14][15][16][17]. The advantage of fractional derivative is that it has nonlocal property. The value of current state depends on both recent values and historical values of objective function. This excellent property is suitable for modeling many financial variable series, mainly because of the fact that the financial and economic variable series always exhibit time-dependent memory effect, such as interest rate, stock price, exchange amount of future, and so on [18][19][20][21]. There are only a few papers on the study of financial and economic behaviors by using fractional derivative models. In [22], a fractional order financial model based on the fractional Chen system is proposed. It involves the macroeconomic variables such as investment, interest, and price index and exhibits the interesting dynamics behavior of them. In [23,24], the chaos dynamics and chaos control have been studied using sliding mode method and feedback control method, respectively. In [25], a delayed fractional-order financial system is proposed and the complex dynamical behaviors of this system are discussed by numerical simulations. A great variety of interesting dynamical behaviors of such a system including single-periodic, multiple-periodic, and chaotic motions are shown. The effect of time delay and fractional order on the chaotic behavior are investigated; it is verified that a proper time delay can enhance or suppress the emergence of chaos.
In most of the recent researches, the dynamics of fractional-order financial and economic systems has been investigated via several mathematical methods (see [22,25] and references therein). However, the physical and natural interpretations of those dynamical phenomena obtained in numerical simulations are still not clear. We have known that the occurrence of chaos in fractional-order financial system depends on both parameters and fractional order, but there is no particular study which concentrates on how to select the proper fractional order and parameters for a concrete group of economic data. For instance, although we have several methods to control the chaos in the fractionalorder financial systems, we have no idea what the economic interpretation of these procedures of controlling chaos is. In practice, the financial and economic data are determined with inner randomness of the real-world issues. Consequently, there should exist a group of corresponding parameters and fractional order, which can be used for depicting the particular financial and economic variables. In this paper, we propose a new nonlinear dynamical model of financial and economic systems by econometric techniques. First of all, we discretize the fractional derivative defined by Jumarie's definition and then derive the estimation of parameters in the model based on the least squared principle. Finally, we will take the macroeconomic data of Japan, for an example, to study the reasonability of the econometric model empirically.
To close this introduction, we would like to make several remarks. First, the study of financial and economic models which are simplified gained much interest in recent years both from integer and fractional calculus. However, most of them deal with the mathematical properties of the models and reveal the dynamic behaviors in numerical simulations, which lead those results, although very beautiful in mathematical point of view, to be very difficult to apply in real finance or economic life. Second, we are not going to employ the existing financial models but construct a new discrete financial model for depicting the real financial and economic data from a particular country, that is, Japan. The suitable parameters and fractional orders are determined by considering the macroeconomic data of it. Finally, based on the novel discrete financial model, we can analyze fitting effectiveness of the empirical model. To the authors' knowledge, this is the first time that the realistic financial and economic data of particular country are modeled by using continuous or discrete fractional-order nonlinear dynamic financial model. The remainder of this paper is organized as follows. In Section 2, the mathematical preliminaries are introduced. In Section 3, the nonlinear dynamic econometric model of financial system with fraction-order derivative is proposed, and the parameters of the nonlinear fractional-order financial model are estimated by least squares principal. In Section 4, the empirical study of macroeconomic data of Japan is shown. Finally, the conclusions are given in Section 5.

Mathematical Preliminaries
In this section, we introduce some preliminaries of fractional derivative. More properties of fractional derivatives could be found in the books and recent papers; for example, see [14,26,27].
Let ( ) : R → R denote a continuous function, and let ℎ > 0 denote a constant discretization span. The fractional difference of order ( ∈ R, 0 < ≤ 1) of ( ) is defined as follows (see [27]): and then its fractional derivative of order is defined by In [27], the Jumarie's modified fractional derivative of order for a continuously differentiable function where < < + 1, = [ ], [ ] denotes the integer part of the real number . Furthermore, if 0 < < 1, then * There are several definitions of fractional derivative in [14], but in this paper we only use the Jumarie's fractional derivative. The advantage of applying Jumarie's fractional derivative is that we can select a small step size ℎ and remove the limit operation in (4); then, the discrete form of fractional derivative can be represented by the classical difference of function, multiplied by some coefficients.
For simplicity in writing, we remove the star " * " in Jumarie's fractional derivative. In what follows, we construct our continuous and discrete financial model by using Jumarie's fractional derivative.

Model Description and Estimation Methodology
There exist several financial models reported in recent years. For instance, the study of investment, interest rate, and price index by using a chaotic fractional Chen system is discussed in [22] as where , , and represent the interest rate, investment, and inflation respectively. The subscript indicates that the variable depends on . Parameters , , and are nonnegative coefficients with economic interpretation. ∈ (0, 1] ( = 1, 2, 3) represent the fractional order of the derivatives. If 1 = 2 = 3 = 1, (5) reduces to the integer-order Chen system.
Instead of considering the same expressions in fractional chaotic Chen system, we assume a more general form of the our financial model as where ( = 1, 2, 3) are the random errors which are assumed to the white noise generally. , , and indicate that the variables , , and depend on time .
Let ( , , , 1 ) = 1 + 11 + 12 + 13 + 14 ( , , , 2 ) = 2 + 21 + 22 + 23 + 24 ℎ ( , , , 3 ) = 3 + 31 + 32 + 33 + 34 then, the model can be rewritten as According to (4), when 0 < < 1, the model (8) can be discretized as We estimate (9) based on empirical data to determine the relationship of these variables. From the form of the model (9), it is easy to find that there do not exist common parameters in three equations of it. Therefore, the above three multivariate regression equations can be estimated separately. To state the technical procedures, we take the first equation as an example. The estimation for the parameters in the other two equations is similar. is the number of sample studied. Similar to the procedure of estimating the multivariate regression by the method of least squares, we can obtain the least squares estimator of the model aŝ0 and = ( 1 , 2 , . . . , −1 ) . The superscript indicates the transposition of matrix and vector.
In what follows, we consider the first regression equation and estimate the parameters ( 1 , 1 ). The corresponding least squares estimation is subjected to According to the evaluation result with 1 = 1, the minimum of the second product part implies that 1 =̂0.
and the minimum of SSR( 1 , 1 ) can be obtained aŝ Equations (18)- (19) are the least squares estimation of 1 and 1 in the first regression equation of system (9), respectively. It is easy to find that the estimator of 1 is not related to the sample observations, and can be computed by numerically. Using the same technique, we can deal with 2 and 3 in system (9) and obtain the optimal estimators of̂2 and̂3. In the next section, we will consider the dynamics of our new model and the prediction results based on the macroeconomic data of Japan.

Empirical Results of Discrete Fractional Financial System: Evidence from Japan
In this section, we present the study of discrete financial system based on the macroeconomic data of Japan.

Data Description.
In financial model (9), the nonlinear dynamic behaviors of interest rate, investment demand, and inflation are studied. This work chooses six-month London interbank offered rate (LIBOR) data to reflect interest rate change in Japan. The total investment percent of GDP is used to measure the investment demand. Average consumer prices percent change rate will be used to reflect the inflation. The annual data starts from year 1980 to 2011. The data about LIBOR, investment percent of GDP, and average consumer prices percent change rate are obtained from EconStats which is organized by IMF.

Empirical
Results. The optimal fractional orders , = 1, 2, 3 with different step sizes Δ = 1, 0.9, 0.8, 0.7, 0.6 are performed in Table 1. We do not consider the case of Δ < 0.5 because of the fact that the fractional order decreases and approaches to zero, which reduces our model to be a linear one but not the fractional financial system. Tables 2-4 show the results about the estimated coefficients, sum squared residuals, and values of statistical tests in the equations of empirical model under the situation of different time steps of discretization. From Table 2, we can find that the sum squared residuals are less than 5 × 10 −3 in the first equation of (9), and the values of the -statistics are less than 5% in the estimation results of different time step steps, which demonstrate that the estimation results based on the empirical data are reasonable statistically. Moreover, the time step sizes of discretization do not influence the values of -test for coefficients and value of -statistics for the empirical model. Finally, the structure of terms at which corresponding coefficients estimated are significant at 5% level in the empirical equation about interest rate includes terms , , , , 2 , constant term, and is independent of the variation of the time step sizes of discretization. Note. * * * , * * , and * denote statistical significance at the 1%, 5%, and 10% levels. SSR is the sum squared residuals. Prob is the P value of F-statistic.  Note. * * * , * * , and * : denote statistical significance at the 1%, 5%, and 10% levels. SSR is the sum squared residuals. Prob is the P value of F-statistic.
The result is consistent with the fact that the estimation of is not related to the sample observations in (19).
In Table 3, we observe that the sum squared residuals are less than 2 × 10 −3 for the empirical equation of the second equation of (9), and the values of the -statistics are 0.0068 which is less than 5%, which also demonstrate that the estimation results about the second equation of (9) are reasonable statistically. The result that the time step sizes of discretization do not influence the values can also be found in Table 3. Moreover, the structure of terms at which corresponding coefficients estimated are significant at 5% level in the empirical equation about interest rate includes terms , , , , 2 , 2 , constant term. Similar to the first empirical equation, and it is not influenced by the variation of the time step sizes of discretization.
In Table 4, we observe that the sum squared residuals are less than 1.5 × 10 −3 for the empirical equation of the third equation of (9), which suggests that the accuracy of fitness about the third empirical equation is satisfying. The values of the -statistics are 0.0000 which is less than 5%, which demonstrate that the estimation of the third equation is reasonable in the statistical sense. Moreover, the structure of terms at which corresponding coefficients estimated are significant at 5% level in the empirical equation about interest rate includes terms , , , 2 , constant term. Both of the structure of terms and values of the coefficients are not changing with the variation of the time step sizes of discretization.    describe the actual data effectively. In Figure 2, the blue line almost passes through all of the circle notations. The property of the figure declares that the actual investment data can be fitted by the empirical equation accurately. In Figure 3, there are only very few circle notations that deviate from the blue line in the beginning and ending parts of the period. Figure 3 suggests that the empirical equation about inflation fits the actual data ideally. Above all, Figures 1-3 suggest that the empirical model estimated by the methodology proposed in the paper depicts the actual data reasonably. In order to analyze the effectiveness of the consequence about the prediction of the empirical model, the multistep predictions of interest rate, investment, and inflation in the case that Δ = 1 are shown in Figures 4-6, respectively. From Figure 4, we can find that the predictions are very close to the actual interest rate from the black line which corresponds to the four-step predictions, and then the bias of predictions from the actual data increases with the steps of prediction increasing, such as purple, yellow, and green line. It suggests that the prediction about interest rate is meaningful in the first four years, and the accuracy of predictions becomes more Year Interest rate Actual data Two-step predition Three-step predition Four-step predition Five-step predition Six-step predition Seven-step predition   frustrating when the time steps of the prediction are more than four time steps. In Figure 5, the bias of the prediction about investment increases as the time steps increase, but the prediction with a few steps still approaches the actual data, such as the black line which is corresponding to the four-step prediction. In Figure 6, we also can find that the accuracy of the prediction of the inflation becomes unsatisfying in the situations that the time step of prediction is not less than four, but the predictions with steps being less than four are close to the actual inflation. The above results suggest that the four-step prediction of the empirical model estimated is meaningful for interest rate, investment, and inflation.

Conclusions
In this paper we propose a new nonlinear dynamic financial econometric model by using the Jumarie's fractional-order Inflation Actual data Two-step predition Three-step predition Four-step predition Five-step predition Six-step predition Seven-step predition derivative. The corresponding discrete financial model is generated by removing the limit operation in the Jumarie's derivative. The model overcomes the problems which cannot be depicted by some simplified nonlinear financial model in literatures, and it provides a feasible technique for describing the actual macroeconomic data of one particular place by nonlinear model. Our model is proved to be reasonable via empirical analysis. Based on the macroeconomic data of Japan, we evaluate the parameters of the financial model. The suitable fractional order for Japan's data is obtained. In the empirical study, we observe that the fractional order has an apparent influence on the dynamics behavior of financial system. With the optimal fractional order, our new fractional financial model can be used to predict the dynamic behavior of financial system of Japan in the coming years reasonably.