`Journal of Applied MathematicsVolume 2012, Article ID 735068, 12 pageshttp://dx.doi.org/10.1155/2012/735068`
Research Article

## Statistical Behavior of a Financial Model by Lattice Fractal Sierpinski Carpet Percolation

Department of Mathematics, Key Laboratory of Communication and Information System, Beijing Jiaotong University, Beijing 100044, China

Received 5 September 2011; Accepted 10 November 2011

Copyright © 2012 Xu Wang and Jun Wang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The lattice fractal Sierpinski carpet and the percolation theory are applied to develop a new random stock price for the financial market. Percolation theory is usually used to describe the behavior of connected clusters in a random graph, and Sierpinski carpet is an infinitely ramified fractal. In this paper, we consider percolation on the Sierpinski carpet lattice, and the corresponding financial price model is given and investigated. Then, we analyze the statistical behaviors of the Hong Kong Hang Seng Index and the simulative data derived from the financial model by comparison.

#### 1. Introduction

Financial fluctuation system is one of complex systems, and the statistical behavior of fluctuation of stock price changes has long been a focus of financial research. With the flourishing research of complex systems, it becomes more and more attractive to find universal rules and principles of these systems and further to answer the origination of financial complex system. Recent research is no longer restricted to the traditional areas but concentrated on the more comprehensive domains, leading to the birth of many burgeoning disciplines through the interaction and amalgamation of mathematics and other fields such as finance, biology, and sociology. For example, the theory of stochastic interacting particle systems (see [16]) recently has been applied to study the behaviors of market fluctuations, see [715]. And the study of financial market prices has been found to exhibit some universal properties similar to those observed in interacting particle systems with a large number of interacting units.

Percolation theory, as a model (in interacting particle systems) for a disordered medium, has brought new understanding and techniques to a broad range of topics in nature and society. First we consider the bond percolation on , that is, for , the distance from to is defined by , where and . By adding edges (or bonds) between all pairs , of points of with , we establish the -dimensional lattice , and we write for the set of the edges. Suppose that each bond of lattice is either open (occupied) with probability or closed (empty) with probability , then connected components of this graph are called open clusters. Let denote the open cluster containing the vertex , and be the probability that the origin belongs to an infinite open cluster. When the intensity increases from zero to one, at some sharp percolation threshold (or critical point) , for the first time, one infinite cluster appears; for all we have exactly one infinite cluster, for all we have no infinite cluster, and at critical value the incipient infinite clusters are supposed to be fractal.

A lattice fractal is a graph which corresponds to a fractal, all of them have a self-similarity, but most of them have no translation invariance, see [1, 1619]. The Sierpinski gasket and the Sierpinski carpet are well-known examples of fractals. The former is a finitely ramified fractal (i.e., it can be disconnected by removing a finite number of points) and the latter is an infinitely ramified fractal. Fractals also have close relations to financial markets [17], electrical conductivity, superconductivity, and mechanical properties of percolating systems, and so forth. In [1], it shows that the Ising model on the lattice Sierpinski carpet does exhibit the phase transition in any dimension, but the Ising model on the lattice Sierpinski gasket has no phase transition in any dimension (because of the character of the finitely ramified fractal). Similar results of phase transitions can be obtained for percolation on the lattice Sierpinski carpet and on the lattice Sierpinski gasket, see [18].

In the present paper, a new method is introduced to model and describe the fluctuations of market prices, namely, we use the lattice fractal Sierpinski carpet percolation to establish a new random market price in a financial market. In this financial model, the local interaction or influence among traders in one stock market is constructed, and a cluster of percolation is used to define the cluster of traders sharing the same opinion about the market. For the comparison, we also consider the most important index of Hong Kong financial market, the Hong Kong Hang Seng Index. We analyze the statistical properties of Hong Kong Hang Seng Index and the simulative data derived from the price model by comparison, which including the sharp peak and the fat-tail distribution for the price changes, the distribution of returns decays with power law in the tails, the price fluctuations are not invariant against time reversal (i.e., they show a forward-backward asymmetry), and so forth. Moreover, the behaviors of long memory and long-range correlation in volatility series of market returns are exhibited.

#### 2. Description of Price Model on Lattice Sierpinski Carpet Percolation

First we give a brief description of percolation on the lattice Sierpinski carpet (for ), which is defined as follows: consider as a graph in the usual sense and set

where . To make the graph more symmetric, let be the union of and its reflections in every coordinate hyperplane. Then we define the lattice Sierpinski carpet as

Similarly to Section 1, we define the corresponding edges set of as . Next we consider random graph (bond percolation) on the lattice , see Figure 1. Let (the intensity value) satisfies , each edge of is declared to be open with probability and closed with probability independently. We denote the product probability by (or ), and define , where is the open cluster containing the origin on , and is the number of vertices in . Let , then percolation on the Sierpinski carpet exhibits the existence of a phase transition, that is, for , for details see [1, 18].

Figure 1: Lattice percolation on lattice Sierpinski carpet.

Next we consider a price model of auctions for a stock in a stock market. Assume that each trader can trade the stock several times at each day , but at most one unit number of the stock at each time. Let denote the daily closing price of th trading day. And let be a subset of , where

and be a random open cluster on . Suppose that this stock consists of ( is large enough) investors, who are located in lattice. And is a random set of the selected traders who receive the information. At the beginning of trading in each day, suppose that the investors receive some news. We define a random variable for these investors, suppose that these investors taking buying positions selling positions , or neutral positions with probability or , respectively. Then these investors send bullish, bearish or neutral signal to the market. According to bond percolation on , investors can affect each other or the news can be spread, which is assumed as the main factor of price fluctuations. For a fixed , let

From the above definitions and mathematical finance theory [2024], we define the stock price at th trading day as

where is the initial stock price at time 0, and (>0) represents the depth function of the market at trading day . Then we have

The formula of the single-period stock logarithmic returns from to is given as follows:

#### 3. Experiment Analysis of Market Return Distribution

In order to make empirical research on the financial price model and an actual stock market by comparison, we select the daily closing prices of Hang Seng Index in the 20-year period from September 3, 1990 to September 3, 2010, the total number of observed data is about 4942. Recent research shows that returns on financial markets are not Gaussian, but exhibit excess kurtosis and fatter tails than the normal distribution, which is usually called the “fat-tail” phenomenon, see [21, 2530]. The general explanation for this phenomenon is thought to be the “herd effect” of investors in the market. The time series of returns by simulating the price model which is developed on the Sierpinski carpet percolation is plotted in Figure 2(a). The returns distributions of Hang Seng Index and the financial model are plotted in Figure 2(b), the part (where the probability is above the 75th or below 25th percentiles of the samples) deviates from the dash line. This implies that the probability distributions of returns deviate from the corresponding normal distributions at the tail parts.

Figure 2: (a) The returns time series of simulation data for the price model with the intensity . (b) The comparison of returns distributions for 20-year period Hang Seng index and the simulative data with .

For further analyzing the character of returns distributions for the simulative data and Hang Seng Index, we make the single-sample Kolmogorov-Smirnov test by the statistical method, the basic statistics of the corresponding returns is displayed in Table 1. The value of two-tail test is 0.000, thus the hypothesis is denied that the distribution of returns follows the Gaussian distribution.

Table 1: The Kolmogorov-Smirnov test.

In this part, we study the properties of skewness and kurtosis on the returns for the simulative data and Hang Seng Index. Next we give the definitions of skewness and kurtosis as follows: where denotes the return of th trading day, is the mean of , is the total number of trading dates, and is the corresponding standard variance. Kurtosis shows the centrality of data, and the skewness shows the symmetry of the data; it is a measure of the “peakedness” of the probability distribution of a real-valued random variable, and the infrequent extreme deviations lead higher kurtosis. Skewness is important because kurtosis is not independent of skewness, and the latter may “induce” the former. It is known that the skewness of standard normal distribution is 0 and the kurtosis is 3. Next we investigate the statistical behaviors of the returns for different intensity values , where the value changes from 0.39 to 0.55 with the interval length 0.01.

Table 2 gives a description of the statistics for 17 group data of the price model. This shows that the distribution of the returns deviates from the Gaussian distribution with the intensity values increasing, and the kurtosis distribution of the returns has a sharper peak, longer and fatter tails for larger . From the definitions in Section 2, is the intensity for the Sierpinski carpet lattice percolation and represents the strength of information spread in the price model. The wider the information spread, the larger the value of is. In the following, we hope to exhibit that the numerical characteristics of simulation results for some intensity are very close to those of the real data. We analyze the probability distributions of the logarithmic returns and the cumulative distributions of the normalized returns for these data in Figure 3, where the intensity values of the model are , , and , respectively.

Table 2: The analysis of the price model for different intensity values .
Figure 3: The plots (a), (b), and (c) are the probability distributions of the logarithmic returns, and the plot (d) is the cumulative distributions of the normalized price returns. The data is selected from Hang Seng Index and from the simulation data with the different values , , and .

#### 4. Long Memory Test of the Model and Hang Seng Index

We analyze the long memory of the returns by using Lo’s modified rescaled range statistic [31]. The long memory is measured by the Hurst exponent , calculated by Lo’s modified rescaled range statistic. For , the series exhibits the long-term persistence, with the maintenance of tendency; for , the series is the antipersistent, presenting reversion to the mean; and for , the series corresponds to a random walk. We consider a sample of series and let denote the sample mean. Then the modified rescaled range statistic, denoted by , is defined by

where

and denote the sample variance and the autocovariance estimators of .

In order to make the statistical inference for the above-modified rescaled range statistics, Lo derived that converges weakly to a random variable , where is the range of a Brownian bridge on the unit interval. Then the corresponding distribution function of is given by

Form this function , we can get test threshold for any level of significance (by examining significant of ), this reflects the long memory behavior for the time series. It is important to select the window wide ; we take the experience value

where is the estimation of first-order autocorrelation coefficient of the time series. Then the Hurst exponent is defined as the limit of the ratio . At the same time, it shows a linear growth trend between modified statistic and sample size , by using regression

With some optimal value, the statistics of returns by the modified statistic is given in Table 3 and Figure 4. Figure 4 also shows the fluctuations of exponent of returns for the price model and Hang Seng Index.

Table 3: Statistics of returns for and .
Figure 4: (a) The plots of modified statistics and the fluctuation of exponent for the price model. (b) The corresponding plots for the actual data from Hang Seng Index.

#### 5. Long-Range Correlation of the Model and Hang Seng Index

In this section, detrended fluctuation analysis (DFA) method is applied on the lattice Sierpinski carpet percolation. The DFA is a technique used to estimate a scaling exponent from the behavior of the average fluctuation of a random variable around its local trend, for the details see [26]. The cumulative deviation of time series is given by

We divide into intervals of nonoverlapping and equal length of time . Then the root mean square fluctuation for all such length interval is defined as

where is the fitting polynomial of the interval. The above definition is repeated for all the divided intervals. There is a power-law relation between and , namely,

The parameter is the scaling exponent or the correlation exponent, which exhibits the long-range correlation of the time series. For , it indicates that the time series is uncorrelated (white noise); for the value , it indicates the anticorrelations; for , the time series has the persistent long-range correlation. According to DFA method and computer simulation, the scaling exponents of the returns of the price model and Hang Seng Index are 0.51755 and 0.51531, respectively, in Figure 5. Although both the exponent values are larger than 0.5, they are very close to 0.5. This shows that there is some strong indication of long-range correlations for the returns.

Figure 5: (a) DFA analysis of the returns for the simulation data with the intensity . (b) DFA analysis of the returns for Hang Seng Index.

#### 6. Conclusion

A new random stock price model is developed by the lattice Sierpinski carpet percolation in the present paper, and a cluster of carpet percolation is applied to describe the cluster of traders sharing the same opinion about the market. The statistical properties of the returns are investigated and analyzed for different intensity values, and the behaviors of long memory and long-range correlation in volatility series are exhibited. Further, Hang Seng Index is also introduced and investigated by comparison; the empirical results show that the price model is accord with the real market to some degree.

#### Acknowledgment

The authors were supported in part by National Natural Science Foundation of China Grants no. 70771006 and no. 10971010, BJTU Foundation no. S11M00010.

#### References

1. M.-F. Chen, From Markov Chains to Non-Equilibrium Particle Systems, World Scientific Publishing, River Edge, NJ, USA, 1992.
2. R. Durrett, Lecture Notes on Particle Systems and Percolation, Wadsworth & Brooks, 1998.
3. G. Grimmett, Percolation, vol. 321, Springer, Berlin, Germany, 2nd edition, 1999.
4. T. M. Liggett, Stochastic Interacting Systems: Contact, Voter and Exclusion Processes, vol. 324, Springer, New York, NY, USA, 1999.
5. T. M. Liggett, Interacting Particle Systems, vol. 276, Springer, New York, NY, USA, 1985.
6. D. Stauffer and A. Aharony, Introduction to Percolation Theory, Taylor & Francis, London, UK, 2001.
7. P. Bak, M. Paczuski, and M. Shubik, “Price variations in a stock market with many agents,” Physica A, vol. 246, no. 3-4, pp. 430–453, 1997.
8. R. Cont and J. P. Bouchaud, “Herd behavior and aggregate fluctuations in financial markets,” Macroeconomic Dynamics, vol. 4, no. 2, pp. 170–196, 2000.
9. G. Iori, “Avalanche dynamics and trading fRICTIon effects on stock market returns,” International Journal of Modern Physics C, vol. 10, no. 6, pp. 1149–1162, 1999.
10. T. Kaizoji, S. Bornholdt, and Y. Fujiwara, “Dynamics of price and trading volume in a spin model of stock markets with heterogeneous agents,” Physica A, vol. 316, no. 1–4, pp. 441–452, 2002.
11. T. Lux and M. Marchesi, “Scaling and criticality in a stochastic multi-agent model of a financial market,” Nature, vol. 397, no. 6719, pp. 498–500, 1999.
12. D. Stauffer and D. Sornette, “Self-organized percolation model for stock market fluctuations,” Physica A, vol. 271, no. 3-4, pp. 496–506, 1999.
13. H. Tanaka, “A percolation model of stock price fluctuations,” Mathematical Economics, no. 1264, pp. 203–218, 2002.
14. J. Wang and S. Deng, “Fluctuations of interface statistical physics models applied to a stock market model,” Nonlinear Analysis: Real World Applications, vol. 9, no. 2, pp. 718–723, 2008.
15. J. Wang, Q. Wang, and J. Shao, “Fluctuations of stock price model by statistical physics systems,” Mathematical and Computer Modelling, vol. 51, no. 5-6, pp. 431–440, 2010.
16. B. B. Mandelbrot, The Fractal Geometry of Nature, W. H. Freeman, San Francisco, Calif, USA, 1982.
17. E. E. Peters, Fractal Market Analysis: Applying Chaos Theory to Investment and Economics, John Wiley & Sons, New York, NY, USA, 1994.
18. M. Shinoda, “Existence of phase transition of percolation on Sierpiński carpet lattices,” Journal of Applied Probability, vol. 39, no. 1, pp. 1–10, 2002.
19. J. Wang, “Supercritical ising model on the lattice fractal—the Sierpinski carpet,” Modern Physics Letters B, vol. 20, no. 8, pp. 409–414, 2006.
20. R. Gaylord and P. Wellin, Computer Simulations with Mathematica: Explorations in the Physical, Biological and Social Science, Springer, New York, NY, USA, 1995.
21. K. Ilinski, Physics of Finance: Gauge Modeling in Non-Equilibrium Pricing, John Wiley, New York, NY, USA, 2001.
22. D. Lamberton and B. Lapeyre, Introduction to Stochastic Calculus Applied to Finance, Chapman & Hall, London, UK, 2000.
23. S. M. Ross, An Introduction to Mathematical Finance, Cambridge University Press, Cambridge, UK, 1999.
24. J. Wang, Stochastic Process and Its Application in Finance, Tsinghua University Press and Beijing Jiaotong University Press, Beijing, China, 2007.
25. F. Black and M. Scholes, “The pricing of options and corporate liabilities,” Journal of Political Economy, vol. 81, pp. 637–654, 1973.
26. P. Grau-Carles, “Long-range power-law correlations in stock returns,” Physica A, vol. 299, no. 3-4, pp. 521–527, 2001.
27. Y. Guo and J. Wang, “Simulation and statistical analysis of market return fluctuation by Zipf method,” Mathematical Problems in Engineering, vol. 2011, Article ID 253523, 13 pages, 2011.
28. M. L'evy and S. Solomon, Microscopic Simulation of Financial Markets, Academic Press, New York, NY, USA, 2000.
29. D. Pirino, “Jump detection and long range dependence,” Physica A, vol. 388, no. 7, pp. 1150–1156, 2009.
30. Z. Zheng, Matlab Programming and the Applications, China Railway Publishing House, Beijing, China, 2003.
31. E. Zivot and J. H. Wang, Modeling Financial Time Series with S-PLUS, Springer, New York, NY, USA, 2006.