Longevity and Mortality Models Andrs M. Villegas white School of - - PowerPoint PPT Presentation

longevity and mortality models
SMART_READER_LITE
LIVE PREVIEW

Longevity and Mortality Models Andrs M. Villegas white School of - - PowerPoint PPT Presentation

Longevity and Mortality Models Andrs M. Villegas white School of Risk and Actuarial Studies, CEPAR, UNSW Sydney white white 19 July 2018, UNSW Sydney CEPAR Workshop Longevity and Long-Term Care Risks and Products Agenda Motivation and


slide-1
SLIDE 1

Longevity and Mortality Models

Andrés M. Villegas white

School of Risk and Actuarial Studies, CEPAR, UNSW Sydney white white 19 July 2018, UNSW Sydney CEPAR Workshop Longevity and Long-Term Care Risks and Products

slide-2
SLIDE 2

Agenda

◮ Motivation and Preliminaries ◮ Single population (discrete) stochastic mortality models ◮ Mortality Improvement rate models ◮ Multipopulation mortality models ◮ Continuous time mortality models ◮ Recent developments and outlook

slide-3
SLIDE 3

Motivation and Preliminaries

slide-4
SLIDE 4

Main demographic trends

◮ Expansion over time ◮ Rectangularisation over time ◮ Increasing trend over time in life expectancy ◮ Downward trend over time in death rates

slide-5
SLIDE 5

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-6
SLIDE 6

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-7
SLIDE 7

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-8
SLIDE 8

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-9
SLIDE 9

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-10
SLIDE 10

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-11
SLIDE 11

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-12
SLIDE 12

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-13
SLIDE 13

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-14
SLIDE 14

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-15
SLIDE 15

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-16
SLIDE 16

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-17
SLIDE 17

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-18
SLIDE 18

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-19
SLIDE 19

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-20
SLIDE 20

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-21
SLIDE 21

Male life table distribution of deaths (dx), England and Wales 1850-2009

Source: Human Mortality Database

slide-22
SLIDE 22

Recent trends in mortality and life expectancy

0.00 0.01 0.02 0.03 0.04 1 9 5 − 1 9 5 5 1 9 5 5 − 1 9 6 1 9 6 − 1 9 6 5 1 9 6 5 − 1 9 7 1 9 7 − 1 9 7 5 1 9 7 5 − 1 9 8 1 9 8 − 1 9 8 5 1 9 8 5 − 1 9 9 1 9 9 − 1 9 9 5 1 9 9 5 − 2 2 − 2 5 2 5 − 2 1 2 1 − 2 1 5

mortality rate

Mortality rate at age 60−64 (Females)

40 50 60 70 80 1 9 5 − 1 9 5 5 1 9 5 5 − 1 9 6 1 9 6 − 1 9 6 5 1 9 6 5 − 1 9 7 1 9 7 − 1 9 7 5 1 9 7 5 − 1 9 8 1 9 8 − 1 9 8 5 1 9 8 5 − 1 9 9 1 9 9 − 1 9 9 5 1 9 9 5 − 2 2 − 2 5 2 5 − 2 1 2 1 − 2 1 5

Life Expectancy (years)

Period Life Expectancy at brith (Females)

Region

AFRICA ASIA Australia EUROPE LATIN AMERICA AND THE CARIBBEAN NORTHERN AMERICA OCEANIA WORLD

United Nations World Population Prospects 2017

◮ Good-news! ◮ Important social and financial implications for governments,

insurers, individuals.

◮ Need to model and project these trends

slide-23
SLIDE 23

IMF assessment of global of three year longevity shock longevity

Advanced economies Emerging economies Source: IMF (2012)

slide-24
SLIDE 24

Mortality forecasting methodologies

A good overview of methodologies is given in the review papers Booth and Tickle (2008), Wong-Fupuy and Haberman (2004), Pitacco (2004) and in book Pitacco et al. (2009)

◮ Expert based ◮ Explanatory

◮ Structural Modelling (Explanatory or Econometric). ◮ Cause of death decomposition

◮ Extrapolation

◮ Trend modelling

slide-25
SLIDE 25

A timeline of “recent” mortality modelling methodologies

slide-26
SLIDE 26

Single population (discrete) stochastic mortality models

Discussion based on: Villegas, A. M., Millossovich, P., & Kaishev, V. K. (2018). StMoMo: An R Package for Stochastic Mortality Modelling. JSS Journal of Statistical Software, 84(3). https://doi.org/10.18637/jss.v084.i03

slide-27
SLIDE 27

Advances in single population mortality modelling

◮ Lee-Carter model (Lee and Carter, 1992)

◮ Add more bilinear age-period components (Renshaw and

Haberman, 2003)

◮ Add a cohort effect (Renshaw and Haberman, 2006)

◮ Two factor CBD model (Cairns et al., 2006)

◮ Add cohort effect, quadratic age term (Cairns et al., 2009) ◮ Combine with features of the Lee-Carter (Plat, 2009b)

◮ Many more models proposed in the literature (e.g. Aro and

Pennanen (2011), O’Hare and Li (2012), Börger et al. (2013), Alai and Sherris (2014))

slide-28
SLIDE 28

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1921)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-29
SLIDE 29

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1925)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-30
SLIDE 30

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1930)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-31
SLIDE 31

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1935)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-32
SLIDE 32

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1940)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-33
SLIDE 33

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1945)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-34
SLIDE 34

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1950)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-35
SLIDE 35

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1955)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-36
SLIDE 36

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1960)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-37
SLIDE 37

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1965)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-38
SLIDE 38

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1970)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-39
SLIDE 39

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1975)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-40
SLIDE 40

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1980)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-41
SLIDE 41

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1985)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-42
SLIDE 42

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1990)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-43
SLIDE 43

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1995)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-44
SLIDE 44

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (2000)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-45
SLIDE 45

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (2005)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-46
SLIDE 46

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (2010)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-47
SLIDE 47

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (2014)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-48
SLIDE 48

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1921−2014)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-49
SLIDE 49

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1921−2014)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-50
SLIDE 50

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1921−2014)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-51
SLIDE 51

Lee-Carter model

20 40 60 80 100 −10 −8 −6 −4 −2

Australia: male death rates (1921−2014)

age log death rates

log µxt = αx +

N

  • i=1

β(i)

x κt(i)

+ β(0)

x γt−x

slide-52
SLIDE 52

Lee-Carter model

20 40 60 80 100 −8 −6 −4 −2

αx vs. x

age 20 40 60 80 100 0.000 0.010

βx

(1) vs. x age 1970 1980 1990 2000 2010 −40 20

κt

(1) vs. t year

slide-53
SLIDE 53

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1961)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-54
SLIDE 54

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1962)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-55
SLIDE 55

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1963)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-56
SLIDE 56

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1964)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-57
SLIDE 57

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1965)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-58
SLIDE 58

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1970)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-59
SLIDE 59

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1980)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-60
SLIDE 60

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (1990)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-61
SLIDE 61

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (2000)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-62
SLIDE 62

Cairns-Blake-Dowd model

logit qxt = κ(1)

t +(x−¯

x)κ(2)

t

50 60 70 80 90 100 −6 −5 −4 −3 −2 −1

EW: death probability male (2010)

age log qx 1960 1970 1980 1990 2000 2010 −3.2 −2.8 −2.4

κt (1)vs t

year 1960 1970 1980 1990 2000 2010 0.095 0.100 0.105

κt (2)vs t

year

slide-63
SLIDE 63

Cohort effects

◮ Statistical evidence both period and

cohort effect have an impact on mortality improvements (Willets 1999, 2004)

◮ Period effects approximate

contemporary factors

◮ General health status of the

population

◮ Healthcare services available ◮ Critical weather conditions

◮ Cohort effects approximate historical

factors

◮ World War II ◮ Diet ◮ Welfare State (in the UK) ◮ Smoking habits

slide-64
SLIDE 64

Age-Period-Cohort models: Classic APC model

◮ A widely used structured used in medicine, psychology and

demography (Hobcraft et al. (1982), Wilmoth (1990)) log µxt = αx + κt + γt−x

◮ No unique set of parameters resulting in optimal fit due to

c = t − x (αx, κt, γt−x) → (αx + φ1 − φ2x, κt + φ2t, γt−x − φ1 − φ2(t − x)) (αx, κt, γt−x) → (αx + c1, κt − c1, γt−x)

◮ Impose constraints

  • t

κt = 0,

  • c

γc = 0,

  • c

cγc = 0

slide-65
SLIDE 65

Age-Period-Cohort models: Identifiability

20 40 60 80 −10 −8 −6 −4 −2 2 αx vs. x age (x)

tilt 1 tilt 2 tilt 3 tilt 4 tilt 5

20 40 60 80 0.005 0.01 0.015 0.02 β(1)

x vs. x

age (x) 20 40 60 80 0.005 0.01 0.015 0.02 β(0)

x vs. x

age (x) 1970 1980 1990 2000 −200 −100 100 κt vs. t calendar year (t) 1880 1900 1920 1940 1960 1980 2000 −400 −200 200 400 ιy vs. y year of birth (y)

slide-66
SLIDE 66

Other stochastic mortality models

◮ Lee-Carter extensions

◮ Add cohorts

log µxt = αx + β(1)

x κ(1) t

+ β(0)

x γt−x

◮ CBD extensions

◮ M6: Add cohorts

logit qxt = ηxt = κ(1)

t

+ (x − ¯ x)κ(2)

t

+ γt−x

◮ M7: Add cohorts and quadratic age effect

logit qxt = ηxt = κ(1)

t

+(x −¯ x)κ(2)

t

+

  • (x − ¯

x)2 − ˆ σ2

x

  • κ(3)

t

+γt−x

◮ Plat model combines the Lee-Carter and the CBD

log µxt = αx + κ(1)

t

+ (¯ x − x)κ(2)

t

+ (¯ x − x)+κ(3)

t

+ γt−x

slide-67
SLIDE 67

Generalised Age-Period-Cohort stochastic mortality models

Recent research has proposed a unifying framework discrete stochastic mortality models

◮ General Age-Period-Cohort model structure Hunt and Blake

(2015a)

◮ Generalised (non-)linear model Currie (2016) ◮ R Implementation of GAPC models Villegas et al. (2018)

slide-68
SLIDE 68

Generalised Age-Period-Cohort stochastic mortality models

  • 1. Random Component:

Dxt ∼ Poisson(E c

xtµxt)

  • r

Dxt ∼ Binomial(Ext, qxt)

  • 2. Systematic Component:

ηxt = αx +

N

  • i=1

β(i)

x κ(i) t

+ β(0)

x γt−x

◮ Lee-Carter type β(i)

x , non-parametric

◮ CBD type β(i)

x

≡ f (i)(x), pre-specified parametric function

  • 3. Link Function:

g

  • E

Dxt

Ext

  • = ηxt

◮ log-Poisson: ηxt = log µxt ◮ logit-Binomial: ηxt = logit qxt

slide-69
SLIDE 69

Generalised Age-Period-Cohort stochastic mortality models

  • 4. Set of parameter constraints:

◮ Need parameters constraints to ensure identifiability

  • 5. Forecasting and simulation

◮ Period indexes: Multivariate random walk with drift

κt = δ + κt−1 + ξκ

t ,

κt =    κ(1)

t

. . . κ(N)

t

   , ξκ

t ∼ N(0, Σ),

◮ Cohort effect: ARIMA(p, q, d) with drift

∆dγc = δ0+φ1∆dγc−1+· · ·+φp∆dγc−p+ǫc+δ1ǫc−1+· · ·+δqǫc−q

GAPC models can be implemented with the R package StMoMo (http://cran.r-project.org/package=StMoMo)

slide-70
SLIDE 70

Other key areas in single population discrete mortality models I

◮ Parameter estimation:

◮ Poisson (Brouhns et al., 2002) ◮ Negative-Binomial (Delwarde et al., 2007a, Li et al. (2009)) ◮ Bayesian and state-space setting (Czado et al., 2005, Pedroza

(2006), Kogure et al. (2009), Fung et al. (2016))

◮ Parameter Smoothing (Delwarde et al., 2007b) and functional

data approach (Hyndman and Ullah, 2007)

◮ Bootstrapping and parameter uncertainty

◮ Semiparametric (Brouhns et al., 2005) ◮ Parametric (Koissi et al., 2006) ◮ Comparison of methods (Renshaw and Haberman, 2008)

◮ Modelling of errors: residual dependence (Debón, 2008, Debón

et al. (2010), Mavros et al. (2017))

◮ Identifiability

slide-71
SLIDE 71

Other key areas in single population discrete mortality models II

◮ Age-Period Models (Hunt and Blake, 2015b) ◮ APC models (Hunt and Blake, 2015c) ◮ Impact on estimation (Hunt and Villegas, 2015)

◮ Modelling and projecting of period and cohort factors

◮ Optimal calibration period (Booth et al., 2002, Denuit2005) ◮ Regime-switching (Milidonis et al., 2011) ◮ Structural changes (Coelho and Nunes, 2011, van Berkum et al.

(2014))

◮ Selection criteria of most appropriate model

◮ Goodness-of-fit (Cairns et al., 2009, Dowd et al. (2010b)) ◮ Backtesting (Dowd et al., 2010a) ◮ Qualitative properties of forecasts (Cairns et al., 2011b) ◮ Overall performance for England and Wales and the USA

(Haberman and Renshaw, 2011).

slide-72
SLIDE 72

Mortality improvement rate models

Discussion based on: Hunt, A., & Villegas, A. M. (2017). Mortality Improvement Rates: Modeling and Parameter Uncertainty. In Living to 100, Society of Actuaries International

  • Symposium. Retrieved from https://www.soa.org/essays-monographs/

2017-living-to-100/2017-living-100-monograph-hunt-villegas-paper.pdf

slide-73
SLIDE 73

Two parallel approaches to modelling and projecting mortality

Mortality Rates Improvement Rates What? x qx,t, µx,t, mx,t 1 −

qx,t qx,t−1 , − ln( mx,t mx,t−1 )

Who?Lee and Carter (1992) xxxxxx Cairns et al. (2006) xxxxxx xx Brouhns et al. (2005) xxxxxxxxxxxxxxx CMIB (1978) xxxxxx CMI (2002, 2009, 2016) SOA (1995, 2012) How? ln mx,t = αx +

N

  • i=1

β(i)

x κ(i) t

+ γt−x qx,T+n = qx,T(1 − Rx)n

slide-74
SLIDE 74

Improvement rates in actuarial practice – AA and BB scales

qx,T+n = qx,T

  • Base table

× (1 − AAx)n

  • Reduction factor
slide-75
SLIDE 75

Improvement rates in actuarial practice – CMI model

Latest CMI projection model uses an APC model on improvement rates (CMI Working paper 90): −∆ ln mx,t = αx + κt + γt−x

slide-76
SLIDE 76

Recent academic interest on improvement rate modelling

◮ Mitchell et al. (2013)

◮ ln

ˆ mx,t+1 ˆ mx,t

  • = αx + βxκt + ǫx,t

◮ Benefits of detrending ◮ Estimation with singular value decomposition

◮ Haberman and Renshaw (2012)

ˆ mx,t−1 − ˆ mx,t 0.5( ˆ mx,t−1 + ˆ mx,t) = ηx,t

◮ Predictor structures ηx,t borrowed from mortality rate modelling ◮ Duality between mortality rate and improvement rate modelling ◮ Estimation using Gaussian model with variable dispersion

◮ Haberman and Renshaw (2013), Plat (2011), Danesi et al.

(2015), Chuliá et al. (2016), Njenga and Sherris (2011)

slide-77
SLIDE 77

Complications with improvement rate modelling

1960 1970 1980 1990 2000 2010 0.02 0.04 0.06

Mortality rate

year Mortality rate at age 70 1960 1970 1980 1990 2000 2010 −0.10 0.00 0.10

Improvement rate

year Improvement rate at age 70 20 40 60 80 100 1e−04 5e−03 5e−01 age Mortality rate in 2011 (log) 20 40 60 80 100 −0.2 0.0 0.2 age Improvement rate in 2011

◮ Patterns are not that clear ◮ Non-standard distribution ◮ Heteroscedasticity ◮ Parameter uncertainty

slide-78
SLIDE 78

Two alterantive approaches to modelling improvement rates

Diagram source: Li et al. (2017)

◮ Route A (Direct): Direct modelling of improvement rates ◮ Route B (Indirect)) Derive improvement rates from mortality rate model

slide-79
SLIDE 79

Route B: Estimation and equivalent mortality rate structure

ln(mx,t) = Ax − αxt +

N

  • i=1

β(i)

x K (i) t

+ Γt−x

1960 1970 1980 1990 2000 2010 −0.8 −0.4 0.0

Kt (1) vs. t

year 1880 1900 1920 1940 0.00 0.15 0.30

Γt−x vs. t−x

cohort

κ(i)

t

= −∆K (i)

t

⇐ ⇒

xγc = −∆Γc −∆ ln mx,t = αx +

N

  • i=1

β(i)

x κ(i) t

+ γt−x

1970 1980 1990 2000 2010 −0.05 0.05

κt (1) vs. t

year 1880 1900 1920 1940 −0.10 0.00

γt−x vs. t−x

cohort

slide-80
SLIDE 80

Route B: Estimation and equivalent mortality rate structure

ln(mx,t) = Ax − αxt +

N

  • i=1

β(i)

x K (i) t

+ Γt−x

1960 1970 1980 1990 2000 2010 −0.8 −0.4 0.0

Kt (1) vs. t

year 1880 1900 1920 1940 0.00 0.15 0.30

Γt−x vs. t−x

cohort

κ(i)

t

= −∆K (i)

t

⇐ ⇒

xγc = −∆Γc −∆ ln mx,t = αx +

N

  • i=1

β(i)

x κ(i) t

+ γt−x

1970 1980 1990 2000 2010 −0.05 0.05

κt (1) vs. t

year 1880 1900 1920 1940 −0.10 0.00

γt−x vs. t−x

cohort

slide-81
SLIDE 81

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

slide-82
SLIDE 82

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t x

slide-83
SLIDE 83

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t x

slide-84
SLIDE 84

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2

slide-85
SLIDE 85

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-86
SLIDE 86

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-87
SLIDE 87

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-88
SLIDE 88

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-89
SLIDE 89

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-90
SLIDE 90

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-91
SLIDE 91

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-92
SLIDE 92

Estimation routes illustration: ηx,t = αx

1960 1970 1980 1990 2000 2010 −6.6 −6.4 −6.2 −6.0

ln m(x,t) vs. t

year Mortality rate at age 40 (log)

Fitted Crude 0.8% 1.0% 1.2% 1.4%

Improvement rate

Estimated αx

Route 2: “Indirect” ln mx,t = Ax − αxt + ǫx,t αx =

T

  • t=0

(t − ¯ t)(ln ˆ mx,t − ln ˆ mx,t)

T

  • t=0

(t − ¯ t)2 Route 1: “Direct” αx = 1 T

T

  • t=1

−∆ ln ˆ mx,t αx = ln ˆ mx,0 − ln ˆ mx,T T x

slide-93
SLIDE 93

Parameter estimates: LC and CBD

Black lines: ‘’Indirect” approach, Red lines: ”Direct” approach LC

20 30 40 50 60 70 80 90 −0.5 0.5 1.0 1.5 2.0

βx (1) vs. x

age 1970 1980 1990 2000 2010 −0.02 0.02 0.04 0.06

κt (1) vs. t

year

CBD

1970 1980 1990 2000 2010 −0.02 0.00 0.02 0.04

κt (1) vs. t

year 1970 1980 1990 2000 2010 −0.002 0.000 0.002

κt (2) vs. t

year

slide-94
SLIDE 94

Key message

◮ Improvement rates are an intuitive and natural way to

interpret mortality data

◮ Compelling reasons for formulating and communicating

projection models in terms of improvement rates

◮ Important differences between approached to fitting

improvement rate models

◮ Considerable parameter uncertainty ◮ Implication for projections and robustness

◮ Compelling reasons for estimating projection models in terms

  • f mortality rates
slide-95
SLIDE 95

Multipopulation (discrete) mortality models

Discussion based on: Villegas, A. M., Haberman, S., Kaishev, V. K., & Millossovich, P. (2017). A Comparative Study of Two-Population Models for the Assessment of Basis Risk in Longevity Hedges. ASTIN Bulletin, 47(03), 631–679. https://doi.org/10.1017/asb.2017.18

slide-96
SLIDE 96

Universe of Multipopulation Models

Extensions of the Lee-Carter

Joint-κ log mi

xt = αi x + βi xκt

Carter and Lee (1992), Li and Hardy (2011), Wilmoth and Valkonen (2001), Delwarde et al. (2006) Three-way Lee-Carter log mi

xt = αi x + βxλiκt

Russolillo et al. (2011) Common Factor log mi

xt = αi x + βxκt

Carter and Lee (1992), Li and Lee (2005), Li and Hardy (2011) Stratified Lee-Carter log mi

xt = αx + αi + βxκt

Butt and Haberman (2009), Debón et al. (2011) Augmented Common Factor log mi

xt = αi x + βxκt + β(i) x κi t

Li and Lee (2005),Li and Hardy (2011) Hyndman et al. (2013), Li (2012) Augmented Common Factor + Cohorts log mi

xt = αi x + βxκt+

N

j=1 β(j,i) x

κ(j,i)

t

+ β(0,i)

x

γi

t−x

Yang et al. (2016) Relative Lee-Carter + Cohorts log mi

xt = αx + β(1) x κt + γt−x+

αi

x + β(2) x κi t

Villegas and Haberman (2014) Co-integrated Lee-Carter log mi

xt = αi x + βi xκi t

Carter and Lee (1992), Li and Hardy (2011), Yang and Wang (2013) Lee-Carter + VAR/VECM log mi

xt = αi x + βxκi t

Zhou et al. (2014) Common Age Effect log mi

xt = αi x + j βj xκ(j,i) t

Kleinow (2015) Bayesian two-population APC log mi

xt = αi x + κi t + γi t−x

Cairns et al. (2011a) Gravity model - Two-population APC log mi

xt = αi x + κi t + γi t−x

Dowd et al. (2011)

Extensions of the CBD model

Two-population M7 logit qi

xt = κ(i,1) t

+ (x − ¯ x)κ(i,2)

t

+ (x − ¯ x)2 − ˆ σ2

x

κ(i,3)

t

+ γi

t−x

Li et al. (2015) Two-population M6 logit qi

xt = κ(i,1) t

+ (x − ¯ x)κ(i,2)

t

+ γi

t−x

Li et al. (2015) Two-population CBD (M5) logit qi

xt = κ(i,1) t

+ (x − ¯ x)κ(i,2)

t

Li et al. (2015)

Other models

Plat Relative model Plat (2009a) Saint model Jarner and Kryger (2011a) Plat + Lee-Carter Wan and Bertschi (2015) Multipopulation GLM Hatzopoulos and Haberman (2013), Ahmadi and Li (2014) Relative P-Splines Biatat and Currie (2010)

slide-97
SLIDE 97

Coherent forecasts for multiple populations

◮ See for example Li and Lee (2005), Hyndman et al. (2013)

Source: Hyndman et al. (2013)

slide-98
SLIDE 98

Quantifying socio-economic differences in mortality

◮ See for example Villegas and Haberman (2014), Cairns et al.

(2016)

Source: Villegas and Haberman (2014)

log mi

xt = αx + αi x + βxλiκt

slide-99
SLIDE 99

Assessing longevity basis risk

◮ See for example Li and Hardy (2011), Villegas et al. (2017), Li

et al. (2018)

Source: Villegas et al. (2017)

slide-100
SLIDE 100

Borrowing strength from a larger population

◮ Derive the trend from the larger population and model the

spread (ratio) between the larger and the smaller population (Jarner and Kryger, 2011b; van Berkum et al., 2017)

Source: van Berkum et al. (2017)

slide-101
SLIDE 101

Other developments in multipopulation modelling

◮ Ensuring consistency between subpopulation and total

population projections (Shang and Hyndman, 2017; Shang and Haberman, 2017)

◮ Modelling of period effects in multipopulation models (Zhou

et al., 2014; Li et al., 2015)

◮ Clustering of mortality trends in multiple-populations (Debón

et al., 2017)

slide-102
SLIDE 102

Recent developments and outlook

slide-103
SLIDE 103

Affine Continuos time mortality models

◮ Satisfy important requirements for financial applications ◮ Continuous time mortality modelling framework for insurance

application Dahl (2004); Biffis (2005)

◮ Affine mortality models

◮ m factor-model Schrager (2006) ◮ Consistent 3-factor model Blackburn and Sherris (2013) ◮ Affine processes and multi-cohort factors Xu et al. (2018) ◮ two population multiple cohorts (Sherris et al., 2018)

slide-104
SLIDE 104

Recent developments

◮ Applications of statistical machine learning techniques to

mortality modelling

◮ Sparse Vector Autoregression (Li and Lu, 2017) ◮ High Dimensional VAR with elastic nets (Guibert et al., 2017) ◮ Gaussian processes (Ludkovski et al., 2016) ◮ Neural networks (Hainaut, 2018) ◮ Random fields (Doukhan et al., 2017)

◮ Scope for applying more of these techniques

◮ Incorporation of additional information (other populations,

macro-economic data, etc)

◮ Evaluation and selection of models

slide-105
SLIDE 105

Outlook

◮ Mortality modelling remains an issue of current interest ◮ Saturation in the traditional Lee-Carter style-approaches ◮ But many topics remain relevant

◮ Modelling of populations of small populations or with scarce

data

◮ Mortality trends at older ages ◮ Understanding of trends in causes of death

◮ Modelling of dependence between ages, cohorts and

populations

◮ Quantification of trends differences between

sociology-economic groups

◮ Incorporation of individual level data ◮ Integration between financial models and mortality models

(continuous time approaches)

slide-106
SLIDE 106

References I

Seyed Saeed Ahmadi and Johnny Siu-Hang Li. Coherent mortality forecasting with generalized linear models: A modified time-transformation approach. Insurance: Mathematics and Economics, 59:194–221, 2014. ISSN 01676687. doi: 10.1016/j.insmatheco.2014.09.007. Daniel H. Alai and Michael Sherris. Rethinking age-period-cohort mortality trend models. Scandinavian Actuarial Journal, (3): 208–227, 2014. Helena Aro and Teemu Pennanen. A user-friendly approach to stochastic mortality modelling. European Actuarial Journal, 1: 151–167, 2011. V.D. Biatat and Iain D Currie. Joint models for classification and comparison of mortality in different countries. In Proceedings of 25rd International Workshop on Statistical Modelling, Glasgow, pages 89–94, 2010.

slide-107
SLIDE 107

References II

Enrico Biffis. Affine processes for dynamic mortality and actuarial

  • valuations. Insurance: Mathematics and Economics, 37(3):

443–468, 2005. ISSN 01676687. doi: 10.1016/j.insmatheco.2005.05.003. Craig Blackburn and Michael Sherris. Consistent dynamic affine mortality models for longevity risk applications. Insurance: Mathematics and Economics, 53(1):64–73, 2013. ISSN 01676687. doi: 10.1016/j.insmatheco.2013.04.007. URL http://linkinghub.elsevier.com/retrieve/pii/S0167668713000632. Heather Booth and Leonie Tickle. Mortality modelling and forecasting: A review of methods. Annals of Actuarial Science, 1 (2):3–43, 2008. Heather Booth, John Maindonald, and Len Smith. Applying Lee-Carter under conditions of variable mortality decline. Population Studies, 56(3):325–36, 2002.

slide-108
SLIDE 108

References III

Matthias Börger, Daniel Fleischer, and Nikita Kuksin. Modeling the mortality trend under modern solvency regimes. ASTIN Bulletin, 44(1):1–38, 2013. Natacha Brouhns, Michel Denuit, and J.K. Vermunt. A Poisson log-bilinear regression approach to the construction of projected

  • lifetables. Insurance: Mathematics and Economics, 31(3):

373–393, 2002. Natacha Brouhns, Michel Denuit, and Ingrid Van Keilegom. Bootstrapping the Poisson log-bilinear model for mortality

  • forecasting. Scandinavian Actuarial Journal, (3):212–224, 2005.

Zoltan Butt and Steven Haberman. Ilc: A collection of R functions for fitting a class of Lee-Carter mortality models using iterative fitting algorithms. Actuarial Research Paper, Cass Business School, 2009.

slide-109
SLIDE 109

References IV

Andrew J. G. Cairns, David Blake, and Kevin Dowd. A two-factor model for stochastic mortality with parameter uncertainty: theory and calibration. Journal of Risk and Insurance, 73(4):687–718, 2006. Andrew J. G. Cairns, David Blake, Kevin Dowd, Guy D. Coughlan,

  • D. Epstein, A. Ong, and I. Balevich. A quantitative comparison of

stochastic mortality models using data from England and Wales and the United States. North American Actuarial Journal, 13(1): 1–35, 2009. Andrew J. G. Cairns, David Blake, Kevin Dowd, and Guy D.

  • Coughlan. Bayesian stochastic mortality modelling for two
  • populations. ASTIN Bulletin, 41:29–59, 2011a.

Andrew J. G. Cairns, David Blake, Kevin Dowd, Guy D. Coughlan, David Epstein, and Marwa Khalaf-Allah. Mortality density forecasts: An analysis of six stochastic mortality models. Insurance: Mathematics and Economics, 48(3):355–367, 2011b.

slide-110
SLIDE 110

References V

Lawrence R. Carter and Ronald D. Lee. Modeling and forecasting US sex differentials in mortality. International Journal of Forecasting, 8(3):393–411, 1992. Helena Chuliá, Montserrat Guillén, and Jorge M. Uribe. Modeling Longevity Risk With Generalized Dynamic Factor Models and Vine-Copulae. ASTIN Bulletin, 46(01):165–190, 2016. ISSN 0515-0361. doi: 10.1017/asb.2015.21. URL http: //www.journals.cambridge.org/abstract{_}S0515036115000215.

  • CMIB. Report n3. Continuous Mortality Investigation Bureau.

Instute of Actuaries and Faculty of Actuaries, 1978. Edviges Coelho and Luis C. Nunes. Forecasting mortality in the event of a structural change. Journal of the Royal Statistical Society: Series A (Statistics in Society), 174(3):713–736, jul 2011. ISSN 09641998. doi: 10.1111/j.1467-985X.2010.00687.x. URL http://doi.wiley.com/10.1111/j.1467-985X.2010.00687.x. Iain D Currie. On fitting generalized linear and non-linear models of

  • mortality. Scandinavian Actuarial Journal, (4):356–383, 2016.
slide-111
SLIDE 111

References VI

Claudia Czado, Antoine Delwarde, and Michel Denuit. Bayesian Poisson log-bilinear mortality projections. 36:260–284, 2005. Mikkel Dahl. Stochastic mortality in life insurance: Market reserves and mortality-linked insurance contracts. Insurance: Mathematics and Economics, 35(1):113–136, 2004. ISSN 01676687. doi: 10.1016/j.insmatheco.2004.05.003. Ivan Luciano Danesi, Steven Haberman, Pietro Millossovich, Ivan Luciano, Steven Haberman, and Pietro Millossovich. Forecasting mortality in subpopulations using Lee-Carter type models: a

  • comparison. Insurance: Mathematics and Economics, 62:151–161,
  • 2015. ISSN 0167-6687. doi: 10.1016/j.insmatheco.2015.03.010.

URL http://dx.doi.org/10.1016/j.insmatheco.2015.03.010. Ana Debón. Modelling residuals dependence in dynamic life tables: A geostatistical approach. Computational Statistics & Data Analysis, 52(6):3128–3147, feb 2008. ISSN 01679473. doi: 10.1016/j.csda.2007.08.006. URL http://linkinghub.elsevier.com/retrieve/pii/S0167947307003040.

slide-112
SLIDE 112

References VII

Ana Debón, F Martínez-Ruiz, and F Montes. A geostatistical approach for dynamic life tables: The effect of mortality on remaining lifetime and annuities. Insurance: Mathematics and Economics, 47(3):327–336, 2010. Ana Debón, Francisco Montes, and Francisco Martínez-Ruiz. Statistical methods to compare mortality for a group with non-divergent populations: an application to Spanish regions. European Actuarial Journal, 1(2):291–308, 2011. ISSN 2190-9733. doi: 10.1007/s13385-011-0043-z. Ana Debón, L. Chaves, S. Haberman, and F. Villa. Characterization

  • f between-group inequality of longevity in European Union
  • countries. Insurance: Mathematics and Economics, 2017. ISSN
  • 01676687. doi: 10.1016/j.insmatheco.2017.05.005. URL

http://linkinghub.elsevier.com/retrieve/pii/S0167668716300920.

slide-113
SLIDE 113

References VIII

Antoine Delwarde, Michel Denuit, Montserrat Guillén, and

  • A. Vidiella-i Anguera. Application of the Poisson log-bilinear

projection model to the G5 mortality experience. Belgian Actuarial Bulletin, 6(1):54–68, 2006. Antoine Delwarde, Michel Denuit, and Partrat Christian. Negative binomial version of the Lee–Carter model for mortality forecasting. Applied Stochastic Models in Business and Industry, 23(5): 385–401, 2007a. ISSN 1524-1904. doi: 10.1002/asmb. Antoine Delwarde, Michel Denuit, and Paul Eilers. Smoothing the Lee-Carter and Poisson log-bilinear models for mortality forecasting: a penalized log-likelihood approach. Statistical Modelling, 7(1):29–48, 2007b. ISSN 1471-082X. doi: 10.1177/1471082X0600700103.

slide-114
SLIDE 114

References IX

  • P. Doukhan, D. Pommeret, J. Rynkiewicz, and Y. Salhi. A class of

random field memory models for mortality forecasting. Insurance: Mathematics and Economics, 77:97–110, 2017. ISSN 01676687. doi: 10.1016/j.insmatheco.2017.08.010. URL http://linkinghub.elsevier.com/retrieve/pii/S0167668716303857. Kevin Dowd, Andrew J. G. Cairns, David Blake, Guy D. Coughlan, David Epstein, and Marwa Khalaf-Allah. Backtesting stochastic mortality models: An ex-post evaluation of multi-period-ahead density forecasts. North American Actuarial Journal, 14(3): 281–298, 2010a. Kevin Dowd, Andrew J. G. Cairns, David Blake, Guy D. Coughlan, David Epstein, and Marwa Khalaf-Allah. Evaluating the goodness

  • f fit of stochastic mortality models. Insurance: Mathematics and

Economics, 47(3):255–265, dec 2010b. ISSN 01676687. doi: 10.1016/j.insmatheco.2010.06.006. URL http://linkinghub.elsevier.com/retrieve/pii/S0167668710000752.

slide-115
SLIDE 115

References X

Kevin Dowd, Andrew J. G. Cairns, David Blake, Guy D. Coughlan, and Marwa Khalaf-allah. A gravity model of mortality rates for two related populations. North American Actuarial Journal, 15(2): 334–356, 2011. Man Chung Fung, Gareth W. Peters, and Pavel V. Shevchenko. A unified approach to mortality modelling using state-space framework: characterisation, identification, estimation and

  • forecasting. Annals of Actuarial Science, (May):46, 2016. ISSN

1748-4995. doi: 10.1017/S1748499517000069. URL http://arxiv.org/abs/1605.09484. Quentin Guibert, Olivier Lopez, Pierrick Piette, Quentin Guibert, Olivier Lopez, Pierrick Piette, Forecasting Mortality, and Rate

  • Improvements. Forecasting Mortality Rate Improvements with a

High-Dimensional VAR. 2017. Steven Haberman and Arthur Renshaw. A comparative study of parametric mortality projection models. Insurance: Mathematics and Economics, 48(1):35–55, 2011.

slide-116
SLIDE 116

References XI

Steven Haberman and Arthur Renshaw. Parametric mortality improvement rate modelling and projecting. Insurance: Mathematics and Economics, 50(3):309–333, 2012. Steven Haberman and Arthur Renshaw. Modelling and projecting mortality improvement rates using a cohort perspective. Insurance: Mathematics and Economics, may 2013. ISSN

  • 01676687. doi: 10.1016/j.insmatheco.2013.04.006. URL

http://linkinghub.elsevier.com/retrieve/pii/S0167668713000620. Donatien Hainaut. a Neural-Network Analyzer for Mortality Forecast. ASTIN Bulletin, (2008):1–28, 2018. ISSN 0515-0361. doi: 10.1017/asb.2017.45. URL https://www.cambridge.org/core/ product/identifier/S0515036117000459/type/journal{_}article.

  • P. Hatzopoulos and Steven Haberman. Common mortality modeling

and coherent forecasts. An empirical analysis of worldwide mortality data. Insurance: Mathematics and Economics, 52(2): 320–337, 2013. ISSN 01676687. doi: 10.1016/j.insmatheco.2012.12.009.

slide-117
SLIDE 117

References XII

John Hobcraft, Jane Menken, and Samuel Preston. Age, period, and cohort effects in demography: a review. Population Index, 48 (1):4–43, 1982. Andrew Hunt and David Blake. On the Structure and Classification

  • f Mortality Models Mortality Models. Pension Institute Working

Paper, 2015a. URL http://www.pensions-institute.org/workingpapers/wp1506.pdf. Andrew Hunt and David Blake. Identifiability in Age / Period Mortality Models. 2015b. Andrew Hunt and David Blake. Identifiability in Age / Period / Cohort Mortality Models. 2015c. Andrew Hunt and Andrés M. Villegas. Robustness and convergence in the Lee-Carter model with cohorts. Insurance: Mathematics and Economics, 64:186–202, 2015. Rob J. Hyndman and S Ullah. Robust forecasting of mortality and fertility rates: A functional data approach. Computational Statistics & Data Analysis, 51(10):4942–4956, 2007.

slide-118
SLIDE 118

References XIII

Rob J. Hyndman, Heather Booth, and Farah Yasmeen. Coherent mortality forecasting: the product-ratio method with functional time series models. Demography, 50(1):261–283, 2013. ISSN 0070-3370. doi: 10.1007/s13524-012-0145-5. S F Jarner and E M Kryger. Modelling adult mortality in small populations: The Saint Model. ASTIN Bulletin, 41(2):377–418, 2011a. S F Jarner and EM Kryger. Modelling adult mortality in small populations: The SAINT model. Astin Bulletin, (August 2013): 377–418, 2011b. doi: 10.2143/AST.41.2.2136982. URL http://demografi.dk/andet/saint.workingpaper.pdf. Torsten Kleinow. A Common Age Effect Model for the Mortality of Multiple Populations. Insurance: Mathematics and Economics, 63:147–152, 2015. ISSN 01676687. doi: 10.1016/j.insmatheco.2015.03.023.

slide-119
SLIDE 119

References XIV

Atsuyuki Kogure, Kenji Kitsukawa, and Y. Kurachi. A Bayesian comparison of models for changing mortalities toward evaluating longevity risk in Japan. Asia-Pacific Journal of Risk and Insurance, 3(2):1–21, 2009. ISSN 2153-3792. doi: 10.2202/2153-3792.1036. Marie-Claire Koissi, A Shapiro, and G Hognas. Evaluating and extending the Lee-Carter model for mortality forecasting: Bootstrap confidence interval. Insurance: Mathematics and Economics, 38(1):1–20, 2006. Ronald D. Lee and Lawrence R. Carter. Modeling and forecasting U.S. mortality. Journal of the American Statistical Association, 87(419):659–671, 1992. Hong Li and Yang Lu. Coherent Forecasting of Mortality Rates: a Sparse Vector-Autoregression Approach. ASTIN Bulletin, pages 1–38, 2017. ISSN 0515-0361. doi: 10.1017/asb.2016.37. URL https://www.cambridge.org/core/product/identifier/ S0515036116000374/type/journal{_}article.

slide-120
SLIDE 120

References XV

Jackie Li. A Poisson common factor model for projecting mortality and life expectancy jointly for females and males. Population Studies, 67(1):111–126, 2012. ISSN 0032-4728. doi: 10.1080/00324728.2012.689316. Jackie Li, Johnny Siu-Hang Li, Chong It Tan, and Leonie Tickle. Assessing basis risk in index-based longevity swap transactions.

  • 2018. ISBN 1748499518000. doi: 10.1017/S1748499518000179.

URL https://www.cambridge.org/core/product/identifier/ S1748499518000179/type/journal{_}article. Johnny Siu-Hang Li and Mary R Hardy. Measuring basis risk in longevity hedges. North American Actuarial Journal, 15(2): 177–200, 2011. Johnny Siu-Hang Li, Mary R Hardy, and Ken Seng Tan. Uncertainty in Mortality Forecasting: An Extension to the Classical Lee-Carter

  • Approach. ASTIN Bulletin, 39(1):137–164, may 2009. ISSN

0515-0361. doi: 10.2143/AST.39.1.2038060. URL http://poj. peeters-leuven.be/content.php?url=article{&}id=2038060.

slide-121
SLIDE 121

References XVI

Johnny Siu-Hang Li, Rui Zhou, and Mary R Hardy. A step-by-step guide to building two-population stochastic mortality models. Insurance: Mathematics and Economics, 63:121–134, 2015. ISSN

  • 01676687. doi: 10.1016/j.insmatheco.2015.03.021.

Johnny Siu-hang Li, Rui Zhou, and Yanxin Liu. Components of Historical Mortality Improvement Volume 1. Society of Actuaries, 1(September), 2017. Nan Li and Ronald D. Lee. Coherent mortality forecasts for a group

  • f populations: An extension of the Lee-Carter method.

Demography, 42(3):575–594, 2005. Michael Ludkovski, James Risk, and Howard Zail. Gaussian Process Models for Mortality Rates and Improvement Factors. SSRN Electronic Journal, 2016. ISSN 1556-5068. doi: 10.2139/ssrn.2831831. URL https://www.ssrn.com/abstract=2831831.

slide-122
SLIDE 122

References XVII

George Mavros, Andrew J. G. Cairns, George Streftaris, and Torsten

  • Kleinow. Stochastic Mortality Modeling: Key Drivers and

Dependent Residuals. North American Actuarial Journal, 21(3): 343–368, 2017. ISSN 10920277. doi: 10.1080/10920277.2017.1286992. URL http://dx.doi.org/10.1080/10920277.2017.1286992. Andreas Milidonis, Yijia Lin, and Samuel H. Cox. Mortality Regimes and Pricing. North American Actuarial Journal, 15(2):266–289,

  • 2011. ISSN 10920277. doi: 10.1080/10920277.2011.10597621.

URL http://www.longevity-risk.org/five/presentations/ Workshop{_}Session7/Discussion7C{_}Bauer.pdf. Daniel Mitchell, Patrick Brockett, Rafael Mendoza-Arriaga, and Kumar Muthuraman. Modeling and forecasting mortality rates. Insurance: Mathematics and Economics, 52(2):275–285, 2013.

slide-123
SLIDE 123

References XVIII

Carolyn Ndigwako Njenga and Michael Sherris. Longevity Risk and the Econometric Analysis of Mortality Trends and Volatility. Asia-Pacific Journal of Risk and Insurance, 5(2), jul 2011. ISSN 2153-3792. doi: 10.2202/2153-3792.1115. URL http://www.bepress.com/apjri/vol5/iss2/2. Colin O’Hare and Youwei Li. Explaining young mortality. Insurance: Mathematics and Economics, 50(1):12–25, 2012. Claudia Pedroza. A Bayesian forecasting model: predicting U.S. male mortality. Biostatistics (Oxford, England), 7(4):530–50, oct

  • 2006. ISSN 1465-4644.

Ermanno Pitacco. Survival models in a dynamic context: a survey. Insurance: Mathematics and Economics, 35(April):279–298, 2004. doi: 10.1016/j.insmatheco.2004.04.001. Ermanno Pitacco, Michel Denuit, Steven Haberman, and Annamaria

  • Olivieri. Modelling longevity dynamics for pensions and annuity
  • business. Oxford University Press, Oxford, 2009.
slide-124
SLIDE 124

References XIX

Richard Plat. Stochastic portfolio specific mortality and the quantification of mortality basis risk. Insurance: Mathematics and Economics, 45(1):123–132, 2009a. ISSN 0167-6687. doi: 10.1016/j.insmatheco.2009.05.002. Richard Plat. On stochastic mortality modeling. Insurance: Mathematics and Economics, 45(3):393–404, 2009b. Richard Plat. One-year Value-at-Risk for longevity and mortality. Insurance: Mathematics and Economics, 49(3):462–470, nov

  • 2011. ISSN 01676687. doi: 10.1016/j.insmatheco.2011.07.002.

URL http://linkinghub.elsevier.com/retrieve/pii/S0167668711000795. Arthur Renshaw and Steven Haberman. Lee-Carter mortality forecasting with age-specific enhancement. Insurance: Mathematics and Economics, 33(2):255–272, 2003. Arthur Renshaw and Steven Haberman. A cohort-based extension to the Lee-Carter model for mortality reduction factors. Insurance: Mathematics and Economics, 38(3):556–570, 2006.

slide-125
SLIDE 125

References XX

Arthur Renshaw and Steven Haberman. On simulation-based approaches to risk measurement in mortality with specific reference to Poisson Lee-Carter modelling. Insurance: Mathematics and Economics, 42(2):797–816, 2008. Maria Russolillo, Giuseppe Giordano, and Steven Haberman. Extending the Lee-Carter model: a three-way decomposition. Scandinavian Actuarial Journal, (2):96–117, 2011. ISSN 0346-1238. doi: 10.1080/03461231003611933. D Schrager. Affine stochastic mortality. Insurance: Mathematics and Economics, 38(1):81–97, feb 2006. ISSN 01676687. doi: 10.1016/j.insmatheco.2005.06.013. URL http://linkinghub.elsevier.com/retrieve/pii/S0167668705000946. Han Lin Shang and Steven Haberman. Grouped multivariate and functional time series forecasting: An application to annuity

  • pricing. Insurance: Mathematics and Economics, 75:166–179,
  • 2017. ISSN 01676687. doi: 10.1016/j.insmatheco.2017.05.007.

URL http://dx.doi.org/10.1016/j.insmatheco.2017.05.007.

slide-126
SLIDE 126

References XXI

Han Lin Shang and Rob J. Hyndman. Grouped Functional Time Series Forecasting: An Application to Age-Specific Mortality

  • Rates. Journal of Computational and Graphical Statistics, 26(2):

330–343, 2017. ISSN 15372715. doi: 10.1080/10618600.2016.1237877. Michael Sherris, Yajing Xu, and Jonathan Ziveyi. Cohort and Value-Based Multi-Country Longevity Risk Management. 2018. Frank van Berkum, Katrien Antonio, and Michel Vellekoop. The impact of multiple structural changes on mortality predictions. Scandinavian Actuarial Journal, 2014. Frank van Berkum, Katrien Antonio, and Michel Vellekoop. a Bayesian Joint Model for Population and Portfolio-Specific

  • Mortality. ASTIN Bulletin, pages 1–33, 2017. ISSN 0515-0361.

doi: 10.1017/asb.2017.17. URL https://www.cambridge.org/core/product/identifier/ S0515036117000174/type/journal{_}article.

slide-127
SLIDE 127

References XXII

Andrés M. Villegas and Steven Haberman. On the modeling and forecasting of socioeconomic mortality differentials: an application to deprivation and mortality in England. North American Actuarial Journal, 18(1):168–193, 2014. ISSN 1092-0277. doi: 10.1080/10920277.2013.866034. Andrés M. Villegas, Steven Haberman, Vladimir K. Kaishev, and Pietro Millossovich. A Comparative Study of Two-Population Models for the Assessment of Basis Risk in Longevity Hedges. ASTIN Bulletin, 47(03):631–679, 2017. ISSN 0515-0361. doi: 10.1017/asb.2017.18. URL https://www.cambridge.org/core/ product/identifier/S0515036117000186/type/journal{_}article. Andrés M. Villegas, Pietro Millossovich, and Vladimir K Kaishev. StMoMo: An R Package for Stochastic Mortality Modelling. Journal of Statistical Software, 84(3), 2018. doi: 10.18637/jss.v084.i03.

slide-128
SLIDE 128

References XXIII

Cheng Wan and Ljudmila Bertschi. Swiss coherent mortality model as a basis for developing longevity de-risking solutions for Swiss pension funds: A practical approach. Insurance: Mathematics and Economics, 63:66–75, 2015. ISSN 01676687. doi: 10.1016/j.insmatheco.2015.03.025. J R Wilmoth. Variation in vital rates by age, period, and cohort. Sociological Methodology, 20:295–335, 1990. ISSN 0081-1750. doi: 10.3109/10826080009148427. JR Wilmoth and T Valkonen. A parametric representation of mortality differentials over age and time. In Fith seminar of EAPS Working Group on Differential in Health, Morbidity and Mortality in Europe, 2001.

  • C. Wong-Fupuy and Steven Haberman. Projecting mortality trends:

recent developments in the United Kingdom and the United

  • States. North American Actuarial Journal, 8(2):56–83, 2004.
slide-129
SLIDE 129

References XXIV

Yajing Xu, Michael Sherris, and Jonathan Ziveyi. Market Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing. 2018. URL https://ssrn.com/abstract=3121520. Bowen Yang, Jackie Li, and Uditha Balasooriya. Cohort extensions

  • f the Poisson common factor model for modelling both genders
  • jointly. Scandinavian Actuarial Journal, (2):93–112, 2016. ISSN

0346-1238. doi: 10.1080/03461238.2014.908411. Sharon S. Yang and Chou-Wen Wang. Pricing and securitization of multi-country longevity risk with mortality dependence. Insurance: Mathematics and Economics, 52(2):157–169, 2013. ISSN

  • 01676687. doi: 10.1016/j.insmatheco.2012.10.004.

Rui Zhou, Yujiao Wang, Kai Kaufhold, Johnny Siu-Hang Li, and Ken Seng Tan. Modeling period effects in multi-population mortality models: applications to Solvency II. North American Actuarial Journal, 18(1):150–167, 2014. ISSN 1092-0277. doi: 10.1080/10920277.2013.872553.