Английская Википедия:Additive process
An additive process, in probability theory, is a cadlag, continuous in probability stochastic process with independent increments. An additive process is the generalization of a Lévy process (a Lévy process is an additive process with identically distributed increments). An example of an additive process that is not a Lévy process is a Brownian motion with a time-dependent drift.Шаблон:Sfn The additive process was introduced by Paul Lévy in 1937.Шаблон:Sfn
There are applications of the additive process in quantitative financeШаблон:Sfn (this family of processes can capture important features of the implied volatility) and in digital image processing.Шаблон:Sfn
Definition
An additive process is a generalization of a Lévy process obtained relaxing the hypothesis of identically distributed increments. Thanks to this feature an additive process can describe more complex phenomenons than a Lévy process.
A stochastic process <math>\{X_t\}_{t \geq 0}</math> on <math>\mathbb R^d</math> such that <math>X_0=0 </math> almost surely is an additive process if it satisfy the following hypothesis:
- It has independent increments.
- It is continuous in probability.Шаблон:Sfn
Main properties
Independent increments
A stochastic process <math>\{X_t\}_{t \geq 0}</math> has independent increments if and only if for any <math>0\leq p<r\leq s<t </math> the random variable <math>X_t-X_s</math> is independent from the random variable <math>X_r-X_p</math>.Шаблон:SfnШаблон:What
Continuity in probability
A stochastic process <math>\{X_t\}_{t \geq 0}</math> is continuous in probability if, and only if, for any <math>t>0</math>
- <math>\lim_{s \to t^-} \Pr \left(\big| X_s- X_t \big| \geq \varepsilon \right) = 0.</math>Шаблон:Sfn
Lévy–Khintchine representation
There is a strong link between additive process and infinitely divisible distributions. An additive process at time <math>t</math> has an infinitely divisible distribution characterized by the generating triplet <math>(\gamma_t, A_t, \nu_t)</math>. <math> \gamma_t</math> is a vector in <math>\mathbb R^d</math>, <math> A_t</math> is a matrix in <math>\mathbb R^{d\times d}</math> and <math>\nu_t</math> is a measure on <math>\mathbb R^d</math> such that <math>\nu_t(\{0\})=0 </math> and <math>\int_{\mathbb R^d}(1\wedge x^2)\nu_t(dx)<\infty</math>. Шаблон:Sfn
<math> \gamma_t</math> is called drift term, <math> A_t</math> covariance matrix and <math>\nu_t</math> Lévy measure. It is possible to write explicitly the additive process characteristic function using the Lévy–Khintchine formula:
- <math>\varphi_X(u)(t) := \operatorname E \left[e^{iu'X_t}\right] = \exp \left(u' \gamma_t i - \frac{1}{2}u' A_t u + \int_{\mathbb R^d} \left(e^{i u' x}-1 -iu'x\mathbf{I}_{|x|<1}\right)\,\nu_t(dx) \right),
</math>
where <math>u</math> is a vector in <math>\mathbb R^d</math> and <math>\mathbf{I_C}</math> is the indicator function of the set <math>C</math>.Шаблон:Sfn
A Lèvy process characteristic function has the same structure but with <math>\gamma_t =t\gamma, \nu_t = t\nu</math> and <math>A_t = At</math> with <math>\gamma</math> a vector in <math>\mathbb R^d</math>, <math>A</math> a positive definite matrix in <math>\mathbb R^{d \times d}</math> and <math>\nu</math> is a measure on <math>\mathbb R^d</math>.Шаблон:Sfn
Existence and uniqueness in law of additive process
The following result together with the Lévy–Khintchine formula characterizes the additive process.
Let <math>\{X_t\}_{t \geq 0}</math> be an additive process on <math>\mathbb R^d</math>. Then, its infinitely divisible distribution is such that:
- For all <math>t</math>, <math>A_t</math> is a positive definite matrix.
- <math>\gamma_0=0, A_0=0, \nu_0=0</math> and for all <math>s, t</math> is such that <math>s<t</math>, <math>A_t-A_s</math> is a positive definite matrix and <math>\nu_t(B)\geq \nu_s(B)</math> for every <math>B</math> in <math>\mathbf{B}(\mathbb R^d)</math>.
- If <math> s\to t</math> <math>\gamma_s\to \gamma_t, A_s \to A_t</math> and <math>\nu_s(B)\to \nu_t(B)</math> every <math>B</math> in <math>\mathbf{B}(\mathbb R^d)</math>, <math>0\not\in B</math>.
Conversely for family of infinitely divisible distributions characterized by a generating triplet <math>(\gamma_t, A_t, \nu_t)</math> that satisfies 1, 2 and 3, it exists an additive process <math>\{X_t\}_{t \geq 0}</math> with this distribution.Шаблон:SfnШаблон:Sfn
Subclass of additive process
Additive Logistic Process
Family of additive processes with generalized logistic distribution. Their 5 parameters characteristic function is
- <math> \operatorname E \left[e^{iuX_t}\right] = \left(\frac{B(\alpha_t+i \sigma_t u, \beta_t- i \sigma u )}{B(\alpha_t, \beta_t)} \right)^{\delta_t} e^{ i \mu_t u}\;\;. </math>
Two subcases of additive logistic process are the symmetric logistic additive process with standard logistic distribution (<math>\alpha_t=1 </math>, <math>\beta_t=1 </math>, <math>\delta_t=1 </math>) and the conjugate-power Dagum additive process with Dagum distribution (<math>\alpha_t=1 </math>, <math>\beta_t=1-\sigma(t) </math>, <math>\alpha_t=1 </math>).
The function <math>\mu_t </math> can always be chosen s.t. the additive process is a martingale.Шаблон:Sfn
Additive Normal Tempered Stable Process
Extension of the Lévy normal tempered stable processes; some well-known Lévy normal tempered stable processes have normal-inverse Gaussian distribution and the variance-gamma distribution. Additive normal tempered stable processesШаблон:Sfn have the same characteristic function of Lévy normal tempered stable processes but with time dependent parameters <math>\sigma_t</math> (the level of the volatility), <math>k_t</math> (the variance of jumps) and <math>\eta_t</math> (linked to the skew):
- <math> \operatorname E \left[e^{iuX_t}\right] = {\cal L}_t \left(iu \left(\frac{1}{2}+\eta_t \right)\sigma_t^2+\frac{u^2\sigma^2_t}{2};\;k_t,\;\alpha \right)e^{iu\varphi_tt}, </math>
where
- <math> \ln {\cal L}_t \left(u;\;k_t,\;\alpha\right) :=
\begin{cases} \displaystyle \frac{t}{k_t} \displaystyle \frac{1-\alpha}{\alpha} \left \{1- \left(1+\frac{u \; k_t}{1-\alpha}\right)^\alpha \right \} & \mbox{if } \; 0< \alpha < 1 \\[4mm] \displaystyle -\frac{t}{k_t} \ln \left(1+u \; k_t\right) & \mbox{if } \; \alpha = 0 \end{cases} </math> The function <math>\varphi_t </math> can always be chosen s.t. the additive process is a martingale.Шаблон:Sfn
Additive Subordinator
A positive non decreasing additive process <math>\{S_t\}_{t \geq 0}</math> with values in <math>\mathbb R</math> is an additive subordinator. An additive subordinator is a semimartingale (thanks to the fact that it is not decreasing) and it is always possible to rewrite its Laplace transform as
- <math>\operatorname E\left[ e^{-u S_t} \right] = \exp\left(u b_t + \int_{\mathbb R^d} (e^{i u x}-1) \nu_t(dx)\right). </math>Шаблон:Sfn
It is possible to use additive subordinator to time-change a Lévy process obtaining a new class of additive processes.Шаблон:Sfn
Sato Process
An additive self-similar process <math>\{Z_t\}_{t \geq 0}</math> is called Sato process.Шаблон:Sfn It is possible to construct a Sato process from a Lévy process <math>\{X_t\}_{t \geq 0}</math> such that <math>Z_t</math> has the same law of <math>t^hX_1</math>.
An example is the variance gamma SSD, the Sato process obtained starting from the variance gamma process.
The characteristic function of the Variance gamma at time <math>t=1</math> is
- <math> \operatorname E \left[e^{iuX_1}\right] = \left(\frac{1}{1-iu\theta\nu+0.5\sigma^2\nu u^2}\right)^{1/\nu}, </math>
where <math>\theta, \nu</math> and <math>\sigma</math> are positive constant.
The characteristic function of the variance gamma SSD is
- <math> \operatorname E \left[e^{iuZ_t}\right] = \left(\frac{1}{1-iut^h\theta\nu+0.5\sigma^2\nu u^2t^2h}\right)^{1/\nu}</math>Шаблон:Sfn
Simulation
Simulation of Additive process is computationally efficient thanks to the independence of increments. The additive process increments can be simulated separately and simulation can also be parallelized.Шаблон:Sfn
Jump simulation
Jump simulation is a generalization to the class of additive processes of the jump simulation technique developed for Lévy processes. The method is based on truncating small jumps below a certain threshold and simulating the finite number of independent jumps. Moreover, Gaussian approximation can be applied to replace small jumps with a diffusive term. It is also possible to use the Ziggurat algorithm to speed up the simulation of jumps.Шаблон:Sfn
Characteristic function inversion
Simulation of Lévy process via characteristic function inversion is a well established technique in the literature.Шаблон:Sfn This technique can be extended to additive processes. The key idea is obtaining an approximation of the cumulative distribution function (CDF) by inverting the characteristic function. The inversion speed is enhanced by the use of the Fast Fourier transform. Once the approximation of the CDF is available is it possible to simulate an additive process increment just by simulating a uniform random variable. The method has similar computational cost as simulating a standard geometric Brownian motion.Шаблон:Sfn
Applications
Quantitative finance
Lévy process is used to model the log-returns of market prices. Unfortunately, the stationarity of the increments does not reproduce correctly market data. A Lévy process fit well call option and put option prices (implied volatility) for a single expiration date but is unable to fit options prices with different maturities (volatility surface). The additive process introduces a deterministic non-stationarity that allows it to fit all expiration dates.Шаблон:Sfn
A four-parameters Sato process (self-similar additive process) can reproduce correctly the volatility surface (3% error on the S&P 500 equity market). This order of magnitude of error is usually obtained using models with 6-10 parameters to fit market data.Шаблон:Sfn A self-similar process correctly describes market data because of its flat skewness and excess kurtosis; empirical studies had observed this behavior in market skewness and excess kurtosis.Шаблон:Sfn Some of the processes that fit option prices with a 3% error are VGSSD, NIGSSD, MXNRSSD obtained from variance gamma process, normal inverse Gaussian process and Meixner process.Шаблон:Sfn
Additive normal tempered stable processes fit accurately equity market data ( error below 0.8% on the S&P 500 equity market) specifically for short maturities. These family of processes reproduces very well also the equity market implied volatility skew. Moreover, an interesting power scaling characteristic arises in calibrated parameters <math>k_t=\bar{k}t^\beta</math> and <math>\eta_t=\bar{\eta}t^\delta</math>. There is statistical evidence that <math>\beta=1</math> and <math>\delta=-1/2</math>.Шаблон:Sfn
Lévy subordination is used to construct new Lévy processes (for example variance gamma process and normal inverse Gaussian process). There is a large number of financial applications of processes constructed by Lévy subordination. An additive process built via additive subordination maintains the analytical tractability of a process built via Lévy subordination but it reflects better the time-inhomogeneus structure of market data.Шаблон:Sfn Additive subordination is applied to the commodity marketШаблон:Sfn and to VIX options.Шаблон:Sfn
Digital image processing
An estimator based on the minimum of an additive process can be applied to image processing. Such estimator aims to distinguish between real signal and noise in the picture pixels.Шаблон:Sfn
References
Sources
- Шаблон:Cite book
- Шаблон:Cite book
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal
- Шаблон:Cite journal