Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives

Similar documents
arxiv: v1 [q-fin.cp] 1 Aug 2015

Pension Risk Management with Funding and Buyout Options

Pricing Pension Buy-ins and Buy-outs 1

Annuity Decisions with Systematic Longevity Risk. Ralph Stevens

ROBUST HEDGING OF LONGEVITY RISK. Andrew Cairns Heriot-Watt University, and The Maxwell Institute, Edinburgh

Longevity risk: past, present and future

3.4 Copula approach for modeling default dependency. Two aspects of modeling the default times of several obligors

Geographical Diversification of life-insurance companies: evidence and diversification rationale

Economathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t

Robust Longevity Risk Management

Evaluating Hedge Effectiveness for Longevity Annuities

The Impact of Natural Hedging on a Life Insurer s Risk Situation

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model

Enhancing Insurer Value Via Reinsurance Optimization

Exam M Fall 2005 PRELIMINARY ANSWER KEY

Stochastic modelling of electricity markets Pricing Forwards and Swaps

Risk analysis of annuity conversion options with a special focus on decomposing risk

Stochastic Processes and Stochastic Calculus - 9 Complete and Incomplete Market Models

A Cohort-Based Value Index for Longevity Risk Management

COMPARING LIFE INSURER LONGEVITY RISK MANAGEMENT STRATEGIES IN A FIRM VALUE MAXIMIZING FRAMEWORK

M5MF6. Advanced Methods in Derivatives Pricing

Comparison of Pricing Approaches for Longevity Markets

Portability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans

Hedging Credit Derivatives in Intensity Based Models

HEDGING LONGEVITY RISK: A FORENSIC, MODEL-BASED ANALYSIS AND DECOMPOSITION OF BASIS RISK

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the

SYSM 6304 Risk and Decision Analysis Lecture 2: Fitting Distributions to Data

Immunization and Hedging of Longevity Risk

θ(t ) = T f(0, T ) + σ2 T

LIBOR models, multi-curve extensions, and the pricing of callable structured derivatives

1. The force of mortality at age x is given by 10 µ(x) = 103 x, 0 x < 103. Compute E(T(81) 2 ]. a. 7. b. 22. c. 23. d. 20

Market Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing

Modeling and Pricing of Variance Swaps for Local Stochastic Volatilities with Delay and Jumps

1.1 Basic Financial Derivatives: Forward Contracts and Options

Interest rate models and Solvency II

Inflation-indexed Swaps and Swaptions

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

Advances in Valuation Adjustments. Topquants Autumn 2015

May 2012 Course MLC Examination, Problem No. 1 For a 2-year select and ultimate mortality model, you are given:

Structural Models of Credit Risk and Some Applications

Advanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives

The Black-Scholes Model

Near-expiration behavior of implied volatility for exponential Lévy models

HEALTH INSURANCE: ACTUARIAL ASPECTS

Implementing an Agent-Based General Equilibrium Model

MATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS

Longevity Risk Management and the Development of a Value-Based Longevity Index

Pricing and Risk Management of guarantees in unit-linked life insurance

Immunization and Hedging of Longevity Risk

1. For a special whole life insurance on (x), payable at the moment of death:

Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing

Ordinary Mixed Life Insurance and Mortality-Linked Insurance Contracts

Pricing theory of financial derivatives

BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS

Exploring Volatility Derivatives: New Advances in Modelling. Bruno Dupire Bloomberg L.P. NY

DEFERRED ANNUITY CONTRACTS UNDER STOCHASTIC MORTALITY AND INTEREST RATES: PRICING AND MODEL RISK ASSESSMENT

Credit Risk using Time Changed Brownian Motions

Modelling, Estimation and Hedging of Longevity Risk

The stochastic calculus

An Overview of Volatility Derivatives and Recent Developments

induced by the Solvency II project

Decomposition of life insurance liabilities into risk factors theory and application to annuity conversion options

A. 11 B. 15 C. 19 D. 23 E. 27. Solution. Let us write s for the policy year. Then the mortality rate during year s is q 30+s 1.

GLWB Guarantees: Hedge E ciency & Longevity Analysis

Basis Risk and Optimal longevity hedging framework for Insurance Company

Managing Systematic Mortality Risk with Group Self Pooling and Annuitisation Schemes

The Black-Scholes Model

November 2012 Course MLC Examination, Problem No. 1 For two lives, (80) and (90), with independent future lifetimes, you are given: k p 80+k

MORNING SESSION. Date: Friday, May 11, 2007 Time: 8:30 a.m. 11:45 a.m. INSTRUCTIONS TO CANDIDATES

(1) Consider a European call option and a European put option on a nondividend-paying stock. You are given:

Optimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University

Counterparty Credit Risk Simulation

4. Black-Scholes Models and PDEs. Math6911 S08, HM Zhu

A Consistent Pricing Model for Index Options and Volatility Derivatives

Immunization and Hedging of Post Retirement Income Annuity Products

The Black-Scholes Model

Practical example of an Economic Scenario Generator

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

Risk management. VaR and Expected Shortfall. Christian Groll. VaR and Expected Shortfall Risk management Christian Groll 1 / 56

Option Pricing Models for European Options


Hedging with Life and General Insurance Products

Asset Pricing Models with Underlying Time-varying Lévy Processes

Solvency, Capital Allocation and Fair Rate of Return in Insurance

IEOR E4602: Quantitative Risk Management

Application of Stochastic Calculus to Price a Quanto Spread

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

An Analytical Approximation for Pricing VWAP Options

Investigation of Dependency between Short Rate and Transition Rate on Pension Buy-outs. Arık, A. 1 Yolcu-Okur, Y. 2 Uğur Ö. 2

Asset-based Estimates for Default Probabilities for Commercial Banks

Saddlepoint Approximation Methods for Pricing. Financial Options on Discrete Realized Variance

Geographical diversification in annuity portfolios

GRANULARITY ADJUSTMENT FOR DYNAMIC MULTIPLE FACTOR MODELS : SYSTEMATIC VS UNSYSTEMATIC RISKS

ACTUARIAL APPLICATIONS OF THE LINEAR HAZARD TRANSFORM

Chapter 2 ( ) Fall 2012

Two Hours. Mathematical formula books and statistical tables are to be provided THE UNIVERSITY OF MANCHESTER. 22 January :00 16:00

Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.

KØBENHAVNS UNIVERSITET (Blok 2, 2011/2012) Naturvidenskabelig kandidateksamen Continuous time finance (FinKont) TIME ALLOWED : 3 hours

Stat 475 Winter 2018

Interest Rate Caps and Vaulation

Transcription:

Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives Simon Man Chung Fung, Katja Ignatieva and Michael Sherris School of Risk & Actuarial Studies University of New South Wales Australia ASTIN, AFIR/ERM and IACA Colloquia Sydney, 23-27 August, 25 Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 / 26

Introduction Provision of Retirement Income Products Examples of retirement income products: life annuities, deferred annuities, variable annuities with guarantees (financial options) Selling them can be risky - long term nature Duration of payments depend on survival of annuitants Idiosyncratic mortality risk can be handled by diversification Systematic mortality risk (longevity risk) is significant Hedging longevity risk using longevity swaps and caps Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 2 / 26

Outline Pricing and hedging analysis of longevity derivatives on a hypothetical life annuity portfolio subject to longevity risk" Propose and calibrate a two-factor Gaussian mortality model Derive analytical pricing formulas for longevity swaps and caps Investigate hedging features of longevity swaps and caps w.r.t. different assumptions on () the market price of longevity risk (2) term to maturity of the hedging instruments (3) portfolio size Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 3 / 26

Longevity risk modelling How to model survival probability? Intensity-based approach: Cox process with stochastic intensity µ. If µ is known then deal with inhomogeneous Poisson process, where probability of k events in [t, T ] is given by ( T t ) k µ(s) ds µ(s) ds P(N T N t = k F T ) = e T t k! Death time is modelled as the first jump time of a Cox process Given only F t, we set k = in Eq. () and use iterated expectation to obtain the expected survival probability ( S x+t (t, T ) = E P e T ) µ x+s (s)ds t Ft (2) Useful for pricing and the corresponding density function of death time can be obtained () Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 4 / 26

Mortality model The mortality intensity µ x (t) (in full: µ x+t (t)) of a cohort aged x at time t = is modelled by µ x (t) = Y (t) + Y 2 (t), (3) where Y (t) is a general trend that is common to all ages, and Y 2 (t) is an age-specific factor, satisfying the following dynamics where dw (t) dw 2 (t) = ρdt dy (t) = α Y (t) dt + σ dw (t) (4) dy 2 (t) = (αx + β)y 2 (t) dt + σe γx dw 2 (t) (5) Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 5 / 26

Survival probability Proposition Under the two-factor mortality model, the expected (T t)-survival probability of a person aged x + t at time t is given by S x+t (t, T ) = Et (e P T ) µ x (s)ds t = e 2 Γ(t,T ) Θ(t,T ) (6) where Θ(t, T ) and Γ(t, T ) are the mean and the variance of the integral µ x (s) ds respectively. T t We will use the fact that the integral T µ t x (s) ds is normally distributed with known mean and variance to derive analytical pricing formulas for longevity derivatives Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 6 / 26

Parameter estimation Figure: Central death rates m(x, t)..4.3 m (x,t).2. 9 8 Age (x) 7 6 97 98 99 Year (t) 2 2 Australian male ages x = 6,..., 95 and years t = 97,...28 Intensity assumed constant over each integer age and calendar year; approximated by central death rates m(x, t) Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 7 / 26

Parameter estimation Figure: Difference of central death rates m(x, t)..5. m (x,t).5.5. 9 8 Age (x) 7 6 97 98 99 Year (t) 2 2 Sample variance: Var( m x ); Model variance: Var( µ x )=(σ 2 + 2σ σρe γx + σ 2 e 2γx ) t Minimizing 9 x=6,65... (Var( µ x σ, σ, γ, ρ) Var( m x )) 2 with respect to the parameters {σ, σ, γ, ρ}. Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 8 / 26

5 5 2 25 3 35 4 3.5 3 2.5 Age x = 65 5th perc. 25th perc. 5th perc. 75th perc. 95th perc. 4 3.5 3 2.5 Age x = 75 5th perc. 25th perc. 5th perc. 75th perc. 95th perc. µ x (t) 2 µ x (t) 2.5.5.5.5 t 5 5 2 25 3 35 t Survival Probability S x (,T).8.6.4.2 Age x = 65 99% CI Mean Survival Probability S x (,T).8.6.4.2 Age x = 75 99% CI Mean 5 5 2 25 3 35 T 5 5 2 25 3 35 T Other parameters are calibrated to the empirical survival curves aged 65 and 75 in 28 by minimizing T x x=65,75 j= ) 2 (Ŝx (, j) S x (, j) (7) Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 9 / 26

Risk-adjusted measure Assuming, under the risk-adjusted measure Q, we have dy (t) = α Y (t) dt + σ d W (t) (8) dy 2 (t) = (αx + β λσe γx ) Y 2 (t) dt + σe γx d W 2 (t) (9) where λ is the market price of longevity risk; λ is estimated using the proposed price of the BNP/EIB longevity bond (Meyricke & Sherris (24)) Survival Probability S 65 (,T).9.8.7.6.5.4.3 λ = λ = 4.5 λ = 8.5 λ = 2.5.2. 5 5 2 25 3 35 Horizon T Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 / 26

(Index-based) Longevity swaps Consider an annuity provider/hedger who has a future liability that depends on the (stochastic) survival probability of a cohort The provider wants certainty for the estimated survival probability S-forward: At maturity T, pays fixed leg, K(T ) (, ), and receives floating leg - the realized survival probability e T µ x (s) ds. The payoff from the S-forward is, assuming the notional amount is, e T Zero price at inception means that e r T E Q (e T µ x (s) ds K(T ) () ) µ x (s) ds K(T ) = () and hence K(T ) = E Q (e T ) µ x (s) ds (2) A longevity swap is a portfolio of S-forwards Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 / 26

The mark-to-market price process F (t) of an S-forward is F (t) = e r(t t) Et (e Q T ) µ x (s) ds K = e r(t t) Et (e Q t µx (s) ds e T ) µ x (s) ds t K = e r(t t) ( Sx (, t) S x+t (t, T ) K ) (3) S x (, t) := e t µx (s) ds Ft is the realized survival probability, which is observable given F t. Let ˆn denotes the number of survivors in [, t] and the initial population of the cohort is n, then we have S x (, t) ˆn n (4) S x+t (t, T ) is the risk-adjusted survival probability; If an analytical expression for S x+t (t, T ) is available then an explicit formula for F (t) is obtained Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 2 / 26

where S t = S x (, t), S t = S x+t (t, T ), d = Γ(t,T ) ( ln {K/( St St )} + 2 Γ(t, T ) ) and Φ( ) denotes the CDF of the standard normal R.V. Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 3 / 26 (Index-based) Longevity caps A longevity cap is a portfolio of caplets - the payoff of a longevity caplet is max {(e T ) } µ x (s) ds K, (5) Similar to an S-forward but is able to capture the upside potential - when survival probability is overestimated The price of the longevity caplet: Cl(t) = e r(t t) Et ((e Q T )) + µ x (s) ds K Proposition Under the two-factor mortality model, the price at time t of a longevity caplet Cl(t), with maturity T and strike K, is given by ( ) Cl(t) = S t S t e r (T t) Φ Γ(t, T ) d Ke r (T t) Φ ( d) (6)

Idea of proof: since the integral I := T µ t x (s) ds is normally distributed, e I is a log-normal random variable The standard deviation Γ(t, T ) of the integral can be interpreted as the volatility of the (risk-adjusted) aggregated longevity risk for an individual aged x + t at time t, for the period from t to T Figure: Caplet price as a function of (left panel) T and K and (right panel) λ where K =.4 and T = 2..5.5.45.4.4 Cl ( ; T, K).3.2 Caplet Price.35..3 5 2 Time to Maturity (T) 25 3.5.4.3 Strike (K).2..25.2 2 4 6 8 2 4 λ Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 4 / 26

Managing longevity risk in a life annuity portfolio Consider a life annuity portfolio consists of n policyholders aged 65 in year 28; Premium of the annuity ($ per year upon survival) is given by ω x a x = B(, T ) S x (, T ; λ) (7) T = For the annuity provider the present value (P.V.) of asset is A = n a x P.V. of (random) liability for policyholder k: τ k L k = B(, T ) (8) T = P.V. of liability for the whole portfolio: L = n k= L k Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 5 / 26

Interested in the discounted surplus distribution per policy where D no = A L D no n (9) 9 8 7 Discounted Surplus Distribution Per Policy n = 2 n = 4 n = 6 n = 8 6 Density 5 4 3 2.5.5 $ Figure: No longevity risk (σ = σ = ); Idiosyncratic mortality risk is reduced as n increases Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 6 / 26

Swap-hedged annuity portfolio: where F swap = n ˆT T = D swap = A L + F swap (2) B(, T ) (e T ) µ x (s) ds S x (, T ) is the random cash flow from the longevity swap; n here acts as the notional amount; ˆT : term to maturity Cap-hedged annuity portfolio: where F cap = n ˆT T = (2) D cap = A L + F cap C cap (22) B(, T ) max {(e T ) } µ x (s) ds S x (, T ), is the random cash flow from the longevity cap and C cap = n ˆT T = (23) Cl (; T, S x (, T )) (24) is the price of the longevity cap Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 7 / 26

Examining hedge features of longevity swaps and caps with respect to (w.r.t.) different assumptions on () the market price of longevity risk λ (2) term to maturity of the hedging instruments ˆT and (3) the portfolio size n Table: Parameters for the base case. λ ˆT (years) n 8.5 3 4, Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 8 / 26

3.5 3 2.5 λ = no hedge swap cap 3.5 3 2.5 λ = 4.5 no hedge swap cap 2 2.5.5.5.5.5.5.5.5.5.5.5.5 3.5 3 2.5 λ = 8.5 no hedge swap cap 3.5 3 2.5 λ = 2.5 no hedge swap cap 2 2.5.5.5.5.5.5.5.5.5.5.5.5 Figure: Effect of the market price of longevity risk λ on the discounted surplus distribution per policy Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 9 / 26

Table: Hedging features of a longevity swap and a cap w.r.t. market price of longevity risk λ. Mean Std.dev. Skewness VaR.99 ES.99 λ = No hedge -.76.3592 -.284 -.922 -.27 Swap-hedged -.89.78 -.99 -.84 -.223 Cap-hedged -.86.254.855 -.393 -.355 λ = 4.5 No hedge.52.3592 -.284 -.766 -.943 Swap-hedged.48.78 -.99 -.73 -.294 Cap-hedged.682.254.855 -.2425 -.2746 λ = 8.5 No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged.24.78 -.99 -.547 -.938 Cap-hedged.25.254.855 -.93 -.2224 λ = 2.5 No hedge.4475.3592 -.284 -.465 -.6476 Swap-hedged.398.78 -.99 -.354 -.744 Cap-hedged.69.254.855 -.489 -.8 Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 2 / 26

3.5 3 2.5 Term = years no hedge swap cap 3.5 3 2.5 Term = 2 years no hedge swap cap 2 2.5.5.5.5.5.5.5.5.5.5.5.5 3.5 3 2.5 Term = 3 years no hedge swap cap 3.5 3 2.5 Term = 4 years no hedge swap cap 2 2.5.5.5.5.5.5.5.5.5.5.5.5 Figure: Effect of the term to maturity ˆT of the hedging instruments on the discounted surplus distribution per policy Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 2 / 26

Table: Hedging features of a longevity swap and cap w.r.t. term to maturity ˆT. Mean Std.dev. Skewness VaR.99 ES.99 ˆT = Years No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged.282.29 -.387 -.577 -.749 Cap-hedged.2893.2989 -.266 -.58 -.7592 ˆT = 2 Years No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged.74.794 -.757 -.3656 -.56 Cap-hedged.2234.23.26 -.387 -.5259 ˆT = 3 Years No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged.24.78 -.99 -.547 -.938 Cap-hedged.25.254.855 -.93 -.2224 ˆT = 4 Years No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged -.9.668.277 -.66 -.869 Cap-hedged.984.999.527 -.99 -.23 Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 22 / 26

Table: Longevity risk reduction R = Var( D ) of a longevity swap and a cap w.r.t. Var( D) different portfolio size (n). n 2, 4, 6, 8, R swap 92.6% 96.% 97.2% 97.7% R cap 64.9% 67.3% 68.% 68.6% Var( D ) variance of discounted surplus for hedged portfolio Var( D) variance of discounted surplus for unhedged portfolio Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 23 / 26

Table: Hedging features of a longevity swap and a cap w.r.t. portfolio size (n). Mean Std.dev. Skewness VaR.99 ES.99 n = 2, No hedge.2973.3646 -.2662 -.636 -.87 Swap-hedged.2.99 -.65 -.22 -.2653 Cap-hedged.2.26.922 -.2432 -.2944 n = 4, No hedge.2978.3592 -.284 -.648 -.7973 Swap-hedged.24.78 -.99 -.547 -.938 Cap-hedged.25.254.855 -.93 -.2224 n = 6, No hedge.2977.3566 -.2786 -.6363 -.8 Swap-hedged.24.594 -.3346 -.259 -.66 Cap-hedged.24.26.59 -.639 -.25 n = 8, No hedge.2982.3554 -.292 -.66 -.7876 Swap-hedged.29.536 -.556 -.9 -.595 Cap-hedged.29.992.66 -.598 -.99 Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 24 / 26

Summary The proposed two-factor Gaussian mortality model is capable of modelling mortality intensities of different ages simultaneously Longevity swaps and caps can be priced analytically under the model; Standard derivation of the integral T µ t x (s) ds plays the role of volatility in longevity derivatives pricing Swap-hedged annuity portfolio is insensitive to the market price of risk λ Longevity swaps and caps are markedly different as hedging instruments when term to maturity ˆT or portfolio size n is large Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 25 / 26

Thank you for your attention Katja Ignatieva ASTIN, AFIR/ERM and IACA 25 Sydney 25 26 / 26