Heath–Jarrow–Morton modelling of longevity bonds and the risk minimization of life insurance portfolios

This paper has two parts. In the first, we apply the Heath–Jarrow–Morton (HJM) methodology to the modelling of longevity bond prices. The idea of using the HJM methodology is not new. We can cite Cairns et al. [Cairns A.J., Blake D., Dowd K, 2006. Pricing death: framework for the valuation and the securitization of mortality risk. Astin Bull., 36 (1), 79–120], Miltersen and Persson [Miltersen K.R., Persson S.A., 2005. Is mortality dead? Stochastic force of mortality determined by arbitrage? Working Paper, University of Bergen] and Bauer [Bauer D., 2006. An arbitrage-free family of longevity bonds. Working Paper, Ulm University]. Unfortunately, none of these papers properly defines the prices of the longevity bonds they are supposed to be studying. Accordingly, the main contribution of this section is to describe a coherent theoretical setting in which we can properly define these longevity bond prices. A second objective of this section is to describe a more realistic longevity bonds market model than in previous papers. In particular, we introduce an additional effect of the actual mortality on the longevity bond prices, that does not appear in the literature. We also study multiple term structures of longevity bonds instead of the usual single term structure. In this framework, we derive a no-arbitrage condition for the longevity bond financial market. We also discuss the links between such HJM based models and the intensity models for longevity bonds such as those of Dahl [Dahl M., 2004. Stochastic mortality in life insurance: Market reserves and mortality-linked insurance contracts, Insurance: Math. Econom. 35 (1) 113–136], Biffis [Biffis E., 2005. Affine processes for dynamic mortality and actuarial valuations. Insurance: Math. Econom. 37, 443–468], Luciano and Vigna [Luciano E. and Vigna E., 2005. Non mean reverting affine processes for stochastic mortality. ICER working paper], Schrager [Schrager D.F., 2006. Affine stochastic mortality. Insurance: Math. Econom. 38, 81–97] and Hainaut and Devolder [Hainaut D., Devolder P., 2007. Mortality modelling with Lévy processes. Insurance: Math. Econom. (in press)], and suggest the standard pricing formula of these intensity models could be extended to more general settings.In the second part of this paper, we study the asset allocation problem of pure endowment and annuity portfolios. In order to solve this problem, we study the “risk-minimizing” strategies of such portfolios, when some but not all longevity bonds are available for trading. In this way, we introduce different basis risks.     

Discrete -convex extremal distributions: Theory and applications

Given a nondegenerate moment space with s fixed moments, explicit formulas for the discrete s-convex extremal distribution have been derived for s=1,2,3 (see [M. Denuit, Cl. Lefèvre, Some new classes of stochastic order relations among arithmetic random variables, with applications in actuarial sciences, Insurance Math. Econom. 20 (1997) 197–214]). If s=4, only the maximal distribution is known (see [M. Denuit, Cl. Lefèvre, M. Mesfioui, On s-convex stochastic extrema for arithmetic risks, Insurance Math. Econom. 25 (1999) 143–155]). This work goes beyond this limitation and proposes a method for deriving explicit expressions for general nonnegative integer s. In particular, we derive explicitly the discrete 4-convex minimal distribution. For illustration, we show how this theory allows one to bound the probability of extinction in a Galton–Watson branching process. The results are also applied to derive bounds for the probability of ruin in the compound binomial and Poisson insurance risk models.     

On the dual risk model with tax payments

In this paper, we study the dual risk process in ruin theory (see e.g. Cramér, H. 1955. Collective Risk Theory: A Survey of the Theory from the Point of View of the Theory of Stochastic Processes. Ab Nordiska Bokhandeln, Stockholm, Takacs, L. 1967. Combinatorial methods in the Theory of Stochastic Processes. Wiley, New York and Avanzi, B., Gerber, H.U., Shiu, E.S.W., 2007. Optimal dividends in the dual model. Insurance: Math. Econom. 41, 111–123) in the presence of tax payments according to a loss-carry forward system. For arbitrary inter-innovation time distributions and exponentially distributed innovation sizes, an expression for the ruin probability with tax is obtained in terms of the ruin probability without taxation. Furthermore, expressions for the Laplace transform of the time to ruin and arbitrary moments of discounted tax payments in terms of passage times of the risk process are determined. Under the assumption that the inter-innovation times are (mixtures of) exponentials, explicit expressions are obtained. Finally, we determine the critical surplus level at which it is optimal for the tax authority to start collecting tax payments.     

The private value of public pensions

As individual retirement savings accounts replace public pensions and defined benefit schemes, more retirees will decumulate using commercial income streams rather than public or corporate annuities. Here we use an approximation to the retirement income problem [Huang, H., Milevsky, M.A., Wang, J., 2004. Ruined moments in your life: How good are the approximations? Insurance: Math. Econom. 34, 421–447] to compute the cost of replicating a public real life annuity (the Australian Age Pension) using commercial decumulation products. We treat the public pension as a phased withdrawal plan, matching insurance and payment features, and back out the stochastic present value of the plan under an arbitrarily small ruin probability. To reproduce the pension payment with 99% certainty, a male retiree needs 3.6 times the current average retirement savings account balance, and a female retiree needs more than 10 times the average female account balance. At 95% certainty, required wealth falls by around 25%. We measure separately the impact of gender, investment strategy, retirement age and management fees on this valuation.     

Static super-replicating strategies for a class of exotic options

In this paper, we investigate static super-replicating strategies for European-type call options written on a weighted sum of asset prices. This class of exotic options includes Asian options and basket options among others. We assume that there exists a market where the plain vanilla options on the different assets are traded and hence their prices can be observed in the market. Both the infinite market case (where prices of the plain vanilla options are available for all strikes) and the finite market case (where only a finite number of plain vanilla option prices are observed) are considered. We prove that the finite market case converges to the infinite market case when the number of observed plain vanilla option prices tends to infinity.We show how to construct a portfolio consisting of the plain vanilla options on the different assets, whose pay-off super-replicates the pay-off of the exotic option. As a consequence, the price of the super-replicating portfolio is an upper bound for the price of the exotic option. The super-hedging strategy is model-free in the sense that it is expressed in terms of the observed option prices on the individual assets, which can be e.g. dividend paying stocks with no explicit dividend process known. This paper is a generalization of the work of Simon et al. [Simon, S., Goovaerts, M., Dhaene, J., 2000. An easy computable upper bound for the price of an arithmetic Asian option. Insurance Math. Econom. 26 (2–3), 175–184] who considered this problem for Asian options in the infinite market case. Laurence and Wang [Laurence, P., Wang, T.H., 2004. What’s a basket worth? Risk Mag. 17, 73–77] and Hobson et al. [Hobson, D., Laurence, P., Wang, T.H., 2005. Static-arbitrage upper bounds for the prices of basket options. Quant. Fin. 5 (4), 329–342] considered this problem for basket options, in the infinite as well as in the finite market case.As opposed to Hobson et al. [Hobson, D., Laurence, P., Wang, T.H., 2005. Static-arbitrage upper bounds for the prices of basket options. Quant. Fin. 5 (4), 329–342] who use Lagrange optimization techniques, the proofs in this paper are based on the theory of integral stochastic orders and on the theory of comonotonic risks.     

To split or not to split: Capital allocation with convex risk measures

Convex risk measures were introduced by Deprez and Gerber [Deprez, O., Gerber, H.U., 1985. On convex principles of premium calculation. Insurance: Math. Econom. 4 (3), 179-189]. Here the problem of allocating risk capital to subportfolios is addressed, when convex risk measures are used. The Aumann-Shapley value is proposed as an appropriate allocation mechanism. Distortion-exponential measures are discussed extensively and explicit capital allocation formulas are obtained for the case that the risk measure belongs to this family. Finally the implications of capital allocation with a convex risk measure for the stability of portfolios are discussed. It is demonstrated that using a convex risk measure for capital allocation can produce an incentive for infinite fragmentation of portfolios.     

Bounds and approximations for sums of dependent log-elliptical random variables

Dhaene, Denuit, Goovaerts, Kaas and Vyncke [Dhaene, J., Denuit, M., Goovaerts, M.J., Kaas, R., Vyncke, D., 2002a. The concept of comonotonicity in actuarial science and finance: theory. Insurance Math. Econom. 31 (1), 3-33; Dhaene, J., Denuit, M., Goovaerts, M.J., Kaas, R., Vyncke, D., 2002b. The concept of comonotonicity in actuarial science and finance: Applications. Insurance Math. Econom. 31 (2), 133-161] have studied convex bounds for a sum of dependent random variables and applied these to sums of log-normal random variables. In particular, they have shown how these convex bounds can be used to derive closed-form approximations for several of the risk measures of such a sum. In this paper we investigate to which extent their general results on convex bounds can also be applied to sums of log-elliptical random variables which incorporate sums of log-normals as a special case. Firstly, we show that unlike the log-normal case, for general sums of log-ellipticals the convex lower bound does no longer result in closed-form approximations for the different risk measures. Secondly, we demonstrate how instead the weaker stop-loss order can be used to derive such closed-form approximations. We also present numerical examples to show the accuracy of the proposed approximations.     

Bounds for Asian basket options

In this paper we propose pricing bounds for European-style discrete arithmetic Asian basket options in a Black and Scholes framework. We start from methods used for basket options and Asian options. First, we use the general approach for deriving upper and lower bounds for stop-loss premia of sums of non-independent random variables as in Kaas et al. [Upper and lower bounds for sums of random variables, Insurance Math. Econom. 27 (2000) 151–168] or Dhaene et al. [The concept of comonotonicity in actuarial science and finance: theory, Insurance Math. Econom. 31(1) (2002) 3–33]. We generalize the methods in Deelstra et al. [Pricing of arithmetic basket options by conditioning, Insurance Math. Econom. 34 (2004) 55–57] and Vanmaele et al. [Bounds for the price of discrete sampled arithmetic Asian options, J. Comput. Appl. Math. 185(1) (2006) 51–90]. Afterwards we show how to derive an analytical closed-form expression for a lower bound in the non-comonotonic case. Finally, we derive upper bounds for Asian basket options by applying techniques as in Thompson [Fast narrow bounds on the value of Asian options, Working Paper, University of Cambridge, 1999] and Lord [Partially exact and bounded approximations for arithmetic Asian options, J. Comput. Finance 10 (2) (2006) 1–52]. Numerical results are included and on the basis of our numerical tests, we explain which method we recommend depending on moneyness and time-to-maturity.     

Worst VaR scenarios: A remark

Theorem 15 of Embrechts et al. [Embrechts, Paul, Höing, Andrea, Puccetti, Giovanni, 2005. Worst VaR scenarios. Insurance: Math. Econom. 37, 115-134] proves that comonotonicity gives rise to the on-average-most-adverse Value-at-Risk scenario for a function of dependent risks, when the marginal distributions are known but the dependence structure between the risks is unknown. This note extends this result to the case where, rather than no information, partial information is available on the dependence structure between the risks. A result of Kaas et al. [Kaas, Rob, Dhaene, Jan, Goovaerts, Marc J., 2000. Upper and lower bounds for sums of random variables. Insurance: Math. Econom. 23, 151-168] is also generalized for this purpose.     

Risk-neutral valuation of participating life insurance contracts in a stochastic interest rate environment

Over the last years, the valuation of life insurance contracts using concepts from financial mathematics has become a popular research area for actuaries as well as financial economists. In particular, several methods have been proposed of how to model and price participating policies, which are characterized by an annual interest rate guarantee and some bonus distribution rules. However, despite the long terms of life insurance products, most valuation models allowing for sophisticated bonus distribution rules and the inclusion of frequently offered options assume a simple Black–Scholes setup and, more specifically, deterministic or even constant interest rates.We present a framework in which participating life insurance contracts including predominant kinds of guarantees and options can be valuated and analyzed in a stochastic interest rate environment. In particular, the different option elements can be priced and analyzed separately. We use Monte Carlo and discretization methods to derive the respective values.The sensitivity of the contract and guarantee values with respect to multiple parameters is studied using the bonus distribution schemes as introduced in [Bauer, D., Kiesel, R., Kling, A., Ruß, J., 2006. Risk-neutral valuation of participating life insurance contracts. Insurance: Math. Econom. 39, 171–183]. Surprisingly, even though the value of the contract as a whole is only moderately affected by the stochasticity of the short rate of interest, the value of the different embedded options is altered considerably in comparison to the value under constant interest rates. Furthermore, using a simplified asset portfolio and empirical parameter estimations, we show that the proportion of stock within the insurer’s asset portfolio substantially affects the value of the contract.     

Some results on the CTE-based capital allocation rule

Tasche [Tasche, D., 1999. Risk contributions and performance measurement. Working paper, Technische Universität München] introduces a capital allocation principle where the capital allocated to each risk unit can be expressed in terms of its contribution to the conditional tail expectation (CTE) of the aggregate risk. Panjer [Panjer, H.H., 2002. Measurement of risk, solvency requirements and allocation of capital within financial conglomerates. Institute of Insurance and Pension Research, University of Waterloo, Research Report 01-15] derives a closed-form expression for this allocation rule in the multivariate normal case. Landsman and Valdez [Landsman, Z., Valdez, E., 2002. Tail conditional expectations for elliptical distributions. North American Actuarial J. 7 (4)] generalize Panjer’s result to the class of multivariate elliptical distributions.In this paper we provide an alternative and simpler proof for the CTE-based allocation formula in the elliptical case. Furthermore, we derive accurate and easy computable closed-form approximations for this allocation formula for sums that involve normal and lognormal risks.     

Neural networks approach for determining total claim amounts in insurance

In this study, we present an approach based on neural networks, as an alternative to the ordinary least squares method, to describe the relation between the dependent and independent variables. It has been suggested to construct a model to describe the relation between dependent and independent variables as an alternative to the ordinary least squares method. A new model, which contains the month and number of payments, is proposed based on real data to determine total claim amounts in insurance as an alternative to the model suggested by Rousseeuw et al. (1984) [Rousseeuw, P., Daniels, B., Leroy, A., 1984. Applying robust regression to insurance. Insurance: Math. Econom. 3, 67–72] in view of an insurer.     

Fractional Poisson process

A fractional non-Markov Poisson stochastic process has been developed based on fractional generalization of the Kolmogorov–Feller equation. We have found the probability of n arrivals by time t for fractional stream of events. The fractional Poisson process captures long-memory effect which results in non-exponential waiting time distribution empirically observed in complex systems. In comparison with the standard Poisson process the developed model includes additional parameter μ. At μ=1 the fractional Poisson becomes the standard Poisson and we reproduce the well known results related to the standard Poisson process.As an application of developed fractional stochastic model we have introduced and elaborated fractional compound Poisson process.     

A locally risk-minimizing hedging strategy for unit-linked life insurance contracts in a Lévy process financial market

In [Riesner, M., 2006. Hedging life insurance contracts in a Lévy process financial market. Insurance Math. Econom. 38, 599–608] the (locally) risk-minimizing hedging strategy for unit-linked life insurance contracts is determined in an incomplete financial market driven by a Lévy process. The considered risky asset is not a martingale under the original measure and therefore, a change of measure to the minimal martingale measure is performed.The goal of this paper is to show that the risk-minimizing hedging strategy under the new martingale measure which is found in the paper cited above is not the locally risk-minimizing strategy under the original measure. Finally, the real locally risk-minimizing strategy is derived and a relationship between the number of risky assets held in the proposed portfolio cited in the above-mentioned paper and the one proposed here is given.     

Optimal dividend strategies for a risk process under force of interest

In the classical Cramér–Lundberg model in risk theory the problem of maximizing the expected cumulated discounted dividend payments until ruin is a widely discussed topic. In the most general case within that framework it is proved [Gerber, H.U., 1968. Entscheidungskriterien fuer den zusammengesetzten Poisson-prozess. Schweiz. Aktuarver. Mitt. 1, 185–227; Azcue, P., Muler, N., 2005. Optimal reinsurance and dividend distribution policies in the Cramér–Lundberg model. Math. Finance 15 (2) 261–308; Schmidli, H., 2008. Stochastic Control in Insurance. Springer] that the optimal dividend strategy is of band type. In the present paper we discuss this maximization problem in a generalized setting including a constant force of interest in the risk model. The value function is identified in the set of viscosity solutions of the associated Hamilton–Jacobi–Bellman equation and the optimal dividend strategy in this risk model with interest is derived, which in the general case is again of band type and for exponential claim sizes collapses to a barrier strategy. Finally, an example is constructed for Erlang(2)-claim sizes, in which the bands for the optimal strategy are explicitly calculated.     

Worst allocations of policy limits and deductibles

In the literature, orderings of optimal allocations of policy limits and deductibles were established with respect to a policyholder’s preference. However, from the viewpoint of an insurer, the orderings are not enough for the purpose of pricing. In this paper, by applying the equivalent utility premium principle, we study worst allocations of policy limits and deductibles for an insurer, which give rise to the maximum fair premiums. Closed-form solutions are derived. Then we present a result concerning the optimality in a general risk-sharing scheme, by which we obtain optimal allocations for policyholders directly from worst allocations for an insurer. Several results in Cheung [Cheung, K.C., 2007. Optimal allocation of policy limits and deductibles. Insurance Math. Econom. 41, 382–391] are generalized here.     

On the regulator–insurer interaction in a structural model

In this paper, we provide a new insight to the previous work of Briys and de Varenne [E. Briys, F. de Varenne, Life insurance in a contingent claim framework: Pricing and regulatory implications, Geneva Papers on Risk and Insurance Theory 19 (1) (1994) 53–72], Grosen and Jørgensen [A. Grosen, P.L. Jørgensen, Life insurance liabilities at market value: An analysis of insolvency risk, bonus policy, and regulatory intervention rules in a barrier option framework, Journal of Risk and Insurance 69 (1) (2002) 63–91] and Chen and Suchanecki [A. Chen, M. Suchanecki, Default risk, bankruptcy procedures and the market value of life insurance liabilities, Insurance: Mathematics and Economics 40 (2007) 231–255]. We show that the particular risk management strategy followed by the insurance company can significantly change the risk exposure of the company, and that it should thus be taken into account by regulators. We first study how the regulator establishes regulation intervention levels in order to control for instance the default probability of the insurance company. This part of the analysis is based on a constant volatility. Given that the insurance company is informed of regulatory rules, we study how results can be significantly different when the insurance company follows a risk management strategy with non-constant volatilities. We thus highlight some limits of the prior literature and believe that the risk management strategy of the company should be taken into account in the estimation of the risk exposure as well as in that of the market value of liabilities.     

Optimal investment and reinsurance of an insurer with model uncertainty

We introduce a novel approach to optimal investment–reinsurance problems of an insurance company facing model uncertainty via a game theoretic approach. The insurance company invests in a capital market index whose dynamics follow a geometric Brownian motion. The risk process of the company is governed by either a compound Poisson process or its diffusion approximation. The company can also transfer a certain proportion of the insurance risk to a reinsurance company by purchasing reinsurance. The optimal investment–reinsurance problems with model uncertainty are formulated as two-player, zero-sum, stochastic differential games between the insurance company and the market. We provide verification theorems for the Hamilton–Jacobi–Bellman–Isaacs (HJBI) solutions to the optimal investment–reinsurance problems and derive closed-form solutions to the problems.     

A note on the Swiss Solvency Test risk measure

In this paper we examine whether the Swiss Solvency Test risk measure is a coherent measure of risk as introduced in Artzner et al. [Artzner, P., Delbaen, F., Eber, J.M., Heath, D., 1999. Coherent measures of risk. Math. Finance 9, 203–228; Artzner, P., Delbaen, F., Eber, J.M., Heath, D., Ku, H., 2004. Coherent multiperiod risk adjusted values and Bellman’s principle. Working Paper. ETH Zurich]. We provide a simple example which shows that it does not satisfy the axiom of monotonicity. We then find, as a monotonic alternative, the greatest coherent risk measure which is majorized by the Swiss Solvency Test risk measure.     

Edgeworth expansion for an estimator of the adjustment coefficient

We establish an Edgeworth expansion for an estimator of the adjustment coefficient R, directly related to the geometric-type estimator for general exponential tail coefficients, proposed in [Brito, M., Freitas, A.C.M., 2003. Limiting behaviour of a geometric-type estimator for tail indices. Insurance Math. Econom. 33, 211-226].Using the first term of the expansion, we construct improved confidence bounds for R. The accuracy of the approximation is illustrated using an example from insurance (cf. [Schultze, J., Steinebach, J., 1996. On least squares estimates of an exponential tail coefficient. Statist. Dec. 14, 353-372]).