Modelling, Estimation and Hedging of Longevity Risk
|
|
- Linda Blair
- 6 years ago
- Views:
Transcription
1 IA BE Summer School 2016, K. Antonio, UvA 1 / 50 Modelling, Estimation and Hedging of Longevity Risk Katrien Antonio KU Leuven and University of Amsterdam IA BE Summer School 2016, Leuven Module II: Fitting a Lee-Carter mortality model
2 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 2 / 50 Parametric models for mortality Ambition: (dates back to the 18th century with Abraham de Moivre, 1725) specify a parametric model for the death rate (m x ), the mortality rate (q x ) or the force of mortality (µ x ). A functional expression is suggested, depending on unknown parameters which should be estimated. When the parametric model is well chosen: very useful and elegant; dimensionality of problem is reduced. If it s not well chosen: completely wrong.
3 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 3 / 50 Parametric models for mortality Gompertz (1825) s law for µ x : µ x = θ 2 θ x 3, with θ 2 > 0 and θ 3 > 1. The law tries to describe the process of human ageing. Makeham (1860) adjusted the Gompertz law: µ x = θ 1 + θ 2 θ x 3, with θ 1 0, θ 2 > 0 and θ 3 > 1.
4 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 4 / 50 Parametric models for mortality The law of Heligman & Pollard (1980) is popular among anglo saxon actuaries (see McDonald et al., 1998). Heligman & Pollard: q x (x+θ 2)θ3 = θ 1 + θ 4 exp ( θ 5 (ln x ln θ 6 ) 2 ) + θ 7 θ p 8. x x Building blocks: - first term = mortality of infants; - second term = accident hump; - third term = Gompertz law for adult mortality.
5 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 5 / 50 Parametric models for mortality Several methods have been proposed to estimate the parameters in such parametric laws for mortality. Namely: - Frère (1968); - Ballegeer (1973); - De Vylder (1975).
6 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 6 / 50 Parametric models for mortality Several methods have been proposed to estimate the parameters in such parametric laws for mortality. Note various parameterizations of mortality laws, e.g. with Makeham µ x = θ 1 + θ 2 θ3 x ( then t p x = exp θ 1 t θ 2θ3 x ) (θ3 t 1) ln θ 3 or tp x = θ t 5 θ (θt 3 1)θx 3 6 or ln p x = θ 1 + θ 7 θ x 3.
7 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 7 / 50 Parametric models for mortality Example of fitting strategy: least-squares method of Ballegeer (1973): - use e.g. Makeham s law with p x = θ 5 θ (θ3 1)θx 3 6, then ln p x = ln (θ 5 ) + (θ 3 1) θ x 3 ln (θ 6); - minimize the following least squares problem x max x=x min (ln(θ 5 ) + (θ 3 1)θ x 3 ln (θ 6 ) ln p x ) 2, - numerical techniques, like Newton-Raphson, are available to solve this optimization problem (see further).
8 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 8 / 50 Parametric models for mortality Example of fitting strategy: ML method of De Vylder (1975): - consider the d x as realizations from a random variable D x following a Binomial law D x Bin(λ x, q x ), - use Maximum Likelihood (ML) to estimate θ 1, θ 3 and θ 7 from x max λ x! (λ x=x x d x )!d x! min dx pλx x qx dx, where ln (p x ) = θ 1 + θ 7 θ3 x. This problem reduces to maximizing x max x=x min ((λ x d x ) ln (p x ) + d x ln (q x )).
9 IA BE Summer School 2016, K. Antonio, UvA Parametric models for mortality 9 / 50 Parametric models for mortality Heligman-Pollard mortality law calibrated on Belgian male q x (from 2012): Belgian 2012 male log death probabilities log q x Observed Infant Accident hump Adult Full Age(x)
10 IA BE Summer School 2016, K. Antonio, UvA Mortality at old ages 10 / 50 Mortality at old ages Very often mortality data at old ages are of bad quality. Many techniques have been developed to extrapolate mortality at old ages i.e. the closing of mortality tables. We cover the approach proposed by Kannistö (1992).
11 IA BE Summer School 2016, K. Antonio, UvA Mortality at old ages 11 / 50 Mortality at old ages Use Kannistö to close the mortality table for old ages, say x {91, 92,..., 120}: - estimate (φ 1, φ 2 ) using the relation (see Doray, 2008) logit(µ x ) = log (φ 1 ) + φ 2 x, with OLS on the ages - say - x {80, 81,..., 90}; - use these ( ˆφ 1, ˆφ 2 ) in the parametric law for µ x as proposed by Kannistö: µ x = φ 1 exp (φ 2 x) 1 + φ 1 exp (φ 2 x), for ages x > 90; - p x follows from the connection: p x = exp with piecewise constant µ x (see further). ( ) 1 0 µ x+r dr = exp ( µ x )
12 IA BE Summer School 2016, K. Antonio, UvA Mortality at old ages 12 / 50 Mortality at old ages Mortality at old ages: Kannistö applied to Belgian male data (from 2013): Belgian 2013 male death rates Kannistö regression for old ages source: ADSEI Belgian males 2013 data, source: ADSEI m x 1e 04 1e 03 1e 02 1e 01 1e+00? logit m x Age(x) Age(x) Closed mortality rate Belgian males 2013 data, source: ADSEI m x 1e 04 1e 03 1e 02 1e 01 1e Age(x)
13 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: set-up 13 / 50 Forecasting mortality: set-up Allowing for future mortality trends (+uncertainty?) is required in a number of actuarial calculations: - pensions, life annuities, long term care ( LTC ) covers, whole life sickness products; - survival probabilities over a long time horizon are required!. To avoid underestimation of the relevant liabilities - use appropriate forecasts of future mortality!
14 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: set-up 14 / 50 Forecasting mortality: set-up Consider mortality as a function of both age (x) and calendar year (or: period) (t). A dynamic mortality model is a function, say Γ(x, t), of x and t. E.g. consider one year probabilities of death, mortality odds, force of mortality,... The projected mortality model is Γ(x, t) t > t, with t the year for which the most recent (reliable) period life table is available.
15 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: set-up 15 / 50 Forecasting mortality: set-up q x,t is the probability that a person who is alive at January 1st of year t and who was born on January 1st of year t x will be death on January 1st of year t + 1. Thus, q x,t is the probability that an individual aged exactly x at exact time t will die between t and t + 1. In scientific literature: - we model q x,t ( mortality rate ) OR µ x,t ( force of mortality ); - we assume (see further) µ x+s (t + s) = µ x (t) for all 0 s < 1, and then q x,t = 1 exp ( µ x,t ). This expression allows straightforward transition from q x,t to µ x,t and vice versa.
16 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: set-up 16 / 50 Forecasting mortality: set-up A mortality/life table can be read as follows: (1) vertical arrangement: q 0 (t), q 1 (t),..., q x (t),... period life tables... t t + 1 t q 0(t) q 0(t + 1) q 0(t + 2) q 1(t) q 1(t + 1) q 1(t + 2) x... q x(t) q x(t + 1) q x(t + 2)... x q x+1(t) q x+1(t + 1) q x+1(t + 2) ω 1... q ω 1(t) q ω 1(t + 1) q ω 1(t + 2)...
17 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: set-up 17 / 50 Forecasting mortality: set-up A mortality/life table can be read as follows: (2) diagonal arrangement: q 0 (t), q 1 (t + 1),..., q x (t + x),... cohort life tables.... t t + 1 t q 0(t) q 0(t + 1) q 0(t + 2) q 1(t) q 1(t + 1) q 1(t + 2) q 2(t) q 2(t + 1) q 2(t + 2) x... q x(t) q x(t + 1) q x(t + 2)... x q x+1(t) q x+1(t + 1) q x+1(t + 2) ω 1... q ω 1(t) q ω 1(t + 1) q ω 1(t + 2)......
18 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: cohort and period life expectancy 18 / 50 Forecasting mortality: period and cohort life expectancy The period life expectancy for an x year old in year t is ex per (t) = 1 exp ( µ x,t) µ x,t + k 1 exp ( µ x+j,t ) 1 exp ( µ x+k,t). µ x+k,t k 1 j=0 The cohort life expectancy for an x year old in year t is ex coh (t) = 1 exp ( µ x,t) µ x,t + k 1 exp ( µ x+j,t+j ) 1 exp ( µ x+k,t+k). µ x+k,t+k k 1 j=0
19 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: cohort and period life expectancy 19 / 50 Forecasting mortality: period and cohort life expectancy Assume piecewise constant µ x (t) then ξ p x (t) = (p x (t)) ξ and: e per x (t) = ξ 0 = k 0 ξp x (t)dξ kp x (t) 1 0 ξp x+k (t)dξ = 1 exp ( µ x(t)) µ x (t) + k 1 exp ( µ x+j (t)) 1 exp ( µ x+k(t)), µ x+k (t) k 1 j=0 where we use (e.g.) 1 0 ξ p x (t)dξ = 1 0 (p x(t)) ξ dξ = px (t) 1 ln p x (t) = 1 exp (µx (t)) µ x (t).
20 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: deterministic techniques 20 / 50 Forecasting mortality: deterministic techniques Using (graphical) extrapolation procedures: - assume trend observed in past years can be graduated (or smoothed); - suppose the trend will continue in future years; - future mortality is estimated by extrapolating the trend itself; - an example: AG projection used in The Netherlands.
21 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: deterministic techniques 21 / 50 Forecasting mortality: deterministic techniques Projection by graphical extrapolation of annual probabilities of death:
22 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: deterministic techniques 22 / 50 Forecasting mortality: deterministic techniques How to use mortality laws in a dynamic context? - The age-pattern of mortality is summarized by (small number of) parameters; - Apply projection procedure to parameters (instead of age-specific probabilities).
23 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: deterministic techniques 23 / 50 Forecasting mortality: deterministic techniques Projection in a parametric framework with mortality laws:
24 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: stochastic models 24 / 50 Forecasting mortality: stochastic models More rigorous approach requires stochastic assumptions about mortality. As a result of the projection procedures we obtain: - point estimates; - interval estimates.
25 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: stochastic models 25 / 50 Forecasting mortality: impact on valuation of life annuities Scenario A: deterministic scenario.
26 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: stochastic models 26 / 50 Forecasting mortality: impact on valuation of life annuities Scenario B: discrete setting.
27 IA BE Summer School 2016, K. Antonio, UvA Forecasting mortality: stochastic models 27 / 50 Forecasting mortality: impact on valuation of life annuities Scenario C: continuous setting.
28 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 28 / 50 Lee Carter model Lee & Carter (1992, JASA) specify for m x,t (the central death rate): ln m x,t = β (1) x + β (2) x κ t + ɛ x,t. The dependence on age is governed by the sequences of β x (1) s and β x (2) s, where x {1,..., X }. The dependence on time by the κ t s where t {1,..., T }. The error terms ɛ x,t (with mean 0 and variance σ 2 ɛ ) reflect influences not captured by the model. Reading list: - Lee-Carter model; - Girosi & King, Understanding the Lee Carter mortality forecasting method, 2007.
29 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 29 / 50 Lee Carter model Identifiability issue? (See Nielsen & Nielsen, 2010.) Consider the set of parameters θ = (β (1) 1,..., β(1) X, β(2) 1,..., β(2) X, κ 1,..., κ T ). These are not identified without additional constraints as for any scalar c and d 0 ln m x,t = β x (1) + β x (2) κ t = (β x (1) β x (2) = (1) (2) β x + β x κ t. c) + β(2) x d {d(κ t + c)} Thus, parametrization θ and θ are equivalent.
30 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 30 / 50 Lee Carter model Solution (from literature on stochastic mortality)? Impose an identification scheme (or: constraints) which solve the estimation problem. For Lee-Carter: (for example) x β (2) x = 1 κ t = 0. t
31 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 31 / 50 Lee Carter model Interpretation of the parameters: (using the set of constraints on previous sheet) (1) β (1) x is β (1) x = 1 T ln m x,t. t (2) β (2) x : d dt ln m x,t = β x (2) d dt κ t, indicates the sensitivity of the logarithm of the death rate at age x to variations in the time index κ t. (3) κ t : evolution over time of mortality ( ) ln m x,t β x (1) x = κ t.
32 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 32 / 50 Lee Carter: ln m x,t 0 log m xt = 3 log m xt Age (x)
33 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 33 / 50 Lee Carter: ln m x,t = β (1) x 0 log m xt = β x (1) 3 log m xt Age (x)
34 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 34 / 50 Lee Carter: ln m x,t = β (1) x + κ t 0 log m xt = β x (1) + κ t 3 log m xt Age (x)
35 IA BE Summer School 2016, K. Antonio, UvA Lee-Carter model: specification 35 / 50 Lee Carter: ln m x,t = β x (1) + β x (2) κ t 0 log m xt = (1) (2) β x + β x κ t 3 log m xt 6 Range Age (x)
36 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 36 / 50 Lee Carter: calibration A least squares approach. Our objective function: ( β (1), β (2), ˆκ) = arg min β (1),β (2),κ ( ) ln m x,t β x (1) β x (2) 2 κ t = arg min β (1),β (2),κ O LS(β (1), β (2), κ). x t We find parameter estimates β (1), β (2), ˆκ by applying - a singular value decomposition; - or Newton Raphson; to the objective function.
37 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 37 / 50 Lee Carter: calibration If the estimates for β x (1), β x (2) and κ t do not satisfy the identifiability constraints: with: β (1) x β (1) x + (2) β x κ κ t ( κ t κ) β (2) β (2) x β (2) x / β (2), - β (2) the sum of the (2) β x s (from SVD or Newton Raphson); - κ the average of the ˆκ t s.
38 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 38 / 50 Lee Carter: SVD calibration Estimation of the parameters β (1) x : β (1) x since t κ t = 0 we get t O LS = 0 ln m x,t = T β x (1) + β x (2) β (1) x = 1 T ln m x,t. t κ t, t
39 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 39 / 50 Lee Carter: SVD calibration Estimation of the parameters β x (2) s and κ t s. Death rates can be combined into: m 1,1... m 1,T M =..... m X,1... m X,T, of dimension X T. Now create the matrix β (1) x Z = ln M ln m 1,1 =. ln m X,1 β (1) 1... ln m 1,T.... β (1) 1 (1) (1) β X... ln m X,T β X.
40 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 40 / 50 Lee Carter: SVD calibration The idea is to approximate Z with a product of a column and a row vector Z β (2) ˆκ, with β (2) = ( β (2) (2) 1,..., β X ) and ˆκ = (ˆκ 1,..., ˆκ T ). Therefore, we have to minimize (z x,t β x (2) κ t ) 2. Solution? Use SVD of Z. x t
41 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 41 / 50 Lee Carter: SVD calibration Find the (2) β x s and ˆκ t s that minimize: Õ LS (β (2), κ) = ( ) z x,t β (2) 2 x κ t. x t Define the square matrices Z Z and ZZ. - Let u1 be the eigenvector corresponding to the largest eigenvalue of Z Z. - Let v 1 be the eigenvector corresponding to the largest eigenvalue of Z Z. The best approximation of Z in the least squares sense is: Z Z = λ 1 v 1 u 1.
42 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 42 / 50 Lee Carter: SVD calibration The parameter estimates then become: β (2) = v 1 j v 1j ˆκ = λ 1 j v 1j u 1, satisfying the Lee Carter constraints, because u 1j = 0, since t z x,t = 0. j
43 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 43 / 50 Lee Carter: SVD calibration SVD only applicable to rectangular matrices. Missing values in the Z matrix will prevent the method from working. This is a rather old fashioned way to estimate parameters in a Lee Carter model (see further).
44 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 44 / 50 Lee Carter: OLS using Newton Raphson Update each parameter in turn using a univariate Newton Raphson recursive scheme. Take derivative of O LS (β (1), β (2), κ) wrt β x (1), β x (2) and κ t ( x, t). Solve a set of equations; each one is of the form f (ξ) = 0 (see next sheet). Starting from ξ (0) the (k + 1)th iteration gives ξ (k+1) from ξ (k) : ξ (k+1) = ξ (k) f (ξ(k) ) f (ξ (k) ). This alternative to SVD does not require a rectangular array of data.
45 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 45 / 50 Lee Carter: OLS using Newton Raphson Recall: O LS (β (1), β (2), κ) = x The system to solve is: t ( ln m x,t β (1) x β (2) x κ t ) 2. 0 = t (ln m x,t β (1) x β (2) x κ t ), for x X ; 0 = x β x (ln m x,t β (1) x β (2) x κ t ), for t T ; 0 = t κ t (ln m x,t β (1) x β (2) x κ t ), for x X.
46 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 46 / 50 Lee Carter: adjustment of ˆκ t Lee & Carter (1992) suggest an additional step: - adjust ˆκ t such that d xt x = x (1) (2) e xt exp ( ˆβ x + ˆβ x ˆκ t ), (or use another criterium, see R demo).
47 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 47 / 50 Lee Carter: calibration strategies Lee Carter model: calibration using SVD and Least Squares with Newton Raphson: Belgium Lee Carter fit Belgium Lee Carter fit β x (1) LS SVD β x (2) LS SVD Age(x) Age(x) Belgium Lee Carter fit κ t LS SVD SVD second stage Time(t)
48 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 48 / 50 Lee Carter: calibration strategies Lee Carter model: calibration using SVD and Least Squares with Newton Raphson: log q x observed q x LS fit Belgium 2012 male death probabilities Lee Carter fit SVD fit SVD second stage Age(x)
49 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 49 / 50 Demographics: computational illustration R scripts are available which: - reproduce the graphs and calculations in this module; - connect to HMD (through the demography package) and download data directly from this source. You will need an internet connection to run these scripts.
50 IA BE Summer School 2016, K. Antonio, UvA Lee Carter: calibration 50 / 50 Demographics: computational illustration R scripts are available which: - calibrate Makeham law (RdemoMakeham.R); - calibrate Heligman & Pollard law (RdemoHPlaw.R); - close for old age mortality using Kannistö (RdemoOldAgesKannisto); - calibrate Lee Carter using SVD and iterative least squares approach (RdemoLeeCarterSVDLS).
Managing Systematic Mortality Risk in Life Annuities: An Application of Longevity Derivatives
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
More informationCoale & Kisker approach
Coale & Kisker approach Often actuaries need to extrapolate mortality at old ages. Many authors impose q120 =1but the latter constraint is not compatible with forces of mortality; here, we impose µ110
More informationAleš Ahčan Darko Medved Ermanno Pitacco Jože Sambt Robert Sraka Ljubljana,
Aleš Ahčan Darko Medved Ermanno Pitacco Jože Sambt Robert Sraka Ljubljana, 11.-12-2011 Mortality data Slovenia Mortality at very old ages Smoothing mortality data Data for forecasting Cohort life tables
More informationRobust Longevity Risk Management
Robust Longevity Risk Management Hong Li a,, Anja De Waegenaere a,b, Bertrand Melenberg a,b a Department of Econometrics and Operations Research, Tilburg University b Netspar Longevity 10 3-4, September,
More informationLongevity risk and stochastic models
Part 1 Longevity risk and stochastic models Wenyu Bai Quantitative Analyst, Redington Partners LLP Rodrigo Leon-Morales Investment Consultant, Redington Partners LLP Muqiu Liu Quantitative Analyst, Redington
More informationLongevity risk: past, present and future
Longevity risk: past, present and future Xiaoming Liu Department of Statistical & Actuarial Sciences Western University Longevity risk: past, present and future Xiaoming Liu Department of Statistical &
More informationSurvival models. F x (t) = Pr[T x t].
2 Survival models 2.1 Summary In this chapter we represent the future lifetime of an individual as a random variable, and show how probabilities of death or survival can be calculated under this framework.
More informationThe Impact of Natural Hedging on a Life Insurer s Risk Situation
The Impact of Natural Hedging on a Life Insurer s Risk Situation Longevity 7 September 2011 Nadine Gatzert and Hannah Wesker Friedrich-Alexander-University of Erlangen-Nürnberg 2 Introduction Motivation
More informationExam M Fall 2005 PRELIMINARY ANSWER KEY
Exam M Fall 005 PRELIMINARY ANSWER KEY Question # Answer Question # Answer 1 C 1 E C B 3 C 3 E 4 D 4 E 5 C 5 C 6 B 6 E 7 A 7 E 8 D 8 D 9 B 9 A 10 A 30 D 11 A 31 A 1 A 3 A 13 D 33 B 14 C 34 C 15 A 35 A
More informationModelling Longevity Dynamics for Pensions and Annuity Business
Modelling Longevity Dynamics for Pensions and Annuity Business Ermanno Pitacco University of Trieste (Italy) Michel Denuit UCL, Louvain-la-Neuve (Belgium) Steven Haberman City University, London (UK) Annamaria
More informationTo apply SP models we need to generate scenarios which represent the uncertainty IN A SENSIBLE WAY, taking into account
Scenario Generation To apply SP models we need to generate scenarios which represent the uncertainty IN A SENSIBLE WAY, taking into account the goal of the model and its structure, the available information,
More informationMartingale Pricing Theory in Discrete-Time and Discrete-Space Models
IEOR E4707: Foundations of Financial Engineering c 206 by Martin Haugh Martingale Pricing Theory in Discrete-Time and Discrete-Space Models These notes develop the theory of martingale pricing in a discrete-time,
More informationDynamic Replication of Non-Maturing Assets and Liabilities
Dynamic Replication of Non-Maturing Assets and Liabilities Michael Schürle Institute for Operations Research and Computational Finance, University of St. Gallen, Bodanstr. 6, CH-9000 St. Gallen, Switzerland
More informationA Markov Chain Approach. To Multi-Risk Strata Mortality Modeling. Dale Borowiak. Department of Statistics University of Akron Akron, Ohio 44325
A Markov Chain Approach To Multi-Risk Strata Mortality Modeling By Dale Borowiak Department of Statistics University of Akron Akron, Ohio 44325 Abstract In general financial and actuarial modeling terminology
More informationMulti-armed bandit problems
Multi-armed bandit problems Stochastic Decision Theory (2WB12) Arnoud den Boer 13 March 2013 Set-up 13 and 14 March: Lectures. 20 and 21 March: Paper presentations (Four groups, 45 min per group). Before
More informationConsistently modeling unisex mortality rates. Dr. Peter Hieber, Longevity 14, University of Ulm, Germany
Consistently modeling unisex mortality rates Dr. Peter Hieber, Longevity 14, 20.09.2018 University of Ulm, Germany Seite 1 Peter Hieber Consistently modeling unisex mortality rates 2018 Motivation European
More informationDefinition 9.1 A point estimate is any function T (X 1,..., X n ) of a random sample. We often write an estimator of the parameter θ as ˆθ.
9 Point estimation 9.1 Rationale behind point estimation When sampling from a population described by a pdf f(x θ) or probability function P [X = x θ] knowledge of θ gives knowledge of the entire population.
More informationMulti-armed bandits in dynamic pricing
Multi-armed bandits in dynamic pricing Arnoud den Boer University of Twente, Centrum Wiskunde & Informatica Amsterdam Lancaster, January 11, 2016 Dynamic pricing A firm sells a product, with abundant inventory,
More informationPricing Pension Buy-ins and Buy-outs 1
Pricing Pension Buy-ins and Buy-outs 1 Tianxiang Shi Department of Finance College of Business Administration University of Nebraska-Lincoln Longevity 10, Santiago, Chile September 3-4, 2014 1 Joint work
More informationThe method of Maximum Likelihood.
Maximum Likelihood The method of Maximum Likelihood. In developing the least squares estimator - no mention of probabilities. Minimize the distance between the predicted linear regression and the observed
More informationProxy Function Fitting: Some Implementation Topics
OCTOBER 2013 ENTERPRISE RISK SOLUTIONS RESEARCH OCTOBER 2013 Proxy Function Fitting: Some Implementation Topics Gavin Conn FFA Moody's Analytics Research Contact Us Americas +1.212.553.1658 clientservices@moodys.com
More informationAnnuities. Lecture: Weeks 8-9. Lecture: Weeks 8-9 (Math 3630) Annuities Fall Valdez 1 / 41
Annuities Lecture: Weeks 8-9 Lecture: Weeks 8-9 (Math 3630) Annuities Fall 2017 - Valdez 1 / 41 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationEC316a: Advanced Scientific Computation, Fall Discrete time, continuous state dynamic models: solution methods
EC316a: Advanced Scientific Computation, Fall 2003 Notes Section 4 Discrete time, continuous state dynamic models: solution methods We consider now solution methods for discrete time models in which decisions
More informationCharacterization of the Optimum
ECO 317 Economics of Uncertainty Fall Term 2009 Notes for lectures 5. Portfolio Allocation with One Riskless, One Risky Asset Characterization of the Optimum Consider a risk-averse, expected-utility-maximizing
More informationJARAMOGI OGINGA ODINGA UNIVERSITY OF SCIENCE AND TECHNOLOGY
OASIS OF KNOWLEDGE JARAMOGI OGINGA ODINGA UNIVERSITY OF SCIENCE AND TECHNOLOGY SCHOOL OF MATHEMATICS AND ACTUARIAL SCIENCE UNIVERSITY EXAMINATION FOR DEGREE OF BACHELOR OF SCIENCE ACTUARIAL 3 RD YEAR 1
More informationWolpin s Model of Fertility Responses to Infant/Child Mortality Economics 623
Wolpin s Model of Fertility Responses to Infant/Child Mortality Economics 623 J.R.Walker March 20, 2012 Suppose that births are biological feasible in the first two periods of a family s life cycle, but
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationAmath 546/Econ 589 Univariate GARCH Models
Amath 546/Econ 589 Univariate GARCH Models Eric Zivot April 24, 2013 Lecture Outline Conditional vs. Unconditional Risk Measures Empirical regularities of asset returns Engle s ARCH model Testing for ARCH
More informationUnderstanding the Death Benefit Switch Option in Universal Life Policies
1 Understanding the Death Benefit Switch Option in Universal Life Policies Nadine Gatzert, University of Erlangen-Nürnberg Gudrun Hoermann, Munich 2 Motivation Universal life policies are the most popular
More informationRisk Management and Time Series
IEOR E4602: Quantitative Risk Management Spring 2016 c 2016 by Martin Haugh Risk Management and Time Series Time series models are often employed in risk management applications. They can be used to estimate
More informationMarket Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing
1/51 Market Price of Longevity Risk for A Multi-Cohort Mortality Model with Application to Longevity Bond Option Pricing Yajing Xu, Michael Sherris and Jonathan Ziveyi School of Risk & Actuarial Studies,
More informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationAn alternative approach for the key assumption of life insurers and pension funds
2018 An alternative approach for the key assumption of life insurers and pension funds EMBEDDING TIME VARYING EXPERIENCE FACTORS IN PROJECTION MORTALITY TABLES AUTHORS: BIANCA MEIJER JANINKE TOL Abstract
More informationPrepared by Ralph Stevens. Presented to the Institute of Actuaries of Australia Biennial Convention April 2011 Sydney
Sustainable Full Retirement Age Policies in an Aging Society: The Impact of Uncertain Longevity Increases on Retirement Age, Remaining Life Expectancy at Retirement, and Pension Liabilities Prepared by
More information2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises
96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with
More informationOutline. 1 Introduction. 2 Algorithms. 3 Examples. Algorithm 1 General coordinate minimization framework. 1: Choose x 0 R n and set k 0.
Outline Coordinate Minimization Daniel P. Robinson Department of Applied Mathematics and Statistics Johns Hopkins University November 27, 208 Introduction 2 Algorithms Cyclic order with exact minimization
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationOptimal liquidation with market parameter shift: a forward approach
Optimal liquidation with market parameter shift: a forward approach (with S. Nadtochiy and T. Zariphopoulou) Haoran Wang Ph.D. candidate University of Texas at Austin ICERM June, 2017 Problem Setup and
More informationM.I.T Fall Practice Problems
M.I.T. 15.450-Fall 2010 Sloan School of Management Professor Leonid Kogan Practice Problems 1. Consider a 3-period model with t = 0, 1, 2, 3. There are a stock and a risk-free asset. The initial stock
More informationA. 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.
Solutions to the Spring 213 Course MLC Examination by Krzysztof Ostaszewski, http://wwwkrzysionet, krzysio@krzysionet Copyright 213 by Krzysztof Ostaszewski All rights reserved No reproduction in any form
More informationSubject CS2A Risk Modelling and Survival Analysis Core Principles
` Subject CS2A Risk Modelling and Survival Analysis Core Principles Syllabus for the 2019 exams 1 June 2018 Copyright in this Core Reading is the property of the Institute and Faculty of Actuaries who
More information. Social Security Actuarial Balance in General Equilibrium. S. İmrohoroğlu (USC) and S. Nishiyama (CBO)
....... Social Security Actuarial Balance in General Equilibrium S. İmrohoroğlu (USC) and S. Nishiyama (CBO) Rapid Aging and Chinese Pension Reform, June 3, 2014 SHUFE, Shanghai ..... The results in this
More informationSteven Heston: Recovering the Variance Premium. Discussion by Jaroslav Borovička November 2017
Steven Heston: Recovering the Variance Premium Discussion by Jaroslav Borovička November 2017 WHAT IS THE RECOVERY PROBLEM? Using observed cross-section(s) of prices (of Arrow Debreu securities), infer
More informationLecture 17: More on Markov Decision Processes. Reinforcement learning
Lecture 17: More on Markov Decision Processes. Reinforcement learning Learning a model: maximum likelihood Learning a value function directly Monte Carlo Temporal-difference (TD) learning COMP-424, Lecture
More informationEvaluating Hedge Effectiveness for Longevity Annuities
Outline Evaluating Hedge Effectiveness for Longevity Annuities Min Ji, Ph.D., FIA, FSA Towson University, Maryland, USA Rui Zhou, Ph.D., FSA University of Manitoba, Canada Longevity 12, Chicago September
More informationLecture 3: Factor models in modern portfolio choice
Lecture 3: Factor models in modern portfolio choice Prof. Massimo Guidolin Portfolio Management Spring 2016 Overview The inputs of portfolio problems Using the single index model Multi-index models Portfolio
More informationVersion A. Problem 1. Let X be the continuous random variable defined by the following pdf: 1 x/2 when 0 x 2, f(x) = 0 otherwise.
Math 224 Q Exam 3A Fall 217 Tues Dec 12 Version A Problem 1. Let X be the continuous random variable defined by the following pdf: { 1 x/2 when x 2, f(x) otherwise. (a) Compute the mean µ E[X]. E[X] x
More informationDynamic Portfolio Choice II
Dynamic Portfolio Choice II Dynamic Programming Leonid Kogan MIT, Sloan 15.450, Fall 2010 c Leonid Kogan ( MIT, Sloan ) Dynamic Portfolio Choice II 15.450, Fall 2010 1 / 35 Outline 1 Introduction to Dynamic
More informationAnnuities. Lecture: Weeks 8-9. Lecture: Weeks 8-9 (Math 3630) Annuities Fall Valdez 1 / 41
Annuities Lecture: Weeks 8-9 Lecture: Weeks 8-9 (Math 3630) Annuities Fall 2017 - Valdez 1 / 41 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationMortality Improvement Rates: Modelling and Parameter Uncertainty
Mortality Improvement Rates: Modelling and Parameter Uncertainty Andrew Hunt a, Andrés M. Villegas b a Pacific Life Re, London, UK b School of Risk and Actuarial Studies and ARC Centre of Excellence in
More informationLife Tables and Selection
Life Tables and Selection Lecture: Weeks 4-5 Lecture: Weeks 4-5 (Math 3630) Life Tables and Selection Fall 2017 - Valdez 1 / 29 Chapter summary Chapter summary What is a life table? also called a mortality
More informationLife Tables and Selection
Life Tables and Selection Lecture: Weeks 4-5 Lecture: Weeks 4-5 (Math 3630) Life Tables and Selection Fall 2018 - Valdez 1 / 29 Chapter summary Chapter summary What is a life table? also called a mortality
More informationMarket Risk Analysis Volume I
Market Risk Analysis Volume I Quantitative Methods in Finance Carol Alexander John Wiley & Sons, Ltd List of Figures List of Tables List of Examples Foreword Preface to Volume I xiii xvi xvii xix xxiii
More information9.1 Principal Component Analysis for Portfolios
Chapter 9 Alpha Trading By the name of the strategies, an alpha trading strategy is to select and trade portfolios so the alpha is maximized. Two important mathematical objects are factor analysis and
More informationOn modelling of electricity spot price
, Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction
More informationA Hybrid Commodity and Interest Rate Market Model
A Hybrid Commodity and Interest Rate Market Model University of Technology, Sydney June 1 Literature A Hybrid Market Model Recall: The basic LIBOR Market Model The cross currency LIBOR Market Model LIBOR
More informationCSCI 1951-G Optimization Methods in Finance Part 07: Portfolio Optimization
CSCI 1951-G Optimization Methods in Finance Part 07: Portfolio Optimization March 9 16, 2018 1 / 19 The portfolio optimization problem How to best allocate our money to n risky assets S 1,..., S n with
More informationGPD-POT and GEV block maxima
Chapter 3 GPD-POT and GEV block maxima This chapter is devoted to the relation between POT models and Block Maxima (BM). We only consider the classical frameworks where POT excesses are assumed to be GPD,
More informationThis page is left blank intentionally
Mortality projections for Belgium, general population over 1950 2015, with corrections for adverse selection Lee-Carter modelling and adaptation to insurance mortality using Reacfin s tools Thanatos and
More informationPROBABILITY AND STATISTICS
Monday, January 12, 2015 1 PROBABILITY AND STATISTICS Zhenyu Ye January 12, 2015 Monday, January 12, 2015 2 References Ch10 of Experiments in Modern Physics by Melissinos. Particle Physics Data Group Review
More information2D penalized spline (continuous-by-continuous interaction)
2D penalized spline (continuous-by-continuous interaction) Two examples (RWC, Section 13.1): Number of scallops caught off Long Island Counts are made at specific coordinates. Incidence of AIDS in Italian
More information1. For a special whole life insurance on (x), payable at the moment of death:
**BEGINNING OF EXAMINATION** 1. For a special whole life insurance on (x), payable at the moment of death: µ () t = 0.05, t > 0 (ii) δ = 0.08 x (iii) (iv) The death benefit at time t is bt 0.06t = e, t
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationAnnuities. Lecture: Weeks Lecture: Weeks 9-11 (Math 3630) Annuities Fall Valdez 1 / 44
Annuities Lecture: Weeks 9-11 Lecture: Weeks 9-11 (Math 3630) Annuities Fall 2017 - Valdez 1 / 44 What are annuities? What are annuities? An annuity is a series of payments that could vary according to:
More informationThe Weibull in R is actually parameterized a fair bit differently from the book. In R, the density for x > 0 is
Weibull in R The Weibull in R is actually parameterized a fair bit differently from the book. In R, the density for x > 0 is f (x) = a b ( x b ) a 1 e (x/b) a This means that a = α in the book s parameterization
More informationIt Takes Two: Why Mortality Trend Modeling is more than modeling one Mortality Trend
It Takes Two: Why Mortality Trend Modeling is more than modeling one Mortality Trend Johannes Schupp Joint work with Matthias Börger and Jochen Russ IAA Life Section Colloquium, Barcelona, 23 th -24 th
More informationChapter 4 - Insurance Benefits
Chapter 4 - Insurance Benefits Section 4.4 - Valuation of Life Insurance Benefits (Subsection 4.4.1) Assume a life insurance policy pays $1 immediately upon the death of a policy holder who takes out the
More informationPakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks
Pakes (1986): Patents as Options: Some Estimates of the Value of Holding European Patent Stocks Spring 2009 Main question: How much are patents worth? Answering this question is important, because it helps
More informationIIntroduction the framework
Author: Frédéric Planchet / Marc Juillard/ Pierre-E. Thérond Extreme disturbances on the drift of anticipated mortality Application to annuity plans 2 IIntroduction the framework We consider now the global
More informationAnnuity Decisions with Systematic Longevity Risk. Ralph Stevens
Annuity Decisions with Systematic Longevity Risk Ralph Stevens Netspar, CentER, Tilburg University The Netherlands Annuity Decisions with Systematic Longevity Risk 1 / 29 Contribution Annuity menu Literature
More informationRetirement Financing: An Optimal Reform Approach. QSPS Summer Workshop 2016 May 19-21
Retirement Financing: An Optimal Reform Approach Roozbeh Hosseini University of Georgia Ali Shourideh Wharton School QSPS Summer Workshop 2016 May 19-21 Roozbeh Hosseini(UGA) 0 of 34 Background and Motivation
More informationCourse information FN3142 Quantitative finance
Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken
More informationChapter 7: Portfolio Theory
Chapter 7: Portfolio Theory 1. Introduction 2. Portfolio Basics 3. The Feasible Set 4. Portfolio Selection Rules 5. The Efficient Frontier 6. Indifference Curves 7. The Two-Asset Portfolio 8. Unrestriceted
More informationDuration Models: Parametric Models
Duration Models: Parametric Models Brad 1 1 Department of Political Science University of California, Davis January 28, 2011 Parametric Models Some Motivation for Parametrics Consider the hazard rate:
More informationInterest Rate Curves Calibration with Monte-Carlo Simulatio
Interest Rate Curves Calibration with Monte-Carlo Simulation 24 june 2008 Participants A. Baena (UCM) Y. Borhani (Univ. of Oxford) E. Leoncini (Univ. of Florence) R. Minguez (UCM) J.M. Nkhaso (UCM) A.
More informationMay 2012 Course MLC Examination, Problem No. 1 For a 2-year select and ultimate mortality model, you are given:
Solutions to the May 2012 Course MLC Examination by Krzysztof Ostaszewski, http://www.krzysio.net, krzysio@krzysio.net Copyright 2012 by Krzysztof Ostaszewski All rights reserved. No reproduction in any
More informationBasis Risk and Optimal longevity hedging framework for Insurance Company
Basis Risk and Optimal longevity hedging framework for Insurance Company Sharon S. Yang National Central University, Taiwan Hong-Chih Huang National Cheng-Chi University, Taiwan Jin-Kuo Jung Actuarial
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More informationFair value of insurance liabilities
Fair value of insurance liabilities A basic example of the assessment of MVM s and replicating portfolio. The following steps will need to be taken to determine the market value of the liabilities: 1.
More informationUnderstanding, Measuring & Managing Longevity Risk. Longevity Modelling Technical Paper
Longevity Modelling Technical Paper Table of Contents Table of Figures and Tables... 4 1.0 Introduction... 6 1.1 The Importance of Understanding Longevity Risk... 6 1.2 Deterministic vs. Stochastic Models...
More informationOptimal investments under dynamic performance critria. Lecture IV
Optimal investments under dynamic performance critria Lecture IV 1 Utility-based measurement of performance 2 Deterministic environment Utility traits u(x, t) : x wealth and t time Monotonicity u x (x,
More informationUnobserved Heterogeneity Revisited
Unobserved Heterogeneity Revisited Robert A. Miller Dynamic Discrete Choice March 2018 Miller (Dynamic Discrete Choice) cemmap 7 March 2018 1 / 24 Distributional Assumptions about the Unobserved Variables
More informationOptimal construction of a fund of funds
Optimal construction of a fund of funds Petri Hilli, Matti Koivu and Teemu Pennanen January 28, 29 Introduction We study the problem of diversifying a given initial capital over a finite number of investment
More information1. 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
1 of 17 1/4/2008 12:01 PM 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 3 c. 23 3 d. 20 3 e. 8 2. Suppose 1 for 0 x 1 s(x) = 1 ex 100 for 1
More informationThe Lost Generation of the Great Recession
The Lost Generation of the Great Recession Sewon Hur University of Pittsburgh January 21, 2016 Introduction What are the distributional consequences of the Great Recession? Introduction What are the distributional
More informationChapter 2 and 3 Exam Prep Questions
1 You are given the following mortality table: q for males q for females 90 020 010 91 02 01 92 030 020 93 040 02 94 00 030 9 060 040 A life insurance company currently has 1000 males insured and 1000
More informationHuman capital formation and macroeconomic performance in an ageing small open economy: Further Results
Human capital formation and macroeconomic performance in an ageing small open economy: Further Results Ben J. Heijdra University of Groningen; Institute for Advanced Studies; CESifo; Netspar Ward E. Romp
More informationModelling Returns: the CER and the CAPM
Modelling Returns: the CER and the CAPM Carlo Favero Favero () Modelling Returns: the CER and the CAPM 1 / 20 Econometric Modelling of Financial Returns Financial data are mostly observational data: they
More informationStatistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography
Statistical and Computational Inverse Problems with Applications Part 5B: Electrical impedance tomography Aku Seppänen Inverse Problems Group Department of Applied Physics University of Eastern Finland
More informationLONGEVITY RISK IN LIVING BENEFITS
Working Paper 23/02 LONGEVITY RISK IN LIVING BENEFITS Ermanno Pitacco Presented at the third annual CeRP conference Developing an Annuity Market in Europe Moncalieri, Turin, 21 22 June 2002 Longevity risk
More informationSimulating Continuous Time Rating Transitions
Bus 864 1 Simulating Continuous Time Rating Transitions Robert A. Jones 17 March 2003 This note describes how to simulate state changes in continuous time Markov chains. An important application to credit
More informationLearning From Data: MLE. Maximum Likelihood Estimators
Learning From Data: MLE Maximum Likelihood Estimators 1 Parameter Estimation Assuming sample x1, x2,..., xn is from a parametric distribution f(x θ), estimate θ. E.g.: Given sample HHTTTTTHTHTTTHH of (possibly
More informationEE641 Digital Image Processing II: Purdue University VISE - October 29,
EE64 Digital Image Processing II: Purdue University VISE - October 9, 004 The EM Algorithm. Suffient Statistics and Exponential Distributions Let p(y θ) be a family of density functions parameterized by
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationSupplemental Online Appendix to Han and Hong, Understanding In-House Transactions in the Real Estate Brokerage Industry
Supplemental Online Appendix to Han and Hong, Understanding In-House Transactions in the Real Estate Brokerage Industry Appendix A: An Agent-Intermediated Search Model Our motivating theoretical framework
More informationBellagos Social Long-Term Care Insurance Program: Evaluation of Sustainability
Bellagos Social Long-Term Care Insurance Program: Evaluation of Sustainability 2018 Society of Actuaries Case Study The Red Hawk (Pi)rates Team Members Jacob Jakubowicz, Jillian Milicki, Steven Rynne,
More informationOnline Appendix (Not intended for Publication): Federal Reserve Credibility and the Term Structure of Interest Rates
Online Appendix Not intended for Publication): Federal Reserve Credibility and the Term Structure of Interest Rates Aeimit Lakdawala Michigan State University Shu Wu University of Kansas August 2017 1
More information