Abstract
We propose a model for an insurance loss index and the claims process of a single insurance company holding a fraction of the total number of contracts that captures both ordinary losses and losses due to catastrophes. In this model we price a catastrophe derivative by the method of utility indifference pricing. The associated stochastic optimization problem is treated by techniques for piecewise deterministic Markov processes. A numerical study illustrates our results.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
Costly natural catastrophes in the recent past (hurricane Andrew in 1992, hurricane Katrina in 2005, the earthquake and tsunami in Japan 2011 resulting in the nuclear disaster at Fukushima, floods in Thailand 2011) all caused severe stress to the (re-)insurance industry. However, these losses are still small relative to losses of the US stock and bond markets. Therefore securitization (i.e. transferring part of the risk to the financial market) is an efficient alternative to reinsuring catastrophe (CAT) losses, cf. [5].
Contracts of this kind are insurance-linked derivatives.Footnote 1 They are usually written on insurance industry catastrophe loss indices, insurer-specific catastrophe losses, or parametric indices based on the physical characteristics of catastrophe events. We focus on the first kind of products; they involve more basis risk, but are less exposed to moral hazard than the others, cf. [4].
Derivatives written on insurance industry catastrophe loss indices were first issued in 1992 by the Chicago Board of Trade; these were futures and later also call- and put spread options written on aggregate CAT-loss indices, cf. [4].
A call spread option is a combination of a call option long and a call option short with a higher strike. Another popular type of catastrophe derivative is the CAT bond. This is a classical bond combined with an option that is triggered by a (predefined) catastrophe event. Note that the buyer of the bond thereby sells the embedded option. The issuer is typically a (re-)insurance company that wants to reinsure parts of its risk exposure on the financial market. In return the investor receives a coupon.
CAT derivatives are interesting for investors who seek to diversify their risk, since they are largely uncorrelated with classical financial instruments.
The challenges in pricing CAT derivatives are that the underlying index is not a traded asset, that they are not liquidly traded themselves and, maybe most of all, the modeling of catastrophe events.
In the following we review the existing literature. For a more detailed literature overview we refer to [17].
Geman and Yor [12] study European vanilla call options written on an insurance loss index, which is modeled by a jump-diffusion. Cox et al. [3] model the aggregate loss of an insurance company by a Poisson process with constant arrival rate of catastrophe events and derive a pricing formula for CAT-puts. Jaimungal and Wang [14] model the aggregate loss by a compound Poisson process to describe the dynamic losses more accurately. Muermann [18] derives the market price of insurance risk from CAT derivative prices in a compound Poisson model. Leobacher and Ngare [15] use the method of utility indifference pricing to price CAT derivatives written on an insurance loss index modeled by a compound Poisson process.
For catastrophe events, the assumption that the resulting claims occur at jump times of a Poisson process as adopted by most previous studies is not beyond justifiable critique. A generalization was proposed in [9], who model an insurance loss index by a doubly stochastic Poisson process (Cox process), i.e. the arrival rate of claims is a stochastic process itself; they price CAT futures in this model. Lin et al. [16] also model the arrival of CAT events by a doubly stochastic Poisson process. See also [11] for no-arbitrage pricing of CAT bonds in this context. Dassios and Jang [6] study the valuation of CAT derivatives by risk neutral valuation, where the underlying is modeled as a Cox process with shot noise intensity.
In this paper, we introduce a novel model for an insurance loss index and for a single insurance portfolio that captures ordinary insurance losses as well as catastrophe losses. We model the ordinary claims in the loss index by a compound Poisson process with constant intensity and we model the arrival of catastrophes by a Poisson process with constant intensity, where a jump triggers another stochastic variable that determines the number of claims in case of a catastrophe.
The claims process of a single insurance company holding a fraction of the total number of contracts is then a dynamic thinning of the process describing the index. Our model has the advantage that the jump height distribution does not need to capture both many small claims and outliers caused by catastrophes, but these outliers are split into many smaller claims.
The dynamic thinning is reached in a very convenient way (by drawing from a uniform distribution on [0, 1]) that we believe to be applicable in many other situations.
Using this model we present a pricing mechanism for CAT derivatives (like CAT spread options). Since the insurance loss index is not a tradable asset, and since the market for CAT derivatives is not liquid, risk neutral valuation is not applicable. Instead we use the method of utility indifference pricing. For this we need a hedging mechanism, which will be an active management of the risk portfolio. The pricing method requires solving an associated stochastic optimization problem.
Our paper extends [15] by a more realistic modeling approach for the insurance loss index and also for the thinning. In our paper a catastrophe event may partly hit the considered insurance company, whereas in [15] a catastrophe event always only affects one company. Their model for the claims process of a single insurance company is a thinning (a change of the intensity) of the Poisson process driving the number of claims, while ours is a dynamic thinning of the claims for each event and thus has a different distribution of the jumps.
The model presented here is technically harder to handle; we provide the mathematical toolkit in this paper. Using this new model instead of a simpler one is justified by our numerical results, which show that the new model has a significant impact on the price of a CAT derivative as it reflects catastrophes more accurately.
We also introduce a way to compute the utility indifference price of the derivative by Fourier techniques. This method also allows to compute the residual risk and the profit-loss distribution and therefore to evaluate coherent risk measures.
The paper is organized as follows. In Sect. 2 we model the insurance loss index and the claims process of a single insurance company. The state process based on which the CAT derivative is priced, is identified as a piecewise deterministic Markov process (PDMP), see [2, 7]. In Sect. 3 we recall the general concept of utility indifference pricing and we solve the associated stochastic optimization problem. In Sect. 4 we show how the utility indifference price and also quantities relevant for risk management can be computed efficiently, and we present a numerical study.
2 The model
Let \((\Omega ,\mathcal{F},\mathbb{P})\) be a probability space carrying all stochastic variables appearing below.
Suppose we have a global claims process C, which keeps track of all property insurance claims in a given country and we consider an insurance company in the same country, so that the index will contain the losses of that particular insurance company among others.
The portfolio income rate consisting of the premium revenues from the risk portfolio is given by a continuous function q of the company’s market share \(\xi \in [0,1]\). The function q is not necessarily linear in \(\xi\), since demand for insurance might depend on the premium the company charges. The wealth process of the insurance company can be controlled by managing the insured portfolio, i.e. by controlling the market share \(\xi\). This allows for optimizing the management strategy for maximizing utility from terminal wealth.
Therefore, we can apply the method of utility indifference pricing for the valuation of CAT-derivatives.
The global claims process is given by
and where \(N^1=(N^1_t)_{t\ge 0}\) and \(N^2=(N^2_t)_{t\ge 0}\) are independent Poisson processes with intensities \(\lambda ^1,\lambda ^2\) and jump times \((\tau ^1_{i})_{i\ge 1},(\tau ^2_{i})_{i\ge 1}\). The jump heights \((Y_{i,j})_{i,j\ge 1}\) are iid random variables representing the damage of, e.g., single houses. The random variables \((\tilde{A}_i)_{i\ge 1}\), \(\tilde{A}_i \in {\mathbb {N}}\backslash \{1\}\), describe the number of claims in case of a catastrophe. The process \(N^1\) describes the occurrence of regular claims whereas a jump of \(N^2\) indicates an accumulation of \(\tilde{A}_1\) claims due to a catastrophe event.
If the insurance company holds the \(\xi\)-th part of the whole risk, it is exposed to the \(\xi\)-th part of the claims. We model this as
The random variables \((U_{i,j})_{i,j\ge 1}\) are iid and \(U_{1,1}\sim \mathcal{U}([0,1])\); they determine whether the company is affected by the corresponding claim or not. For fixed \(\xi\) this is a thinning of the original process, cf. [19, Section 3.12.1].
We assume independence of \(N^1,N^2,(Y_{i,j})_{i,j\ge 1},(\tilde{A}_i)_{i\ge 1},(U_{i,j})_{i,j\ge 1}\).
It is possible to write (1) as a single sum by adapting the jump intensity and the distribution of the \(\tilde{A}_i\), which we will do to ease the notation in the following. Note that the jump height distribution does not need to be adapted so that we do not loose the favourable properties for modeling catastrophe events. Let \(L=(L_t)_{t\ge 0}\) be a Poisson process with intensity \(\lambda =\lambda ^1+\lambda ^2\) and jump times \((\tau _{i})_{i\ge 1}\) and let the number of claims per jump of L be denoted by \((A_i)_{i\ge 1}\) with
We can write the insurance loss index as
Both \(\lambda\) and the distribution of \(A_1\) are chosen such that (1) and (4) are equivalent.
The claims process of the insurance company holding the \(\xi\)-th part of the risk becomes
Denote by \(a_k:=\mathbb{P} (A_1=k)\). For all \(s\in [0,1]\) the generating function of \(A_1\) is given by \(G_{A_1}(s):=\sum _{k=1}^\infty a_{k} s^k\), where
Assumption 2.1
We assume that
-
\({\mathbb {E}}(e^{\eta Y_{1,1}})<\infty\);
-
\(\limsup _{k\rightarrow \infty }a_{k+1}/a_k<1/{\mathbb {E}}(e^{\eta Y_{1,1}})\).
Assumption 2.1 implies that the convergence radius of the generating function \(G_{A_1}\) is greater than \({\mathbb {E}}(e^{\eta Y_{1,1}})\) and hence
Note that if \({\mathbb {E}}\left( e^{\eta Y_{1,1}A_1}\right) <\infty\), also \({\mathbb {E}}\left( {\mathbb {E}}(e^{\eta Y_{1,1}})^{A_1}\right) <\infty\) by Jensen’s inequality.
In contrast to a model where the claims process is a simple compound Poisson process, here assuming the existence of exponential moments of the claim size distribution is not a great restriction, since we model catastrophes as an accumulation of small claims rather than one big claim.
The dynamics of the wealth process \(X^\xi =(X^\xi _t)_{t\ge 0}\) of the insurance company with initial wealth x is given by:
2.1 PDMP characterization
The two-dimensional process \((C,X^\xi )\) is a PDMP in the sense of [7]. We also refer to [2, Chapter 8] or [1]. Our PDMP has the following characteristics:
-
state space \({\mathbb {R}}_0^+ \times {\mathbb {R}}\);
-
control space [0, 1];
-
deterministic flow \(d(C_t,X^\xi _t)=(0, q(\xi _t)) dt\) between jumps;
-
jump intensity \(\lambda\);
-
jump kernel Q,
$$\begin{aligned} Q(B|(c,x),\xi )= \sum _{k=0}^\infty a_k Q_k(B|(c,x),\xi ), \end{aligned}$$where
$$\begin{aligned} Q_k(B|(c,x),\xi )= \sum _{\mathcal{K}\subseteq \{1,\dots ,k\}} \xi ^{|\mathcal{K}|} \left( 1-\xi \right) ^{k-|\mathcal{K}|} \mathbb{P} \left( \left( \sum _{j=1}^k Y_{1,j},\sum _{j\in \mathcal{K}} Y_{1,j} \right) \in B-(c,x) \right) , \end{aligned}$$and where we use the notation \(B-(c,x)=\{(b_1-c,b_2-x):(b_1,b_2)\in B\}\);
-
zero running reward rate;
-
zero discount rate.
Denoting by \(\tau\) the time of a jump of the PDMP and by \((C_\tau ,X_\tau )\) the state immediately after that jump, we define the set of bounded Markov controls \(\mathcal{M}_b\) as the set of all measurable functions assigning to given input data \((\tau ,C_\tau ,X_\tau )\) a control until the next jump, i.e.
3 Utility indifference pricing
The method of utility indifference pricing for the valuation of derivatives in incomplete markets has been introduced in [13]. It relies on the fact that even if the derivative cannot be replicated, it may still be the case that much of its variation can be hedged.
In [8] utility indifference pricing is used to price structured catastrophe bonds. However, there is a difference in modeling the hedging possibility. In our setup this is done via managing the insured portfolio. The main idea is that the loss in the portfolio of a single insurance company is necessarily correlated with the insurance loss index. The introduction of the derivative has therefore an influence on the pricing policy of the insurance company.
We will first explain the notion of utility indifference pricing and then apply it to our problem.
Assume the investor has a utility function u and initial wealth x. Define \(J(x,\ell ):=\sup _{X_T}{\mathbb {E}}(u(X_T+ \ell \psi ))\), where the supremum is taken over all possible wealths \(X_T\) that can be generated from x. The random variable \(\psi\) is the payment from a European claim with expiry T, and \(\ell\) is the number of claims that are bought.
The utility indifference bid price \(p^b(\ell )\) is the price at which the investor has the same utility whether she pays nothing and does not receive the claim \(\psi\), or she pays \(p^b(\ell )\) now and receives \(\ell\) units of the claim \(\psi\) at time T. Therefore, \(p^b(\ell )\) is the largest amount of money the investor is willing to pay for buying \(\ell\) units of the claim \(\psi\); it solves \(J(x-p^b(\ell ),\ell )=J(x,0)\).
The utility indifference ask price \(p^a(\ell )\) is the smallest amount of money the investor is willing to accept for selling \(\ell\) units of the claim \(\psi\); it solves \(J(x+p^a(\ell ),-\ell )=J(x,0)\).
The two prices are related via \(p^b(\ell )=-p^a(-\ell )\). With this in mind we can define the utility indifference price \(p:=p^b(1)\).
Assumption 3.1
-
The insurance company has exponential utility \(u(x)=-\exp (-\eta x)\), \(\eta >0\).
-
\(X_T\) is of the form \(x+\Gamma ^\xi _T\) for some control \(\xi\) and \(\Gamma ^\xi _T\) does not depend on the initial wealth x.
In that case
(provided that the arguments in the logarithms are finite), and hence p does not depend on the initial wealth x.
Note that exponential utility is a natural choice for insurance companies as often such a utility function is used to calculate insurance premia. As an example where exponential utility is used in a stochastic optimal control framework in an insurance context, see [10].
3.1 The stochastic optimization problem
We apply the concept of utility indifference pricing to the model presented in Sect. 2. Our aim is to price a derivative written on the total claims process C with payoff \(\psi (C_T)\), where \(\psi\) is a continuous and bounded function on \({\mathbb {R}}_0^+\).
Example 3.1
We are specifically interested in CAT (spread) options, i.e.
with cap L and strike \(0<K<L\). The option is in the money, if c exceeds K, and the payoff is bounded by \(L-K\).
Note that the main task in pricing CAT bonds also lies in pricing the embedded spread option, since for exponential utility the price of a CAT bond is the sum of a spread option price and a bond price.
We maximize the expected utility from terminal wealth. The corresponding value function is defined by
Since \(\psi\) is bounded we have that for \(\xi \equiv 0\), \({\mathbb {E}}(u(X^\xi _T+\psi {(C_T)})|C_t=c,X^\xi _t=x)>-\infty\) for all t, c, x, and hence \(V(t,c,x)>-\infty\). V is bounded from above since u is bounded. Therefore, V is well-defined.
For \(v:[0,T] \times {\mathbb {R}}_0^+ \times {\mathbb {R}}\longrightarrow {\mathbb {R}}\) bounded and measurable the generator of the jump process is defined by
The Hamilton–Jacobi–Bellman (HJB) equation corresponding to optimization problem (8) is
We make the ansatz \(v(t,c,x)=u(x)\exp (-\eta w(t,c))\) to obtain a backward equation which is independent of the initial wealth x. This yields
Defining
and using that v is negative, we obtain the backward equation for w:
Lemma 3.2
Let W be such that \(V(t,c,x)=u(x)\exp (-\eta W(t,c,x))\). Then W is bounded by \(\Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty\).
Proof
We have \(V(t,c,x)=u(x)\exp (-\eta W(t,c,x))\), i.e.
Denote by \(X^0\) the process \(X^\xi\) with \(\xi \equiv 0\). Then
and
Thus \(|W(t,c,x)|\le \Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty\). \(\square\)
3.2 Verification result
We show that the solution of the HJB equation (9) solves the optimization problem (8). For this we apply results from stochastic control theory for PDMPs; more precisely, a slight variation of the verification theorem [2, Theorem 8.2.8]. For this we recall two definitions from [1].
Definition 3.3
A measurable function \(b:{\mathbb {R}}_0^+ \times {\mathbb {R}}\longrightarrow {\mathbb {R}}_0^+\) is called a bounding function for our piecewise deterministic Markov decision model, if there exist constants \(c_u,c_Q,c_{\text {flow}}\ge 0\) such that for all \((c,x)\in {\mathbb {R}}_0^+ \times {\mathbb {R}}\)
-
i.
\(|u(x+\psi (c))|\le c_u b(c,x)\);
-
ii.
\(\int b(\tilde{c},\tilde{x})Q(d\tilde{c}\times d\tilde{x}|(c,x),\xi )\le c_Q b(c,x)\) for all \((c,x)\in {\mathbb {R}}_0^+ \times {\mathbb {R}}\), \(\xi \in [0,1]\);
-
iii.
\(b(c,x+\int _0^T \int _0^1 q(\xi )r_s(d\xi )ds)\le c_{\text {flow}}b(c,x)\) for all \(r\in \mathcal{R}\).
Here \(\mathcal{R}\) is the space of relaxed policies, i.e. of measurable maps \({\mathbb {R}}_0^+\longrightarrow \mathcal{P}([0,1])\), where \(\mathcal{P}([0,1])\) is the space of all probability measures on the Borel \(\sigma\)-algebra on [0, 1].
Definition 3.4
Let \(b:{\mathbb {R}}_0^+\times {\mathbb {R}}\longrightarrow {\mathbb {R}}_0^+\) be a bounding function and \(\gamma >0\) fixed. Define the Banach space \({\mathbb {B}}_{b,\gamma }:=\{v:[0,T] \times {\mathbb {R}}_0^+ \times {\mathbb {R}}\longrightarrow {\mathbb {R}}_0^+: v \;\text { measurable and }\;\Vert v\Vert _b<\infty \},\) with the norm
Theorem 3.5
Let b be a bounding function for our piecewise deterministic Markov decision model with \({\mathbb {E}}(|b(C_T,X^\xi _T)|\big |C_t=c,X_t=x)<\infty\) for all \(\xi ,t,c,x\). Let \(v\in \mathcal{C}^{1,0,1}([0,T] \times {\mathbb {R}}_0^+ \times {\mathbb {R}})\cap {\mathbb {B}}_{b,\gamma }\) be a solution of the HJB Eq. (9) and let \(\alpha ^*\) be a maximizer for (9), leading to the state process \((C,X^{\xi ^*})\).
Then \(v=V\) and \(\xi ^*=\alpha ^*(t,C_{t-},X^{\xi ^*}_{t-})\) is an optimal feedback-type Markov policy.
Remark 3.6
In the statement of [2, Theorem 8.2.8] there is another condition required, namely that \(\alpha _b<1\) for a constant \(\alpha _b\) depending on b, Q and the arbitrary \(\gamma\) from Definition 3.4. But it is argued in [1] that for finite horizon problems \(\gamma\) can always be chosen large enough to satisfy \(\alpha _b<1\).
For proving Theorem 3.5, we first need to prove existence of a bounding function.
Lemma 3.7
The function b defined by \(b(c,x):=\exp (\eta |x|)\) is a bounding function for our piecewise deterministic Markov decision model.
Proof
We need to check the conditions given in Definition 3.3.
-
i.
\(u(x+\psi (c))=-\exp (-\eta (x+\psi (c))\) such that \(|u(x+\psi (c))|=\exp (-\eta (x+\psi (c))\le \exp (\eta \Vert \psi \Vert _\infty )b(c,x)\).
-
ii.
\(\int b(\tilde{c},\tilde{x})Q(d\tilde{c}\times d\tilde{x}|(c,x),\xi ) =\int \exp (\eta |\tilde{x}|)Q(d\tilde{c}\times d\tilde{x}|(c,x),\xi ) ={\mathbb {E}}(\exp (\eta |x-Z_1^{\xi }|)) \le b(c,x){\mathbb {E}}(\exp (\eta Z_1^{\xi }))\).
-
iii.
\(b\left( c,x+\int _0^T\int _0^1 q(\xi )r_s(d\xi )ds\right) =\exp (\eta |x+\int _0^T\int _0^1 q(\xi )r_s(d\xi )ds|) \le \exp (\eta \Vert q\Vert _{\infty } T)b(c,x)\).
\(\square\)
Now we need to show that the backward equation (11) has a solution and hence also (9) has a solution.
Define \(\mathcal{H}\) on \(\mathcal{C}_b({\mathbb {R}}_0^+)\) by
We show that if \(\varphi \in \mathcal{C}_b({\mathbb {R}}_0^+)\), then \(\mathcal{H}\varphi \in \mathcal{C}_b({\mathbb {R}}_0^+)\) and that \(\mathcal{H}\) is locally Lipschitz. For this we write \(\mathcal{H}=g \circ h \circ f\) and show that g, h, f are locally Lipschitz and g is \(\mathcal{C}_b({\mathbb {R}}_0^+)\)-valued.
Lemma 3.8
For \(\sigma \in \mathcal{C}_b(R_0^+)\) the mapping \(\xi \mapsto {\mathbb {E}}(\sigma (Z_1)\exp (\eta Z_1^{\xi }) )+\lambda /\eta\) is a power series in \(\xi\). Its coefficients are of the form \(h_k(\sigma )={\mathbb {E}}(\delta _k \sigma (Z_1))\) for non-negative random variables \(\delta _k\) with \(\sum _{k=0}^\infty {\mathbb {E}}(\delta _k)<\infty\) that do not dependent on \(\xi\) and \(\sigma\). The power series converges uniformly on [0, 1].
Proof
Let F be the distribution function of \(Y_{1,1}\). Then
Expanding the above expression yields the first and the second claim of the lemma. Setting \(\sigma \equiv 1\), we get \(\sum _{k=0}^\infty {\mathbb {E}}(\delta _k)<\infty\). Setting \(\xi =1\) gives
The right-hand side is finite by Assumption 2.1 and by the assumption of the lemma. \(\square\)
Define the function-space
and let \(h:\mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}_0^+)\longrightarrow \Lambda\) be defined by \(h_k(\phi )(c)={\mathbb {E}}(\delta _k \phi (c,Z_1))\), with \(\delta _k\) as in Lemma 3.8. Hence, for any \(\phi \in \mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}_0^+)\), \({\mathbb {E}}\left( \phi (c,Z_1)\exp (\eta Z_1^\xi )\right) =\sum _{k=0}^\infty h_k(\phi )(c) \xi ^k\) for every \(c\in {\mathbb {R}}_0^+\).
Lemma 3.9
The function h is a bounded linear operator.
Proof
We need to prove that for every \(k\in {\mathbb {N}}_0\) the mapping \(h_k\) is a bounded linear operator \(\mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}_0^+)\longrightarrow \mathcal{C}_b({\mathbb {R}}_0^+)\).
Let \(\phi \in \mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}_0^+)\). We show that the mapping \(c\mapsto {\mathbb {E}}(\delta _k \phi (c,Z_1))\) is continuous and bounded on \({\mathbb {R}}_0^+\). Let \(c_n \rightarrow c\) in \({\mathbb {R}}_0^+\). Then \(\phi (c_n,z) \rightarrow \phi (c,z)\) for all \(z\in {\mathbb {R}}_0^+\). The sequence \((\phi (c_n,\cdot )\delta _k)_{n\ge 0}\) is dominated by \(\Vert \phi \Vert _\infty \delta _k\), which is integrable. Hence, \({\mathbb {E}}(\phi (c_n,Z_1)\delta _k)\rightarrow {\mathbb {E}}(\phi (c,Z_1)\delta _k)\) by the dominated convergence theorem. Thus \(h_k(\phi )\) is continuous. Moreover, \(h_k(\phi )\) is bounded, since \(\left| {\mathbb {E}}(\phi (c,Z_1)\delta _k)\right| \le \Vert \phi \Vert _\infty {\mathbb {E}}(\delta _k)\).
For \(\phi \in \mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}_0^+)\) it holds that \(\Vert h\Vert _\Lambda \le \sum _{k=0}^\infty \Vert h_k(\phi )\Vert _\infty =\sum _{k=0}^\infty \sup _c|{\mathbb {E}}(\delta _k\phi (c,Z_1))| \le \Vert \phi \Vert _\infty \sum _{k=0}^\infty {\mathbb {E}}(\delta _k)\). Thus h is bounded by Lemma 3.8. \(\square\)
For a sequence \((x_k)_{k\ge 0}\) in \({\mathbb {R}}\) define the function \(\tilde{g}\) by
Lemma 3.10
Let \(\tilde{g}\) be defined as in (13). Then
-
1.
\(\tilde{g}\) is defined on \(\ell ^1\) and it is bounded on every norm-bounded subset of \(\ell ^1\);
-
2.
\(\tilde{g}\) is convex;
-
3.
\(\tilde{g}\) is Lipschitz on every norm-bounded subset of \(\ell ^1\).
Proof
For \(x=(x_k)_{k\ge 0}\in \ell ^1\) and \(\xi \in [0,1]\) we have \(|\sum _{k=0}^\infty x_k \xi ^k|\le \Vert x\Vert _1\). The function \(\xi \mapsto \sum _{k=0}^\infty x_k \xi ^k\), \(\xi \in [0,1]\) is well-defined and continuous as a uniform limit of continuous functions on [0, 1]. Since q is also continuous, the first statement follows.
The proof of the second statement is straightforward.
Following [20] we use the convexity of \(\tilde{g}\) to show that \(\tilde{g}\) is Lipschitz on \(\{x\in \ell ^1:\Vert x\Vert \le R\}\). Let \(\Vert x\Vert ,\Vert y\Vert \le R\) and define \(z:=y+\frac{R}{\Vert y-x\Vert }(y-x)\). It holds that \(\Vert z-y\Vert =R\) and hence \(\Vert z\Vert \le 2R\). By the definition of z we have that \(y=\beta z + (1-\beta )x\), where \(\beta =\Vert y-x\Vert /(\Vert y-x\Vert +R)\). Since \(\tilde{g}\) is convex, \(\tilde{g}(y)\le \beta \tilde{g}(z)+(1-\beta )\tilde{g}(x)\) and hence \(\Vert \tilde{g}(y)-\tilde{g}(x)\Vert = \Vert \beta (\tilde{g}(z)-\tilde{g}(x))\Vert \le 2\beta \sup _{\Vert z\Vert \le 2R}|\tilde{g}(z)| \le 2\beta c \le (2C/R) \Vert y-x\Vert\) for some constant \(c>0\), since \(\tilde{g}\) is bounded on \(\{z\in \ell ^1:\Vert z\Vert \le 2R\}\). \(\square\)
For \(\phi \in \Lambda\) let \(\phi (c):=(\phi _k(c))_{k\ge 0}\) and define the function g by \(g(\phi )(c)=\tilde{g}(\phi (c))\).
Lemma 3.11
The function g is \(\mathcal{C}_b({\mathbb {R}}_0^+)\) -valued and locally Lipschitz.
Proof
Let \(\phi =(\phi _k)_{k\ge 0}\in \Lambda\). Let \(c_n\rightarrow c\) in \({\mathbb {R}}_0^+\) and let \(\varepsilon >0\). There exists \(k_0\in {\mathbb {N}}_0\) such that \(\sum _{k\ge k_0} \Vert \phi _k\Vert _\infty <\varepsilon /4\), and for n large enough \(\sum _{k=0}^{k_0} |\phi _k(c_n)-\phi _k(c)|<\varepsilon /2\). Thus, \(\sum _{k=0}^\infty |\phi _k(c_n)-\phi _k(c)| \le \sum _{k=0}^{k_0} |\phi _k(c_n)-\phi _k(c)|+\sum _{k=k_0+1}^{\infty } |\phi _k(c_n)-\phi _k(c)|< \varepsilon\).
The claim that g is locally Lipschitz follows from Lemma 3.10: let \(R>0\) and let \(\phi ^1,\phi ^2\in \Lambda\) with \(\Vert \phi ^1\Vert _\Lambda \le R\) and \(\Vert \phi ^2\Vert _\Lambda \le R\). \(\tilde{g}\) is Lipschitz on the ball with radius R in \(\ell ^1\). Denote the corresponding Lipschitz constant by \(L_R\). Then \(\phi ^1(c),\phi ^2(c)\) lie in the ball with radius R in \(\ell ^1\). Hence, \(\Vert g(\phi ^1)-g(\phi ^2)\Vert _\infty =\sup _c|\tilde{g}(\phi ^1(c))-\tilde{g}(\phi ^2(c))|\le L_R \Vert \phi ^1(c)-\phi ^2(c)\Vert _1 \le L_R\Vert \phi ^1-\phi ^1\Vert _\Lambda\). \(\square\)
Finally, define \(f:\mathcal{C}_b({\mathbb {R}}_0^+) \longrightarrow \mathcal{C}_b({\mathbb {R}}_0^+ \times {\mathbb {R}}^+)\), \(f(w)(c,z):=\exp (-\eta (w(c+z)-w(c))\) and note that f is locally Lipschitz.
Lemma 3.12
Let \(\mathcal{H}\) be defined as in (12). If \(\varphi \in \mathcal{C}_b({\mathbb {R}}_0^+)\) , then \(\mathcal{H}\varphi \in \mathcal{C}_b({\mathbb {R}}_0^+)\) and \(\mathcal{H}\) is locally Lipschitz.
Proof
We have \(\mathcal{H}=g \circ h \circ f\). The first claim follows from Lemma 3.11. Further, \(\mathcal{H}\) is locally Lipschitz as a concatenation of locally Lipschitz functions. The latter follows from Lemmas 3.9 and 3.11. \(\square\)
Now we prove that (11) has a unique maximal local solution.
Lemma 3.13
Let \(\psi \in \mathcal{C}_b({\mathbb {R}})\). Then the backward equation (11) has a unique maximal local solution.
Proof
The backward equation (11) is an initial value problem with \(\mathcal{C}_b({\mathbb {R}}_0^+)\)-valued solution:
By Lemma 3.12, \(\mathcal{H}\) is \(\mathcal{C}_b({\mathbb {R}}_0^+)\)-valued and locally Lipschitz. In particular, \(\mathcal{H}\) is Lipschitz on the ball with radius \(2(\Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty )\). From the Picard-Lindelöf theorem on existence and uniqueness of solutions of ordinary differential equations we get existence and uniqueness of a maximal local solution of (14), i.e. there exists \(\varepsilon >0\) and a solution \(\varphi\) of (14) on \([T-\varepsilon ,T]\) with \(\Vert \varphi (t)\Vert _\infty \le 2 (\Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty )\) for all \(t\in [T-\varepsilon ,T]\). We may choose \(\varepsilon\) maximal such that \(\varepsilon =T\) or \(\Vert \varphi (T-\varepsilon )\Vert _\infty =2 (\Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty )\).
The function \(w:[T-\varepsilon , T]\times {\mathbb {R}}_0^+\longrightarrow {\mathbb {R}}\) defined by \(w(t,c)=\varphi (t)(c)\) is the unique maximal local solution of (11). \(\square\)
Proof (Proof of Theorem 3.5)
Since for every \(\varphi \in \mathcal{C}_b({\mathbb {R}}_0^+)\) and \(c\in {\mathbb {R}}_0^+\) the function \(\xi \mapsto q(\xi ) +(\tilde{\mathcal{A}}^\xi w)(c)\) is continuous on [0, 1] by Lemma 3.8, there exists a maximizer for (9). By Lemmas 3.7 and 3.13 the assumptions of Theorem 3.5 are satisfied on \([T-\varepsilon ,T]\times {\mathbb {R}}_0^+ \times {\mathbb {R}}\longrightarrow {\mathbb {R}}\), where \(\varepsilon\) is as in the proof of Lemma 3.13. Along the lines of the proof of [2, Theorem 8.2.8] it can be shown that \(v:[T-\varepsilon ,T]\times {\mathbb {R}}_0^+ \times {\mathbb {R}}\longrightarrow {\mathbb {R}}\) with \(v(t,c,x)=u(x)\exp (-\eta w(t,c))\) solves the optimization problem (8) for \(t\in [T-\varepsilon ,T]\), i.e. \(v(t,\cdot ,\cdot )=V(t,\cdot ,\cdot )\) for \(t\in [T-\varepsilon ,T]\).
By Lemma 3.2 it holds that \(|w(t,c)|=|-1/\eta \log (v(t,c,x)/u(x))| =|-1/\eta \log (V(t,c,x)/u(x))|=|W(t,c,x)|\le \Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty\) for \(t\in [T-\varepsilon ,T]\). Therefore, \(\Vert w(T-\varepsilon ,.)\Vert _\infty <2( \Vert q\Vert _{\infty } T+\Vert \psi \Vert _\infty )\) and hence \(\varepsilon =T\). Thus w solves (11) on the whole of \([0,T]\times {\mathbb {R}}_0^+\) and therefore v solves the optimization problem (8) in the whole of \([0,T] \times {\mathbb {R}}_0^+ \times {\mathbb {R}}\). \(\square\)
3.3 Utility indifference price
With the solution w of (11) we can compute the utility indifference price p of a derivative with payoff \(\psi\). Given the value of the index c and the amount of wealth x at time t, the maximum expected utility of terminal wealth can be written as
The corresponding value with no derivative bought is given by
where \(w^0\) is the solution of Eq. (11) with \(\psi \equiv 0\). Therefore (7) simplifies to
In particular, p does not depend on x and we omit that parameter from p henceforth. The function \(w^0\) does not depend on c. So \(w^0\) is the solution of an ordinary differential equation.
Lemma 3.14
Let \(w^0\) solve (11) with terminal condition \(w^0(T,c)\equiv 0\). Then
Proof
Since \(w^0\) does not depend on c, the backward equation (11) becomes
We calculate the expected value in (16):
Integration yields the claimed solution. \(\square\)
From (15) we can derive a backward equation for p. Since \(w^0\) does not depend on the second variable, we have
and it holds that \(p_t(t,c)=w_t(t,c)-w^0_t(t,c)=w_t(t,c)-\bar{w}\), where
Hence
4 Computations
In this section, we present a convenient numerical method for computing the expected value in (17) for the case where the distribution of \(Y_{1,1}\) has a smooth density. Denote by \(f^{*k}\) the k-fold convolution of a function f with itself, i.e. \(f^{*2}=f*f\) and \(f^{*(k+1)}=f*f^{*k}\). Let \(\hat{F}\) denote the Fourier transform. It holds that \(\hat{F}(f^{*k})=\hat{F}(f)^k\). The following lemma gives an efficient method for computing \({\mathbb {E}}(\exp (-\eta w(t,c+Z_1))\exp (\eta Z_1^{\xi }))\).
Lemma 4.1
Assume that the distribution of \(Y_{1,1}\) has a piecewise continuous density \(\mu\). Denote by \(\tilde{\mu }(z)=\exp (\eta z) \mu (z)\). Let \(\sigma\) be measurable and bounded.
Then it holds that
where \(\mu _\xi =\hat{F}^{-1}\left( G_{A_1}(\hat{F}(\xi \tilde{\mu }+(1-\xi )\mu ))\right)\).
Proof
Denote \(\bar{Y}_{k_1,k_2}:=\sum _{j=k_1}^{k_2} Y_{1,j}\). We have
where we used that \(\tilde{\mu }^{*k}(z)=\exp (\eta z) \mu ^{*k}(z)\), which can be seen by induction. The first claim now follows by linearity of the convolution. Further,
With this,
Since \((\xi \tilde{\mu }+(1-\xi )\mu )\) is piecewise continuous, \((\xi \tilde{\mu }+(1-\xi )\mu ))^{*k}\) is continuous for \(k\ge 2\) and hence \(\sum _{k=1} a_k\left( \xi \tilde{\mu }+(1-\xi )\mu \right) ^{*k}\) is piecewise continuous.
Thus \(\sum _{k=1}^\infty a_k\left( \xi \tilde{\mu }+(1-\xi )\mu \right) ^{*k}=\hat{F}^{-1}\left( \hat{F}\left( \sum _{k=1}^\infty a_k\left( \xi \tilde{\mu }+(1-\xi )\mu \right) ^{*k}\right) \right)\) a.e. on \({\mathbb {R}}_0^+\). \(\square\)
4.1 Numerical experiments
Our aim is to price a CAT (spread) option, i.e. \(\psi (C_T)=\max (0,\min (C_T-K,L-K))\).
The function q governing the company’s market share is chosen as in [15]. We assume that there are M clients in the market who potentially contribute to the claims process. Let a be the fair annual premium for one client, i.e. \({\mathbb {E}}(C_1)=M a\). The annual premium for one contract therefore has to be greater or equal than a, since otherwise the insurance company will make an almost sure loss in the long run. The premium the insurance company charges for a claim is \(a(1+\theta )\) with \(\theta >0\). Furthermore, the company faces an exogenously given demand curve d for insurance. It is continuous, decreasing in \(\theta\), and satisfies \(d(\theta )=M\) for \(\theta \le 0\) and \(d(\theta )=0\) for \(\theta \ge m\), i.e. the company gets to insure the whole risk, if it does not charge any risk loading (\(\theta =0\)) and it gets 0 contracts, if the risk-loading exceeds some fixed number \(m>0\). With this \(q(\xi )=\xi a (1+\theta (\xi ))\), where \(\theta (\xi )=d^{-1}(\xi M)\). For our numerical example we choose
The model parameters in our example are given by \(M=10^4\), \(m=2\), \(K=10^7\), \(\eta =10^{-6}\), \(T=1\) year. The number of jumps in case of a catastrophe is Poisson distributed, \((\tilde{A}_1-2) \sim \text {Poisson}(40)\); the distribution of \(A_1\) then follows from (3). The claim size distribution is a Gamma distribution, \(Y_{1,1}\sim \text {Gamma}(10,5000)\).
We want to study two effects: the effect that holding a derivative has on the risk loading (which depends on the optimal market share \(\xi ^*\)) and its change over time, and the effect of our model [the clustered claims (CC) model] on the utility indifference price of the derivative and on the risk loading in comparison to the model where the claims process is a simple compound Poisson process as, e.g., in [15] [the single claim (SC) model].
For the SC-model we choose \(\lambda _1=100, \lambda _2=0\).
In order to be able to compare the two models we adapt \(\lambda _1,\lambda _2\) such that the expected annual claim size per contract a stays constant, yielding \(\lambda _1=69, \lambda _2=1\).
Figure 1 shows the utility indifference price p of the CAT spread option in dependence of the value c of the claims process. We see that the price increases in c. As time increases the expected number of claims within the remaining time decreases, and hence also the price decreases; for \(t\rightarrow T\) the prices converges to the payoff. Further, we observe that in the SC-model the price is always lower than in the CC-model, since the latter more accurately accounts for a clustering of claims.
Figure 2 shows the risk loading corresponding to the optimal market share \(\xi ^*\) in dependence of the value c of the claims process. For small c the risk loading is the same as in the case of no derivative held, since the probability that the derivative has a positive payoff is small. For \(c>L\) a further increase of c does not change the payoff and hence the situation is the same as for holding no derivative.
As time increases the probability that the payoff of the derivative grows in c and hence compensates losses during the remaining time decreases, but also the expected number of claims before T decreases. An interesting observation is that for small c the first effect dominates and hence the risk loading increases, whereas for large c the latter effect dominates and hence the risk loading decreases.
In the CC-model the risk loading is in general higher than in the SC-model, since we imposed risk aversion. The risk loading decreases significantly when a derivative is bought.
The effect of holding a derivative is higher in the CC-model; the optimal average risk loading decreases by approximately \(31.9\%\) (compared to \(2.7\%\) in the SC-model).
We observe the same effects when comparing our CC-model to the SC-model for the example of [15], where a bounded claim size distribution was used.
4.2 Concluding remarks
The introduction of a derivative serves as an effective alternative to classical reinsurance and leads to significantly smaller insurance premia. The CC-model introduced in this paper has a significant impact on the price of a CAT derivative as well as on the optimal average risk loading. It reflects catastrophes more accurately.
4.2.1 Risk management
In this section we compute the profit-loss distribution and the residual risk of an insurance company holding a CAT derivative. The former is useful for the derivation of coherent risk measures, the latter quantifies the efficiency of the hedge.
4.2.2 Profit-loss distribution
The profit-loss distribution is the distribution of the (optimally controlled) wealth \(\rho :=X_T^{\xi ^*}+\psi (C_T)-p\) in case the company holds a CAT derivative.
Let \(\xi ^*\) be the optimal control. For \(\varsigma \in \mathbb {C}\) define
The right-hand side in (18) can be interpreted as (two-sided) Laplace transform \(\hat{L}\) of the profit-loss distribution at \(\varsigma\), \(\hat{L}(\varsigma )={\mathbb {E}}(\exp (-\varsigma \rho )|C_t=c,X^{\xi ^*}_t=x)\). We can compute \(V_{\varsigma }\) by making the ansatz \(V_{\varsigma }(t,x,c)=u(x)e^{-\varsigma W_{\varsigma }(t,c)}\), where \(W_\varsigma\) solves the backward equation
with corresponding \(\tilde{\mathcal{A}}^{\xi ^*}_\varsigma\). By solving (19) for different values of \(\varsigma\), we get the density \(\nu\) of the profit-loss distribution by inverting the Laplace transform:
4.2.3 Residual risk
The numerical experiments in Sect. 4.1 showed that the optimal market share \(\xi ^*\) of an insurance company that holds a CAT derivative is higher than for a company that does not (\(\xi ^0\)). The change in the strategy \(\xi ^*-\xi ^0\) when an insurance company buys a CAT derivative, is also the strategy used for hedging the derivative itself.
The (buyer’s) risk of derivative is \(\psi (C_T)-p\); the residual risk, i.e. the remaining risk after hedging is given by
The density of (20) can be computed in the same way as the density of the profit-loss distribution above.
Notes
Details on currently listed insurance-linked derivatives can be found at http://www.artemis.bm/deal_directory.
References
Bäuerle N, Rieder U (2010) Optimal control of piecewise deterministic markov processes with finite time horizon. In: Piunovskiy A (ed) Modern trends in controlled stochastic processes: theory and applications. Luniver Press, Frome, pp 123–143
Bäuerle N, Rieder U (2011) Markov decision processes with applications in finance. Universitext. Springer, Berlin
Cox SH, Fairchild JR, Pedersen HW (2004) Valuation of structured risk management products. Insur Math Econ 34(2):259–272
JD Cummins (2006) Should the government provide insurance for catastrophes? Fed Reserve Bank St. Louis Rev 88(4):337–380
Cummins JD, Lalonde D, Phillips RD (2004) The basis risk of catastrophic-loss index securities. J Financ Econ 71(1):77–111
Dassios A, Jang J-W (2003) Pricing of catastrophe reinsurance and derivatives using the cox process with shot noise intensity. Financ Stoch 7:73–95
Davis MHA (1993) Markov models and optimization. Monographs on statistics and applied probability. Chapman & Hall, London
Egami M, Young VR (2008) Indifference prices of structured catastrophe (CAT) bonds. Insur Math Econ 42(2):771–778
Embrechts P, Meister S (1997) Pricing insurance derivatives: the case of CAT futures. In: Securization of insurance risk: the 1995 Bowles symposium, pp 16–26
Fernández B, Hernández-Hernández D, Meda A, Saavedra P (2008) An optimal investment strategy with maximal risk aversion and its ruin probability. Math Methods Oper Res 68:159–179
Fuita T, Ishimura N, Tanake D (2008) An arbitrage approach to the pricing of catastrophe options involving the cox process. Hitotsubashi J Econ 49:67–74
Geman H, Yor M (1997) Stochastic time changes in catastrophe option pricing. Insur Math Econ 21(3):185–193
Hodges S, Neuberger A (1989) Optimal replication of contingent claim under transaction costs. Rev Futures Mark 8:222–239
Jaimungal S, Wang T (2006) Catastrophe options with stochastic interest rates and compound poisson losses. Insur Math Econ 38(3):469–483
Leobacher G, Ngare P (2016) Utility indifference pricing of derivatives written on industrial loss indexes. J Comput Appl Math 300:68–82
Lin S, Chang C, Powers MR (2009) The valuation of contingent capital with catastrophe risks. Insur Math Econ 45(1):65–73
Muermann A (2004) Catastrophe derivatives. In: Teugels JL, Sundt B (eds) Encyclopedia of actuarial science, vol 1. Wiley, Chichester, pp 231–236
Muermann A (2008) Market price of risk implied by catastrophe derivatives. N Am Actuar J 12(3):221–227
Olofsson P, Andersson M (2012) Probability, statistics, and stochastic processes, 2nd edn. Wiley, Hoboken
Roberts AW, Varberg DE (1974) Another proof that convex functions are locally Lipschitz. Am Math Mon 81:1014–1016
Acknowledgements
Open access funding provided by Austrian Science Fund (FWF). A. Eichler is supported by the Austrian Science Fund (FWF): Project P21196. G. Leobacher is supported by the Austrian Science Fund (FWF): Project F5508-N26, which is part of the Special Research Program “Quasi-Monte Carlo Methods: Theory and Applications” and by the Austrian Science Fund (FWF): Project P21196. This paper was written while G. Leobacher was member of the Department of Financial Methematics and Applied Number Theory, Johannes Kepler University Linz, Altenbergerstraße 69, 4040 Linz, Austria. M. Szölgyenyi is supported by the Vienna Science and Technology Fund (WWTF): Project MA14-031.
Author information
Authors and Affiliations
Corresponding author
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
About this article
Cite this article
Eichler , A., Leobacher, G. & Szölgyenyi, M. Utility indifference pricing of insurance catastrophe derivatives. Eur. Actuar. J. 7, 515–534 (2017). https://doi.org/10.1007/s13385-017-0154-2
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s13385-017-0154-2
Keywords
- Insurance mathematics
- Catastrophe derivatives
- Utility indifference pricing
- Modeling catastrophe losses
- Piecewise deterministic Markov process