This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Analysis of intra-day fluctuations in the Mexican financial market index

-Léster Alfonsoa, Danahe E. Garcia-Ramirezb, Ricardo Mansillac, César A. Terrero-Escalanted a Universidad Autónoma de la Ciudad de México,
C.P. 09790, México D.F., México
b Departamento de Astronomía,
DCNE-CGT, Universidad de Guanajuato,
C.P. 36023, Guanajuato, México.
c Centro de Investigaciones Interdisciplinarias en Ciencias y Humanidades,
Universidad Nacional Autónoma de México,
Ciudad Universitaria, C.P. 04510, México D.F., México.
dFacultad de Ciencias, Universidad de Colima,
Bernal Díaz del Castillo 340, Col. Villas San Sebastián,
C.P. 28045, Colima, Colima, México.
(August 18, 2025)
Abstract

In this paper, a statistical analysis of high frequency fluctuations of the IPC, the Mexican Stock Market Index, is presented. A sample of tick–to–tick data covering the period from January 1999 to December 2002 was analyzed, as well as several other sets obtained using temporal aggregation. Our results indicates that the highest frequency is not useful to understand the Mexican market because almost two thirds of the information corresponds to inactivity. For the frequency where fluctuations start to be relevant, the IPC data does not follows any α\alpha-stable distribution, including the Gaussian, perhaps because of the presence of autocorrelations. For a long range of lower-frequencies, but still in the intra-day regime, fluctuations can be described as a truncated Lévy flight, while for frequencies above two-days, a Gaussian distribution yields the best fit. Thought these results are consistent with other previously reported for several markets, there are significant differences in the details of the corresponding descriptions.

keywords:
stock markets , high frequency fluctuations distribution , tail behavior , autocorrelations
MSC:
91B82 , 60E07 , 62D05
journal: Physica A

1 Introduction

Is a fundamental assumption of the economic theory of markets that the financial markets are efficient in the sense of the community of economic agents being able to discount all the possibilities of arbitrage, incorporating with it all the relevant information for the prices formation, so that it is impossible to design an (always) winning strategy for investment[1, 2]. 111For more details on the background and definitions of the concepts we use in this paper see book [2]. Mathematically, this is described as a process, where a future increase or decrease of the current price is the result of a random event. Underlying such a random walk is a binomial distribution, which after a very large number of steps (price changes), converges to a normal, Gaussian, distribution. In turn, normal distributions are common in the description of systems in equilibrium, a stable state, small fluctuations around which decay exponentially with time. However, it seems that the financial market cannot actually be modeled as this kind of system; data for most financial indexes around the world are statistically described by probability distributions that exhibit a large skewness or kurtosis, relative to that of a normal distribution. It may means that strong perturbations are not, de facto, necessary for the market entering a critical state; recurrent significant deviations of economic variables from their average values could be just the cumulative and unavoidable result of many small-scale processes and interactions occurring within the market system, continuously subjected to a net external action like, for instance, publicly available announcements of annual earnings, stock splits, companies profits forecasts, new securities or, even more, the punctual action of agents having preferential access to restricted or confidential information.

In this regard, for years now particularly good fits to financial data have been obtained using Lévy-stable distributions [3]. Nevertheless, it is obvious that, strictly speaking, the market can neither be such a random process. First, it is unreasonable to expect data from any real (as opposed to hypothetical) process to have an infinite variance. Secondly, it is reasonable to expect data from financial transactions at a given time to have some memory of the previous transactions, i.e., the elements of the corresponding time series should not be really independent each from the others. This is easily noted, for instance, during steeped variations of the stock returns due to panic or euphoria.

Commonly there are two causes of autocorrelation in this kind of data, irregular sampling and causality determined by market microstructure. In agreement with the ‘efficient market hypothesis’, it is intuitive not to expect persistent serial dependence in price changes, otherwise it would be used to influence the market. Consistently with this, autocorrelation is often neutralized after homogenization of the corresponding series by time averaging of subsamples. This typically results in sets with information correlated with those of daily or lower frequency data. It is for these last sets that fitting to stable distributions have already yielded remarkably good results [4, 5, 6, 7, 8]. More precisely, it has been verified that the actual density distribution behind the data could be a so-called Lévy-truncated [9, 10] which, after sequential convolution corresponding to the sum of the values over increasing time intervals, also converges to a Gaussian distribution. This convergence is ultra-slow, therefore the truncation still allows for a relatively high probability of extreme fluctuations which, as it was already noted, seems to be a distinctive feature of most financial markets.

Despite the success in describing markets activity as a Lévy (truncated) flight, it is widely acknowledged that the description will not be complete unless high-frequency data is also brought into the study. On one hand, most market models are based on a variety of hypothesis regarding the long-memory features of volatility which are difficult to extract from daily or lower-frequency data, but can be observed in intraday data. For instance, as it was already mentioned, crisis in finances may not necessarily be linked only to strong perturbations over several days, but also to the cumulative effect of intraday weaker perturbations. On the other hand, the quality of risk analysis of investment depends on the accuracy of measurement of ex post volatility and of forecast evaluation. There is evidence of the improvement in both directions thanks to the availability of high-frequency data from liquid financial markets such as the foreign-exchange, bond or equity-index markets afford [11]. Last, but not least, the analysis of the impact of sample size on the probability estimation of extreme events yields that a large volume of data is needed in order to determine the actual stable distribution corresponding to a given asymptotic scaling[12, 8]. This is the kind of volume that commonly characterizes the sets of high frequency data.

It must be pointed out that market microstructure and short term interactions become relevant while analyzing intra-day fluctuations and, since they depend on local socio-economic factors, the stylized facts of the market as a complex system can be difficult to extract. Therefore, it is very important to study different realizations of high frequency financial data in order to discriminate local and universal properties of the market dynamics. It is with this intention that in this paper we report a first step into analyzing the tick-to-tick data of the Indice de Precios y Cotizaciones (IPC), the main benchmark stock index in Mexico. It has been previously shown in Ref.[8] that it can be dismissed that a normal distribution relies under the data of daily closure values of the IPC. On the other hand, the null hypothesis that it comes from α\alpha-stable Lévy distribution cannot be rejected at the 5% significance level. This implies that the daily data can safely being considered as independent and identically distributed, characterized by an infinite variance. Our aim here is to study how this conclusion changes when tick-to-tick data is analyzed.

The paper is organized as follows. In the next section we provide a brief review of Lévy distributions, in particular of the stable and Lévy truncated distributions which we use for fitting the IPC fluctuations. The description of this data, as well as the sets derived from them and used for our analysis is presented in section 3. Next in sections 4 and 5 we proceed to present the analysis of the probability distribution and serial dependence of the actual fluctuations of the IPC data. We devote the last section to discuss our results and state the main conclusions.

2 Lévy distributions

As it was mentioned in the introduction, Lévy distributions are among the more frequently used to fit data from complex processes. Particularly, it has been found to be an excellent fit to the distribution of stock returns as well as of other financial time series. In this section we will briefly review the definition and properties of the stable and truncated Lévy distributions.

2.1 Stable distributions

While studying the behavior of sums of independent random variables Paul Lévy [3] introduced a skew distribution specified by scale γ\gamma, exponent α\alpha, skewness parameter β\beta and a location parameter μ\mu. Since the analytical form of the Levy stable distribution is known only for a few cases, they are generally specified by their characteristic function. The most popular parameterization is defined by Samorodnitsky and Taqqu [13] with the characteristic function:

ϕ(t)={exp(γ|t|[1+iβ2πsign(t)ln(|t|)+iμt]),if α=1.exp(γα|t|α[1iβtan(πα2)sign(t)+iμt]),otherwise.\phi(t)=\begin{cases}\exp{\left(-\gamma|t|\left[1+i\beta\frac{2}{\pi}\rm{sign}(t)\ln(|t|)+i\mu t\right]\right)},&\text{if $\alpha=1$}.\\ \exp{\left(-\gamma^{\alpha}|t|^{\alpha}\left[1-i\beta\tan\left(\frac{\pi\alpha}{2}\right){\rm sign}(t)+i\mu t\right]\right)},&\text{otherwise}.\end{cases} (1)

where sign(t)sign(t) stands for the sign of tt. Then, the probability density function is calculated from it with the inverse Fourier transform in the form:

f(x;α,β,γ,μ)=12π+ϕ(t)eitx𝑑t.f(x;\alpha,\beta,\gamma,\mu)=\frac{1}{2\pi}\int^{+\infty}_{-\infty}\phi(t){\rm e}^{-itx}dt\,. (2)

Lévy distributions are characterized by the property of being stable under convolution, i.e, the sum of two independent and identically Lévy-distributed random variables, is also Lévy distributed with the same stability index α\alpha. The stability parameter α\alpha lies in the interval (0,2](0,2]. Small α\alpha represents a sharp peak but heavy tails which asymptotically decay as power laws with exponent (α+1)-(\alpha+1). For the normal distribution α=2\alpha=2. For symmetric distribution (like the normal distribution), the skewness parameter β=0\beta=0. The skewness parameter must lie in the range [1,1][-1,1]. When β=+1,1\beta=+1,-1, one tail vanishes completely. The parameter γ\gamma lies in the interval (0,)(0,\infty), while the location parameter μ\mu is in (,+)(-\infty,+\infty).

The asymptotic behavior of the Lévy distributions is described by the expression

f(x;α)|x|1α.f(x;\alpha)\approx|x|^{-1-\alpha}\,. (3)

Hence, the variance of the Levy stable distributions is infinite for all α<2\alpha<2.

2.2 Truncated Lévy distributions

As mentioned in the previous subsection, α\alpha-stable Levy distributions have infinite variance, hence, they have power-law tails that decay too slowly. Therefore, the fit of empirical data by Lévy stable distributions will usually overestimate the probability of extreme events. In particular, real prices fluctuations have finite variance, so their distribution decays slower than a Gaussian, but faster than a Levy-stable distribution, with the tails better described by an exponential law than by a power law. The truncated Lévy flight (TLF) was proposed by Mantegna and Stanley [9] to overcome this problem, and can be defined as a stochastic process with finite variance and scaling relations in a large, but finite interval. They defined the truncated Lévy flight distribution as:

P(x)={0if x>lcPL(x)if lxl0if x<l,P(x)=\left\{\begin{array}[]{lll}0&\mbox{if }x>l\\ cP_{L}(x)&\mbox{if }-l\leq x\leq l\\ 0&\mbox{if }x<-l\,,\end{array}\right. (4)

where PL(x)P_{L}(x) is a symmetric Lévy distribution. As the TLF has a finite variance, with sequential convolution it will converge to a Gaussian process, but the convergence is very slowly, as was demonstrated by Mantegna and Stanley [9]. However, the cutoff in the tail given by (4) is abrupt. This problem was solved by Koponen [10], who introduced an infinitely divisible TLF with an exponential cutoff with the characteristic function:

logϕ(t)=cαcos(πα2)[(t2+λ2)α2cos(αarctan|t|λλα)],\log\phi(t)=-\frac{c^{\alpha}}{\cos(\frac{\pi\alpha}{2})}\left[\left(t^{2}+\lambda^{2}\right)^{\frac{\alpha}{2}}\cos\left(\alpha\arctan\frac{|t|}{\lambda}-\lambda^{\alpha}\right)\right]\,, (5)

where cc is the scaling factor, α\alpha is the stability index, and λ\lambda is the cutoff parameter. The Lévy α\alpha-stable law is restored by setting to zero the cutoff parameter. For small values of xx, the truncated Levy density described by the characteristic function (5), behaves like a Lévy-stable law of index α\alpha [14]. It was used by Matacz [15] to describe the behavior of the Australian All Ordinaries Index.

3 Data sets from IPC values

For our analysis we used the IPC value over the period January 1999-December 2002 which comprises 43214274321427 transactions. Taking into account that the Mexican trading day is of six and a half hours, this give us a mean time between transactions of 5.25.2 seconds. With the aim of analyzing the statistics of nontrivial index fluctuations, the repeated values (the intervals where no change in the returns was recorded) were removed and the set was reduced to N=1164256N=1164256 elements with a mean value over the period of 6209.3926209.392, a variance of 744334.3744334.3 and an excess of kurtosis of 1537.4421537.442. These ticks, YkY_{k}, are now irregularly sampled, with a mean time between fluctuations of 19.3 seconds. This set is plotted in figure 1.

Refer to caption
Figure 1: IPC actual fluctuations. The horizontal red line corresponds to the mean value over the period.

The actual fluctuations are then taken as the difference of the natural logarithm of the YkY_{k},

Sk=lnYk+1lnYk,fork=1,2,,N.S_{k}=\ln Y_{k+1}-\ln Y_{k}\,,\quad{\rm for}\,k=1,2,\dots,N\,. (6)

For this study we also used sets corresponding to the convolution of the density distributions of high-frequency data {Sk}\{S_{k}\}, i.e., sets obtained after summing for different values of NconvN_{conv},

SjNconv=k=1+(j1)×Nconvj×NconvSk,forj=1,2,,N¯,S_{j}^{N_{conv}}=\sum_{k=1+(j-1)\times N_{conv}}^{j\times N_{conv}}S_{k}\,,\quad{\rm for}\,j=1,2,\dots,\bar{N}\,, (7)

where N¯\bar{N} is the multiple of NconvN_{conv} closer to NN.

For completeness we analyzed also a set of closing values of the IPC for the same period downloaded from the Yahoo Finance website.

4 Probability distribution of IPC fluctuations

We started by analyzing the non-convoluted data. The corresponding distribution is presented in figure 2, where a logarithmic scale is used for the vertical axis and the horizontal axis has been rescaled dividing by the standard deviation. The plots of best fits to a Gaussian (narrow blue curve below the data points) and a Lévy (wide red curve above the data points) distributions are also shown.222For the analysis of stable distributions, we used the library developed by J. Nolan [16].

Refer to caption
Figure 2: Distributions for Nconv=0N_{conv}=0. In the horizontal axis ZZ stands for the fluctuations data or the values of the corresponding fits.

It can be observed that the probability distribution function for this data does not correspond to a normal distribution, but it neither does to a Lévy one. This is confirmed by the results of the Kolmogorov-Smirnov test shown in table 4 for the best fit of data to an α\alpha-stable distribution.

Table 1: Results of the Kolmogorov-Smirnov Goodness of Fit Test (K-S test).
NconvN_{conv} α\alpha β\beta γ\gamma δ\delta K-S Statistics p-value
0 1.25651.2565 0.0024-0.0024 0.37960.3796 0.00140.0014 0.01790.0179 0.00.0 Yes
10 1.57881.5788 0.00560.0056 2.31382.3138 0.0051-0.0051 0.00970.0097 0.00.0 Yes
20 1.61071.6107 0.00890.0089 3.84973.8497 0.0114-0.0114 0.00760.0076 0.00260.0026 Yes
30 1.62961.6296 0.00720.0072 5.24205.2420 0.0075-0.0075 0.00810.0081 0.01150.0115 Yes
40 1.65391.6539 0.01400.0140 6.63266.6326 0.03553-0.03553 0.00990.0099 0.00640.0064 Yes
50 1.65231.6523 0.23320.2332 7.82167.8216 0.0497-0.0497 0.01080.0108 0.00840.0084 Yes
60 1.65941.6594 0.01450.0145 8.93918.9391 0.0492-0.0492 0.01000.0100 0.04040.0404 Yes
70 1.66471.6647 0.0053-0.0053 10.079310.0793 0.03590.0359 0.00780.0078 0.26020.2602 No
80 1.66821.6682 0.00100.0010 11.189311.1893 0.03390.0339 0.00890.0089 0.20220.2022 No
90 1.66881.6688 0.01410.0141 12.170512.1705 0.01030.0103 0.01020.0102 0.13650.1365 No
100 1.68451.6845 0.00290.0029 13.294213.2942 0.03290.0329 0.00990.0099 0.20530.2053 No
110 1.69161.6916 0.00800.0080 14.263514.2635 0.0547-0.0547 0.01230.0123 0.08190.0819 No
120 1.68871.6887 0.00490.0049 15.209615.2096 0.0066-0.0066 0.01140.0114 0.16070.1607 No
130 1.68261.6826 0.0291-0.0291 16.027016.0270 0.23030.2303 0.01310.0131 0.09030.0903 No
140 1.68741.6874 0.0211-0.0211 16.909416.9094 0.26230.2623 0.01150.0115 0.21630.2163 No
150 1.69091.6909 0.0027-0.0027 17.876217.8762 0.11310.1131 0.01070.0107 0.33440.3344 No
1200 1.86931.8693 0.3670-0.3670 71.896471.8964 6.88226.8822 0.01500.0150 0.53280.5328 No
2500 1.90651.9065 0.6018-0.6018 108.4828108.4828 12.090312.0903 0.03010.0301 0.78600.7860 No
2700 2.00002.0000 0.79730.7973 122.6070122.6070 5.01045.0104 0.03090.0309 0.79850.7985 No

From this table we can also observe that as NconvN_{conv} reaches a value around 7070, the test cannot longer reject the hypothesis of the probability density function for this data being an α\alpha-stable distribution. As it is shown in fig.3,

Refer to caption
Figure 3: Results of the Kolmogorov-Smirnov Goodness of Fit Test versus NconvN_{conv}.

Lévy scaling holds over a long range of values of NconvN_{conv}. For instance, in fig.4 is presented the data and the best fits for a normal and an α\alpha-stable distribution for Nconv=100N_{conv}=100.

Refer to caption
Figure 4: Distributions for Nconv=100N_{conv}=100. In the horizontal axis ZZ stands for the fluctuations data or the values of the corresponding fits.

The values of α\alpha keep steadily increasing as NconvN_{conv} is also increased. As a matter of fact, as it can be seen in table 4 and it is represented in Fig.5,

Refer to caption
Figure 5: The stable coefficient α\alpha as function of NconvN_{conv}.

the value of the stable coefficient slowly converges to 22, while the convolution involves larger blocks of data. That is, for instance, the value of α\alpha for Nconv=2700N_{conv}=2700. For this case the corresponding statistics are presented in the last row of table 4 and the distributions are plotted in Fig.6.

Refer to caption
Figure 6: Distributions for Nconv=2700N_{conv}=2700. In the horizontal axis ZZ stands for the fluctuations data or the values of the corresponding fits.

As we can see, for Nconv=2700N_{conv}=2700 a normal distribution and the corresponding Lévy distribution (with α=2\alpha=2) give both a very good fit to the convoluted data. This is a strong evidence that, indeed, the variance of the data is finite.

4.1 Convergence to α=2\alpha=2

Note in figure 5 that α2\alpha\approx 2 for NconvN_{conv} around 20002000 too. Also in table 4 it can be seen that the convergence to α=2\alpha=2 is not just slow, as has been noted previously [9], but it is also non-uniform. We believe that this is an effect of the finite number of elements in the sample. To verify that, we simulated data by truncating sets generated using the stable library by Nolan [16] and following expression (4). For the three sets we used the same parameters, but the length of the series are one, two and five millions elements respectively. In fig.7 is plotted how the stable coefficient α\alpha for each original set (without truncation) evolves with convolution.

Refer to caption
Figure 7: The stable coefficient α\alpha as function of NconvN_{conv}, without truncation.

It is noticeable the high quality of the simulated data, since in each case α\alpha converges to a value significantly different from 22 and, the larger the size of the set, the closer this value gets to the α\alpha used for generating the set.

The truncation was done by erasing out the elements of a given set with absolute value greater than nstd×σn_{std}\times\sigma, where σ\sigma is the standard deviation of the corresponding data. It implies that, the smaller nstdn_{std}, the fewer the elements remaining in the truncated set. There is a range of nstdn_{std} when the simulation works. For high nstdn_{std} (i.e., small truncation), the length of the series is not large enough for observing the convergence to α=2\alpha=2, thought it seems to converge to a value significantly different from the value used for generating the set, i.e., the case without convolution presented in fig.7. An example is given in fig.8.

Refer to caption
Figure 8: The stable coefficient α\alpha as function of NconvN_{conv} for nstd=50n_{std}=50.

On the other hand, for low nstdn_{std} (i.e., large truncation), the Kolmogorov-Smirnov test rejects that the data correspond to a stable distribution. Neither it is a normal distribution, but it converges very fast to α=2\alpha=2. An example is given now in fig.9.

5 Refer to caption

Figure 9: The stable coefficient α\alpha as function of NconvN_{conv} for nstd=10n_{std}=10.

Therefore, reasonable truncation can be performed, such that it can be obtained a series that for low NconvN_{conv} the test cannot reject them to be stable-distributed, but with convolution converges to α=2\alpha=2. The corresponding example is given in fig.10.

Refer to caption
Figure 10: The stable coefficient α\alpha as function of NconvN_{conv} for nstd=35n_{std}=35.

From these simulations it can be observed that the convergence is not only slow, but also non uniform. However, the larger the amount of elements in a given set, the smaller the size of the irregularities. If considering the whole population, (a condition for the generalization of the Central Limit Theorem), the convergence can be expected to be uniform.

4.2 Lévy to Gaussian crossover

Since our sample is finite and this affect the estimation of the parameter α\alpha, we further analyzed the transition from Lévy to Gaussian regime by following the procedure proposed in Ref.[17]. We study the behavior of the excess kurtosis for convoluted samples, from Nconv=0N_{conv}=0 to Nconv=3000N_{conv}=3000. The excess kurtosis

k(Siμ)4(Siμ)223,k\equiv\frac{\langle(S_{i}-\mu)^{4}\rangle}{\langle(S_{i}-\mu)^{2}\rangle^{2}}-3\,, (8)

gives a statistical measure of the heaviness of the tail of a distribution with mean value μ\mu. A normal process shows zero excess kurtosis for the population, while it is positive for leptokurtic distributions like Lévy-stable distributions.

The results we obtained are presented in fig. 11.

Refer to caption
Figure 11: Excess kurtosis for the samples as a function of NconvN_{conv}. The circle encloses the Lévy to Gaussian crossover.

It can be seen that the transition between Lévy and Gaussian regimes occurs for values of NconvN_{conv} between 2700 and 3000. This confirms the result discussed above for the fit of the convoluted data to a stable Lévy distribution (see table 4), that gives a value of α=2\alpha=2 for Nconv=2700N_{conv}=2700. This way, the crossover time can be set equal to Nc=2700N_{c}=2700. Recalling that, during the analyzed period (from January 1999 to December 2002), the average time between successive fluctuations is close to 20 seconds and that for the Mexican market one trading day is equal to 6.5 hours, we find that the Levy-Gaussian crossover is approximately 2.3 trading days.

5 Serial dependence in the IPC data

From the previous section we concluded that for convolutions below Nconv=70N_{conv}=70 data does not fit neither a Gaussian nor any Lévy distribution. A common hypothesis for both cases is the data being independent and identically distributed. In figures 12 and 13 we present the results for the analysis of autocorrelation for the set {Sk}\{S_{k}\}.

Refer to caption
Refer to caption
Figure 12: Autocorrelations for Nconv=0N_{conv}=0 (left) and Nconv=10N_{conv}=10 (right)
Refer to caption
Refer to caption
Figure 13: Autocorrelations for Nconv=50N_{conv}=50 (left) and Nconv=150N_{conv}=150 (right)

In principle, the vertical axis would cover values from 1-1 (full anti-correlation) to 11 (full correlation) and the values in the horizontal axis stand for the lag δ\delta, i.e., denote the correlation between SkS_{k} and Sk+δS_{k+\delta}. The blue dashed lines in these figures indicate approximate limits of correlation coefficients expected under a null hypothesis of uncorrelated data. We successfully tested these limits using the simulations of truncated Lévy distributions described in subsection 4.1. As it can be observed the fluctuations exhibit positive autocorrelations, which are correspondingly diluted after convoluting the series. It suggests a mild serial dependence between fluctuations within an interval of about 4848 minutes.

In figures 14

Refer to caption
Refer to caption
Figure 14: Autocorrelations for the absolute value of the log returns for Nconv=0N_{conv}=0 (left) and Nconv=150N_{conv}=150 (right)

are shown the autocorrelations for the series of absolute values of SkS_{k} for Nconv=0N_{conv}=0 (left) and Nconv=150N_{conv}=150 (right). This is a measure of volatility and it exhibits long range serial dependence, a fact consistent with findings reported for other markets [2, 7].

Finally, we present in figure 15

Refer to caption
Refer to caption
Figure 15: Autocorrelations for the log returns of the dayly closing (left) and the corresponding absolute value (right).

the autocorrelations for the data and absolute value of the data of the daily closing set. As it can be observed these last results are totally consistent with those obtained for the tick–to–tick set.

6 Discussion

Our results suggest that the statistical description of the Mexican market strongly depends on the time scale of interest. For the analyzed period, around 73% of the IPC tick–to–tick data, sampled every 5.25.2 seconds, shows no activity. Fluctuations occur, in average, every 19.319.3 seconds, but the underlying probability is unlikely to follow a Lévy-stable distribution, including the Gaussian. Our analysis shows that one of the reasons could be that these data are not independent. The fluctuations exhibit mild positive autocorrelations that persist for about 4848 minutes before falling below the level of noise. This is, for instance, twice the value reported in Ref.[2] for the S&PS\&P 500500 index, sampled at a 1 min time scale, and also sometimes reported for various asset returns [7]. The autocorrelation in the IPC are diluted by convolution, i.e., the aggregation of elements of the set of fluctuations in blocks of length NconvN_{conv}. This is characteristic of random walks with short memory. Conversely, in the case of a deterministic process with noise, even if autocorrelations are initially hidden, they surface and get more noticeable with convolutions. Therefore, this serial dependence seems to reflect less the internal mechanics of the market (due to the law of supply and demand) than the (complex, noisy) external action on it.

It is worthy to note that sets obtained by convolution of the IPC tick–to–tick data can be safely used to describe the activity at lower frequencies. We tested this by comparing, for instance, the fit for Nconv=1200N_{conv}=1200, which corresponds to one trading day, with the daily data from Yahoo finance for the same period (from January 1999 to December 2002). From table 4 we see that α=1.8693\alpha=1.8693, while for the closure data we obtained 1.86461.8646. Moreover, the analysis of serial dependence for the closure data (and their absolute value) presented in figures 15 is totally consistent with those obtained for the intra-day set described in the paragraph above.

Taking all of this into account, the IPC fluctuations can be described by Lévy-stable distributions in the wide range of sampling interval from 20 seconds up to two trading days. After that time, data seems to obey normal probability distributions. We would like to note that this value for the Lévy to Gaussian crossover also differs from results obtained by other authors. For example, Mantegna and Stanley [5], for the S&P500S\&P500 (during the six year period from January 1984 to December 1989), estimated the crossover time to be of the order of one month; Matacz [15], for the Australian All Ordinaries share market index for the period 1993–1997, found that the crossover time is approximately 19 trading days, and in Cuoto Miranda and Riera [17], a crossover of approximately 20 days was found for the Sao Paulo Stock Exchange Index in Brazil (IBOVESPA), during the 15 years period 1986–2000.

The differences outlined here with respect to the stylized facts found for other markets seems to indicate that in the given period the behavior of the Mexican market was atypical: the effects induced by external factors (for instance, different kinds of expected announcements and unexpected action of privileged agents) led the dynamics for periods of about an hour. As this effect vanishes, significant (as compared to normal) fluctuations of the index from its average value were likely within an interval of a couple of trading days, but if analyzed over longer intervals, data describes a usual random walk, i.e., with the size of the fluctuations decaying exponentially.

Nevertheless, there are reasons for expecting the statistical description of the IPC to also depend on the period analyzed, and the variation from period to period not just being given by the presence of critical events. As mentioned before, in Ref.[8] we analyzed the daily closure data for the IPC covering the period from 04/09/200004/09/2000 to 04/09/201004/09/2010. For this whole period we obtained α=1.64\alpha=1.64, which is significantly lower than the value of 1.86461.8646 we obtained for the daily data over the period 01/199901/1999 to 12/200212/2002 studied in this paper. Both these values are still larger than, for example, the one reported by Mantegna and Stanley [5] for the S&P500S\&P500 about a decade earlier.

7 Acknowledgments

This research was supported by the Sistema Nacional de Investigadores (México). It was also partially funded by a grant from the Consejo Nacional de Ciencia y Tecnología of Mexico (SEP–CONACyT) CB-284482. D. G-R. thanks the support of the Facultad de Ciencias, Universidad de Colima, where most of her contribution was done.

References

  • [1] E. F. Fama, J. Finance, 46 1575–613 (1991)
  • [2] R. N. Mantegna and H. E. Stanley, An introduction to econophysics: correlations and complexity in finance., (Cambridge University Press, New York, (1999))
  • [3] P. Lévy, Calcul des Probabilites, (Paris : Gauthier-Villars et Cie, (1925))
  • [4] B. B. Mandelbrot, Journal of Business 36, 394 (1963).
  • [5] R. N. Mantegna and H. E. Stanley, Nature 376 46 (1995).
  • [6] P. Gopikrishnan, V. Plerou, L. A. N. Amaral, M. Meyer, and H. E. Stanley, Phys. Rev. E 60, 5305–5316 (1999)
  • [7] Cont, R., Quantitative Finance, 1: 223-236, (2001)
  • [8] L. Alfonso, R. Mansilla, and C.A. Terrero-Escalante, Physica A: Statistical Mechanics and its Applications 391, 2990–2996 (2012).
  • [9] R. N. Mantegna and H. E. Stanley, Physical Review Letters, 73 22, 2946 (1994).
  • [10] I. Koponen, Phys. Rev. E 52, 1197 (1995)
  • [11] T. G. Andersen, Journal of Business and Economic Statistics 18, no. 2 (2000): 146-53. doi:10.2307/1392552.
  • [12] R. Weron, International Journal of Modern Physics C 12 2 (2001).
  • [13] G. Samorodnitsky and M. S. Taqqu, Stable Non-Gaussian Random Processes. Stochastic Models with Infinite Variance., (Chapman and Hall, New York - London (1994))
  • [14] R. Cont, M. Potters, J. P. Bouchaud, Scale Invariance and Beyond. (Centre de Physique des Houches, vol 7. Springer, Berlin, Heidelberg (1997))
  • [15] A. Matacz, International Journal of Theoretical and Applied Finance, 3 (01), 143-160 (2000).
  • [16] http://fs2.american.edu/jpnolan/www/stable/stable.html
  • [17] L. C. Miranda, and R. Riera,  Physica A: Statistical Mechanics and its Applications, 297 (3-4), 509-520 (2001).