This module implements the fractal approach to understanding financial markets that was pioneered by Mandelbrot. In particular, it implements the multifractal random walk model of asset returns as developed by Bacry, Kozhemyak, and Muzy, 2006, Continuous cascade models for asset returns and many other papers by Bacry et al. See http://www.cmap.polytechnique.fr/~bacry/ftpPapers.html
See also Mandelbrot’s The Misbehavior of Markets for a motivated introduction to the general idea of using a self-similar approach to modeling asset returns.
One of the main goals of this implementation is that everything is highly optimized and ready for real world high performance simulation work.
AUTHOR:
Returns the partial sums of a fractional Gaussian noise simulation with the same input parameters.
INPUT:
OUTPUT:
List of n time series.
EXAMPLES:
sage: set_random_seed(0)
sage: finance.fractional_brownian_motion_simulation(0.8,0.1,8,1)
[[-0.0754, 0.1874, 0.2735, 0.5059, 0.6824, 0.6267, 0.6465, 0.6289]]
sage: set_random_seed(0)
sage: finance.fractional_brownian_motion_simulation(0.8,0.01,8,1)
[[-0.0239, 0.0593, 0.0865, 0.1600, 0.2158, 0.1982, 0.2044, 0.1989]]
sage: finance.fractional_brownian_motion_simulation(0.8,0.01,8,2)
[[-0.0167, 0.0342, 0.0261, 0.0856, 0.1735, 0.2541, 0.1409, 0.1692],
[0.0244, -0.0153, 0.0125, -0.0363, 0.0764, 0.1009, 0.1598, 0.2133]]
Return n simulations with N steps each of fractional Gaussian noise with Hurst parameter H and innovations variance sigma2.
INPUT:
OUTPUT:
List of n time series.
EXAMPLES:
We simulate a fractional Gaussian noise:
sage: set_random_seed(0)
sage: finance.fractional_gaussian_noise_simulation(0.8,1,10,2)
[[-0.1157, 0.7025, 0.4949, 0.3324, 0.7110, 0.7248, -0.4048, 0.3103, -0.3465, 0.2964],
[-0.5981, -0.6932, 0.5947, -0.9995, -0.7726, -0.9070, -1.3538, -1.2221, -0.0290, 1.0077]]
The sums define a fractional Brownian motion process:
sage: set_random_seed(0)
sage: finance.fractional_gaussian_noise_simulation(0.8,1,10,1)[0].sums()
[-0.1157, 0.5868, 1.0818, 1.4142, 2.1252, 2.8500, 2.4452, 2.7555, 2.4090, 2.7054]
ALGORITHM:
See Simulating a Class of Stationary Gaussian Processes using the Davies-Harte Algorithm, with Application to Long Meoryy Processes, 2000, Peter F. Craigmile for a discussion and references for why the algorithm we give – which uses the stationary_gaussian_simulation() function.
Return a list of n simulations of a multifractal random walk using the log-normal cascade model of Bacry-Kozhemyak-Muzy 2008. This walk can be interpreted as the sequence of logarithms of a price series.
INPUT:
OUTPUT:
List of time series.
EXAMPLES:
sage: set_random_seed(0)
sage: a = finance.multifractal_cascade_random_walk_simulation(3770,0.02,0.01,0.01,10,3)
sage: a
[[-0.0096, 0.0025, 0.0066, 0.0016, 0.0078, 0.0051, 0.0047, -0.0013, 0.0003, -0.0043],
[0.0003, 0.0035, 0.0257, 0.0358, 0.0377, 0.0563, 0.0661, 0.0746, 0.0749, 0.0689],
[-0.0120, -0.0116, 0.0043, 0.0078, 0.0115, 0.0018, 0.0085, 0.0005, 0.0012, 0.0060]]
The corresponding price series:
sage: a[0].exp()
[0.9905, 1.0025, 1.0067, 1.0016, 1.0078, 1.0051, 1.0047, 0.9987, 1.0003, 0.9957]
MORE DETAILS:
The random walk has n-th step , where is gaussian white noise of variance and is renormalized gaussian magnitude, which is given by a stationary gaussian simulation associated to a certain autocovariance sequence. See Bacry, Kozhemyak, Muzy, 2006, Continuous cascade models for asset returns for details.
Implementation of the Davies-Harte algorithm which given an autocovariance sequence (ACVS) s and an integer N, simulates N steps of the corresponding stationary Gaussian process with mean 0. We assume that a certain Fourier transform associated to s is nonnegative; if it isn’t, this algorithm fails with a NotImplementedError.
INPUT:
OUTPUT:
A list of n time series.
EXAMPLES:
We define an autocovariance sequence:
sage: N = 2^15
sage: s = [1/math.sqrt(k+1) for k in [0..N]]
sage: s[:5]
[1.0, 0.70710678118654746, 0.57735026918962584, 0.5, 0.44721359549995793]
We run the simulation:
sage: set_random_seed(0)
sage: sim = finance.stationary_gaussian_simulation(s, N)[0]
Note that indeed the autocovariance sequence approximates s well:
sage: [sim.autocovariance(i) for i in [0..4]]
[0.98665816086255..., 0.69201577095377..., 0.56234006792017..., 0.48647965409871..., 0.43667043322102...]
Warning
If you were to do the above computation with a small value of N, then the autocovariance sequence would not approximate s very well.
REFERENCES:
This is a standard algorithm that is described in several papers. It is summarized nicely with many applications at the beginning of Simulating a Class of Stationary Gaussian Processes Using the Davies-Harte Algorithm, with Application to Long Memory Processes, 2000, Peter F. Craigmile, which is easily found as a free PDF via a Google search. This paper also generalizes the algorithm to the case when all elements of s are nonpositive.
The book Wavelet Methods for Time Series Analysis by Percival and Walden also describes this algorithm, but has a typo in that they put a instead of a certain sum. That book describes exactly how to use Fourier transform. The description is in Section 7.8. Note that these pages are missing from the Google Books version of the book, but are in the Amazon.com preview of the book.