Derivation of moment and covariance formulas for Affine Jump Diffusion processes.
Project description
ajdmom
Description
The ajdmom
package represents an advancement in the analysis of Affine Jump Diffusion (AJD) models. It enables the automated derivation of closed-form moment and covariance formulas of any order, which is facilitated by a recursive moment-computing equation developed by the authors. The package advances the analytical capabilities of AJD models, surpassing the limitations of closed-form conditional Characteristic Function (CF). This advancement constitutes the core contribution of our work.
Currently, ajdmom
supports computations of moments, central moments, and covariances for the Heston Stochastic Volatility model and its three AJD extensions. Moreover, the package can compute partial derivatives of these quantities with respect to model parameters.
The moments and covariances obtained through ajdmom
have far-reaching implications for multiple domains, including financial modelling, simulation and parameter estimation. For simulations, these moments can establish the underlying probability distributions, leading to significant reductions in computational time when contrasted with conventional numerical CF inversion techniques. In parameter estimation, the moments serve to formulate explicit moment estimators while the likelihood functions are not analytically solvable. Consequently, ajdmom
has the potential to become an essential instrument for researchers and practitioners demanding comprehensive AJD model analysis.
Simple Usage
To get the formula for the first moment $\mathbb{E}[y_n]$ for the Heston Stochastic Volatility model ( $y_n$ denotes the return over the nth interval of length $h$ ), run the following code snippet:
from ajdmom import mdl_1fsv # mdl_1fsv -> mdl_1fsvj, mdl_2fsv, mdl_2fsvj
from pprint import pprint
m1 = mdl_1fsv.moment_y(1) # 1 in moment_y(1) -> 2,3,4...
# moment_y() -> cmoment_y() : central moment
# dpoly(m1, wrt), wrt = 'k','theta',... : partial derivative
msg = "which is a Poly with attribute keyfor = \n{}"
print("moment_y(1) = "); pprint(m1); print(msg.format(m1.keyfor))
which produces:
moment_y(1) =
{(0, 1, 0, 0, 1, 0, 0, 0): Fraction(-1, 2),
(0, 1, 0, 1, 0, 0, 0, 0): Fraction(1, 1)}
which is a Poly with attribute keyfor =
('e^{-kh}', 'h', 'k^{-}', 'mu', 'theta', 'sigma_v', 'rho', 'sqrt(1-rho^2)')
Within the produced results, the two principal key-value pairs, namely (0,1,0,0,1,0,0,0): Fraction(-1,2) and (0,1,0,1,0,0,0,0): Fraction(1,1), correspond to the following expressions:
$$ -\frac{1}{2}\times e^{-0kh}h^1k^{-0}\mu^0\theta^1\sigma_v^0\rho^0\left(\sqrt{1-\rho^2}\right)^0, $$
$$ 1\times e^{-0kh}h^1k^{-0}\mu^1\theta^0\sigma_v^0\rho^0\left(\sqrt{1-\rho^2}\right)^0, $$
respectively. The summation of these terms yields the first moment of the One-Factor SV model: $\mathbb{E}[y_n] = (\mu-\theta/2)h$. This demonstrates that the ajdmom
package successfully encapsulates the model's dynamics into a computationally manipulable form, specifically leveraging a custom dictionary data structure, referred to as Poly
, to encode the moment's expression.
Documentation
The documentation is hosted on http://www.yyschools.com/ajdmom/
Ongoing Development
This code is being developed on an on-going basis at the author's Github site.
Support
For support in using this software, submit an issue.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.