Modern Stochastics: Theory and Applications logo


  • Help
Login Register

  1. Home
  2. Issues
  3. Volume 5, Issue 1 (2018)
  4. Cliquet option pricing with Meixner proc ...

Cliquet option pricing with Meixner processes
Volume 5, Issue 1 (2018), pp. 81–97
Markus Hess  

Authors

 
Placeholder
https://doi.org/10.15559/18-VMSTA96
Pub. online: 12 February 2018      Type: Research Article      Open accessOpen Access

Received
27 September 2017
Revised
2 January 2018
Accepted
21 January 2018
Published
12 February 2018

Abstract

We investigate the pricing of cliquet options in a geometric Meixner model. The considered option is of monthly sum cap style while the underlying stock price model is driven by a pure-jump Meixner–Lévy process yielding Meixner distributed log-returns. In this setting, we infer semi-analytic expressions for the cliquet option price by using the probability distribution function of the driving Meixner–Lévy process and by an application of Fourier transform techniques. In an introductory section, we compile various facts on the Meixner distribution and the related class of Meixner–Lévy processes. We also propose a customized measure change preserving the Meixner distribution of any Meixner process.

1 Introduction

Cliquet option based contracts constitute a customized subclass of equity indexed annuities. The underlying options commonly are of monthly sum cap style paying a credited yield based on the sum of monthly-capped rates associated with some reference stock index. In this regard, cliquet type investments belong to the class of path-dependent exotic options. In [15] cliquet options are regarded as “the height of fashion in the world of equity derivatives”. In the literature, there are different pricing approaches for cliquet options involving e.g. partial differential equations (see [15]), Monte Carlo techniques (see [2]), numerical recursive algorithms related to inverse Laplace transforms (see [9]) and analytical computation methods (see [3, 7, 8]). The present article belongs to the last category.
The aim of the present paper is to provide analytical pricing formulas for globally-floored locally-capped cliquet options with multiple resetting times where the underlying reference stock index is driven by a pure-jump time-homogeneous Meixner–Lévy process. In this setup, we derive cliquet option price formulas under two different approaches: once by using the distribution function of the driving Meixner–Lévy process and once by applying Fourier transform techniques (as proposed in [8]). All in all, the present article can be seen as an accompanying (but to a large degree self-contained) paper to [8], as it presents a specific application of the results derived in [8] to the class of Meixner–Lévy processes.
The paper is organized as follows: In Section 2 we compile facts on the Meixner distribution and the related class of stochastic Meixner–Lévy processes. In Section 3 we introduce a geometric pure-jump stock price model driven by a Meixner–Lévy process. In Section 3.1 we establish a customized structure preserving measure change from the risk-neutral to the physical probability measure. Section 4 is dedicated to the pricing of cliquet options. We obtain semi-analytic expressions for the cliquet option price by using the probability distribution function of the driving Meixner–Lévy process in Section 4.1 and by an application of Fourier transform techniques in Section 4.2. In Section 5 we draw the conclusions.

2 A review of Meixner processes

Let $(\Omega ,\mathbb{F},(\mathcal{F}_{t})_{t\in [0,T]},\mathbb{Q})$ be a filtered probability space satisfying the usual hypotheses, i.e. $\mathcal{F}_{t}=\mathcal{F}_{t+}:=\cap _{s>t}\mathcal{F}_{s}$ constitutes a right-continuous filtration and $\mathbb{F}$ denotes the sigma-algebra augmented by all $\mathbb{Q}$-null sets (cf. p. 3 in [10]). Here, $\mathbb{Q}$ is a risk-neutral probability measure and $0<T<\infty $ denotes a finite time horizon. In the following, we compile various facts on the Meixner distribution and Meixner–Lévy processes from [1, 6, 12, 13] and [14].
A real-valued, càdlàg, pure-jump, time-homogeneous Lévy process $M=(M_{t})_{t\in [0,T]}$ (with independent and stationary increments) satisfying $M_{0}=0$ is called Meixner (-Lévy) process with scaling parameter $\alpha >0$, shape/skewness parameter $\beta \in (-\pi ,\pi )$, peakedness parameter $\delta >0$ and location parameter $\mu \in \mathbb{R}$, if $M_{t}$ possesses the Lévy–Itô decomposition
(2.1)
\[ M_{t}=\theta t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{Q}}(s,z)\]
where $\mathbb{R}_{0}:=\mathbb{R}\setminus \{0\}$, the drift parameter
(2.2)
\[ \theta :=\mu +\delta \alpha \tan (\beta /2)\]
is a real-valued constant and the $\mathbb{Q}$-compensated Poisson random measure (PRM) is given by
(2.3)
\[ d{\tilde{N}}^{\mathbb{Q}}(s,z):=dN(s,z)-d{\nu }^{\mathbb{Q}}(z)ds\]
with positive and finite Meixner-type Lévy measure
(2.4)
\[ d{\nu }^{\mathbb{Q}}(z):=\delta \frac{{e}^{\beta z/\alpha }}{z\sinh (\pi z/\alpha )}dz\]
(cf. [12, 13], Eq. (3) in [6]) satisfying ${\nu }^{\mathbb{Q}}(\{0\})=0$ and
\[ \int _{\mathbb{R}_{0}}\big(1\wedge {z}^{2}\big)d{\nu }^{\mathbb{Q}}(z)<\infty .\]
We denote the Lévy triplet of $M_{t}$ by $(\theta ,0,{\nu }^{\mathbb{Q}})$. (Note that this notation is not entirely consistent with [6, 8].) We recall that $M_{t}$ possesses moments of all orders (cf. Section 5.3.10 in [13]). Evidently, $M_{t}$ has no Brownian motion part. Since
\[ \int _{\mathbb{R}_{0}}|z|d{\nu }^{\mathbb{Q}}(z)=\infty \]
the process $M_{t}$ possesses infinite variation (cf. Section 5.3.10 in [13]). We write for any fixed $t\in [0,T]$
\[ M_{t}\sim \mathcal{M}(\alpha ,\beta ,\delta t,\mu t)\]
(cf. Section 3.6 in [1]) and say that M is Meixner distributed under $\mathbb{Q}$ with parameters α, β, δ and μ. From (2.1) and (2.2) we instantly receive the mean value
(2.5)
\[ \mathbb{E}_{\mathbb{Q}}[M_{t}]=\theta t=\mu t+\delta t\alpha \tan (\beta /2)\]
standing in accordance with Eq. (11) in [6]. The variance, skewness and kurtosis of $M_{t}$ are respectively given by
\[\begin{array}{r@{\hskip0pt}l}\displaystyle \mathbb{Var}_{\mathbb{Q}}[M_{t}]& \displaystyle =\frac{\delta t}{2}\frac{{\alpha }^{2}}{{\cos }^{2}(\beta /2)},\hspace{2em}\mathbb{S}_{\mathbb{Q}}[M_{t}]=\sqrt{2/(\delta t)}\sin (\beta /2),\\{} \displaystyle \mathbb{K}_{\mathbb{Q}}[M_{t}]& \displaystyle =3+\frac{2-\cos (\beta )}{\delta t}\end{array}\]
(cf. Table 6 in [1]). Furthermore, for all $x\in \mathbb{R}$ and $t\in [0,T]$ the real-valued probability density function (pdf) of $M_{t}$ under $\mathbb{Q}$ reads as
(2.6)
\[ f_{M_{t}}(x):=\frac{{(2\cos (\beta /2))}^{2\delta t}}{2\pi \alpha \varGamma (2\delta t)}{e}^{\beta (x-\mu t)/\alpha }{\bigg|\varGamma \bigg(\delta t+i\frac{x-\mu t}{\alpha }\bigg)\bigg|}^{2}\]
(cf. [1, 12], Eq. (4) in [6]) wherein
\[ \varGamma (\zeta ):={\int _{0}^{\infty }}{u}^{\zeta -1}{e}^{-u}du\]
denotes the gamma function which is defined for all $\zeta \in \mathbb{C}$ with $Re(\zeta )>0$. Taking the definition of the gamma function and Euler’s formula into account, we get
\[\begin{array}{r@{\hskip0pt}l}\displaystyle \varGamma \bigg(\delta t+i\frac{x-\mu t}{\alpha }\bigg)& \displaystyle ={\int _{{0}^{+}}^{\infty }}{u}^{\delta t-1}{e}^{-u}\cos \bigg(\frac{x-\mu t}{\alpha }\ln u\bigg)du\\{} & \displaystyle \hspace{1em}+i{\int _{{0}^{+}}^{\infty }}{u}^{\delta t-1}{e}^{-u}\sin \bigg(\frac{x-\mu t}{\alpha }\ln u\bigg)du\end{array}\]
which implies
\[\begin{array}{r@{\hskip0pt}l}\displaystyle {\bigg|\varGamma \bigg(\delta t+i\frac{x-\mu t}{\alpha }\bigg)\bigg|}^{2}& \displaystyle ={\Bigg({\int _{{0}^{+}}^{\infty }}{u}^{\delta t-1}{e}^{-u}\cos \bigg(\frac{x-\mu t}{\alpha }\ln u\bigg)du\Bigg)}^{2}\\{} & \displaystyle \hspace{1em}+{\Bigg({\int _{{0}^{+}}^{\infty }}{u}^{\delta t-1}{e}^{-u}\sin \bigg(\frac{x-\mu t}{\alpha }\ln u\bigg)du\Bigg)}^{2}.\end{array}\]
Note that the latter object appears in (2.6). The cumulative distribution function (cdf) of $M_{t}$ does not possess a closed form representation but it can be computed numerically. Further on, the characteristic function of $M_{t}$ can be computed by the Lévy–Khinchin formula (see e.g. [4, 5, 11, 13]) due to
(2.7)
\[ \phi _{M_{t}}(u):=\mathbb{E}_{\mathbb{Q}}\big[{e}^{iuM_{t}}\big]={e}^{\psi (u)t}\]
with ${i}^{2}=-1$, $u\in \mathbb{R}$, $t\in [0,T]$ and a characteristic exponent
(2.8)
\[ \psi (u):=iu\bigg[\mu +\delta \alpha \tan \bigg(\frac{\beta }{2}\bigg)\bigg]+\delta \int _{\mathbb{R}_{0}}\frac{{e}^{iuz}-1-iuz}{z}\frac{{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz.\]
Moreover, let us define the Fourier transform, respectively inverse Fourier transform, of a deterministic function $q\in {\mathcal{L}}^{1}(\mathbb{R})$ via
\[ \hat{q}(y):=\int _{\mathbb{R}}q(x){e}^{iyx}dx,q(x)=\frac{1}{2\pi }\int _{\mathbb{R}}\hat{q}(y){e}^{-iyx}dy.\]
Then for all $u\in \mathbb{R}$ and $t\in [0,T]$ we receive the well-known relationship
\[ \phi _{M_{t}}(u)=\hat{f}_{M_{t}}(u)\]
where $\hat{f}_{M_{t}}$ denotes the Fourier transform of the density function $f_{M_{t}}$ defined in (2.6). An application of the inverse Fourier transform yields
\[ f_{M_{t}}(x)=\frac{1}{2\pi }\int _{\mathbb{R}}{e}^{\psi (u)t-iux}du\]
thanks to (2.7). On the other hand, from Eq. (1) in [6] we know that
(2.9)
\[ \phi _{M_{t}}(u)={e}^{iu\mu t}{\bigg(\frac{\cos (\beta /2)}{\cosh ((\alpha u-i\beta )/2)}\bigg)}^{2\delta t}\]
where $u\in \mathbb{R}$ and $t\in [0,T]$. Taking the logarithm in (2.7) and (2.9), we finally deduce
(2.10)
\[ \psi (u)=iu\mu +2\delta \bigg[\ln \cos \bigg(\frac{\beta }{2}\bigg)-\ln \cosh \bigg(\frac{\alpha u-i\beta }{2}\bigg)\bigg].\]
Further on, for the Meixner distribution the following properties are well-known (cf. [14], Section 5.3.10 in [13], Section 3.6 in [1], Corollary 1 in [6]).
Lemma 2.1.
  • (a) If $X\sim \mathcal{M}(\alpha ,\beta ,\delta ,\mu )$, then $cX+m\sim \mathcal{M}(c\alpha ,\beta ,\delta ,c\mu +m)$ with constants $c>0$ and $m\in \mathbb{R}$.
  • (b) If $X_{1}\sim \mathcal{M}(\alpha ,\beta ,\delta _{1},\mu _{1})$ and $X_{2}\sim \mathcal{M}(\alpha ,\beta ,\delta _{2},\mu _{2})$ are independent random variables, then $X_{1}+X_{2}\sim \mathcal{M}(\alpha ,\beta ,\delta _{1}+\delta _{2},\mu _{1}+\mu _{2})$.
  • (c) The characteristic function $\phi _{X}(u;\alpha ,\beta ,\delta ,\mu )$ of a Meixner distributed random variable $X\sim \mathcal{M}(\alpha ,\beta ,\delta ,\mu )$ satisfies
    \[ \phi _{X}(u;\alpha ,\beta ,\delta ,\mu )=\phi _{X}{(u;\alpha ,\beta ,\delta /n,\mu /n)}^{n}\]
    for arbitrary $n\in \mathbb{N}$ such that the Meixner distribution is infinitely divisible.

3 A stock price model driven by a Meixner process

Let $t\in [0,T]$ and define the stochastic stock price process $S_{t}$ via
(3.1)
\[ S_{t}:=S_{0}{e}^{M_{t}+bt}\]
with deterministic initial value $S_{0}$, a constant $b\in \mathbb{R}$ and a real-valued Meixner process
\[ M_{t}=\theta t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{Q}}(s,z)\]
such as introduced in (2.1)–(2.4). Here, the constant b provides some additional degree of freedom which is introduced in order to ensure the arbitrage-freeness of the stock price model. More details on this topic will be given below. Verify that (3.1) belongs to the same model class (geometric Lévy models) as (2.2)–(2.3) in [8]. We next introduce the historical filtration
\[ \mathcal{F}_{t}:=\sigma \{S_{u}:0\le u\le t\}=\sigma \{M_{u}:0\le u\le t\}.\]
Using Itô’s formula, we obtain the stochastic differential equation (SDE)
\[ \frac{dS_{t}}{S_{t-}}=\bigg(\theta +b+\int _{\mathbb{R}_{0}}\big[{e}^{z}-1-z\big]d{\nu }^{\mathbb{Q}}(z)\bigg)dt+\int _{\mathbb{R}_{0}}\big[{e}^{z}-1\big]d{\tilde{N}}^{\mathbb{Q}}(t,z)\]
under $\mathbb{Q}$. Let us further define the discounted stock price via
\[ \hat{S}_{t}:=\frac{S_{t}}{B_{t}}\]
where $S_{t}$ is such as defined in (3.1) and $B_{t}:={e}^{rt}$ is the value of a bank account with normalized initial capital $B_{0}=1$ and risk-less interest rate $r>0$. Due to (3.1) we find
\[ \hat{S}_{t}=S_{0}{e}^{M_{t}+(b-r)t}\]
while Itô’s formula yields the following SDE under $\mathbb{Q}$
\[ \frac{d\hat{S}_{t}}{\hat{S}_{t-}}=\bigg(\theta +b-r+\int _{\mathbb{R}_{0}}\big[{e}^{z}-1-z\big]d{\nu }^{\mathbb{Q}}(z)\bigg)dt+\int _{\mathbb{R}_{0}}\big[{e}^{z}-1\big]d{\tilde{N}}^{\mathbb{Q}}(t,z).\]
In accordance to no-arbitrage theory, the discounted stock price $\hat{S}$ must form a martingale under the risk-neutral probability measure $\mathbb{Q}$. For this reason, we require the drift restriction
(3.2)
\[ b=r-\theta -\int _{\mathbb{R}_{0}}\big[{e}^{z}-1-z\big]d{\nu }^{\mathbb{Q}}(z).\]
With this particular choice of the coefficient b, we deduce
\[ \frac{dS_{t}}{S_{t-}}=rdt+\int _{\mathbb{R}_{0}}\big[{e}^{z}-1\big]d{\tilde{N}}^{\mathbb{Q}}(t,z)\]
under $\mathbb{Q}$. Combining (3.1) and (3.2), we receive
\[ S_{t}=S_{0}{e}^{rt}\exp \Bigg\{{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{Q}}(s,z)-{\int _{0}^{t}}\int _{\mathbb{R}_{0}}\big[{e}^{z}-1-z\big]d{\nu }^{\mathbb{Q}}(z)ds\Bigg\}\]
where the last factor on the right hand side constitutes a Doléans-Dade exponential which again shows the $\mathbb{Q}$-martingale property of the discounted stock price process $\hat{S}_{t}=S_{t}{e}^{-rt}$. Moreover, taking (2.2), (2.4), (2.8) and (2.10) into account, Eq. (3.2) can be expressed as
(3.3)
\[ b=r-\psi (-i)=r-\mu -2\delta \ln \bigg(\frac{\cos (\beta /2)}{\cos ((\alpha +\beta )/2)}\bigg).\]
Unless otherwise stated, from now on we assume that the constant $b\in \mathbb{R}$ appearing in (3.1) is such as given in (3.3). Though constituting an admissible choice, taking $b=0$ in (3.1)–(3.3) might be too restrictive in practical applications. In the following, we investigate the log-returns related to our model (3.1). For an arbitrary time step $\varDelta >0$ and $t\le T-\varDelta $ we obtain
\[ \ln \bigg(\frac{S_{t+\varDelta }}{S_{t}}\bigg)\cong M_{\varDelta }+b\varDelta \sim \mathcal{M}\big(\alpha ,\beta ,\delta \varDelta ,(\mu +b)\varDelta \big)\]
by Lemma 2.1 (a). Here, the symbol ≅ denotes equality in distribution. Hence, in our stock price model (3.1) the log-returns are Meixner distributed. We stress that in [12] it was shown that the Meixner distribution fits empirical financial log-returns very well. Furthermore, for $n\in \mathbb{N}$ we introduce the time partition $\mathcal{P}:=\{0<t_{0}<t_{1}<\cdots <t_{n}\le T\}$ and define the return/revenue process associated with the period $[t_{k-1},t_{k}]$ via
\[ R_{k}:=\frac{S_{t_{k}}-S_{t_{k-1}}}{S_{t_{k-1}}}\]
where $k\in \{1,\dots ,n\}$. A substitution of (3.1) into the latter equation yields
(3.4)
\[ R_{k}={e}^{Y_{t_{k}}-Y_{t_{k-1}}}-1\]
where
\[ Y_{t}:=M_{t}+bt\sim \mathcal{M}\big(\alpha ,\beta ,\delta t,(\mu +b)t\big)\]
is a Meixner–Lévy process. Taking (2.1)–(2.4) and (3.3) into account, we get
(3.5)
\[ Y_{t}=\gamma t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zdN(s,z)\]
where
\[ \gamma :=r+\delta \bigg[\alpha \tan \bigg(\frac{\beta }{2}\bigg)-2\ln \bigg(\frac{\cos (\beta /2)}{\cos ((\alpha +\beta )/2)}\bigg)-\int _{\mathbb{R}_{0}}\frac{{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz\bigg]\]
is a real-valued constant. Recall that the Meixner–Lévy process Y given in (3.5) above just is a special case of the more general Lévy process X defined in Eq. (2.3) in [8]. For this reason, the cliquet option pricing results derived in [8] simultaneously apply in our current Meixner modeling case. More details on this topic are given in Section 4 below. Also note that $R_{1},\dots ,R_{n}$ are $\mathbb{Q}$-independent random variables and that $R_{k}>-1$ $\mathbb{Q}$-almost surely for all k. Since Y is a Lévy process under $\mathbb{Q}$, we observe $Y_{t_{k}}-Y_{t_{k-1}}\cong Y_{\tau }$ (stationary increments) where $\tau :=t_{k}-t_{k-1}$ (equidistant partition). Here, the symbol ≅ denotes equality in distribution. For the sake of notational simplicity, we always work under the assumption of equidistant time points in the following, unless otherwise stated. Taking (3.4) into account, we obtain the subsequent relationship between the cumulative distribution functions of $R_{k}$ and $Y_{\tau }$
(3.6)
\[ \mathbb{Q}(R_{k}\le \xi )=\mathbb{Q}\big(Y_{\tau }\le \ln (1+\xi )\big)\]
where $\xi >-1$ is an arbitrary real-valued constant.

3.1 A structure preserving measure change to the physical probability measure

Recall that we worked under the risk-neutral probability measure $\mathbb{Q}$ in the previous sections. Since log-returns of financial assets are commonly observed under the physical measure $\mathbb{P}$ (instead of under $\mathbb{Q}$), we establish a measure change from $\mathbb{Q}$ to $\mathbb{P}$ in the sequel. In this context, we have to pay special attention to the so-called structure preserving property of the measure change, as the log-returns under $\mathbb{P}$ shall again follow a Meixner distribution. In other words, the Meixner process $M_{t}$ introduced in (2.1) under $\mathbb{Q}$ shall also be a Meixner process under $\mathbb{P}$. First of all, for $t\in [0,T]$ we define the Radon–Nikodym density process
\[ \varLambda _{t}:=\frac{d\mathbb{P}}{d\mathbb{Q}}\bigg|_{\mathcal{F}_{t}}:=\exp \Bigg\{{\int _{0}^{t}}\int _{\mathbb{R}_{0}}h(z)d{\tilde{N}}^{\mathbb{Q}}(s,z)-{\int _{0}^{t}}\int _{\mathbb{R}_{0}}\big[{e}^{h(z)}-1-h(z)\big]d{\nu }^{\mathbb{Q}}(z)ds\Bigg\}\]
where the $\mathbb{Q}$-compensated PRM ${\tilde{N}}^{\mathbb{Q}}$ and the corresponding Lévy measure ${\nu }^{\mathbb{Q}}$ are such as defined in (2.3), respectively (2.4), while $h(z)$ is a time-independent deterministic function on $\mathbb{R}_{0}$. Recall that we may write
\[ \varLambda _{t}=\frac{{e}^{L_{t}}}{\mathbb{E}_{\mathbb{Q}}[{e}^{L_{t}}]}\]
with a local $\mathbb{Q}$-martingale process
\[ L_{t}:={\int _{0}^{t}}\int _{\mathbb{R}_{0}}h(z)d{\tilde{N}}^{\mathbb{Q}}(s,z)\]
such that the density process Λ is detected to be of Esscher transform type. Note that Λ is a discontinuous Doléans-Dade exponential which constitutes a local martingale under $\mathbb{Q}$ satisfying the SDE
\[ d\varLambda _{t}=\varLambda _{t-}\int _{\mathbb{R}_{0}}\big[{e}^{h(z)}-1\big]d{\tilde{N}}^{\mathbb{Q}}(t,z).\]
In accordance to Theorem 12.21 in [5], we further impose the Novikov condition
\[ \mathbb{E}_{\mathbb{Q}}\Bigg[\exp \Bigg\{{\int _{0}^{t}}\int _{\mathbb{R}_{0}}\big[1-{e}^{h(z)}+h(z){e}^{h(z)}\big]d{\nu }^{\mathbb{Q}}(z)ds\Bigg\}\Bigg]<\infty \]
for all $t\in [0,T]$. Then it holds $\mathbb{E}_{\mathbb{Q}}[\varLambda _{t}]\equiv 1$ for all $t\in [0,T]$ such that Λ constitutes a true $\mathbb{Q}$-martingale. Hence, we may apply Girsanov’s theorem stating that
(3.7)
\[ d{\tilde{N}}^{\mathbb{P}}(s,z):=dN(s,z)-d{\nu }^{\mathbb{P}}(z)ds\]
constitutes the $\mathbb{P}$-compensated Poisson random measure with Lévy measure
(3.8)
\[ d{\nu }^{\mathbb{P}}(z):={e}^{h(z)}d{\nu }^{\mathbb{Q}}(z).\]
Note that the Novikov condition is equivalent to requiring that
\[ \int _{\mathbb{R}_{0}}\big[1-{e}^{h(z)}+h(z){e}^{h(z)}\big]d{\nu }^{\mathbb{Q}}(z)<\infty \]
since h and ${\nu }^{\mathbb{Q}}$ both are deterministic. A combination of (2.1), (2.3), (3.7) and (3.8) yields the following Lévy–Itô decomposition
(3.9)
\[ M_{t}=\bigg(\theta +\int _{\mathbb{R}_{0}}z\big[{e}^{h(z)}-1\big]d{\nu }^{\mathbb{Q}}(z)\bigg)t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{P}}(s,z)\]
under $\mathbb{P}$ where θ, ${\nu }^{\mathbb{Q}}$ and ${\tilde{N}}^{\mathbb{P}}$ are such as defined in (2.2), (2.4) and (3.7), respectively. The remaining challenge now consists in finding an appropriate function $h(z)$ which, firstly, fulfills the Novikov condition, secondly, guarantees that ${\nu }^{\mathbb{P}}$ in (3.8) constitutes a Lévy measure of Meixner-type and, thirdly, ensures that $M_{t}$ in (3.9) is a Meixner–Lévy process. In this regard, we propose to work with the specification
(3.10)
\[ h(z):=\frac{{\beta }^{\ast }-\beta }{\alpha }z\]
from now on. Herein, the constant skewness parameter ${\beta }^{\ast }\in (-\pi ,\pi )$ satisfies ${\beta }^{\ast }\ne \beta $ while $\alpha >0$ is the scaling parameter introduced above. Note that taking ${\beta }^{\ast }=\beta $ would imply $h(z)\equiv 0$ and hence, $\mathbb{P}=\mathbb{Q}$. Combining (3.8) with (2.4) and (3.10), we deduce
(3.11)
\[ d{\nu }^{\mathbb{P}}(z)=\delta \frac{{e}^{{\beta }^{\ast }z/\alpha }}{z\sinh (\pi z/\alpha )}dz\]
which constitutes a Meixner-type Lévy measure with parameters α, ${\beta }^{\ast }$ and δ [recall Eq. (2.4)]. Moreover, with respect to (3.8) and (3.10), we obtain
\[ \int _{\mathbb{R}_{0}}\big[1-{e}^{h(z)}+h(z){e}^{h(z)}\big]d{\nu }^{\mathbb{Q}}(z)={\nu }^{\mathbb{Q}}(\mathbb{R}_{0})-{\nu }^{\mathbb{P}}(\mathbb{R}_{0})+\frac{{\beta }^{\ast }-\beta }{\alpha }\int _{\mathbb{R}_{0}}zd{\nu }^{\mathbb{P}}(z)\]
which is finite, because the Lévy measures ${\nu }^{\mathbb{Q}}$ and ${\nu }^{\mathbb{P}}$ are finite. Thus, the function $h(z)$ defined in (3.10) indeed fulfills the Novikov condition. Further on, we take (2.2), (2.4), (3.9) and (3.10) into account and receive
(3.12)
\[ \mathbb{E}_{\mathbb{P}}[M_{t}]={\theta }^{\ast }t\]
with drift parameter
\[ {\theta }^{\ast }:={\mu }^{\ast }+\delta \alpha \tan \bigg(\frac{{\beta }^{\ast }}{2}\bigg)\]
and a constant and real-valued location parameter
(3.13)
\[ {\mu }^{\ast }:=\mu +\delta \bigg[\alpha \tan \bigg(\frac{\beta }{2}\bigg)-\alpha \tan \bigg(\frac{{\beta }^{\ast }}{2}\bigg)+\int _{\mathbb{R}_{0}}\frac{{e}^{{\beta }^{\ast }z/\alpha }-{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz\bigg].\]
Note in passing that (3.12) possesses the same structure as (2.5). Also verify that
(3.14)
\[ {\theta }^{\ast }-\theta =\delta \int _{\mathbb{R}_{0}}\frac{{e}^{{\beta }^{\ast }z/\alpha }-{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz=\int _{\mathbb{R}_{0}}z\big[d{\nu }^{\mathbb{P}}(z)-d{\nu }^{\mathbb{Q}}(z)\big]\]
due to (2.2), (2.4), (3.11), (3.12) and (3.13). All in all, combining (3.9) with (2.4), (3.10) and (3.14), we conclude that
(3.15)
\[ M_{t}={\theta }^{\ast }t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{P}}(s,z)\]
which constitutes a Meixner–Lévy process under $\mathbb{P}$ with distribution
(3.16)
\[ M_{t}\sim \mathcal{M}\big(\alpha ,{\beta }^{\ast },\delta t,{\mu }^{\ast }t\big).\]
For this reason, we call the recently introduced measure change structure preserving. Recall that in Section 2 under $\mathbb{Q}$ we observed $M_{t}\sim \mathcal{M}(\alpha ,\beta ,\delta t,\mu t)$ on the other hand. Thus, the proposed measure change does neither affect the scaling parameter α nor the peakedness parameter δ whereas both the skewness parameter β and the location parameter μ are changed. Moreover, the Lévy triplet of the process M claimed in (3.15) is given by $({\theta }^{\ast },0,{\nu }^{\mathbb{P}})$. In analogy to the result provided in the sequel of (3.3), we remark that under $\mathbb{P}$ it holds
\[ \ln \bigg(\frac{S_{t+\varDelta }}{S_{t}}\bigg)\cong M_{\varDelta }+b\varDelta \sim \mathcal{M}\big(\alpha ,{\beta }^{\ast },\delta \varDelta ,\big({\mu }^{\ast }+b\big)\varDelta \big)\]
due to Lemma 2.1 (a). Here, $\varDelta >0$ is a constant and b is such as given in (3.3). Hence, if we specify the Radon–Nikodym function $h(z)$ like in (3.10), then the log-returns again are Meixner distributed under the real-world probability measure $\mathbb{P}$. Further note that
\[ \theta t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{Q}}(s,z)=M_{t}={\theta }^{\ast }t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{P}}(s,z)\]
holds $\mathbb{P}$- respectively $\mathbb{Q}$-almost surely for all $t\in [0,T]$ thanks to (2.1) and (3.15). We obtain the following expressions for the variance, skewness and kurtosis of $M_{t}$ under $\mathbb{P}$
\[\begin{array}{r@{\hskip0pt}l}\displaystyle \mathbb{Var}_{\mathbb{P}}[M_{t}]& \displaystyle =\frac{\delta t}{2}\frac{{\alpha }^{2}}{{\cos }^{2}({\beta }^{\ast }/2)},\hspace{2em}\mathbb{S}_{\mathbb{P}}[M_{t}]=\sqrt{2/(\delta t)}\sin \big({\beta }^{\ast }/2\big),\\{} \displaystyle \mathbb{K}_{\mathbb{P}}[M_{t}]& \displaystyle =3+\frac{2-\cos ({\beta }^{\ast })}{\delta t}\end{array}\]
while under $\mathbb{P}$ the density and characteristic function of $M_{t}$ are such as given in (2.6) and (2.9) but with μ and β therein replaced by ${\mu }^{\ast }$ and ${\beta }^{\ast }$, respectively.

3.1.1 A generalized structure preserving measure change

In this section, we present a generalized structure preserving measure change from the risk-neutral to the physical probability measure. Recall that the measure change proposed above only affects the skewness parameter β and the location parameter μ whereas both the scaling parameter α and the peakedness parameter δ remain untouched. From a practical point of view, this fact might be regarded as an advantage, as there is no need to recalibrate the parameters α and δ when changing from the risk-neutral to the physical probability measure. Conversely, the described feature might likewise cause some difficulties when it comes to calibrating under $\mathbb{P}$, since the values of the parameters α and δ have to be the same as under $\mathbb{Q}$ yielding some loss of flexibility. To avoid this disadvantage, we now propose a generalized measure change which affects each of the four parameters of the Meixner distribution. For this purpose, we presently require that the Meixner–Lévy measure under $\mathbb{P}$ is of the form
(3.17)
\[ d{\nu }^{\mathbb{P}}(z)={\delta }^{\ast }\frac{{e}^{{\beta }^{\ast }z/{\alpha }^{\ast }}}{z\sinh (\pi z/{\alpha }^{\ast })}dz\]
[cf. Equation (2.4)] with new parameters ${\alpha }^{\ast }>0$, ${\beta }^{\ast }\in (-\pi ,\pi )$ and ${\delta }^{\ast }>0$ which are different from α, β and δ introduced previously under $\mathbb{Q}$. Following this approach, we are led to the equality
\[ {e}^{h(z)}=\frac{{\delta }^{\ast }\sinh (\pi z/\alpha )}{\delta \sinh (\pi z/{\alpha }^{\ast })}\exp \bigg\{\bigg(\frac{{\beta }^{\ast }}{{\alpha }^{\ast }}-\frac{\beta }{\alpha }\bigg)z\bigg\}\]
due to (3.17), (3.8) and (2.4). Taking the logarithm in the latter equation, we receive
(3.18)
\[ h(z)=\bigg(\frac{{\beta }^{\ast }}{{\alpha }^{\ast }}-\frac{\beta }{\alpha }\bigg)z+\ln \bigg(\frac{{\delta }^{\ast }\sinh (\pi z/\alpha )}{\delta \sinh (\pi z/{\alpha }^{\ast })}\bigg)\]
which corresponds to (3.10) above. Note that (3.18) is well-defined for all $z\in \mathbb{R}_{0}$ and that we obtain (3.10), if we take ${\alpha }^{\ast }=\alpha $ and ${\delta }^{\ast }=\delta $ in (3.18). Hence, (3.10) is a special case of (3.18). Moreover, if we take ${\alpha }^{\ast }=\alpha $ and ${\beta }^{\ast }=\beta $ in (3.18), then we get $h(z)=\ln {\delta }^{\ast }-\ln \delta $ which is constant and independent of z. We summarize our findings in the following proposition.
Proposition 3.1.
Consider the measure change from the risk-neutral to the physical probability measure with Radon–Nikodym density process Λ such as defined at the beginning of Section 3.1. Then the new Lévy measure ${\nu }^{\mathbb{P}}$ under $\mathbb{P}$ is of Meixner-type again, if and only if the Radon–Nikodym function $h(z)$ is of the form (3.18).
In the sequel, we investigate the distributional properties of the corresponding Meixner–Lévy process under $\mathbb{P}$ related to the Radon–Nikodym function $h(z)$ given in (3.18). A substitution of (2.2), (2.4) and (3.18) into (3.9) yields the following Lévy–Itô decomposition under $\mathbb{P}$
(3.19)
\[ M_{t}=\overline{\theta }t+{\int _{0}^{t}}\int _{\mathbb{R}_{0}}zd{\tilde{N}}^{\mathbb{P}}(s,z)\]
with deterministic and real-valued drift parameter
\[ \overline{\theta }:=\mu +\delta \alpha \tan \bigg(\frac{\beta }{2}\bigg)+{\delta }^{\ast }\int _{\mathbb{R}_{0}}\frac{{e}^{{\beta }^{\ast }z/{\alpha }^{\ast }}}{\sinh (\pi z/{\alpha }^{\ast })}dz-\delta \int _{\mathbb{R}_{0}}\frac{{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz.\]
The Meixner–Lévy process $M_{t}$ given in (3.19) possesses the Lévy triplet $(\overline{\theta },0,{\nu }^{\mathbb{P}})$ where ${\nu }^{\mathbb{P}}$ is such as claimed in (3.17). In the next step, we require that $\overline{\theta }$ is of the form (2.2), i.e.
\[ \overline{\theta }=\overline{\mu }+{\delta }^{\ast }{\alpha }^{\ast }\tan \big({\beta }^{\ast }/2\big)\]
with some new location parameter $\overline{\mu }\in \mathbb{R}$. Following this onset, we deduce
\[ \overline{\mu }\hspace{0.1667em}=\hspace{0.1667em}\mu +\delta \alpha \tan \bigg(\frac{\beta }{2}\bigg)-{\delta }^{\ast }{\alpha }^{\ast }\tan \bigg(\frac{{\beta }^{\ast }}{2}\bigg)+{\delta }^{\ast }\int _{\mathbb{R}_{0}}\frac{{e}^{{\beta }^{\ast }z/{\alpha }^{\ast }}}{\sinh (\pi z/{\alpha }^{\ast })}dz-\delta \int _{\mathbb{R}_{0}}\frac{{e}^{\beta z/\alpha }}{\sinh (\pi z/\alpha )}dz.\]
Hence, if the measure change from $\mathbb{Q}$ to $\mathbb{P}$ is performed with the Radon–Nikodym function $h(z)$ defined in (3.18), then the corresponding Meixner–Lévy process $M_{t}$ again is Meixner distributed under $\mathbb{P}$ with parameters
(3.20)
\[ M_{t}\sim \mathcal{M}\big({\alpha }^{\ast },{\beta }^{\ast },{\delta }^{\ast }t,\overline{\mu }t\big).\]
If we compare (3.20) with (3.16), we see that in the generalized measure change related to (3.18) each of the four parameters of the Meixner–Lévy process M is affected.

4 Cliquet option pricing in a geometric Meixner model

This section is devoted to the pricing of cliquet options in the Meixner stock price model presented in Chapter 3. Since the Meixner process Y in (3.5) above just is a special case of the more general Lévy process X defined in Eq. (2.3) in [8], the cliquet option pricing results derived in [8] simultaneously apply to our present Meixner–Lévy modeling case. The details are worked out in the remainder of the current section. Parallel to [8] and Eq. (1.1) in [3], we consider a monthly sum cap style cliquet option with payoff
\[ H_{T}=K+K\max \Bigg\{g,{\sum \limits_{k=1}^{n}}\min \{c,R_{k}\}\Bigg\}\]
where T is the maturity time, K denotes the notional (i.e. the initial investment), g is the guaranteed rate at maturity, $c\ge 0$ is the local cap and $R_{k}$ is the return process given in (3.4). Recall that the payoff $H_{T}$ is globally-floored by the constant $K(1+g)$ and locally-capped by c. By a case distinction, we get
\[ H_{T}=K\max \Bigg\{1+g,1+{\sum \limits_{k=1}^{n}}\min \{c,R_{k}\}\Bigg\}=K\Bigg(1+g+\max \Bigg\{0,{\sum \limits_{k=1}^{n}}Z_{k}\Bigg\}\Bigg)\]
where for all $k\in \{1,\dots ,n\}$ the appearing objects
(4.1)
\[ Z_{k}:=\min \{c,R_{k}\}-g/n\]
are independent and identically distributed (i.i.d.) random variables. Note that $R_{k}$ is $\mathcal{F}_{t_{k}}$-measurable such that $H_{T}$ is $\mathcal{F}_{t_{n}}$-measurable. Since $t_{n}\le T$, it holds $\mathcal{F}_{t_{n}}\subseteq \mathcal{F}_{T}$ such that $H_{T}$ constitutes an $\mathcal{F}_{T}$-measurable claim. As before, let us denote the constant interest rate by $r>0$. Then the price at time $t\le T$ of a cliquet option with payoff $H_{T}$ at maturity T is given by the discounted risk-neutral conditional expectation of the payoff, i.e.
\[ C_{t}={e}^{-r(T-t)}\mathbb{E}_{\mathbb{Q}}(H_{T}\mid \mathcal{F}_{t}).\]
Combining the latter equations, we obtain
(4.2)
\[ C_{0}=K{e}^{-rT}\Bigg(1+g+\mathbb{E}_{\mathbb{Q}}\Bigg[\max \Bigg\{0,{\sum \limits_{k=1}^{n}}Z_{k}\Bigg\}\Bigg]\Bigg)\]
which shows that the considered cliquet option with payoff $H_{T}$ essentially is a plain-vanilla call option with strike zero written on the basket-style underlying ${\sum _{k=1}^{n}}Z_{k}$.
Proposition 4.1 (Cliquet option price).
Let $k\in \{1,\dots ,n\}$ and consider the independent and identically distributed random variables $Z_{k}=\min \{c,R_{k}\}-g/n$ where $c\ge 0$ is the local cap, $R_{k}$ is the return process given in (3.4) and g is the guaranteed rate at maturity. Denote the maturity time by T, the notional by K and the risk-less interest rate by r. Then the price at time zero of a cliquet option with payoff $H_{T}$ can be represented as
(4.3)
\[ C_{0}=K{e}^{-rT}\Bigg(1+g+\frac{n}{2}\mathbb{E}_{\mathbb{Q}}[Z_{1}]+\frac{1}{\pi }{\int _{{0}^{+}}^{\infty }}\frac{1-Re(\phi _{Z}(x))}{{x}^{2}}dx\Bigg)\]
where $Re$ denotes the real part and the characteristic function $\phi _{Z}(x)$ is defined via
(4.4)
\[ \phi _{Z}(x):={\prod \limits_{k=1}^{n}}\phi _{Z_{k}}(x)={\prod \limits_{k=1}^{n}}\mathbb{E}_{\mathbb{Q}}\big[{e}^{ixZ_{k}}\big]={\big(\phi _{Z_{1}}(x)\big)}^{n}={\big(\mathbb{E}_{\mathbb{Q}}\big[{e}^{ixZ_{1}}\big]\big)}^{n}.\]
Proof.
See the proof of Prop. 3.1 in [3], respectively of Prop. 3.1 in [8].  □
In the subsequent sections, we derive explicit expressions for $\phi _{Z}(x)$ and $\mathbb{E}_{\mathbb{Q}}[Z_{1}]$ appearing in the pricing formula (4.3). As before, we stick to the presumption of equidistant resetting times and set $\tau =t_{k}-t_{k-1}$ for all $k\in \{1,\dots ,n\}$ in the following.

4.1 Cliquet option pricing with distribution functions

Let us first apply a method involving probability distribution functions (cf. [3] and Section 3.1 in [8]). We initially investigate the treatment of $\phi _{Z}(x)$ defined in (4.4).
Proposition 4.2.
Let $Y_{\tau }\sim \mathcal{M}(\alpha ,\beta ,\delta \tau ,(\mu +b)\tau )$ and suppose that $Z_{k}=\min \{c,R_{k}\}-g/n$ where $k\in \{1,\dots ,n\}$. Then the characteristic function of $Z_{k}$ under $\mathbb{Q}$ can be represented as
(4.5)
\[\begin{array}{r@{\hskip0pt}l}\displaystyle \phi _{Z_{k}}(x)& \displaystyle ={e}^{-ix(1+g/n)}\Bigg({e}^{ix(1+c)}{\int _{1+c}^{\infty }}f_{Y_{\tau }}(u)du\\{} & \displaystyle \hspace{1em}+{\int _{-\infty }^{0}}f_{Y_{\tau }}(u)du+{\int _{0}^{1+c}}{e}^{ixu}f_{Y_{\tau }}(u)du\Bigg)\end{array}\]
where
(4.6)
\[ f_{Y_{\tau }}(u)=\frac{{(2\cos (\beta /2))}^{2\delta \tau }}{2\pi \alpha \varGamma (2\delta \tau )}{e}^{\beta (u-(\mu +b)\tau )/\alpha }{\bigg|\varGamma \bigg(\delta \tau +i\frac{u-(\mu +b)\tau }{\alpha }\bigg)\bigg|}^{2}\]
constitutes the probability density function of the Meixner–Lévy process Y given in (3.5) and b is the real-valued constant claimed in (3.3).
Proof.
By similar arguments as in the proof of Prop. 3.2 in [8], we obtain
(4.7)
\[ \phi _{Z_{k}}(x)={e}^{-ix(1+g/n)}\Bigg({e}^{ix(1+c)}-ix{\int _{0}^{1+c}}{e}^{ixw}\mathbb{Q}(R_{k}\le w-1)dw\Bigg).\]
Using (3.6) and the definition of the distribution function, we get for the last integral in (4.7)
\[ {\int _{0}^{1+c}}{e}^{ixw}\mathbb{Q}(R_{k}\le w-1)dw={\int _{0}^{1+c}}{\int _{-\infty }^{\ln (w)}}{e}^{ixw}f_{Y_{\tau }}(u)dudw\]
where $Y_{\tau }\sim \mathcal{M}(\alpha ,\beta ,\delta \tau ,(\mu +b)\tau )$ is the Meixner–Lévy process given in (3.5) and $f_{Y_{\tau }}(u)$ constitutes the probability density function of $Y_{\tau }$ under $\mathbb{Q}$ claimed in (4.6). Applying Fubini’s theorem and hereafter splitting up the resulting outer integral, we deduce
\[\begin{array}{r@{\hskip0pt}l}& \displaystyle {\int _{0}^{1+c}}{e}^{ixw}\mathbb{Q}(R_{k}\le w-1)dw\\{} & \displaystyle \hspace{1em}={\int _{-\infty }^{0}}{\int _{0}^{1+c}}{e}^{ixw}f_{Y_{\tau }}(u)dwdu+{\int _{0}^{1+c}}{\int _{u}^{1+c}}{e}^{ixw}f_{Y_{\tau }}(u)dwdu.\end{array}\]
We next compute the emerging $dw$-integrals and finally substitute the resulting expression into (4.7) which yields (4.5).  □
If we insert (4.5) into (4.4), we receive a representation for the characteristic function $\phi _{Z}(x)$. Let us proceed with the computation of $\mathbb{E}_{\mathbb{Q}}[Z_{k}]$.
Proposition 4.3.
Suppose that $Z_{k}=\min \{c,R_{k}\}-g/n$ where $k\in \{1,\dots ,n\}$. Then the first moment of $Z_{k}$ under $\mathbb{Q}$ is given by
(4.8)
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=-1-\frac{g}{n}+(1+c){\int _{1+c}^{\infty }}f_{Y_{\tau }}(u)du+{\int _{0}^{1+c}}uf_{Y_{\tau }}(u)du\]
where $f_{Y_{\tau }}(u)$ is the probability density function of $Y_{\tau }$ under $\mathbb{Q}$ given in (4.6).
Proof.
In accordance to Prop. 2.4 in [4], we have
(4.9)
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=\frac{1}{i}\frac{\partial }{\partial x}\big(\phi _{Z_{k}}(x)\big)\big|_{x=0}.\]
A substitution of (4.5) into (4.9) instantly yields (4.8).  □
As mentioned in Section 2, we recall that the cumulative distribution function (cdf) of the Meixner-Lévy process $Y_{t}$, i.e.
\[ F_{Y_{t}}(x):={\int _{-\infty }^{x}}f_{Y_{t}}(u)du\]
does not possess a closed form representation, but it can be computed efficiently with numerical methods. Also note that all integrals appearing in (4.5) and (4.8) are finite, since $f_{Y_{\tau }}(\boldsymbol{\cdot })$ constitutes a probability density function while the Meixner–Lévy process $Y_{\tau }$ possesses moments of all orders (cf. [13]).

4.2 Cliquet option pricing with Fourier transform techniques

There is an alternative method to derive expressions for $\mathbb{E}_{\mathbb{Q}}[Z_{k}]$, $\phi _{Z}(x)$ and $C_{0}$ involving Fourier transforms and the Lévy–Khinchin formula. In the following, we present this method which has firstly been proposed in [8] in a cliquet option pricing context.
Proposition 4.4.
Let $Y_{t}\sim \mathcal{M}(\alpha ,\beta ,\delta t,(\mu +b)t)$ be the Meixner–Lévy process considered in (3.5). Suppose that $Z_{k}=\min \{c,R_{k}\}-g/n$ where $k\in \{1,\dots ,n\}$ and let $\vartheta >0$ be a finite real-valued dampening parameter. Then the first moment of $Z_{k}$ under $\mathbb{Q}$ can be represented as
(4.10)
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=c-\frac{g}{n}-\frac{1}{2\pi }\int _{\mathbb{R}}\frac{{(c+1)}^{1+\vartheta +iy}}{(\vartheta +iy)(1+\vartheta +iy)}\phi _{Y_{\tau }}(i\vartheta -y)dy\]
where the characteristic function $\phi _{Y_{\tau }}$ is given by
(4.11)
\[ \phi _{Y_{\tau }}(i\vartheta -y)={e}^{-(\vartheta +iy)(\mu +b)\tau }{\bigg(\frac{\cos (\beta /2)}{\cosh ((i(\alpha \vartheta -\beta )-\alpha y)/2)}\bigg)}^{2\delta \tau }.\]
Proof.
The proof follows the same lines as the proof of Prop. 3.4 in [8]. From (4.1) and the equality
\[ \min \{c,R_{k}\}=c-{[c-R_{k}]}^{+}\]
we deduce
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=c-g/n-\mathbb{E}_{\mathbb{Q}}\big[{(c-R_{k})}^{+}\big].\]
Taking (3.4) into account, we next receive
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=c-g/n-\mathbb{E}_{\mathbb{Q}}\big[{\big(c+1-{e}^{Y_{\tau }}\big)}^{+}\big]\]
where $\tau =t_{k}-t_{k-1}$ and Y is the real-valued Meixner–Lévy process given in (3.5). With a finite and real-valued dampening parameter $\vartheta >0$ we define the function
\[ \varphi (u):={e}^{\vartheta u}{\big(c+1-{e}^{u}\big)}^{+}.\]
Since $\varphi \in {\mathcal{L}}^{1}(\mathbb{R})$, its Fourier transform exists and reads as
\[ \hat{\varphi }(y)=\frac{{(c+1)}^{1+\vartheta +iy}}{(\vartheta +iy)(1+\vartheta +iy)}.\]
Using the inverse Fourier transform along with Fubini’s theorem, we get
\[ \mathbb{E}_{\mathbb{Q}}\big[{\big(c+1-{e}^{Y_{\tau }}\big)}^{+}\big]=\mathbb{E}_{\mathbb{Q}}\big[{e}^{-\vartheta Y_{\tau }}\varphi (Y_{\tau })\big]=\frac{1}{2\pi }\int _{\mathbb{R}}\hat{\varphi }(y)\mathbb{E}_{\mathbb{Q}}\big[{e}^{-(\vartheta +iy)Y_{\tau }}\big]dy\]
which implies (4.10). The expression for the characteristic function $\phi _{Y_{\tau }}$ given in (4.11) can directly be obtained by virtue of (2.9).  □
Our argumentation in the proof of Proposition 4.4 motivates the following considerations.
Proposition 4.5.
Let $Y_{t}\sim \mathcal{M}(\alpha ,\beta ,\delta t,(\mu +b)t)$ be the Meixner–Lévy process presented in (3.5). Suppose that $Z_{k}=\min \{c,R_{k}\}-g/n$ with $k\in \{1,\dots ,n\}$ and $c\ge 0$. Then the characteristic function of $Z_{k}$ under $\mathbb{Q}$ reads as
(4.12)
\[ \phi _{Z_{k}}(x)={e}^{-ixg/n}\Bigg({e}^{ixc}+{\int _{-\infty }^{\ln (1+c)}}\big[{e}^{ix({e}^{u}-1)}-{e}^{ixc}\big]f_{Y_{\tau }}(u)du\Bigg)\]
where the probability density function $f_{Y_{\tau }}$ of $Y_{\tau }$ under $\mathbb{Q}$ is such as given in (4.6).
Proof.
Similar computations as in the proof of Prop. 3.5 in [8] yield (4.12).  □
There is an alternative method involving (4.9) to derive an expression for $\mathbb{E}_{\mathbb{Q}}[Z_{k}]$ which is presented in the following.
Corollary 4.6.
In the setup of Proposition 4.5, we receive the representation
(4.13)
\[ \mathbb{E}_{\mathbb{Q}}[Z_{k}]=c-\frac{g}{n}+{\int _{-\infty }^{\ln (1+c)}}\big[{e}^{u}-1-c\big]f_{Y_{\tau }}(u)du.\]
Proof.
The claimed representation immediately follows from Eq. (3.16) in [8].  □
Inspired by the Fourier transform techniques applied in the proof of Proposition 4.4, we now focus on the derivation of an alternative representation for the cliquet option price $C_{0}$ given in (4.2). The corresponding result reads as follows.
Theorem 4.7 (Fourier transform cliquet option price).
Let $k\in \{1,\dots ,n\}$ and consider the independent and identically distributed random variables $Z_{k}\hspace{0.1667em}=\hspace{0.1667em}\min \{c,R_{k}\}-g/n$ where $c\ge 0$ is the local cap, g is the guaranteed rate at maturity and $R_{k}$ is the return process defined in (3.4). For $n\in \mathbb{N}$ we set $\varrho :=nc-g$ and denote the maturity time by T, the notional by K and the riskless interest rate by r. Let $Y_{t}\sim \mathcal{M}(\alpha ,\beta ,\delta t,(\mu +b)t)$ be the Meixner–Lévy process given in (3.5). Then the price at time zero of a cliquet option paying
\[ H_{T}=K\Bigg(1+g+\max \Bigg\{0,{\sum \limits_{k=1}^{n}}Z_{k}\Bigg\}\Bigg)\]
at maturity can be represented as
(4.14)
\[\begin{array}{r@{\hskip0pt}l}\displaystyle C_{0}& \displaystyle =K{e}^{-rT}\Bigg[1+g+{\int _{{0}^{+}}^{\infty }}\frac{1+iy\varrho -{e}^{iy\varrho }}{2\pi {y}^{2}}\\{} & \displaystyle \hspace{1em}\times {\Bigg(1+{\int _{-\infty }^{\ln (1+c)}}\big[{e}^{iy({e}^{u}-1-c)}-1\big]f_{Y_{\tau }}(u)du\Bigg)}^{n}dy\Bigg]\end{array}\]
where $f_{Y_{\tau }}(u)$ constitutes the probability density function claimed in (4.6).
Proof.
The proof of Theorem 3.7 in [8] here applies equally, if we replace $f_{X_{\tau }}$ therein by $f_{Y_{\tau }}$.  □

5 Conclusions

In this paper, we investigated the pricing of a monthly sum cap style cliquet option with underlying stock price modeled by a geometric pure-jump Meixner–Lévy process. In Section 2, we compiled various facts on the Meixner distribution and the related class of stochastic Meixner–Lévy processes. In Section 3, we introduced a stock price model driven by a Meixner–Lévy process and established a customized structure preserving measure change from the risk-neutral to the physical probability measure. Moreover, we obtained semi-analytic expressions for the cliquet option price by using the probability distribution function of the driving Meixner–Lévy process in Section 4.1 and by an application of Fourier transform techniques in Section 4.2. To read more on cliquet option pricing in a jump-diffusion Lévy model, the reader is referred to the accompanying article [8].

References

[1] 
Albrecher, H., Ladoucette, S., Schoutens, W.: A generic one-factor Lévy model for pricing synthetic CDOs. In: Advances in Mathematical Finance, pp. 259–277. Birkhäuser, (2007). MR2359372. https://doi.org/10.1007/978-0-8176-4545-8_14
[2] 
Bernard, C., Boyle, P., Gornall, W.: Locally-capped Contracts and the Retail Investor. Journal of Derivatives 18(4), 72–88 (2011). https://doi.org/10.3905/jod.2011.18.4.072
[3] 
Bernard, C., Li, W.: Pricing and Hedging of Cliquet Options and Locally-capped Contracts. SIAM Journal on Financial Mathematics 4, 353–371 (2013). MR3038023. https://doi.org/10.1137/100818157
[4] 
Cont, R., Tankov, P.: Financial Modeling with Jump Processes, 1st edn. Chapman & Hall/CRC, (2004). MR2042661
[5] 
Di Nunno, G., Øksendal, B., Proske, F.: Malliavin Calculus for Lévy Processes with Applications to Finance, 1st edn. Springer, (2009). MR2460554
[6] 
Grigelionis, B.: Processes of Meixner Type. Lithuanian Mathematical Journal 39(1), 33–41 (1999). MR1711971
[7] 
Haifeng, Y., Jianqi, Y., Limin, L.: Pricing Cliquet Options in Jump-diffusion Models. Stochastic Models 21, 875–884 (2005). MR2179304. https://doi.org/10.1080/15326340500294587
[8] 
Hess, M.: Cliquet Option Pricing in a Jump-Diffusion Lévy Model (2017). SSRN working paper. https://ssrn.com/abstract=2979296
[9] 
Iseger P., den Oldenkamp, E.: Cliquet options: Pricing and Greeks in deterministic and stochastic volatility models (2005). SSRN working paper. https://ssrn.com/abstract=1013510
[10] 
Protter, P.: Stochastic Integration and Differential Equations, 2nd edn. Springer, (2005). MR2273672
[11] 
Sato, K.: Lévy Processes and Infinitely Divisible Distributions. Cambridge studies in advanced mathematics, vol. 68 (1999). MR1739520
[12] 
Schoutens, W.: The Meixner Process: Theory and Applications in Finance, Eurandom Report 2001-004, Eindhoven, pp. 1–24 (2002)
[13] 
Schoutens, W.: Lévy Processes in Finance: Pricing Financial Derivatives. John Wiley & Sons, Ltd., (2003)
[14] 
Schoutens, W., Teugels, J.: Lévy Processes, Polynomials and Martingales. Stochastic Models 14(1–2), 335–349 (1998). MR1617536
[15] 
Wilmott, P.: Cliquet Options and Volatility Models, Wilmott magazine, December (2002)
Exit Reading PDF XML


Table of contents
  • 1 Introduction
  • 2 A review of Meixner processes
  • 3 A stock price model driven by a Meixner process
  • 4 Cliquet option pricing in a geometric Meixner model
  • 5 Conclusions
  • References

Export citation

Copy and paste formatted citation
Placeholder

Download citation in file


Share


RSS

MSTA

MSTA

  • Online ISSN: 2351-6054
  • Print ISSN: 2351-6046
  • Copyright © 2018 VTeX

About

  • About journal
  • Indexed in
  • Editors-in-Chief

For contributors

  • Submit
  • OA Policy
  • Become a Peer-reviewer

Contact us

  • ejournals-vmsta@vtex.lt
  • Mokslininkų 2A
  • LT-08412 Vilnius
  • Lithuania
Powered by PubliMill  •  Privacy policy