On the Importance of Dispersion Modeling for Claims Reserving: - - PowerPoint PPT Presentation

on the importance of dispersion modeling for claims
SMART_READER_LITE
LIVE PREVIEW

On the Importance of Dispersion Modeling for Claims Reserving: - - PowerPoint PPT Presentation

On the Importance of Dispersion Modeling for Claims Reserving: Application of the Double GLM Theory Danal Davidov under the supervision of Jean-Philippe Boucher Dpartement de mathmatiques Universit du Qubec Montral August 2009


slide-1
SLIDE 1

On the Importance of Dispersion Modeling for Claims Reserving: Application of the Double GLM Theory

Danaïl Davidov under the supervision of Jean-Philippe Boucher

Département de mathématiques Université du Québec à Montréal

August 2009

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 1 / 37

slide-2
SLIDE 2

Table of Contents

1

Introduction

2

Tweedie GLM Model Definition Log-likelihood Overdispersion

3

Variance Modeling Variance in regressions Dispersion Models Double GLMs

4

Swiss Motor Industry Example Analysis of results

5

Conclusion

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 2 / 37

slide-3
SLIDE 3

Introduction

Definition

A loss reserve is a provision for an insurer’s liability for claims

Notes

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 3 / 37

slide-4
SLIDE 4

Introduction

Definition

A loss reserve is a provision for an insurer’s liability for claims A stochastic model uses random variables in a regression framework

Notes

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 3 / 37

slide-5
SLIDE 5

Introduction

Definition

A loss reserve is a provision for an insurer’s liability for claims A stochastic model uses random variables in a regression framework

Notes

Claims reserves models presented here use GLM theory as introduced in England, Verrall 2002

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 3 / 37

slide-6
SLIDE 6

Model Definition

Notations

Ci,j Incremental payments

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 4 / 37

slide-7
SLIDE 7

Model Definition

Notations

Ci,j Incremental payments wi,j Exposure

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 4 / 37

slide-8
SLIDE 8

Model Definition

Notations

Ci,j Incremental payments wi,j Exposure ri,j Incremental number of payments

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 4 / 37

slide-9
SLIDE 9

Model Definition

Notations

Ci,j Incremental payments wi,j Exposure ri,j Incremental number of payments Yi,j Normalized incremental payments

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 4 / 37

slide-10
SLIDE 10

Model Definition

Notations

Ci,j Incremental payments wi,j Exposure ri,j Incremental number of payments Yi,j Normalized incremental payments Yi,j = Ci,j

wi,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 4 / 37

slide-11
SLIDE 11

Model Definition

Hypotheses

Ci,j is a compound Poisson-Gamma distribution

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 5 / 37

slide-12
SLIDE 12

Model Definition

Hypotheses

Ci,j is a compound Poisson-Gamma distribution Frequency ∼ Poisson with mean ϑi,jwi,j and variance ϑi,jwi,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 5 / 37

slide-13
SLIDE 13

Model Definition

Hypotheses

Ci,j is a compound Poisson-Gamma distribution Frequency ∼ Poisson with mean ϑi,jwi,j and variance ϑi,jwi,j Severity ∼ Gamma with mean τi,j and variance ντ 2

i,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 5 / 37

slide-14
SLIDE 14

Model Definition

Hypotheses

Ci,j is a compound Poisson-Gamma distribution Frequency ∼ Poisson with mean ϑi,jwi,j and variance ϑi,jwi,j Severity ∼ Gamma with mean τi,j and variance ντ 2

i,j

Using the following parametrisation p = ν + 2 ν + 1 , p ∈ (1, 2) µ = ϑτ φ = ϑ1−pτ 2−p (2 − p)

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 5 / 37

slide-15
SLIDE 15

Model Definition

Tweedie Model

Yi,j ∼ Tweedie(µi,j, p, φ, wi,j)

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 6 / 37

slide-16
SLIDE 16

Model Definition

Tweedie Model

Yi,j ∼ Tweedie(µi,j, p, φ, wi,j) µi,j = eXβ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 6 / 37

slide-17
SLIDE 17

Model Definition

Tweedie Model

Yi,j ∼ Tweedie(µi,j, p, φ, wi,j) µi,j = eXβ E[Yi,j] = µi,j , Var[Yi,j] =

φ wi,j µp i,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 6 / 37

slide-18
SLIDE 18

Log-likelihood function

Tweedie Model

l =

  • i,j

ri,j log

  • (wi,j/φ)ν+1yν

i,j

(p − 1)ν(2 − p)

  • − log
  • ri,j!Γ
  • ri,jν
  • yi,j
  • + wi,j

φ

  • yi,j

µ1−p

i,j

1 − p − µ2−p

i,j

2 − p

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 7 / 37

slide-19
SLIDE 19

Parameters

FIGURE: Parameter main influence

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 8 / 37

slide-20
SLIDE 20

Parameter p

Tweedie Model

Can be estimated only when the number of payments is known

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 9 / 37

slide-21
SLIDE 21

Parameter p

Tweedie Model

Can be estimated only when the number of payments is known Otherwise, it’s supposed fixed and known

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 9 / 37

slide-22
SLIDE 22

Parameter p

Tweedie Model

Can be estimated only when the number of payments is known Otherwise, it’s supposed fixed and known p and φ need both to be estimated at the same time

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 9 / 37

slide-23
SLIDE 23

Parameter φ

Optimizing φ using the likelihood principle

  • φp

= −

i,j wi,j

  • yi,j

µ1−p

i,j

1−p − µ2−p

i,j

2−p

  • (1 + ν)

i,j ri,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 10 / 37

slide-24
SLIDE 24

Parameter φ

FIGURE: Optimizing p using the likelihood principle for φ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 11 / 37

slide-25
SLIDE 25

Parameter φ

Optimizing φ using the deviance principle

  • φp

=

  • i,j

2 N − Q

  • yi,j

y1−p

i,j

− µ1−p

i,j

1 − p − y2−p

i,j

− µ2−p

i,j

2 − p

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 12 / 37

slide-26
SLIDE 26

Parameter φ

FIGURE: Optimizing p using the deviance principle for φ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 13 / 37

slide-27
SLIDE 27

Parameter wi,j

Note

The exposure has been incorporated in the begining, within the initial hypothesis

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 14 / 37

slide-28
SLIDE 28

Parameter wi,j

Note

The exposure has been incorporated in the begining, within the initial hypothesis Different ways of incorporating the exposure within the initial hypothesis would lead to different models

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 14 / 37

slide-29
SLIDE 29

Frequency vs Severity

Compound Poisson Model

Y =

N

  • k=1

Xi

Case 1 Case 2 Case 3 E[N] 10 20 10 Var[N] 10 20 10 E[X] 10 10 20 Var[X] 100 100 400 E[Y] 100 200 200 Var[Y] 2000 4000 8000 TABLE: Mean and variance of total costs for various situations

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 15 / 37

slide-30
SLIDE 30

Frequency vs Severity

Typical situation in a long-tail business

Decreasing average frequency

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 16 / 37

slide-31
SLIDE 31

Frequency vs Severity

Typical situation in a long-tail business

Decreasing average frequency Increasing average severity

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 16 / 37

slide-32
SLIDE 32

Frequency vs Severity

Typical situation in a long-tail business

Decreasing average frequency Increasing average severity Increasing variance in the severity

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 16 / 37

slide-33
SLIDE 33

Impact of the Distribution

FIGURE: Fitted curve for Normal, Poisson and Gamma models

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 17 / 37

slide-34
SLIDE 34

Dispersion Models

Model Definition

Yi,j ∼ Tweedie(µi,j, p, φi,j, wi,j)

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 18 / 37

slide-35
SLIDE 35

Dispersion Models

Model Definition

Yi,j ∼ Tweedie(µi,j, p, φi,j, wi,j) µi,j = eXβ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 18 / 37

slide-36
SLIDE 36

Dispersion Models

Model Definition

Yi,j ∼ Tweedie(µi,j, p, φi,j, wi,j) µi,j = eXβ E[Yi,j] = µi,j , Var[Yi,j] = φi,jµp

i,j

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 18 / 37

slide-37
SLIDE 37

Dispersion Models

Model Definition

Yi,j ∼ Tweedie(µi,j, p, φi,j, wi,j) µi,j = eXβ E[Yi,j] = µi,j , Var[Yi,j] = φi,jµp

i,j

φi,j = eVγ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 18 / 37

slide-38
SLIDE 38

Dispersion Models

Log-likelihood

l =

  • i,j

ri,j log

  • (wi,j/φi,j)ν+1yν

i,j

(p − 1)ν(2 − p)

  • − log
  • ri,j!Γ
  • ri,jν
  • yi,j
  • + wi,j

φi,j

  • yi,j

µ1−p

i,j

1 − p − µ2−p

i,j

2 − p

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 19 / 37

slide-39
SLIDE 39

Dispersion Models

Notes

The p parameter is optimized at the same time as all other parameters using implicitly the likelihood principle

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 20 / 37

slide-40
SLIDE 40

Dispersion Models

Notes

The p parameter is optimized at the same time as all other parameters using implicitly the likelihood principle Accident years do not have a significant impact on the dispersion parameter

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 20 / 37

slide-41
SLIDE 41

Dispersion Models

Notes

The p parameter is optimized at the same time as all other parameters using implicitly the likelihood principle Accident years do not have a significant impact on the dispersion parameter Due to lack of data in the last column, the model was build so that the last two columns have the same dispersion parameter

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 20 / 37

slide-42
SLIDE 42

Dispersion Models

Notes

The p parameter is optimized at the same time as all other parameters using implicitly the likelihood principle Accident years do not have a significant impact on the dispersion parameter Due to lack of data in the last column, the model was build so that the last two columns have the same dispersion parameter Possibility to incorporate trends in the dispersion parameter by using the Hoerl’s curve parametrisation

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 20 / 37

slide-43
SLIDE 43

Double GLMs

Algorithm

1

Start with the initial exposure and find the normalized incremental payments

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 21 / 37

slide-44
SLIDE 44

Double GLMs

Algorithm

1

Start with the initial exposure and find the normalized incremental payments

2

Find the deviance between fitted and observed values

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 21 / 37

slide-45
SLIDE 45

Double GLMs

Algorithm

1

Start with the initial exposure and find the normalized incremental payments

2

Find the deviance between fitted and observed values

3

Model the deviance

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 21 / 37

slide-46
SLIDE 46

Double GLMs

Algorithm

1

Start with the initial exposure and find the normalized incremental payments

2

Find the deviance between fitted and observed values

3

Model the deviance

4

Establish the new exposure and start over

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 21 / 37

slide-47
SLIDE 47

Double GLMs

FIGURE: Inter-relationship between the two sub-models

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 22 / 37

slide-48
SLIDE 48

Iterative Weighted Least Squares

Algorithm additional specifications

Analogous to Fisher’s weighted scoring method for optimization

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 23 / 37

slide-49
SLIDE 49

Iterative Weighted Least Squares

Algorithm additional specifications

Analogous to Fisher’s weighted scoring method for optimization IWLS implicitly uses the deviance principle for estimating φ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 23 / 37

slide-50
SLIDE 50

Restricted Maximum Likelihood

Notes

Maximum likelihood estimators are biased downwards when the number of estimators is large compared to the number of

  • bservations
  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 24 / 37

slide-51
SLIDE 51

Restricted Maximum Likelihood

Notes

Maximum likelihood estimators are biased downwards when the number of estimators is large compared to the number of

  • bservations

REML produces estimators which are approximately and sometimes exactly unbiased

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 24 / 37

slide-52
SLIDE 52

Restricted Maximum Likelihood

Notes

Maximum likelihood estimators are biased downwards when the number of estimators is large compared to the number of

  • bservations

REML produces estimators which are approximately and sometimes exactly unbiased Approximately maximizes the penalized log-likelihood l ∗p (y; γ; p) = l(y; βγ; γ; p) + 1 2 log

  • X TWX
  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 24 / 37

slide-53
SLIDE 53

Optimizing p

Algorithm

1

Suppose p fixed and known

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 25 / 37

slide-54
SLIDE 54

Optimizing p

Algorithm

1

Suppose p fixed and known

2

Evaluate all the other parameters using the DGLM

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 25 / 37

slide-55
SLIDE 55

Optimizing p

Algorithm

1

Suppose p fixed and known

2

Evaluate all the other parameters using the DGLM

3

Evaluate the penalized log-likelihood

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 25 / 37

slide-56
SLIDE 56

Optimizing p

Algorithm

1

Suppose p fixed and known

2

Evaluate all the other parameters using the DGLM

3

Evaluate the penalized log-likelihood

4

Start over with different values for p and compare the penalized log-likelihood

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 25 / 37

slide-57
SLIDE 57

Reserve Variability

Mean Square Error of Prediction

Dispersion Models : overdispersion included implicitly in the parameter covariance matrix

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 26 / 37

slide-58
SLIDE 58

Reserve Variability

Mean Square Error of Prediction

Dispersion Models : overdispersion included implicitly in the parameter covariance matrix GLMs : overdispersion is included manually in the parameter covariance matrix

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 26 / 37

slide-59
SLIDE 59

Data Analyzed

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 27 / 37

slide-60
SLIDE 60

Data Analyzed

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 28 / 37

slide-61
SLIDE 61

Data Analyzed

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 29 / 37

slide-62
SLIDE 62

Optimizing p

FIGURE: Restricted log-likelihood for various p in a DGLM

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 30 / 37

slide-63
SLIDE 63

Analysis of Results

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 31 / 37

slide-64
SLIDE 64

Analysis of Results

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 32 / 37

slide-65
SLIDE 65

Analysis of Results

FIGURE: Adjusting the most deviant observations has a bigger influence on the deviance principle

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 33 / 37

slide-66
SLIDE 66

Analysis of Results

FIGURE: Contribution of each cell to the dispersion. p = 1.1741 , wi,j ≡ 1

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 34 / 37

slide-67
SLIDE 67

Conclusion

Further Discussion

Lack of observations and abundance of parameters is a hostile environment for DGLMs

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 35 / 37

slide-68
SLIDE 68

Conclusion

Further Discussion

Lack of observations and abundance of parameters is a hostile environment for DGLMs The deviance cannot be estimated in the last column for a DGLM and is hence ignored when estimating φ

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 35 / 37

slide-69
SLIDE 69

Conclusion

Further Discussion

Lack of observations and abundance of parameters is a hostile environment for DGLMs The deviance cannot be estimated in the last column for a DGLM and is hence ignored when estimating φ The algorithm does not converge for p fixed when there are more than 7 parameters

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 35 / 37

slide-70
SLIDE 70

Conclusion

Further Discussion

Lack of observations and abundance of parameters is a hostile environment for DGLMs The deviance cannot be estimated in the last column for a DGLM and is hence ignored when estimating φ The algorithm does not converge for p fixed when there are more than 7 parameters The algorithm does converge for p fixed for one parameter, but p cannot be optimized

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 35 / 37

slide-71
SLIDE 71

Conclusion

Further Discussion

Lack of observations and abundance of parameters is a hostile environment for DGLMs The deviance cannot be estimated in the last column for a DGLM and is hence ignored when estimating φ The algorithm does not converge for p fixed when there are more than 7 parameters The algorithm does converge for p fixed for one parameter, but p cannot be optimized The "bounds" of convergence need to be explored further

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 35 / 37

slide-72
SLIDE 72

References

England, PD and Verrall, RJ (2002). Stochastic Claims Reserving in General Insurance. Institute of Actuaries and Faculty of Actuaries Smyth, G. and Jorgensen, B. (2002). Fitting Tweedie’s Compound Poisson Model to Insurance Claims Data : Dispersion Modelling. ASTIN Bulletin, 1, 143–157 Wüthrich, MV (2003). Claims Reserving Using Tweedie’s Compound Poisson Model. ASTIN Bulletin, 331–346

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 36 / 37

slide-73
SLIDE 73

The end

  • D. Davidov (UQAM)

Loss Reserves Variance Modeling August 2009 37 / 37