Derivatives pricing using signature payoffs

Imanol Perez Arribas

Introduction

This offers a wide range of applications in various fields, including finance. Suppose we have some (possibly high-dimensional) price paths, and we want to study the effects of an unknown real-valued continuous function on price paths. If all we know about the function is its image on a finite set of such price paths, we could infer more knowledge about it by taking a look to its linear effects on the signature of the price paths. We can then use this linear representation of the unkown function in terms of the signature to make predictions about the value of the function on unseen price paths.

The objective of this paper is to leverage this linear representation of functions on paths to quickly price arbitrary financial derivatives. We do so by introducing in Section 4 a family of derivatives called signature payoffs. Similarly to polynomials, these signature payoffs form an algebra: the sum of two signature payoff is a signature payoff; the product of two signature payoffs is a signature payoff. Moreover, signature payoffs are fast to price and can approximate continuous payoffs – i.e. continuous real-valued functions defined on price paths – arbitrarily well.

In Section 6 we applied the methodology to approximate prices of European options, American options, Asian options, lookback options and variance swaps using signature payoffs, and the signature approach proved to be very accurate at pricing them.

Notation

The signature of a path

We shall begin by defining the space of dd-dimensional paths with bounded pp-variation.

where the suppremum is taken over all partitions of [0,T][0,T]. The space of continuous paths of bounded pp-variation is then defined as

Let XX be a path for which the signature is well defined. The first term of the signature is given by

In other words, the first term of the signature is equal to increment of the path over the period [0,T][0,T]. The tensor X2X^{2}, on the other hand, is just the Levy area of the path: the signed area between the path and the chord that joins the initial point X0X_{0} and the ending point XTX_{T}. Higher order terms of the signature capture other aspects of the path, but their geometric intuition becomes less clear.

Signatures are unique representations of paths up to tree-like equivalences () that are parametrisation invariant ([9, Lemma 2.12]). The reduction of the infinite dimensional group of reparametrisations makes signatures a concise way of representating paths, so that it makes sense to try to describe functions on paths by taking a look to their effects on path signatures. This is the idea we followed in Section 4, where we define signature payoffs as payoffs on price paths that act through their signatures in a linear way. As we will see, this class of payoffs is rich enough to approximate continuous payoffs to arbitrary accuracy.

Similarly to the case of the polynomial map discussed in the introduction, the expected signature of a stochastic process turns out to contain enough information about the law of the process to completely determine it ().

Before defining signature payoffs, we will introduce an augmention of a one-dimensional path.

Therefore, the augmention of a path is a 33-dimensional path. A payoff function will then be defined as a real-valued mapping on augmented paths.

The payoff of a lookback call option with floating strike is defined as P(X^):=max⁡0≤t≤TXt−XT\mathscr{P}(\widehat{X}):=\max_{0\leq t\leq T}X_{t}-X_{T} for all X^∈Ap([0,T])\widehat{X}\in\mathcal{A}^{p}([0,T]).

Notice that in these two examples, the payoff function is continuous. This is true for many other derivatives such as American options, Asian options, variance swaps, etc.

Signature payoffs

The assumption that S(X^)S(\widehat{X}) is a pp-geometric rough path is used to ensure uniqueness of the signature, as per . If XX has bounded variation (i.e. p=1p=1) then the signature, defined as Riemann–Stieltjes integrals, will always be a geometric rough path. If XX is a semimartingale, we need to define S(X^)S(\widehat{X}) using Stratonovich integrals to make sure that the signature of X^\widehat{X} is a geometric rough path.

It is important to have examples of payoffs that satisfy the conditions of Theorem 4.2. Most financial derivatives, both vanilla and exotic options, are included in the framework of Theorem 4.2. Some examples include European options, American options, Asian options, lookback options and variance swaps. Notice that barrier options, for example, are not included in this setting in principle because the discontinuity on the barrier implies that they do not satisfy the continuity hypothesis of Theorem 4.2. However, one could fix this by smoothening barrier payoffs – for example, with barrier bending.

A major application of Theorem 4.2 is that as we have found a family of derivatives that approximates continuous payoffs, we could price a given derivative or a basket of derivatives by approximation using signature payoffs. As we will see in Section 5, signature payoffs can be quickly priced. Therefore, pricing a derivative by approximation with signatures becomes especially interesting when pricing the derivative is computationally expensive. In Section 6 we will see that this approach to pricing arbitrary continuous payoffs is extremely accurate.

Pricing signature payoffs

However, since the signature payoffs were defined to be linear functions on signatures, we have:

Suppose now that one has a basket of derivatives {P1,…,PN}\{\mathscr{P}_{1},\ldots,\mathscr{P}_{N}\}. These payoffs could be expensive to price, so that pricing the entire basket would take a considerable amount of time. Using Theorem 4.2, one could use the following procedure to price the basket:

with rr and σ\sigma the interest rate and volatility respectively, which are assumed to be constant, and WW a standard Brownian motion. The augmented path X^t=(t,Xt,tTX0)\widehat{X}_{t}=\left(t,X_{t},\frac{t}{T}X_{0}\right) will then be a diffusion process satisfying the SDE

and the initial condition is X^0=(0,0,X0)\widehat{X}_{0}=(0,0,X_{0}).

It turns out that the function Φ\Phi satisfies a parabolic partial differential equation, as shown in [7, Theorem 4.7].

Let XX be a dd-dimensional Itô diffusion process

with WW a kk-dimensional Brownian motion and

Moreover, Φn\Phi_{n} satisfies the initial condition Φn(0,⋅)=0\Phi_{n}(0,\cdot)=0 for n≥1n\geq 1, and Φ0≡1\Phi_{0}\equiv 1.

Using the theorem above, we may now state and prove a theorem that gives an explicit recurrent relation for the expected signature of the Black–Scholes model.

Define En(t):=exp⁡((nr+n(n−1)σ22)t)E_{n}(t):=\exp\left(\left(nr+n(n-1)\frac{\sigma^{2}}{2}\right)t\right) for 0≤t≤T0\leq t\leq T, and given a word II with alphabet {1,2,3}\{1,2,3\} let α(I)\alpha(I) denote the number of appearences of the letters 2 or 3 in II.

Moreover, the first order term of FF is given by F1(t)=te1+X0(exp⁡(rt)−1)e2+X0te3/TF_{1}(t)=te_{1}+X_{0}(\exp(rt)-1)e_{2}+X_{0}te_{3}/T, and for n≥2n\geq 2 the nth term of FF is characterised by the following:

π(1,i2,…,in)(F(t))=∫0tEα((1,i2,…,in))(s)π(i2,…,in)(F(t−s))ds\pi^{(1,i_{2},\ldots,i_{n})}(F(t))=\int_{0}^{t}E_{\alpha((1,i_{2},\ldots,i_{n}))}(s)\pi^{(i_{2},\ldots,i_{n})}(F(t-s))ds.

π(2,1,i3,…,in)(F(t))=(r+σ2α((1,i3,…,in)))∫0tEα((2,1,i3,…,in))(s)π(1,i3,…,in)(F(t−s))ds\pi^{(2,1,i_{3},\ldots,i_{n})}(F(t))=(r+\sigma^{2}\alpha((1,i_{3},\ldots,i_{n})))\int_{0}^{t}E_{\alpha((2,1,i_{3},\ldots,i_{n}))}(s)\pi^{(1,i_{3},\ldots,i_{n})}(F(t-s))ds.

π(2,2,i3,…,in)(F(t))=(r+σ2α((2,i3,…,in)))∫0tEα((2,2,i3,…,in))π(2,i3,…,in)(F(t−s))ds\pi^{(2,2,i_{3},\ldots,i_{n})}(F(t))=(r+\sigma^{2}\alpha((2,i_{3},\ldots,i_{n})))\int_{0}^{t}E_{\alpha((2,2,i_{3},\ldots,i_{n}))}\pi^{(2,i_{3},\ldots,i_{n})}(F(t-s))ds +12σ2∫0tEα((2,2,i3,…,in))(s)π(i3,…,in)(F(t−s))ds+\frac{1}{2}\sigma^{2}\int_{0}^{t}E_{\alpha((2,2,i_{3},\ldots,i_{n}))}(s)\pi^{(i_{3},\ldots,i_{n})}(F(t-s))ds.

π(3,i2,…,in)(F(t))=1Tπ(1,i2,…,in)(F(t))\pi^{(3,i_{2},\ldots,i_{n})}(F(t))=\frac{1}{T}\pi^{(1,i_{2},\ldots,i_{n})}(F(t)).

Assume XX follows the dynamics given in (3). Then, XX is of the form

where AA is the infinitesimal generator of X^\widehat{X}. Moreover, Φ0≡1\Phi_{0}\equiv 1 and Φn(0,⋅)=0\Phi_{n}(0,\cdot)=0.

Set fn(t,x)=(e1+rxe2+X0e3/T)⊗Φn−1(t,x)+σ2x2e2⊗∂xΦn−1(t,x)+12σ2x2e2⊗e2⊗Φn−2(t,x)f_{n}(t,x)=(e_{1}+rxe_{2}+X_{0}e_{3}/T)\otimes\Phi_{n-1}(t,x)+\sigma^{2}x^{2}e_{2}\otimes\partial_{x}\Phi_{n-1}(t,x)+\dfrac{1}{2}\sigma^{2}x^{2}e_{2}\otimes e_{2}\otimes\Phi_{n-2}(t,x) for n≥2n\geq 2, and f1(t,x)=e1+rxe2+X0e3/Tf_{1}(t,x)=e_{1}+rxe_{2}+X_{0}e_{3}/T. By Feynman-Kac formula, Φn\Phi_{n} will then be given by

Therefore, for I=(1)I=(1), I=(2)I=(2) or I=(3)I=(3), we have πI(Φ(t,X0))=X0α(I)πI(F(t))\pi^{I}(\Phi(t,X_{0}))=X_{0}^{\alpha(I)}\pi^{I}(F(t)) with π(1)(F(t)):=t\pi^{(1)}(F(t)):=t, π(2)(F(t)):=exp⁡(rt)−1\pi^{(2)}(F(t)):=\exp(rt)-1 and π(3)(F(t)):=t/T\pi^{(3)}(F(t)):=t/T.

Assume now that πI(Φ(t,X0))=X0α(I)πI(F(t))\pi^{I}(\Phi(t,X_{0}))=X_{0}^{\alpha(I)}\pi^{I}(F(t)) holds for all words II of length strictly less than nn. We will show that it also holds for words of length nn, with n≥2n\geq 2. We will distinguish four steps.

Case 1: I=(1,i2,…,in)I=(1,i_{2},\ldots,i_{n}). In this case π(1,i2,…,in)(fn(t,x))=π(i2,…,in)(Φn−1(t,x))=xα((i2,…,in))π(i2,…,in)(F(t))\pi^{(1,i_{2},\ldots,i_{n})}(f_{n}(t,x))=\pi^{(i_{2},\ldots,i_{n})}(\Phi_{n-1}(t,x))=\\ x^{\alpha((i_{2},\ldots,i_{n}))}\pi^{(i_{2},\ldots,i_{n})}(F(t)), so that π(1,i2,…,in)(Φn(t,X0))=X0α((1,i2,…,in))π(1,i2,…,in)(F(t))\pi^{(1,i_{2},\ldots,i_{n})}(\Phi_{n}(t,X_{0}))=X_{0}^{\alpha((1,i_{2},\ldots,i_{n}))}\pi^{(1,i_{2},\ldots,i_{n})}(F(t)) with π(1,i2,…,in)(F(t)):=∫0tEα((1,i2,…,in))(s)π(i2,…,in)(F(t−s))ds\pi^{(1,i_{2},\ldots,i_{n})}(F(t)):=\int_{0}^{t}E_{\alpha((1,i_{2},\ldots,i_{n}))}(s)\pi^{(i_{2},\ldots,i_{n})}(F(t-s))ds.

Case 2: I=(2,1,i3,…,in)I=(2,1,i_{3},\ldots,i_{n}). We have

Hence, π(2,1,i3,…,in)(Φn(t,X0))=X0α((2,1,i3,…,in))π(2,1,i3,…,in)(F(t))\pi^{(2,1,i_{3},\ldots,i_{n})}(\Phi_{n}(t,X_{0}))=X_{0}^{\alpha((2,1,i_{3},\ldots,i_{n}))}\pi^{(2,1,i_{3},\ldots,i_{n})}(F(t)) with

Case 3: I=(2,2,i3,…,in)I=(2,2,i_{3},\ldots,i_{n}). We have

Thus, π(2,2,i3,…,in)(Φn(t,X0))=X0α((2,2,i3,…,in))π(2,2,i3,…,in)(F(t))\pi^{(2,2,i_{3},\ldots,i_{n})}(\Phi_{n}(t,X_{0}))=X_{0}^{\alpha((2,2,i_{3},\ldots,i_{n}))}\pi^{(2,2,i_{3},\ldots,i_{n})}(F(t)) with+

Case 4: I=(3,i2,…,in)I=(3,i_{2},\ldots,i_{n}). From the definition of π(3,i2,…,in)(S(X^))\pi^{(3,i_{2},\ldots,i_{n})}(S(\widehat{X})), we have:

Theorem 5.2 provides a way to explicitly find the expected signature. Therefore, we have now found a way of explicitly pricing under the Black–Scholes model any signature payoff. Since the terms in FF can be iteratively computed up to a given order in a fast way, we can price signature payoffs inexpensively:

which corresponds to the well-known fair value of forward contracts.

In this section we studied what the fair value of a signature payoff is, under the Black–Scholes model, and we found an explict closed formula. A natural question would be whether a similar procedure can be followed for more complex models. Assuming that the price path follows some diffusion process, such as a local volatility model, one should be able to apply Theorem 5.1 again in order to get a similar explicit expresion for the fair value.

Numerical experiments

We implemented the proposed approach of pricing using signature payoffs to compute the fair price of different derivatives. We considered European call options with moneyness 99%99\%, American put options with moneyness 99%99\%, Asian options with moneyness 102%102\%, lookback options and variance swaps with strike 20%20\%, all with 1 year maturity. All of these derivatives satisfy the conditions of Theorem 4.2.

For obvious computational reasons, we had to truncate signatures so that instead of considering linear functionals on the full signature, we considered linear functionals on the truncated signature of order n≥1n\geq 1. For this experiments, we fixed n=4n=4, which for a 3-dimensional augmented path produces a linear functional of dimension 1+3+32+33+34=1211+3+3^{2}+3^{3}+3^{4}=121.

The linear functional was estimated applying linear regression against the truncated signature using a dataset of simulated market conditions following the Black–Scholes model (3). A dataset of 100100 different market conditions was used for this task – a tiny dataset for machine learning standards. We then priced the derivatives using the approximated signature payoff in an out-of-sample set of 100 market conditions, and we then compared the price we obtained with the corresponding real prices. As we see in Figure 1, the accuracy is remarkable: we obtained an R2R^{2} higher than 0.99999 in all the derivatives we considered.

Conclusion

In this paper we introduce signature payoffs, a family of derivatives that pay to the holder an amount dependent on the signature of the price process. This dependence on the signature is given by a linear functional, which makes pricing them relatively easy – the task of computing the fair value of signature payoffs is reduced to the task of computing an expected signature. As we have seen, this expected signature is easy to compute if the price process is assumed to follow a specific model. In Section 5.1 we studied the particular case of a Black–Scholes model, but the computations shown in that section can easily be extended to other models. Moreover, one could easily extend the framework to price signature payoffs for dividend-paying underlyings by suitably modifying the definition of the augmention of a path (Definition 3) to incorporate information about the dividends. A similar approach can be followed to price multi-asset signature payoffs.

The power of signature payoffs comes from the capability of signatures to approximate continuous functions on paths. In Theorem 4.2 we show that signature payoffs can approximate arbitrary continuous functions to arbitrary accuracy. This makes pricing entire baskets of derivatives quick: using an accurate representation as a signature payoff of each derivative in the basket, which can be precomputed, we can use Section 5 to approximate the fair value of each derivative in the basket by the corresponding fair values of the signature payoffs. As we saw in Section 6, this approximation turns out to be remarkably accurate: we achieved an R2R^{2} higher than 0.999990.99999 in all the payoffs we considered.

Disclaimer

Opinions and estimates constitute our judgement as of the date of this Material, are for informational purposes only and are subject to change without notice. This Material is not the product of J.P. Morgans Research Department and therefore, has not been prepared in accordance with legal requirements to promote the independence of research, including but not limited to, the prohibition on the dealing ahead of the dissemination of investment research. This Material is not intended as research, a recommendation, advice, offer or solicitation for the purchase or sale of any financial product or service, or to be used in any way for evaluating the merits of participating in any transaction. It is not a research report and is not intended as such. Past performance is not indicative of future results. Please consult your own advisors regarding legal, tax, accounting or any other aspects including suitability implications for your particular circumstances. J.P. Morgan disclaims any responsibility or liability whatsoever for the quality, accuracy or completeness of the information herein, and for any reliance on, or use of this material in any way.

Important disclosures at: www.jpmorgan.com/disclosures.

References