Fun with numeraires!

Author

Quasar

Published

November 28, 2024

Introduction

A proficiency in the change-of-measure technique is useful to the working quant. An excellent summary of the important results is the note Girsanov, Numeraires and all that, by Andrew Lesniewski. In this post, I would like to derive relevant results and then we can enjoy pricing some payoffs together!

Girsanov Theorem

Theorem 1 (Girsanov Theorem) Let (WP,t0) be a P standard brownian motion on (Ω,F,P) and let ϕ be any adapted process. Choose fixed T and defined the process L on [0,T] by:

(1)dLt=ϕtLtdWtP

(2)L0=1

that is:

(3)Lt=exp(0tϕsdWsP0tϕs2ds)

Assume that:

(4)EP[LT]=1

and define the new probability measure Q on Ft by:

(5)LT=dQdP

Then,

(6)dWtP=dWtQ+ϕtdt

where dWtQ is a Q-standard brownian motion.

Proof.

Our claim is that, under the Q measure, the increments (WtQWsQ) are normally distributed with mean 0 and variance (ts). We start with the special case s=0. Using moment generating functions, it is enough to show that:

It is straightforward to derive using Ito’s lemma. Let f(x)=lnx. Then, fx=1x, fxx=1x2.

d(lnLt)=1LtdLt121Lt2(dLt)2=1LtϕtLtdWtP12Lt2ϕt2Lt2dt=ϕtdWtP12ϕt2dtLt=exp(0tϕsdWsP120tϕs2ds)

To prove our main result, we will now use the MGF of the increments. For nN and (tj,jn) a partition of [0,T], with tn=T, I will show that :

(7)EQ[exp(j=0n1λj(Wtj+1QWtjQ))]=exp[j=0n1λj2(tj+1tj)]

This proves that the increments are the ones of standard brownian motion.

Let (Ftj,jn) be the filtrations of the Brownian motion at the time of the partition. The proof is by successively conditioning from tn1 to t1. We have:

EQ[exp(j=0n1λj(Wtj+1QWtjQ))]=EP[EP[Ltnexp(j=0n1λj(Wtj+1QWtjQ))|Ftn1]]=EP[j=0n2λj(Wtj+1QWtjQ)EP[Ltnexp(λn1(WtnQWtn1Q))|Ftn1]]=EP[j=0n2λj(Wtj+1QWtjQ)EP[exp(tn1tnϕsdWsQ+12tn1tnϕs2ds)exp(λn1tn1tndWsQ)|Ftn1]]=EP[j=0n2λj(Wtj+1QWtjQ)EP[exp(tn1tn(ϕs+λn1)dWsQ+12tn1tnϕs2ds)|Ftn1]]=EP[j=0n2λj(Wtj+1QWtjQ)EP[exp(tn1tn(ϕs+λn1)(dWsPθds)+12tn1tnϕs2ds)|Ftn1]]=EP[j=0n2λj(Wtj+1QWtjQ)exp(12tn1tnϕs2+tn1tnϕs(ϕs+λn1)ds) exp(12tn1tn(ϕs+λn1)2ds)]{ XdWsP  is a NP(0,E[X2]ds) gaussian random variable.} =EP[j=0n2λj(Wtj+1QWtjQ)exp(tn1tn(12ϕs2λn1ϕs+12ϕs2+λn1ϕs+λn12)ds)]=EP[j=0n2λj(Wtj+1QWtjQ)exp(λn1tn1tnds)]=exp(λn1(tntn1)EP[j=0n2λj(Wtj+1QWtjQ)]

Here, I used the fact that Ltn1 is Ftn1 measurable. I can now condition on Ftn2 down to Ft1 and proceed as above to obtain the desired result.

The process ϕt is called the Girsanov kernel.

What is a numeraire?

As Shreve puts it, a numeraire is the unit of account in which other assets are denominated. In practice, we tend to choose numeraires that simply the payoff expression.

Any strictly positive (non-dividend paying) price process can be chosen as a numeraire. A numeraire must be a tradable asset.

Consider a unit of stock worth St. It can be used as numeraire, because the price process ertSt (assume a constant short rate) is a martingale under risk-neutral measure QM. Powers of the stock price Stα cannot be used as numeraires, because their discounted values are not martingales under the risk-neutral measure. Clearly, set the short rate r=0, then EQM[ST2](EQM[ST])2=S02 by the Jensen’s inequality.

The price-process Vt of a derivative contract that pays VT=ST2 is a martingale under Q and can be used as a numeraire.

Consider the price of a contract that pays a unit sum 1 at maturity T. This instrument is the zero-coupon bond. Its an observable and tradable asset. Its price process P(t,T)=EQM[1/MT] can be used as a numeraire. QT is called the T-forward measure.

Abstract Bayes Formula

Theorem 2 (Abstract Bayes Formula) Let (Ω,F,P) be a probability space and let Q be any other probability measure on it. By the Radon-Nikodym theorem, L=dQdP, L0 with EP[L]=1. Then we have:

(8)EQ[X|G]=EP[LX|G]EP[L|G]

Proof.

By the definition of conditional expectations, recall that if W is any G-measurable random variable, then the conditional expectation must satisfy the relationship:

E[WX]=E[WE[X|G]]

It is sufficient to prove that:

EP[X|G]EQ[L|G]=EP[LX|G]

Let G be an arbitrary set in G. We have:

GEQ[X|G]EP[L|G]dP=GEP[LEQ[X|G]|G]dP{EQ[X|G] is G-measurable }=GLEQ[X|G]dP=GdQdPEQ[X|G]dP=GEQ[X|G]dQ=GXdQ

Also, we have:

GEP[LX|G]dP=GLXdP=GdQdPXdP=GXdQ

Hence, proved.

Note that, the filtration G is the same irrespective of what probability measure we construct on Ω.

Martingale property

Proposition 1 Assume that there exists a numeraire M and a probability measure QM, such that the price of any traded asset X (without intermediate payments) relative to M is a martingale under QM.That is:

XtMt=EQM{XTMT|Ft}

Let Nt be an arbitrary numeraire. Then, there exists a probability measure QN such that the price of X normalized by N is a martingale under QN.

XtNt=EQN{XTNT|Ft}

Moreover, the Radon-Nikodym derivative defining the measure QN is given by:

dQNdQM=NT/N0MT/M0

Proof.

We have:

X0=M0EQM[XTMT]

Imposing the simple fact that, the price of the derivative contract should be the same, even if we switch numeraires from M to N, we should have:

X0=N0EQN[XTNT]

Thus,

N0EQN[XTNT]=M0EQM[XTMT]NTN0×N0EQN[XTNT]=NTN0×M0 EQM[XTMT]{Multiplying both sides by NTN0}EQN[XT]=EQM[NT/N0MT/M0XT]

But, we know that:

EQN[XT]=EQM[dQNdQMXT]

Consequently, our candidate for the Radon-Nikodym derivative should be:

LT=dQNdQM=NT/N0MT/M0

Further (Xt/Nt) is a martingale under QN. Its easy to see that:

EQN[XTNT|Ft]=EQM[LTXTNT|Ft]EQM[LT|Ft]{ Abstract bayes formula }=EQM[LTXTNT|Ft]Lt=EQM[NTN0M0MTXTNT|Ft]NtN0M0Mt=EQM[NTN0M0MTXTNT|Ft]NtN0M0Mt=MtNtEQM[XTMT]=MtNtXtMt=XtNt

Since we determined the relevant likelihood process, it is easy to find the Girsanov Kernel.

Drift transformation under change of numeraire

Suppose we are interested in the dynamics of the stochastic process (Xt,t0). Under QM measure, its dynamics reads:

(9)dX(t)=μXQM(t)dt+cX(t)dWtQM

I supressed μXQM(t,Xt) as μXQM(t) for brevity.

Under the QN measure, its dynamics reads:

(10)dX(t)=μXQN(t)dt+cX(t)dWtQN

Remember that the diffusion coefficients in these equations are unaffected by the change of measure! We assume that QM is associated with the numeraire M(t) whose dynamics is given by:

dM(t)=μM(t)dt+cM(t)dWQM

and that the numeraire N has QM dynamics:

dN(t)=μN(t)dt+cN(t)dWQN

According to the Girsanov theorem, the likelihood process L(t) accompanying this change of measure is a martingale under the measure QM measure and satisfies the stochastic differential equation:

dLt=L(t)θ(t)dWtQM

Explicitly, the likelihood process L(t) is given by the stochastic exponential of the martingale 0tθsdWsQM:

L(t)=exp(0tθsdWsQM120tθs2ds)

On the other hand, from , we have:

Lt=Nt/N0Mt/M0

Differentiating using Ito’s lemma, we have:

dLt=M0N0d(NtMt)=M0N0(NtMt2dMt+1MtdNt+122NtMt3(dMt)21Mt2(dMtdNt))=M0N0(NtMt2(μM(t)dt+cM(t)dWtQM)+1Mt(μN(t)dt+cN(t)dWtQM)+NtMt3cM2(t)dt1Mt2cM(t)cN(t)dt

But since Lt is driftless, we can ignore the dt terms (whatever they are, they are bound to cancel out) and only look at the diffusion coefficient. So, we can write:

dLt=M0N0(NtMt2cM(t)+1MtcN(t))dWtQM=Nt/N0Mt/M0(cN(t)NtcM(t)Mt)dWtQM=Lt(cN(t)NtcM(t)Mt)dWtQM

Comparing this, we can infer that:

θt=cN(t)NtcM(t)Mt

Since we can write:

(11)dXt=μXP(t)dt+cX(t)dWP(t)=μXQ(t)+cX(t)dWQ(t)dWP(t)=μXQ(t)μXP(t)cX(t)+dWQ(t)

Using , we conclude that the change of drift accompanying a change of numeraire is given by:

(12)μXQ(t)μXP(t)=cX(t)(cM(t)M(t)cN(t)N(t))

Examples of numeraires

The basic component of an interest rate model is an instantaneous forward rate process f(t,s). Its value is the future instantaneous interest rate at a future time s, that is the rate for the infinitesimally short term [s,s+ds] observed at time ts.

A zero-coupon bond settling at time T0 and maturing at time T>T0 is the process:

P(t,T0,T)=exp(T0Tf(t,s)ds)

for tT0. In other words, it is the time T0 value of 1$ (without the risk of default) at T, as observed at time tT0. Its current value is given by:

Bank-account numeraire

The bank-account(money-market account) numeraire is simply the value of 1$ deposited in a bank and accruing the (credit-riskless) instantaneous interest rate. In reality, the bank credits interest to the account daily, but this can very well be approximated to a continous process. The associated stochastic price process M(t) is given by:

M(t)=exp(0tr(s)ds)

Here, the spot rate r(t) is the instantaneous forward observed at the time it settles. That is,

r(t)=f(t,t)

Forward numeraire

A zero-coupon bond(ZCB) is a simple contract with unit payoff 1$ at maturity T. By the risk-neutral valuation formula:

P(t,t,T)=V(t)=M(t)EQM[1M(T)]

So, a T-maturity ZCB is a tradable asset and its price P(t,T) can be used as a numeraire. The associated measure is called the T-foward measure QT.

The term(e.g. 3 months) forward rates for settlement at T0 and maturity at T are defined by the equation:

P(t,T0,T)=11+δF(t,T0,T)

where δ is the day-count fraction for the period [T0,T]. Re-arranging, we have:

F(t,T0,T)=1δP(t,T,T)P(t,T0,T)P(t,T0,T)F(t,T0,T)P(t,T0,T)=1δ(P(t,T,T)P(t,T0,T))

Clearly, it is a multiple of a difference P(t,T,T) and P(t,T0,T) normalized by T-maturity zero coupon bond price P(t,T0,T). So, the forward iBOR-rate must be a martingale under the T-forward measure QT.

Proposition 2 (Forward rates are QT expectations of future spot rates.) Any simply compounded forward rate spanning a time interval ending in T is a martingale under the T-forward measure.

EQT[F(t;S,T)|Fu]=F(u;S,T)

for each 0utS<T. In particular, the forward rate spanning the interval [S,T] is the QT expectation of the future simply-compounded spot rate at time S for the maturity T.

(13)EQT[L(S,T)|Ft]=F(t;S,T)

The expected value of any future instantaneous spot interest rate, under the corresponding forward measure, is equal to the related instantaneuous forward rate. That is, ETrT|Ft=f(t,T)

for each 0tT.

Proof.

We have:

EQT[rT|Ft]=1P(t,T)E[rTP(t,T)|Ft]=1P(t,T)E[rTetTr(s)ds|Ft]=1P(t,T)E[TetTr(s)ds|Ft]=1P(t,T)TE[etTr(s)ds|Ft]=1P(t,T)TP(t,T)=f(t,T)

Pricing an IRS

Consider a forward starting interest rate swap(IRS) which settles in T0 and matures in TN years from now. An IRS is a transaction between two counterparties who exchange interest rate payments on an agreed notional principal.

On a vanilla swap, a fixed-coupon interest payments are exchanged for floating rate payments. For the sake of simplicity, we assume that the payment dates on the fixed and floating leg of the swap are the same, and that the floating rate is the same as the discounting rate. The former of these assumptions is a minor simplification, made to lighten up the notation only. The latter is an important simplification, as the basis between the floating rate and the discounting rate may exhibit a complex dynamics.

Let S be the fixed-rate on the swap. By the risk-neutral valuation formula, the fixed leg value at time t can be expressed as:

Vfixed(t)=Ni=1NEQM[er(Tit)Sτ(Ti1,Ti)|Ft]=SNi=1NEQ[er(Tit)τ(Ti1,Ti)|Ft]=SN(i=1NP(t,Ti)τ(Ti1,Ti))

The floating leg can be written as:

Vfloat(t)=Ni=1NP(t,Ti)ETi[L(T,Ti1,Ti)|Ft]τ(Ti1,Ti)=Ni=1NP(t,Ti)L(t,Ti1,Ti)τ(Ti1,Ti)=Ni=1NP(t,Ti)1τ(Ti1,Ti)(P(t,Ti1)P(t,Ti)1)τ(Ti1,Ti)=Ni=1N(P(t,Ti)P(t,Ti1))=NP(t,T0)+NP(t,TN)

where the expectations are under the Ti-forward measure. Note that, I used the fact that the iBOR-rates are martingales under the forward measure.

The par-swap rate S is the fixed-rate which renders the value of the swap zero at the contract start date t.

Vfix+Vfloating=0SN(i=1NP(t,Ti)τ(Ti1,Ti))=NP(t,T0)+NP(t,TN)S(t,T0:N)=P(t,T0)+P(t,TN)(i=1NP(t,Ti)τ(Ti1,Ti))

The Annuity Measure

The annuity is an asset that pays 1$ on each coupon payment day of the swap, accrued according to the swap’s day count convention.

A(t,T0:N)=(i=1NP(t,Ti)τ(Ti1,Ti))

Since, it is a portfolio of zero coupon bonds, it is a tradable asset and its price A(t,T0:N) can be used as numeraire. This is called the Annuity numeraire and the measure QT0:N associated with this numeraire is called the (forward) swap measure. The annuity numeraire arises as the natural numeraire when valuing swaptions. It is the mechanism that allows us to link the swaption as an option on a swap to the option on the corresponding swap rate.

The forward swap rate S(t,T0:N) is a martingale in the annuity measure QT0:N.

Pricing the payoff V(T)=ST(STK)+

Under the risk-neutral measure QM associated with the money-market account numeraire M(t), the stock price S(t) evolves according to:

(14)dSt=rStdt+σStdWQM

which has the solution:

(15)St=S0exp[(rσ22)t+σWQM(t)]

By the risk-neutral pricing formula, we have:

(16)V(0)=EQM[V(T)M(T)]=EQM[STMT(STK)+]

By the change-of-measure formula, if N is any other numeraire with associated probability measure QN, we know that:

(17)EQN[V(T)]=EQM[dQNdQMV(T)]

We switch to the stock numeraire St. The Radon-Nikodym derivative LT=dQS/dQM is simply:

(18)LT=dQSdQM=S(T)/S(0)M(T)/M(0)=1S(0)S(T)M(T)

The dynamics of Lt under the QM measure is given by:

dLt=ϕtLtdWQM

where the Girsanov kernel ϕt=σ. Then, the Girsanov transformation is:

(19)WTQM=WTQS+σT

So, the QS dynamics of the stock price is:

(20)ST=S0exp[(rσ22)T+σ(WTQS+σT)]=S0exp[(r+σ22)T+σ(WTQS)]

We can develop the expression in to be:

(21)V(0)=S0EQM[1S0STMT(STK)+]

Applying the change-of-measure (abstract Baye’s formula),we have:

(22)V(0)=S0EQS[(STK)+]

Now, we have:

(23)QS(ST>K)=QS(S0exp[(r+σ22)T+σ(WTQS)]>K)=QS(S0exp[(r+σ22)T+σ(TZ)]>K)=QS(logS0K+(r+σ22)T>σTZ)

where we subbed TZ=WTQS. Recall, the standard normals Z and Z have the same distribution by symmetry.

Define

(24)d=logS0K+(r+σ22)TσT

Then, since ZNQS(0,1):

(25)QS(ST>K)=QS(Z<d)=Φ(d)

So, we can expand the expectation in . The indicator random variable 1{ST>K} is 1 for all points Z<d. So, the limits of integration will be to d.

V(0)=S0[dSTdQSKddQS]=S0[dSTfZQS(z)dzKdfZQS(z)dz]=S0[12πdS0exp[(r+σ22)TσTz]exp(z22)dzKΦ(d)]

Completing the square, we have:

V(0)=S0[12πdS0exp[12(z2+2σTz+(σT)2)]exp[(r+σ2)T]dzKΦ(d)]=S0[S0e(r+σ2)T2πdexp(12(z+σT)2)dzKΦ(d)]

Let u=z+σT. And define:

d+=d+σT=logS0K+(r+32σ2)σT

We have the closed-form formula:

V(0)=S02e(r+σ2)TΦ(d+)KS0Φ(d)

References