Ambit field

From HandWiki

In mathematics, an ambit field is a d-dimensional random field describing the stochastic properties of a given system. The input is in general a d-dimensional vector (e.g. d-dimensional space or (1-dimensional) time and (d − 1)-dimensional space) assigning a real value to each of the points in the field. In its most general form, the ambit field, [math]\displaystyle{ Y }[/math], is defined by a constant plus a stochastic integral, where the integration is done with respect to a Lévy basis, plus a smooth term given by an ordinary Lebesgue integral. The integrations are done over so-called ambit sets, which is used to model the sphere of influence (hence the name, ambit, Latin for "sphere of influence" or "boundary") which affect a given point.

The use and development of ambit fields is motivated by the need of flexible stochastic models to describe turbulence[1] and the evolution of electricity prices[2] for use in e.g. risk management and derivative pricing. It was pioneered by Ole E. Barndorff-Nielsen and Jürgen Schmiegel to model turbulence and tumour growth.[1]

Note, that this article will use notation that includes time as a dimension, i.e. we consider (d − 1)-dimensional space together with 1-dimensional time. The theory and notation easily carries over to d-dimensional space (either including time herin or in a setting involving no time at all).

Intuition and motivation

In stochastic analysis, the usual way to model a random process, or field, is done by specifying the dynamics of the process through a stochastic (partial) differential equation (SPDE). It is known, that solutions of (partial) differential equations can in some cases be given as an integral of a Green's function convolved with another function – if the differential equation is stochastic, i.e. contaminated by random noise (e.g. white noise) the corresponding solution would be a stochastic integral of the Green's function. This fact motivates the reason for modelling the field of interest directly through a stochastic integral, taking a similar form as a solution through a Green's Function, instead of first specifying a SPDE and then trying to find a solution to this. This provides a very flexible and general framework for modelling a variety of phenomena.[2]

Definition

A tempo-spatial ambit field, [math]\displaystyle{ Y }[/math], is a random field in space-time [math]\displaystyle{ \chi \times \mathbb{R} }[/math] taking values in [math]\displaystyle{ \mathbb{R} }[/math]. Let [math]\displaystyle{ \mu \in \mathbb{R}, A_t (x), B_t (x) }[/math] be ambit sets in [math]\displaystyle{ \chi \times \mathbb{R}_{+}, g, q }[/math] deterministic kernel functions, [math]\displaystyle{ a }[/math] a stochastic function, [math]\displaystyle{ \sigma \geq 0 }[/math] a stochastic field (called the energy dissipation field in turbulence and volatility in finance) and [math]\displaystyle{ L }[/math] a Lévy basis. Now, the ambit field [math]\displaystyle{ Y }[/math] is

[math]\displaystyle{ Y_t(x) = \mu + \int_{A_{t}(x)} g(\eta,s,x,t)\sigma_{s}(\eta) L(d\eta,ds) + \int_{B_t(x)} q(\eta,s,x,t)a_{s}(\eta) \, d\eta \, ds }[/math]

Ambit sets

In the above, the ambit sets [math]\displaystyle{ A_t(x) }[/math] and [math]\displaystyle{ B_t(x) }[/math] describe the sphere of influence for a given point in space-time. I.e. at a given point, [math]\displaystyle{ (t,x) \in \chi \times \mathbb{R} }[/math] the sets [math]\displaystyle{ A_t(x) }[/math] and [math]\displaystyle{ B_t(x) }[/math] are the points in space-time which affect the value of the ambit field at [math]\displaystyle{ (t,x), Y_t(x) }[/math]. When time is considered as one of the dimensions, the sets are often taken to only include time-coordinates which are at or prior to the current time, t, so as to preserve causality of the field (i.e. a given point in space-time can only be affected by events that happened prior to time [math]\displaystyle{ t }[/math] and can thus not be affected by the future).

The ambit sets can be of a variety of forms and when using ambit fields for modelling purposes, the choice of ambit sets should be made in a way that captures the desired properties (e.g. stylized facts) of the system considered in the best possible way. In this sense, the sets can be used to make a particular model fit the data as closely as possible and thus provides a very flexible – yet general – way of specifying the model.

Ambit process

Often, the object of interest is not the ambit field itself, but instead a process taking a particular path through the field. Such a process is called an ambit process. As an example such a process can represent the price of a particular financial object – e.g. the price of a forward contract for a certain time and point in space, space representing things such as time to delivery, spot price, period of delivery etc.[2] This motivates the following definition:

Let the ambit field, Y, be given as above and consider a curve in space-time [math]\displaystyle{ \tau(\theta) = (x(\theta), t(\theta)) \in \chi \times \mathbb{R} }[/math]. An ambit process is defined as the value of the field along the curve, i.e.

[math]\displaystyle{ X_\theta = Y_{t(\theta)}(x(\theta)) }[/math]

Stochastic intermittency/volatility

The energy dissipation field/volatility, [math]\displaystyle{ \sigma }[/math], is, in general, stochastic (called intermittency in the context of turbulence), and can be modelled as a stochastic variable or field. Particularly, it may itself be modelled by another ambit field, i.e.

[math]\displaystyle{ \sigma^2_t(x) = \int_{C_t(x)} h(\eta,s,x,t) \tilde{L}(d\eta,ds) }[/math]

where [math]\displaystyle{ \tilde{L} }[/math] is a non-negative Lévy basis.

Integration with respect to a Lévy basis

The stochastic integral, [math]\displaystyle{ \int_{A_{t}(x)} g(\eta,s,x,t)\sigma_{s}(\eta) L(d\eta,ds) }[/math], in the definition of the ambit process is an integral of a stochastic field (the integrand) over Lévy basis (the integrator), and is thus more complicated than the usual stochastic Itô-integral. A new theory of integration was provided by Walsh (1987)[3] where integration is done with respect to random fields and this theory can be extended to integration with respect to so-called Lévy bases,[4] which is the main building block of the ambit field.

Definition of Lévy basis

A family [math]\displaystyle{ (\Lambda(A) : A \in \mathbb{B}_b(S)) }[/math] of random vectors in [math]\displaystyle{ \mathbb{R}^d }[/math] is called a Lévy basis on [math]\displaystyle{ S }[/math] if:

1. The law of [math]\displaystyle{ \Lambda(A) }[/math] is infinitely divisible for all [math]\displaystyle{ A \in \mathbb{B}_{b}(S) }[/math].
2. If [math]\displaystyle{ A_1 , A_2, \ldots, A_n \in \mathbb{B}_b(S) }[/math] are disjoint, then [math]\displaystyle{ \Lambda(A_1), \Lambda(A_2),\ldots, \Lambda(A_n) }[/math] are independent.
3. If [math]\displaystyle{ A_1 , A_2, \ldots \in \mathbb{B}_b(S) }[/math] are disjoint with [math]\displaystyle{ \bigcup_{i=1}^\infty A_i \in \mathbb{B}_b(S) }[/math], then
[math]\displaystyle{ \Lambda(\bigcup_{i=1}^\infty A_i) = \sum_{i=1}^\infty \Lambda(A_i) }[/math], a.s.

where the convergence on the right hand side of 3. is a.s.

Note that properties 2. and 3. define an independently scattered random measure.

A stationary example

In some data (e.g. commodity prices) there is often found a stationary component, which a good model should be able to capture. The ambit field can be made stationary in a straightforward way. Consider the ambit field [math]\displaystyle{ Y }[/math], defined as

[math]\displaystyle{ Y_t = \mu + \int_{A_{t}(x)} g(\eta,t-s,x)\sigma_{s}(\eta) L(d\eta,ds) + \int_{B_{t}(x)} q(\eta,t-s,x)a_{s}(\eta) \, d\eta \, ds }[/math]

where the ambit sets, [math]\displaystyle{ A_{t}(x), B_{t}(x) }[/math] are of the form [math]\displaystyle{ A_{t}(x) = A + (x,t) }[/math] where the time-coordinates of [math]\displaystyle{ A }[/math] are negative (same for [math]\displaystyle{ B }[/math]). Furthermore, we take [math]\displaystyle{ g(\eta,t,x) = q(\eta,t,x) = 0 }[/math] for [math]\displaystyle{ t \leq 0 }[/math] and that [math]\displaystyle{ \sigma }[/math] and [math]\displaystyle{ a }[/math] are also stationary random variables/fields. In particular, we can take [math]\displaystyle{ \sigma }[/math] to be a stationary ambit field itself:

[math]\displaystyle{ \sigma^2_{t}(x) = \int_{C_{t}(x)} h(\eta,t-s,x) \tilde{L}(d\eta,ds) }[/math]

where [math]\displaystyle{ \tilde{L} }[/math] is a non-negative Lévy basis and [math]\displaystyle{ h }[/math] is a positive function.

References

  1. 1.0 1.1 Barndorff-Nielsen, O. E., Schmiegel, J. "Ambit processes; with applications to turbulence and tumour growth", Research report, Thiele Centre, December 2005
  2. 2.0 2.1 2.2 Barndorff-Nielsen, O. E., Benth, F. E., and Veraart, A., "Modelling electricity forward markets by ambit fields", CREATES research center, 2010
  3. Walsh, J., "An introduction to stochastic partial differential equations", Lecture Notes in Mathematics, 1986
  4. Barndorff–Nielsen, O. E., Benth, F. E., and Veraart, A., "Ambit processes and stochastic partial differential equations", CREATES research center, 2010

External links