# Optimal Reinsurance: A Risk Sharing Approach

^{1}

^{2}

^{3}

^{*}

*Keywords:*optimal reinsurance; general risk measure; risk sharing; systemic risk

Next Article in Journal / Special Issue

University Carlos III of Madrid. CL. Madrid 126. 28903 Getafe, Madrid, Spain

University of Castilla la Mancha Avda. Real Fábrica de Seda, s/n. 45600 Talavera, Toledo, Spain

University Complutense of Madrid. Department of Actuarial and Financial Economics. Somosaguas-Campus. 28223 Pozuelo de Alarcón, Madrid, Spain

Author to whom correspondence should be addressed.

Received: 11 June 2013 / Revised: 21 July 2013 / Accepted: 24 July 2013 / Published: 5 August 2013

(This article belongs to the Special Issue Systemic Risk and Reinsurance)

This paper proposes risk sharing strategies, which allow insurers to cooperate and diversify non-systemic risk. We deal with both deviation measures and coherent risk measures and provide general mathematical methods applying to optimize them all. Numerical examples are given in order to illustrate how efficiently the non-systemic risk can be diversified and how effective the presented mathematical tools may be. It is also illustrated how the existence of huge disasters may lead to wrong solutions of our optimal risk sharing problem, in the sense that the involved risk measure could ignore the existence of a non-null probability of "global ruin" after the design of the optimal risk sharing strategy. To overcome this caveat, one can use more conservative risk measures. The stability in the large of the optimal sharing plan guarantees that "the global ruin caveat" may be also addressed and solved with the presented methods.

The optimal reinsurance problem is a classic topic in Actuarial Mathematics. A common approach is to minimize some measure of the first insurer risk after reinsurance. Seminal papers by Borch [1] and Arrow [2] used the variance as the risk measure and proved that the stop-loss reinsurance minimizes the retained risk if premiums are calculated with the Expected Value Premium Principle ($EVPP$). The subsequent research followed similar ideas and tried to take into account more general risk measures and premium principles, which may give optimal contracts other than stop-loss. Recently, Gajec and Zagrodny [3] considered more general symmetric and even asymmetric risk functions, such as the absolute deviation and the truncated variance of the retained loss, under the standard deviation premium principle. Kaluszka [4] studied reinsurance contracts with many convex premium principles (exponential, semi-deviation and semi-variance, Dutch, distortion, etc.). Other well known financial risk measures, such as the Value at Risk ($VaR$) or the Conditional Value at Risk ($CVaR$, also called $AVaR$, $TVaR$, $CTE$,Expected Shortfall, etc.), are also being considered. For example, Kaluszka [4] uses the $CVaR$ as a premium principle, and Cai and Tan [5] calculate the optimal retention for a stop-loss reinsurance by considering the $VaR$ and the $CVaR$, under the $EVPP$. Important extensions are also presented in Cai et al. [6] and Bernard and Tian [7]. Balbás et al. [8] seem to be the first authors proving that the stop-loss reinsurance is still optimal for the $EVPP$ and the $CVaR$, though the reinsurance budget does not have to be always saturated. An interesting summary of all of these findings may be found in Centeno and Simoes [9], though there are still many open problems, and more recent papers may be found (Seng et al. [10], Cui et al. [11], Chi and Tanb [12], etc.).

It is not so usual in the literature to include the reinsurer point of view. The main reason is that the final decision is made by the ceding company. However, the reinsurer must also accept the selected strategy. A third approach considers the objectives of both companies and proposes a “reciprocal reinsurance”. A very recent study may be found in Cai et al. [13], though there are more analyses involving both optimization methods and the theory of cooperative games.

This paper deals with “reciprocal reinsurance contracts” involving a general number of companies. In other words, we propose a risk sharing approach in order to diversify the risk as much as possible, so as to make the “global market risk” (or systemic risk, in this paper) as close as possible to the total sum of partial risks. Risks are measured in a very general setting, since Coherent Risk Measures (Artzner et al. [14]), Deviation Measures (Rockafellar et al. [15] and Expectation Bounded Risk Measures (Rockafellar et al. [15]) may be used.

The paper outline is as follows. Section 2 will be devoted to introducing the main notations and the framework. In Section 3, we will present an “Optimal Risk Sharing Problem” ($ORSP$) in such a way that every insurer attempts to conserve the size of his portfolio (market share) and simultaneously reduce his risk. As said above, it may be done by means of Game Theory or Vector Optimization, and we have selected the second approach, due to the mathematical problems that the use of risk measures may cause. Thus, we propose a two stage $ORSP$, such that, first of all (stage 1), every insurer computes his ideal value (minimum risk with a non-lower market share). As usual, in Vector Optimization, the ideal values compose the Ideal Point (Nakayama et al. [16]), and a vector optimization problem (second stage) enables the insurers to look for a risk sharing plan respecting the ideal point and every market share, as well as approaching the systemic risk as much as possible.

The proposed problems are not easy to solve in practice, since it requires the minimization of general risk functions, which are barely differentiable. Standard mathematical methods do not apply, and Section 4 is devoted to extending a methodology of Balbás et al. [8,17]. We will see that Linear Programming linked methods will be sufficient to overcome the mathematical difficulties, despite the fact that we may deal with both discrete and continuous random variables and the risk measure is quite general.

Section 5 and Section 6 present illustrative numerical examples. We only try to clarify how the developed mathematical methods allow us to solve, in practice, the proposed optimization problems in a simple manner. Simultaneously, the examples will illustrate the practical properties of the optimal risk sharing strategy and a caveat that might arise from time to time. The existence of huge disasters might lead to wrong solutions of our $ORSP$, in the sense that the involved risk measure could ignore the existence of a non-null probability of “global ruin” ( i.e., the ruin of every involved company) after the design of the optimal risk sharing strategy. To overcome this caveat, one can use more conservative risk measures. The stability in the large of the optimal sharing plan guarantees that “this global ruin caveat” may be also addressed with the presented methods.

The last section of the paper summarizes the most important conclusions.

Consider the probability space, $\left(\Omega ,\mathcal{F},\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{P}\right)$, composed of the set of “states of the world” Ω, the $\sigma -$algebra $\mathcal{F}$ and the probability measure, $\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{P}$. Denote by $\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)$ the mathematical expectation of every $\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{R}-$valued random variable y, defined on Ω. Let $p\in \left(1,\infty \right)$ and denote by ${L}^{p}$ the Banach space of random variables, y, on Ω, such that $\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\left|y\right|}^{p}\right)<\infty $, endowed with the norm:
for every $y\in {L}^{p}$. According to the Riesz Representation Theorem, ${L}^{q}$ is the dual space of ${L}^{p}$, where $q\in \left(1,\infty \right)$ is characterized by $1/p+1/q=1$.

$${\u2225y\u2225}_{p}={\left(\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\left|y\right|}^{p}\right)\right)}^{1/p}$$

Let $\left[0,T\right]$ be a time interval. From an intuitive point of view, one can interpret that $y\in {L}^{p}$ may represent claims at T for some arbitrary insurer.

Consider n insurance companies, whose final (within $\left[0,T\right]$) claims will be represented by the non-negative random variables, ${y}_{1},{y}_{2},....{y}_{n}\in {L}^{p}$, respectively. In order to simplify notations, let us assume that all of them deal with the Expected Value Premium Principle with the same loading rate, γ. Thus, the random final (at T) of the wealth of the ${j}^{th}-$company will be $\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)-{y}_{j}$, $j=1,2,...,n$.

In order to simplify the framework, assume that all of the insurers deal with the same risk measure:
Consider the sub-gradient of ρ:
We will assume that ${\Delta}_{\rho}$ is convex and $\sigma \left({L}^{q},{L}^{p}\right)-$compact, and:
holds for every $y\in {L}^{p}$. Furthermore, we will also assume the existence of ${\tilde{E}}_{\rho}\ge 0$, such that the constant random variable, $z={\tilde{E}}_{\rho}$, is in ${\Delta}_{\rho}$ and:
Summarizing, we have:

$$\rho :{L}^{p}\u27f6\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{R}$$

$${\Delta}_{\rho}=\left\{z\in {L}^{q};-\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(yz\right)\le \rho \left(y\right),\phantom{\rule{4.pt}{0ex}}\forall y\in {L}^{p}\right\}\subset {L}^{q}.$$

$$\rho \left(y\right)=Max\phantom{\rule{4.pt}{0ex}}\left\{-\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(yz\right):z\in {\Delta}_{\rho}\right\}$$

$${\Delta}_{\rho}\subset \left\{z\in {L}^{q};\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(z\right)={\tilde{E}}_{\rho}\right\}.$$

Assumption 1 is not at all restrictive, since it is satisfied by every expectation bounded risk measure (Rockafellar et al. [15]) with ${\tilde{E}}_{\rho}=1$ and by every deviation measure (Rockafellar et al. [15]) with ${\tilde{E}}_{\rho}=0$. Examples of expectation bounded risk measures are the Conditional Value at Risk ($CVaR$) and the Weighted Conditional Value at Risk ($WCVaR$), amongst many others. Examples of deviation measures are, amongst others, the classical $p-$deviation:
or the upside and downside $p-$semi-deviation:
and:
If ${\tilde{E}}_{\rho}=1$, then it is easy to see that ρ is also coherent in the sense of Artzner et al. [14], if and only if:
Further details may be also found in Balbás et al. [8, 17].

$${\sigma}_{p}\left(y\right)={\left[\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\left|\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right|}^{p}\right)\right]}^{1/p},$$

$${\sigma}_{p}^{+}\left(y\right)={\left[\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\left|Max\phantom{\rule{4.pt}{0ex}}\left\{y-\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right),0\right\}\right|}^{p}\right)\right]}^{1/p}$$

$${\sigma}_{p}^{-}\left(y\right)={\left[\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\left|Max\phantom{\rule{4.pt}{0ex}}\left\{\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y,0\right\}\right|}^{p}\right)\right]}^{1/p}.$$

$${\Delta}_{\rho}\subset {L}_{+}^{q}=\left\{z\in {L}^{q};\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{P}\left(z\ge 0\right)=1\right\}.$$

Under the framework above, we can consider the risk of the ${j}^{th}-$company, given by:
$j=1,2,...,n$, and the systemic risk, given by:
where:
denotes “global claims”. Since Assumption 1 implies that ρ is sub-additive, we have that:

$$\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)-{y}_{j}\right)$$

$$\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)$$

$$y=\sum _{j=1}^{n}{y}_{j}\in {L}^{p}$$

$$\begin{array}{c}\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)=\rho \left(\left(1+\gamma \right)\sum _{j=1}^{n}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)-\sum _{j=1}^{n}{y}_{i}\right)\\ \le \sum _{j=1}^{n}\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)-{y}_{j}\right).\end{array}$$

Expression (6) may suggest that every company can reduce risk and simultaneously improve expected profits. Actually, companies might attempt to modify their portfolios in such a manner that (6) almost becomes an equality. This purpose may be addressed with two stages.

In the first stage, every insurer can compute his ideal risk level, i.e., the minimum risk associated with the expected wealth of his portfolio. The optimization problem for the ${j}^{th}-$company becomes:
${x}_{j}\in {L}^{p}$ being the decision variable.

$$\left\{\begin{array}{c}Min\phantom{\rule{4.pt}{0ex}}\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)-{x}_{j}\right)\hfill \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)\ge \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)\hfill \\ 0\le {x}_{j}\le y\hfill \end{array}\right.$$

Fix an ideal point,
i.e., ${I}_{j}$ is not higher than the optimal value attained by (7), $j=1,2,...,n$. In the second stage, all the companies, such that:
could collaborate, so as to improve their particular risk. Without loss of generality, we can assume that (8) holds for $j=1,2,...,n$, so the insurers could look for a risk sharing plan solving the optimization problem:
$x=\left({x}_{1},{x}_{2},...,{x}_{n}\right)\in {\left({L}^{p}\right)}^{n}$ being the decision variable and ${w}_{j}>0$ being a weight associated with the ${j}^{th}-$insurer, $j=1,2,...,n$. Every company is guaranteeing that his portfolio size (or market share) will not decrease, and the weights $w={\left({w}_{j}\right)}_{j=1}^{n}$ must be chosen in such a manner that:
$j=1,2,...,n$ must hold. Thus, every company will be approaching its ideal risk level, and the solution:
of (9) will be a Pareto optimum, in the sense that there is not any (9)−feasible allocation,
,such that:
$j=1,2,...,n$, with one strict inequality at least (Nakayama et al. [16]).

$$I=\left({I}_{1},{I}_{2},...,{I}_{n}\right)$$

$${I}_{j}<{\rho}_{j}=\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)-{y}_{j}\right)$$

$$\left\{\begin{array}{cc}Min\phantom{\rule{4.pt}{0ex}}\sum _{j=1}^{n}{w}_{j}\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)-{x}_{j}\right)\hfill & \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)\ge \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right),\hfill & j=1,2,...,n\hfill \\ \sum _{j=1}^{n}{x}_{j}\le y\hfill & \\ 0\le {x}_{j},\hfill & j=1,2,...,n\hfill \end{array}\right.$$

$${I}_{j}<\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)-{x}_{j}\right)\le {\rho}_{j},$$

$${x}^{*}=\left({x}_{1}^{*},{x}_{2}^{*},...,{x}_{n}^{*}\right)\in {\left({L}^{p}\right)}^{n}$$

$$x=\left({x}_{1},{x}_{2},...,{x}_{n}\right)\in {\left({L}^{p}\right)}^{n}$$

$$\rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)-{x}_{j}\right)\le \rho \left(\left(1+\gamma \right)\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}^{*}\right)-{x}_{j}^{*}\right),$$

Notice that the constraints of (9) lead to:
and:
$j=1,2,...,n$, so, if necessary, we can slightly modify the restrictions of (9).

$$\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)=\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)$$

$$\sum _{j=1}^{n}{x}_{j}=y$$

Both optimization problems, (7) and (9), involve the risk function, ρ, which is non-differentiable in general, and therefore, the standard mathematical methods do not apply any more. In order to overcome this caveat, several authors have proposed equivalent linear problems that apply for discrete random variables and particular interesting examples, such as the $CVaR$ or the absolute deviation (Konno et al. [18], Mansini et al. [19], etc.). Balbás et al. [8] dealt with a particular optimal reinsurance problem and gave a new linear dual problem characterizing the primal solutions and applied in a much more general setting, since the probability space, $\left(\Omega ,\mathcal{F},\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{P}\right)$, does not have to be discrete, and the fulfillment of Assumption 1 is the unique hypothesis about the risk measure, ρ. This section is devoted to pointing out how the approach of Balbás et al. [8] may be significantly extended, and both (7) and (9) have a linear dual problem characterizing their solutions. This is important from a computational point of view, since it will allow us to find, in practice, the optimal allocation, $x=\left({x}_{1},{x}_{2},...,{x}_{n}\right)\in {\left({L}^{p}\right)}^{n}$.

With regard to Problem (7), following Balbás et al. [8], and bearing in mind (2) and (3), there is an equivalent linear problem with infinitely many constraints, namely:
$\left(\theta ,{x}_{j}\right)\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{R}\times {L}^{p}$ being the decision variable. Indeed, bearing in mind Expression (2), one can easily prove the equivalence between (7) and (13) with quite parallel arguments to those in Balbás et al. [8]. Problem (13) shows an advantage with respect to Problem (7), because it is linear, but the drawback is provoked by the existence of infinitely many constrains (there is one constraint per every element, ${z}_{j}\in {\Delta}_{\rho}$). Nevertheless, the dual problem of (13) overcomes this caveat, since it remains linear and may be given as follows:
$\left({\lambda}_{j},{z}_{j},{\nu}_{j}\right)\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{R}\times {L}^{q}\times {L}^{q}$ being the decision variable. We can take ${I}_{j}$ as the optimal value of the linear dual problem above, so as to define the ideal point, I.

$$\left\{\begin{array}{cc}Min\phantom{\rule{4.pt}{0ex}}\theta \hfill & \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}{z}_{j}\right)-\left(1+\gamma \right){\tilde{E}}_{\rho}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)\le \theta ,\hfill & \forall {z}_{j}\in {\Delta}_{\rho}\hfill \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)\ge \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right),\hfill & \\ 0\le {x}_{j}\le y\hfill \end{array}\right.$$

$$\left\{\begin{array}{c}Max\phantom{\rule{4.pt}{0ex}}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right){\lambda}_{j}-\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y{\nu}_{j}\right)\hfill \\ {z}_{j}+{\nu}_{j}\ge \left(1+\gamma \right){\tilde{E}}_{\rho}+{\lambda}_{j}\hfill \\ {\lambda}_{j}\ge 0,\phantom{\rule{4.pt}{0ex}}{z}_{j}\in {\Delta}_{\rho},\phantom{\rule{4.pt}{0ex}}{\nu}_{j}\in {L}^{q},\phantom{\rule{4.pt}{0ex}}{\nu}_{j}\ge 0\hfill \end{array}\right.$$

With respect to Problem (9), the equivalent linear problem with infinitely many constraints becomes:
$\left({\left({\theta}_{j}\right)}_{j=1}^{n},{\left({x}_{j}\right)}_{j=1}^{n}\right)\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}{\mathsf{R}}^{n}\times {\left({L}^{p}\right)}^{n}$ being the decision variable, while the linear dual problem is:
being the decision variable.

$$\left\{\begin{array}{cc}Min\phantom{\rule{4.pt}{0ex}}\sum _{j=1}^{n}{\theta}_{j}\hfill & \\ {w}_{j}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}{z}_{j}\right)-{w}_{j}\left(1+\gamma \right){\tilde{E}}_{\rho}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)-{\theta}_{j}\le 0,\hfill & \forall {z}_{j}\in {\Delta}_{\rho},\phantom{\rule{4.pt}{0ex}}j=1,2,...,n\hfill \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}\right)\ge \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right),\hfill & j=1,2,...,n\hfill \\ \sum _{j=1}^{n}{x}_{j}\le y\hfill & \\ 0\le {x}_{j},\hfill & j=1,2,...,n\hfill \end{array}\right.$$

$$\left\{\begin{array}{cc}Max\phantom{\rule{4.pt}{0ex}}\sum _{j=1}^{n}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right){\lambda}_{j}-\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\nu \right)\hfill & \\ -\left(1+\gamma \right){w}_{j}{\tilde{E}}_{\rho}+{w}_{j}{z}_{j}-{\lambda}_{j}+\nu \ge 0,\hfill & j=1,2,...,n\hfill \\ {\lambda}_{j}\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{R},\phantom{\rule{4.pt}{0ex}}{\lambda}_{j}\ge \mathsf{0},\phantom{\rule{4.pt}{0ex}}\nu \in {L}^{q},\nu \ge 0,\phantom{\rule{4.pt}{0ex}}{z}_{j}\in {\Delta}_{\rho}\hfill & j=1,2,...,n\hfill \end{array}\right.$$

$$\left({\left({\lambda}_{j}\right)}_{j=1}^{n},{\left({z}_{j}\right)}_{j=1}^{n},\nu \right)\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}{\mathsf{R}}^{n}\times {L}^{q}\times {L}^{q}$$

We are dealing with infinite-dimensional Banach spaces, so the existence of a duality gap between (9) and (15) might hold, i.e., both problems might attain different optimal values (Luenberger, [20]). However, similar methods to those in Balbás et al. [8] enable us to prove the duality gap absence, and (15) totally characterizes the optimal allocation, ${\left({x}_{j}^{*}\right)}_{j=1}^{n}$, solving (9). Thus, in practice, one can solve the dual problem, which is linear and, therefore, easy to solve, even if the probability space, $\left(\Omega ,\mathcal{F},\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{P}\right)$, is complex and (15) is an infinite-dimensional problem (Anderson and Nash, [21]), and then, one can find the optimal allocation, ${\left({x}_{j}^{*}\right)}_{j=1}^{n}$, by applying the Karush-Kuhn-Tucker conditions below.

Suppose that ${\left({x}_{j}^{*}\right)}_{j=1}^{n}$ is (9)-feasible and:
is (15)-feasible. Then, they solve (9) and (15), if and only if the following Karush-Kuhn-Tucker-like conditions:
hold. ☐

$$\left({\left({\lambda}_{j}^{*}\right)}_{j=1}^{n},{\left({z}_{j}^{*}\right)}_{j=1}^{n},{\nu}^{*}\right)$$

$$\left\{\begin{array}{cc}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}^{*}{z}_{j}\right)\le \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}^{*}{z}_{j}^{*}\right),\hfill & \forall {z}_{j}\in {\Delta}_{\rho},\phantom{\rule{4.pt}{0ex}}j=1,2,...,n\hfill \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}^{*}\left(-\left(1+\gamma \right){w}_{j}{\tilde{E}}_{\rho}+{w}_{j}{z}_{j}^{*}-{\lambda}_{j}^{*}+{\nu}^{*}\right)\right)=0,\hfill & \phantom{\rule{4.pt}{0ex}}j=1,2,...,n\hfill \\ \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{j}^{*}-{y}_{j}\right)=0,\hfill & j=1,2,...,n\hfill \\ \sum _{j=1}^{n}{x}_{j}^{*}=y\hfill \end{array}\right.$$

Let us deal with a simple numerical example in order to illustrate how the theorem above allows us to solve Problems (7) and (9) and obtain the optimal allocation, ${\left({x}_{j}^{*}\right)}_{j=1}^{n}$. This numerical experiment will be used in Section 6 to point out how the optimal allocation may increase the “probability global of ruin”.

Consider three companies and a set of states composed of five scenarios, $\Omega =\left\{{\omega}_{1},{\omega}_{2},{\omega}_{3},{\omega}_{4},{\omega}_{5}\right\}$. Suppose that the probability of all of them equals $0.2$. Matrix:
gives the random cost for annual claims, the ${j}^{th}-$row being associated with the ${j}^{th}-$company, $j=1,2,3$. Suppose, finally, that $\gamma =0.05$ and that ρ is the $CVaR$ with the confidence level, $60\%$ (so ${\tilde{E}}_{\rho}=1$). Though higher levels of confidence are usually recommended, the selected one, $60\%$, will enable us to illustrate many effects we are interested in. Furthermore, the $CVaR$ is consistent with the second order stochastic dominance (and the standard utility functions) for every level of confidence (Ogryczak and Ruszczynski, [22]).

$${C}_{y}=\left(\begin{array}{ccccc}2,& 4,& 6,& 8,& 10\\ 10,& 8,& 6,& 4,& 2\\ 1,& 1,& 1,& 1,& 1000\end{array}\right)$$

It is easy to verify that:
and the systemic risk, $\rho \left(y\right)=289.06$, is strictly lower than the sum of the partial risks, ${\rho}_{1}+{\rho}_{2}+{\rho}_{3}=295.06$. The three companies may attempt to remove the difference (six monetary units) by means of a risk sharing plan.

$$\left(\begin{array}{cccc}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{1}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{2}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{3}\right)=200.8,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{1}+{y}_{2}+{y}_{3}\right)=212.8\\ {\rho}_{1}=2.7,& {\rho}_{2}=2.7,& {\rho}_{3}=289.66,& \begin{array}{c}\rho \left(1.05\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)=289.06\hfill \\ {\rho}_{1}+{\rho}_{2}+{\rho}_{3}=295.06\hfill \end{array}\end{array}\right)$$

Bearing in mind that (Rockafellar et al., [15])
Problem (14) becomes:
$\left(\lambda ,{\left({\nu}_{i}\right)}_{i=1}^{5},{\left({z}_{i}\right)}_{i=1}^{5}\right)\in \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}{\mathsf{R}}^{11}$ being the decision variable and $\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)$ equaling six for $j=1,2$ or $200.8$ for $j=3$. By solving these simple linear problems, we can obtain the optimal value of (7) for the three involved insurers, and the ideal point becomes
As can be seen, the ideal point would significantly improve the risk level of the three companies, while their expected profit would remain the same. The first and second company could reach a negative $CVaR$, and the risk reduction could equal $\frac{2.7+0.65}{2.7}=124.07\%$, whereas this percentage would become $\frac{289.66-271.66}{289.66}\simeq 6.21\%$ for the third insurer. In this situation, a risk sharing plan could be interesting for them all, but recall that the ideal point above is not reachable (actually, the systemic risk, $\rho \left(1.05\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)=289.06$, can never be lower than the sum of partial risks; see (6)). Besides, a significant fall of the third risk could provoke a positive increment of the rest of the ones, so the weights, ${\left({w}_{j}\right)}_{j=1}^{3}$, should be selected, so as to prevent this situation. Suppose that the three companies choose ${\left({w}_{j}\right)}_{j=1}^{3}=\left(10,10,1\right)$. Then, Problem (9) will lead to a Pareto optimum, and Problem (15) will become:
with
being the decision variable and
Solving this problem with standard linear optimization methods and using Conditions (17), we get the optimal allocation:
and the summary (19) becomes:
Thus, the equality:
shows that the three companies can totally diversify the non-systemic risk in this case. Notice that the ideal risk level is not reached by any company.

$${\Delta}_{\rho}=\left\{z=\left({z}_{1},...,{z}_{5}\right);\phantom{\rule{4.pt}{0ex}}0\le {z}_{i}\le 2.5\phantom{\rule{4.pt}{0ex}}and{\phantom{\rule{4.pt}{0ex}}}_{i=1}^{5}{z}_{i}=5\right\},$$

$$\left\{\begin{array}{cc}Max\phantom{\rule{4.pt}{0ex}}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({y}_{j}\right)\lambda -0.{2}_{i=1}^{5}{y}_{i}{\nu}_{i}\hfill & \\ -\lambda +{z}_{i}+{\nu}_{i}\ge 1.05,\hfill & i=1,2,...,5\hfill \\ {}_{i=1}^{5}{z}_{i}=5\hfill & \\ 0\le {z}_{i}\le 2.5,\hfill & i=1,2,...,5\hfill \\ 0\le \lambda ,\phantom{\rule{4.pt}{0ex}}0\le {\nu}_{i},\hfill & i=1,2,...,5\hfill \end{array}\right.$$

$$I=\left(-0.65,-0.65,271.66\right).$$

$$\left\{\begin{array}{cc}Max\phantom{\rule{4.pt}{0ex}}6\left({\lambda}_{1}+{\lambda}_{2}\right)+200.8{\lambda}_{3}-0.2\sum _{i=1}^{5}{y}_{i}{\nu}_{i}\hfill & \\ -{\lambda}_{j}+{w}_{j}{z}_{j,i}+{\nu}_{i}\ge 1,05{w}_{j},\hfill & \begin{array}{c}j=1,2,3\hfill \\ i=1,...,5\hfill \end{array}\hfill \\ {}_{i=1}^{5}{z}_{j,i}=5,\hfill & j=1,2,3\hfill \\ {\lambda}_{j}\ge \mathsf{0},\phantom{\rule{4.pt}{0ex}}0\le {z}_{j,i}\le 2.5,\phantom{\rule{4.pt}{0ex}}0\le {\nu}_{i},\hfill & \begin{array}{c}j=1,2,3\hfill \\ i=1,...,5\hfill \end{array}\hfill \end{array}\right.$$

$$\left({\left({\lambda}_{j}\right)}_{j=1}^{3},{\left({\left({z}_{j,i}\right)}_{j=1}^{3}\right)}_{i=1}^{5},{\left({\nu}_{i}\right)}_{i=1}^{5}\right)$$

$${\left({y}_{i}\right)}_{i=1}^{5}=\left(13,13,13,13,1012\right).$$

$${C}_{x}=\left(\begin{array}{ccccc}6,& 6,& 6,& 6,& 6\\ 6,& 6,& 6,& 6,& 6\\ 1,& 1,& 1,& 1,& 1000\end{array}\right)$$

$$\left(\begin{array}{cccc}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{1}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{2}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{3}\right)=200.8,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({x}_{1}+{x}_{2}+{x}_{3}\right)=212.8\\ {\rho}_{1}=-0.3,& {\rho}_{2}=-0.3,& {\rho}_{3}=289.66,& \begin{array}{c}\rho \left(1.05\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)=289.06\hfill \\ {\rho}_{1}+{\rho}_{2}+{\rho}_{3}=289.06\hfill \end{array}\end{array}\right).$$

$$\rho \left(1.05\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)={\rho}_{1}+{\rho}_{2}+{\rho}_{3}$$

Consider the numerical example above, but suppose that the second stage is addressed with the alternative weights, ${\left({\tilde{w}}_{j}\right)}_{j=1}^{3}=\left(1,1,1\right)$. Then, it is easy to repeat the process and get the new allocation:
Thus, (19) and (21) become:
Obviously, this new solution would not be accepted by the first and second insurers, since their particular risk increases from $2.7$ to $8.7$. As indicated in Section 4, the weights in (9) must be selected in such a manner that (10) holds. However, we could also provide other numerical examples, such that (10) would hold, and the caveat below would also apply.

$${C}_{\tilde{x}}=\left(\begin{array}{ccccc}0,& 0,& 0,& 0,& 30\\ 0,& 0,& 0,& 0,& 30\\ 13,& 13,& 13,& 13,& 952\end{array}\right).$$

$$\left(\begin{array}{cccc}\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\tilde{x}}_{1}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\tilde{x}}_{2}\right)=6,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\tilde{x}}_{3}\right)=200.8,& \mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left({\tilde{x}}_{1}+{\tilde{x}}_{2}+{\tilde{x}}_{3}\right)=212.8\\ {\rho}_{1}=8.7,& {\rho}_{2}=8.7,& {\rho}_{3}=271.66,& \begin{array}{c}\rho \left(1.05\mathsf{I}\phantom{\rule{-1.5pt}{0ex}}\mathsf{E}\left(y\right)-y\right)=289.06\hfill \\ {\rho}_{1}+{\rho}_{2}+{\rho}_{3}=289.06\hfill \end{array}\end{array}\right).$$

Indeed, the risk sharing strategy, ${C}_{\tilde{x}}$, shows the existence of serious caveats that cannot be ignored when designing optimal reinsurance problems. Expression (21) implies that ${C}_{\tilde{x}}$ also satisfies (22), and therefore, both ${C}_{x}$ and ${C}_{\tilde{x}}$ diversify the non-systemic risk. There is only systemic risk after both reinsurance strategies. However, the major difference is caused by the state of nature, ${\omega}_{5}\in \Omega $. If ${\omega}_{5}$ comes out, then the three companies will be facing significant capital losses. In other words, if one compares Strategies ${C}_{y}$ and ${C}_{\tilde{x}}$, the probability of ruin for the three insurers is clearly higher for ${C}_{\tilde{x}}$, despite the fact that this risk sharing plan eliminates the non-systemic risk, while the initial one, ${C}_{y}$, does not.

In general, the existence of very negative catastrophes implying high capital losses generates loss-distributions very asymmetric and with very heavy tails, as well as very closely correlated risks in the insurance industry. In a risk sharing plan, it is very important to diversify the risk provoked by very severe disasters, which can only be done by dealing with very conservative risk measures, related to very risk averse decision makers. Otherwise, more aggressive risk measures will not detect this contagion effect and might indicate a global diversification of the non-systemic risk and, simultaneously, lead to risk sharing strategies, making it grow the ruin probability of the whole system.

The analysis of Balbás et al. [17] may be easily adapted to the setting of this paper. The implication is that, under weak conditions, the optimal risk sharing strategy remains stable as the risk measure becomes more and more conservative. In fact, there is a limit in the large of this strategy. Thus, the method proposed in Section 3 must be complemented with a new analysis verifying the ruin probability of the involved companies. If the result is not adequate, then the optimal strategy must be reached with a more risk averse risk measure. The stability of the optimal strategy in the large will guarantee that “the probability of global ruin” may be controlled, even when facing heavy tails and high correlations. The example in Section 5 shows that the methodology of Section 3 usually leads to successful solutions (${C}_{x}$ in (20) totally diversifies the non-systemic risk and does not make it grow any probability of ruin), but strategies, such as ${C}_{\tilde{x}}$, should be discarded if they were obtained.

The optimal reinsurance problem is a classical topic in Actuarial Mathematics. The usual viewpoint only considers the ceding company objective, though there are other approaches, taking into account the reinsurer opinion, too. This paper proposes a “reciprocal reinsurance” involving a general number of companies. The main purpose is the diversification of the non-systemic risk conserving every market share.

We have addressed the objective above by means of two stage Vector Optimization Problems. In the first step, every insurer computes his ideal risk level, and the second step provides an optimal risk sharing plan that integrates the objectives of all of the involved insurers and respects ideal points and market shares.

The usual mathematical methods do not apply to solve the proposed optimization problems, due to the lack of differentiability, so we have provided specific mathematical tools permitting us to give explicit solutions of the presented problems. These tools apply in a very general framework, since we can deal with both discrete and continuous probability spaces, and the risk measure assumptions are quite weak.

Numerical examples have shown how the proposed tools apply in practice, as well as the properties of the reached optimal risk sharing strategy. Furthermore, the examples have shown that under some particular conditions, the selected risk sharing plan might provoke the existence of scenarios that are very negative for all of the involved companies, making it grow the “probability of global bankruptcy”, with respect to this probability value, before the reinsurance contract. If so, the risk sharing plan should be modified, and a good way to do that is to choose a more conservative risk measure. The stability in the large of the optimal sharing plan guarantees that this “global ruin caveat” may be always solved. ☐

Research partially supported by “Comunidad Autónoma de Madrid” (Spain, Grant $S2009/ESP\u20131594$) and “MEyC” (Spain, Grants $ECO2009\u201314457\u2013C04$ and $ECO2012\u201339031\u2013C02\u201301$). The usual caveat applies.

- K. Borch. “An attempt to determine the optimum amount of stop loss reinsurance.” Transactions of the 16th International Congress of Actuaries I 1 (1960): 597–610. [Google Scholar]
- K.J. Arrow. “Uncertainty and the welfare of medical care.” Am. Econ. Rev. 53 (1963): 941–973. [Google Scholar]
- L. Gajec, and D. Zagrodny. “Optimal reinsurance under general risk measures.” Insur. Math. Econ. 34 (2004): 227–240. [Google Scholar] [CrossRef]
- M. Kaluszka. “Optimal reinsurance under convex principles of premium calculation.” Insur. Math. Econ. 36 (2005): 375–398. [Google Scholar] [CrossRef]
- J. Cai, and K.S. Tan. “Optimal retention for a stop loss reinsurance under the VaR and CTE risk measures.” ASTIN Bull. 37 (2007): 93–112. [Google Scholar]
- J. Cai, K.S. Tan, C. Weng, and Y. Zhang. “Optimal reinsurance under VaR and CTE risk measures.” Insur. Math. Econ. 43 (2008): 185–196. [Google Scholar] [CrossRef]
- C. Bernard, and W. Tian. “Optimal reinsurance arrangements under tail risk measures.” J. Risk Insur. 76 (2009): 709–725. [Google Scholar] [CrossRef]
- A. Balbás, B. Balbás, and A. Heras. “Optimal reinsurance with general risk measures.” Insur. Math. Econ. 44 (2009): 374–384. [Google Scholar] [CrossRef]
- M.L. Centeno, and O. Simoes. “Optimal reinsurance.” RACSAM 103 (2009): 387–405. [Google Scholar] [CrossRef]
- K.T. Seng, C. Wenga, and Y. Zhang. “Optimality of general reinsurance contracts under CTE risk measure.” Insur. Math. Econ. 49 (2011): 175–187. [Google Scholar]
- W. Cui, J. Yang, and L. Wu. “Optimal reinsurance minimizing the distortion risk measure under general reinsurance premium principles.” Insur. Math. Econ. 53 (2013): 74–85. [Google Scholar] [CrossRef]
- Y. Chi, and K.S. Tanb. “Optimal reinsurance with general premium principles.” Insur. Math. Econ. 52 (2013): 180–189. [Google Scholar] [CrossRef]
- J. Cai, Y. Fang, Z. Li, and G.E. Willmot. “Optimal reciprocal reinsurance treaties under the joint survival probability and the joint profitable probability.” J. Risk Insur. 80 (2012): 145–168. [Google Scholar] [CrossRef]
- P. Artzner, F. Delbaen, J.M. Eber, and D. Heath. “Coherent measures of risk.” Math. Finance 9 (1999): 203–228. [Google Scholar] [CrossRef]
- R.T. Rockafellar, S. Uryasev, and M. Zabarankin. “Generalized deviations in risk analysis.” Financ. Stoch. 10 (2006): 51–74. [Google Scholar] [CrossRef]
- H. Nakayama, Y. Sawaragi, and T. Tanino. Theory of Multiobjective Optimization. Waltham, MA, USA: Academic Press, 1985. [Google Scholar]
- A. Balbás, B. Balbás, and A. Heras. “Stable solutions for optimal reinsurance problems involving risk measures.” Eur. J. Oper. Res. 214 (2011): 796–804. [Google Scholar] [CrossRef]
- H. Konno, K. Akishino, and R. Yamamoto. “Optimization of a long-short portfolio under non-convex transaction costs.” Comput. Optim. Appl. 32 (2005): 115–132. [Google Scholar] [CrossRef]
- R. Mansini, W. Ogryczak, and M.G. Speranza. “Conditional value at risk and related linear programming models for portfolio optimization.” Ann. Oper. Res. 152 (2007): 227–256. [Google Scholar] [CrossRef]
- D.G. Luenberger. Optimization by Vector Spaces Methods. Hoboken, NJ, USA: John Wiley & Sons, 1969. [Google Scholar]
- E.J. Anderson, and P. Nash. Linear Programming in Infinite-Dimensional Spaces. Hoboken, NJ, USA: John Wiley & Sons, 1987. [Google Scholar]
- W. Ogryczak, and A. Ruszczynski. “Dual stochastic dominance and related mean risk models.” SIAM J. Optimiz. 13 (2002): 60–78. [Google Scholar] [CrossRef]

© 2013 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).