Exponential tilting explained

Exponential Tilting (ET), Exponential Twisting, or Exponential Change of Measure (ECM) is a distribution shifting technique used in many parts of mathematics.The different exponential tiltings of a random variable

X is known as the natural exponential family of X .

Exponential Tilting is used in Monte Carlo Estimation for rare-event simulation, and rejection and importance sampling in particular.In mathematical finance ^[1] Exponential Tilting is also known as Esscher tilting (or the Esscher transform), and often combined with indirect Edgeworth approximation and is used in such contexts as insurance futures pricing.^[2]

The earliest formalization of Exponential Tilting is often attributed to Esscher^[3] with its use in importance sampling being attributed to David Siegmund.^[4]

Overview

Given a random variable

X with probability distribution P , density f , and moment generating function (MGF) M _{X(\theta)} = E [e ^{\theta} ] < inf ty , the exponentially tilted measure P _{\theta} is defined as follows:P _{\theta(X} \in dx) = (E [e ^{\theta} I [X \in dx]])/(M _{X(\theta)}) = e ^{\theta} P (X\in dx),

where

\kappa(\theta) is the cumulant generating function (CGF) defined as \kappa(\theta) = log E [e ^{\theta} ] = log M _{X(\theta).}

We call

P _{\theta(X} \in dx) = f _{\theta(x)}

the

\theta -tilted density of X . It satisfies f _{\theta(x)} \propto e ^{\theta} f(x) .

The exponential tilting of a random vector

X has an analogous definition:P _{\theta} (X\in dx) = e ^{\theta ^{T} x - \kappa(\theta)} P (X\in dx),

where

\kappa(\theta) = log E [\exp\{\theta ^{T} X\}] .

Example

The exponentially tilted measure in many cases has the same parametric form as that of

X . One-dimensional examples include the normal distribution, the exponential distribution, the binomial distribution and the Poisson distribution.

For example, in the case of the normal distribution,

N( \mu, \sigma ^{2)} the tilted density f _{\theta(x)} is the N( \mu + \theta \sigma ^{2,} \sigma ^{2)} density. The table below provides more examples of tilted density.

Original distribution^[5] ^[6]	θ-Tilted distribution
Gamma (\alpha, \beta)	Gamma (\alpha, \beta - \theta)
Binomial (n, p)	Binomial \left(n, (pe ^{\theta})/(1- p + pe ^{\theta}) \right)
Poisson ( λ )	Poisson ( λ e ^{\theta} )
Exponential ( λ )	Exponential ( λ - \theta)
l{N}(\mu, \sigma ^{2)}	l{N}(\mu + \theta\sigma ^{2,} \sigma ^{2)}
l{N}(\mu, \Sigma)	l{N}(\mu + \Sigma\theta, \Sigma)
\chi ^{2(\kappa)}	Gamma \left( (\kappa)/(2) , (2)/(1-2 \theta) \right)

For some distributions, however, the exponentially tilted distribution does not belong to the same parametric family as

f . An example of this is the Pareto distribution with f(x) = \alpha /( 1 + x) ^{\alpha,} x > 0 , where f _{\theta(x)} is well defined for \theta < 0 but is not a standard distribution. In such examples, the random variable generation may not always be straightforward.^[7]

In statistical mechanics, the energy of a system in equilibrium with a heat bath has the Boltzmann distribution:

P (E\in dE) \propto e ^{-\beta} dE , where \beta is the inverse temperature. Exponential tilting then corresponds to changing the temperature: P _{\theta} (E\in dE) \propto e ^{-(\beta} dE .

Similarly, the energy and particle number of a system in equilibrium with a heat and particle bath has the grand canonical distribution:

P ((N, E)\in (dN, dE) ) \propto e ^{\beta} dNdE , where \mu is the chemical potential. Exponential tilting then corresponds to changing both the temperature and the chemical potential.

Advantages

In many cases, the tilted distribution belongs to the same parametric family as the original. This is particularly true when the original density belongs to the exponential family of distribution. This simplifies random variable generation during Monte-Carlo simulations. Exponential tilting may still be useful if this is not the case, though normalization must be possible and additional sampling algorithms may be needed.

In addition, there exists a simple relationship between the original and tilted CGF,

\kappa _{\theta(η)} = log ( E ^{η X} _{\theta[e} ]) = \kappa(\theta + η ) - \kappa(\theta).

We can see this by observing that

F _{\theta(x)} = \int\limits ^{x\exp\{\theta} _{inf ty} y - \kappa(\theta)\}f(y)dy.

Thus,

\begin{align} \kappa _{\theta} ( η ) & = log \int e ^{η} dF _{\theta} (x) \ & = log \int e ^{η} e ^{\theta} dF(x) \\ & = log E [e ^{(η+\theta)X-\kappa(\theta)} ] \\ & = log (e ^{\kappa(η+\theta)-\kappa(\theta)} ) \\ & = \kappa( η+ \theta) - \kappa(\theta) \end{align} .

Clearly, this relationship allows for easy calculation of the CGF of the tilted distribution and thus the distributions moments. Moreover, it results in a simple form of the likelihood ratio. Specifically,

\ell = (d P)/(d P _{\theta}) = (f(x))/(f _{\theta} (x)) = e ^{-} .

Properties

• If

\kappa( η ) =log E [\exp( η X)] is the CGF of X , then the CGF of the \theta -tilted X is \kappa _{\theta(η)} = \kappa(\theta + η ) - \kappa(\theta).

This means that the

i -th cumulant of the tilted X is \kappa ^{(i)} (\theta) . In particular, the expectation of the tilted distribution is E _{\theta[X]=\tfrac{d}{dη}\kappa} _{\theta(η)|} _{η=0} = \kappa'(\theta) .

The variance of the tilted distribution is

Var ^{2 }{d η} _{\theta[X] = \tfrac{d} ^{2 }\kappa} ^{\theta( η )|} _{η=0} = \kappa''(\theta) .

• Repeated tilting is additive. That is, tilting first by

\theta _{1} and then \theta _{2} is the same as tilting once by \theta _{1+\theta} _{2} .

• If

X is the sum of independent, but not necessarily identical random variables X _{1,} X _{2,} ... , then the \theta -tilted distribution of X is the sum of X _{1,} X _{2,} ... each \theta -tilted individually.

• If

D _{KL(P} \parallel P _{\theta)=E} \left[ log \tfrac{P}{P _{\theta}\right]}

between the tilted distribution

P _{\theta} and the original distribution P of X .

• Similarly, since

E _{\theta} [X] = \kappa'(\theta) , we have the Kullback-Leibler divergence asD _{KL(P} _{\theta} \parallel P) = E _{\theta} \left[ log \tfrac{P _{\theta}{P}} \right] = \theta \kappa'(\theta) - \kappa(\theta) .

Applications

Rare-event simulation

The exponential tilting of

X , assuming it exists, supplies a family of distributions that can be used as proposal distributions for acceptance-rejection sampling or importance distributions for importance sampling. One common application is sampling from a distribution conditional on a sub-region of the domain, i.e. X|X\in A . With an appropriate choice of \theta , sampling from P _{\theta} can meaningfully reduce the required amount of sampling or the variance of an estimator.

Saddlepoint approximation

The saddlepoint approximation method is a density approximation methodology often used for the distribution of sums and averages of independent, identically distributed random variables that employs Edgeworth series, but which generally performs better at extreme values. From the definition of the natural exponential family, it follows that

f _{\theta} (\bar{x}) = f(\bar{x})\exp\{n(\theta \bar{x} - \kappa(\theta))\} .

Applying the Edgeworth expansion for

f _{\theta} (\bar{x}) , we havef _{\theta} (\bar{x}) = \psi(z)( Var [\bar{X}]) ^{-1/2} \left\{ 1 + (\rho _{3(\theta)h} _{3(z)})/(6) + (\rho _{4(\theta)h} _{4(z)})/(24) ... \right\},

where

\psi(z) is the standard normal density of z = (\bar{x)/(-) \kappa _{\bar{x} } ' (\theta)}, \rho _{n(\theta)} = \kappa ^{(n)} (\theta)\{\kappa ''(\theta) ^{n/2} \} ,

and

h _{n} are the hermite polynomials.

When considering values of

\bar{x} progressively farther from the center of the distribution, |z| → inf ty and the h _{n(z)} terms become unbounded. However, for each value of \bar{x} , we can choose \theta such that \kappa '(\theta) = \bar{x}.

This value of

\theta is referred to as the saddle-point, and the above expansion is always evaluated at the expectation of the tilted distribution. This choice of \theta leads to the final representation of the approximation given byf(\bar{x}) ≈ \left( (n)/(2 \pi\kappa ''(\theta)) \right) ^{1/2} \exp\{n(\kappa(\theta) - \theta\bar{x})\}. ^{[8] [9]}

Rejection sampling

Using the tilted distribution

P _{\theta} as the proposal, the rejection sampling algorithm prescribes sampling from f _{\theta(x)} and accepting with probability (1)/(c) \exp( - \theta x + \kappa(\theta)),

where

c = \sup\limits _{x\in} (d P)/(d P _{\theta}) (x).

That is, a uniformly distributed random variable

p \sim Unif ( 0 , 1 ) is generated, and the sample from f _{\theta(x)} is accepted ifp \leq (1)/(c) \exp( - \theta x + \kappa(\theta)).

Importance sampling

Applying the exponentially tilted distribution as the importance distribution yields the equation

E (h(X)) = E _{\theta} [\ell(X)h(X)] ,

where

\ell(X) = (d P)/(d P _{\theta})

is the likelihood function. So, one samples from

f _{\theta} to estimate the probability under the importance distribution P (dX) and then multiplies it by the likelihood ratio. Moreover, we have the variance given byVar (X) = E [(\ell(X)h(X) ^{2]} .

Example

Assume independent and identically distributed

\{X _{i\}} such that \kappa(\theta) < inf ty . In order to estimate P (X _{1} + … + X _{n} > c) , we can employ importance sampling by taking h(X) = I (\sum ^{n} _{i =1} X _{i} > c) .

The constant

c can be rewritten as na for some other constant a . Then, P (\sum ^{n} _{i =1} X _{i} > na) = E _{\theta _{a}} \left[\exp\{ - \theta _{a\sum} ^{n} ^{i =1} X _{i} + n\kappa(\theta _{a)\}I(\sum} ^{n} ^{i =1} X _{i} > na) \right] ,

where

\theta _{a} denotes the \theta defined by the saddle-point equation \kappa '(\theta _{a)} = a .

Stochastic processes

Given the tilting of a normal R.V., it is intuitive that the exponential tilting of

X _{t} , a Brownian motion with drift \mu and variance \sigma ^{2} , is a Brownian motion with drift \mu + \theta\sigma ^{2} and variance \sigma ^{2} . Thus, any Brownian motion with drift under P can be thought of as a Brownian motion without drift under P _{\theta ^{*}} . To observe this, consider the process X _{t} = B _{t} + \mu _{t} . f(X _{t)} = f _{\theta ^{*}} (X _{t) (d P)/(d P _{\theta ^{*}})} = f(B _{t)\exp\{\mu} B _{T} - (1)/(2) \mu ^{2T\}} . The likelihood ratio term, \exp\{\mu B _{T} - (1)/(2) \mu ^{2} T\} , is a martingale and commonly denoted M _{T} . Thus, a Brownian motion with drift process (as well as many other continuous processes adapted to the Brownian filtration) is a P _{\theta ^{*}} -martingale.^{[10] [11]}

Stochastic Differential Equations

dX(t) = \mu(t)dt + \sigma(t) dB(t) : dX _{\theta} (t) = \mu _{\theta} (t) dt + \sigma(t) dB(t) , where \mu _{\theta} (t) = \mu(t) + \theta\sigma(t) . Girsanov's Formula states the likelihood ratio (d P)/(d P _{\theta}) = \exp\{ - \int\limits ^{T (\mu _{\theta} (t) - \mu(t))/(\sigma ^{2(t)})} _{0} dB(t) + \int\limits ^{T( (\sigma ^{2(t)})/(2)} _{0} )dt\} . Therefore, Girsanov's Formula can be used to implement importance sampling for certain SDEs.

Tilting can also be useful for simulating a process

X(t) via rejection sampling of the SDE dX(t) = \mu(X(t))dt + dB(t) . We may focus on the SDE since we know that X(t) can be written \int\limits ^{t} _{0} dX(t) + X( 0 ) . As previously stated, a Brownian motion with drift can be tilted to a Brownian motion without drift. Therefore, we choose P _{proposal} = P _{\theta ^{*}} . The likelihood ratio (d P _{\theta ^{*}})/(d P) (dX(s): 0 \leq s \leq t) = \prod\limits _{\tau\geq} \exp\{\mu(X(\tau))dX(\tau) - (\mu(X(\tau)) ^{2})/(2) \}dt = \exp\{\int\limits ^{t\mu(X(\tau))dX(\tau)} _{0} - \int\limits ^{t (\mu(X(s)) ^{2})/(2)} _{0} \}dt . This likelihood ratio will be denoted M(t) . To ensure this is a true likelihood ratio, it must be shown that E [M(t)] = 1 . Assuming this condition holds, it can be shown that f _{X(t)} (y) = f ^{\theta ^{*}} _{X(t)} (y) E _{\theta ^{*}} [M(t)|X(t) = y] . So, rejection sampling prescribes that one samples from a standard Brownian motion and accept with probability (f _{X(t)} (y))/(f ^{\theta ^{*}} _{X(t)} (y)) (1)/(c) = (1)/(c) E _{\theta ^{*}} [M(t)|X(t) = y] .

Choice of tilting parameter

Siegmund's algorithm

Assume i.i.d. X's with light tailed distribution and

E [X] > 0 . In order to estimate \psi(c) = P (\tau(c) < inf ty) where \tau(c) = inf \{t:\sum\limits ^{t} _{i =1} X _{i>} c\} , when c is large and hence \psi(c) small, the algorithm uses exponential tilting to derive the importance distribution. The algorithm is used in many aspects, such as sequential tests,^[12] G/G/1 queue waiting times, and \psi is used as the probability of ultimate ruin in ruin theory. In this context, it is logical to ensure that P _{\theta(\tau(c)} < inf ty) = 1 . The criterion \theta > \theta _{0} , where \theta _{0} is s.t. \kappa'(\theta _{0)} = 0 achieves this. Siegmund's algorithm uses \theta = \theta ^{*} , if it exists, where \theta ^{*} is defined in the following way: \kappa(\theta ^{*)} = 0 . It has been shown that \theta ^{*} is the only tilting parameter producing bounded relative error (\underset{x → inf ty}{\lim\sup} (Var I _{A(x)})/(P A(x) ^{2}) < inf ty ).^[13]

Black-Box algorithms

We can only see the input and output of a black box, without knowing its structure. The algorithm is to use only minimal information on its structure. When we generate random numbers, the output may not be within the same common parametric class, such as normal or exponential distributions. An automated way may be used to perform ECM. Let

X _{1,} X _{2,...} be i.i.d. r.v.’s with distribution G ; for simplicity we assume X\geq 0 . Define ak{F} _{n} = \sigma(X _{1,...,X} _{n,U} _{1,...,} U _{n)} , where U _{1,} U _{2} , . . . are independent (0, 1) uniforms. A randomized stopping time for X _{1,} X _{2} , . . . is then a stopping time w.r.t. the filtration \{ ak{F} _{n\}} , . . . Let further ak{G} be a class of distributions G on [ 0 , inf ty) with k _{G} = \int ^{inf ty} _{0} e ^{\theta} G(dx) < inf ty and define G _{\theta} by (dG _{\theta})/(dG(x)) = e ^{\theta x - k _{G}} . We define a black-box algorithm for ECM for the given \theta and the given class ak{G} of distributions as a pair of a randomized stopping time \tau and an ak{F} _{\tau-} measurable r.v. Z such that Z is distributed according to G _{\theta} for any G \in ak{G} . Formally, we write this as P _{G} (Z<x) = G _{\theta} (x) for all x . In other words, the rules of the game are that the algorithm may usesimulated values from G and additional uniforms to produce an r.v. from G _{\theta} .^[14]

See also

• Importance sampling
• Rejection sampling
• Monte Carlo method
• Exponential family
• Esscher transform

Notes and References

1. H.U. Gerber & E.S.W. Shiu. 1994. Option pricing by Esscher transforms. Transactions of the Society of Actuaries. 46. 99–191.
2. Book: Cruz, Marcelo. Fundamental Aspects of Operational Risk and Insurance Analytics. Wiley. 2015. 978-1-118-11839-9 . 784–796.
3. Book: Butler, Ronald. Saddlepoint Approximations with Applications. limited. Cambridge University Press. 2007. 9780521872508. 156.
4. Siegmund. D. . 1976 . Importance Sampling in the Monte Carlo Study of Sequential Tests . . 4 . 4 . 673–684. 10.1214/aos/1176343541 . free.
5. Book: Asmussen Soren & Glynn Peter. Stochastic Simulation. Springer. 2007. 978-0-387-30679-7. 130.
6. Fuh. Cheng-Der. Teng. Huei-Wen. Wang. Ren-Her. Efficient Importance Sampling for Rare Event Simulation with Applications. 2013.
7. Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 164–167.
8. Book: Butler, Ronald. Saddlepoint Approximations with Applications. limited. Cambridge University Press. 2007. 9780521872508. 156–157.
9. Book: Seeber, G.U.H.. Advances in GLIM and Statistical Modelling. Springer. 1992. 978-0-387-97873-4 . 195–200.
10. Book: Asmussen Soren & Glynn Peter. Stochastic Simulation. Springer. 2007. 978-0-387-30679-7. 407.
11. Book: Steele, J. Michael. Stochastic Calculus and Financial Applications. limited. Springer. 2001. 978-1-4419-2862-7 . 213–229.
12. D. Siegmund (1985) Sequential Analysis. Springer-Verlag
13. Book: Asmussen Soren & Glynn Peter, Peter. Stochastic Simulation. Springer. 2007. 978-0-387-30679-7. 164–167.
14. Asmussen, Soren & Glynn, Peter (2007). Stochastic Simulation. Springer. pp. 416–420.