Development Pattern and Prediction Error for the Stochastic Bornhuetter-Ferguson Claims Reserving Method

Similar documents
Model Uncertainty within the Tweedie Exponential Dispersion Family

The Insurance Risk in the SST and in Solvency II: Modeling and Parameter Estimators

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger

Composite Hypotheses testing

STATS 306B: Unsupervised Learning Spring Lecture 10 April 30

Computing MLE Bias Empirically

ENG 8801/ Special Topics in Computer Engineering: Pattern Recognition. Memorial University of Newfoundland Pattern Recognition

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Statistics and Probability Theory in Civil, Surveying and Environmental Engineering

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Chapter 7 Generalized and Weighted Least Squares Estimation. In this method, the deviation between the observed and expected values of

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

Rockefeller College University at Albany

Logistic regression models 1/12

8/25/17. Data Modeling. Data Modeling. Data Modeling. Patrice Koehl Department of Biological Sciences National University of Singapore

Solutions Homework 4 March 5, 2018

Stochastic Claims Reserving under Consideration of Various Different Sources of Information

Predictive Analytics : QM901.1x Prof U Dinesh Kumar, IIMB. All Rights Reserved, Indian Institute of Management Bangalore

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

Fall 2012 Analysis of Experimental Measurements B. Eisenstein/rev. S. Errede

Population Design in Nonlinear Mixed Effects Multiple Response Models: extension of PFIM and evaluation by simulation with NONMEM and MONOLIX

See Book Chapter 11 2 nd Edition (Chapter 10 1 st Edition)

Negative Binomial Regression

Parameters Estimation of the Modified Weibull Distribution Based on Type I Censored Samples

Linear Approximation with Regularization and Moving Least Squares

Chain Ladder Method: Bayesian Bootstrap versus Classical Bootstrap

MLE and Bayesian Estimation. Jie Tang Department of Computer Science & Technology Tsinghua University 2012

GARETH W. PETERS, PAVEL V. SHEVCHENKO AND MARIO V. WÜTHRICH ABSTRACT

Lecture 3 Stat102, Spring 2007

4DVAR, according to the name, is a four-dimensional variational method.

9. Binary Dependent Variables

Hidden Markov Models

Effects of Ignoring Correlations When Computing Sample Chi-Square. John W. Fowler February 26, 2012

Maximum Likelihood Estimation of Binary Dependent Variables Models: Probit and Logit. 1. General Formulation of Binary Dependent Variables Models

Tracking with Kalman Filter

STK4080/9080 Survival and event history analysis

Asymptotic Efficiencies of the MLE Based on Bivariate Record Values from Bivariate Normal Distribution

/ n ) are compared. The logic is: if the two

Time-Varying Systems and Computations Lecture 6

Chapter 11: Simple Linear Regression and Correlation

The conjugate prior to a Bernoulli is. A) Bernoulli B) Gaussian C) Beta D) none of the above

Chapter 2 - The Simple Linear Regression Model S =0. e i is a random error. S β2 β. This is a minimization problem. Solution is a calculus exercise.

e i is a random error

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

UNIVERSITY OF TORONTO Faculty of Arts and Science. December 2005 Examinations STA437H1F/STA1005HF. Duration - 3 hours

MACHINE APPLIED MACHINE LEARNING LEARNING. Gaussian Mixture Regression

Chapter 14 Simple Linear Regression

where I = (n x n) diagonal identity matrix with diagonal elements = 1 and off-diagonal elements = 0; and σ 2 e = variance of (Y X).

Maximum Likelihood Estimation

x i1 =1 for all i (the constant ).

STAT 511 FINAL EXAM NAME Spring 2001

STAT 3008 Applied Regression Analysis

Laboratory 1c: Method of Least Squares

A Robust Method for Calculating the Correlation Coefficient

7: Estimation: The Least-Squares Method

Uncertainty as the Overlap of Alternate Conditional Distributions

Radar Trackers. Study Guide. All chapters, problems, examples and page numbers refer to Applied Optimal Estimation, A. Gelb, Ed.

First Year Examination Department of Statistics, University of Florida

Andreas C. Drichoutis Agriculural University of Athens. Abstract

Properties of Least Squares

1. Inference on Regression Parameters a. Finding Mean, s.d and covariance amongst estimates. 2. Confidence Intervals and Working Hotelling Bands

Chapter 3. Two-Variable Regression Model: The Problem of Estimation

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

SIO 224. m(r) =(ρ(r),k s (r),µ(r))

As is less than , there is insufficient evidence to reject H 0 at the 5% level. The data may be modelled by Po(2).

Logistic Regression. CAP 5610: Machine Learning Instructor: Guo-Jun QI

Bias-correction under a semi-parametric model for small area estimation

Statistics for Economics & Business

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

2.3 Nilpotent endomorphisms

Statistics for Managers Using Microsoft Excel/SPSS Chapter 13 The Simple Linear Regression Model and Correlation

WORKSHOP ON THE ANALYSIS OF NEURAL DATA 2001 MARINE BIOLOGICAL LABORATORY WOODS HOLE, MASSACHUSETTS

Econ Statistical Properties of the OLS estimator. Sanjaya DeSilva

MATH 829: Introduction to Data Mining and Analysis The EM algorithm (part 2)

3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. CDS Mphil Econometrics Vijayamohan. 3-Mar-14. CDS M Phil Econometrics.

Econometrics of Panel Data

Efficient nonresponse weighting adjustment using estimated response probability

Stat 543 Exam 2 Spring 2016

Number of cases Number of factors Number of covariates Number of levels of factor i. Value of the dependent variable for case k

On Outlier Robust Small Area Mean Estimate Based on Prediction of Empirical Distribution Function

Continuous vs. Discrete Goods

Unit 10: Simple Linear Regression and Correlation

Lecture 2: Prelude to the big shrink

Finite Mixture Models and Expectation Maximization. Most slides are from: Dr. Mario Figueiredo, Dr. Anil Jain and Dr. Rong Jin

Which estimator of the dispersion parameter for the Gamma family generalized linear models is to be chosen?

Stat 543 Exam 2 Spring 2016

Continuous Time Markov Chain

Estimation of the Probability of Success Based on Communication History

Explaining the Stein Paradox

Department of Quantitative Methods & Information Systems. Time Series and Their Components QMIS 320. Chapter 6

Diagnostics in Poisson Regression. Models - Residual Analysis

3.1 ML and Empirical Distribution

Random Partitions of Samples

[The following data appear in Wooldridge Q2.3.] The table below contains the ACT score and college GPA for eight college students.

Chapter 14: Logit and Probit Models for Categorical Response Variables

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

MIMA Group. Chapter 2 Bayesian Decision Theory. School of Computer Science and Technology, Shandong University. Xin-Shun SDU

β0 + β1xi and want to estimate the unknown

Computation of Higher Order Moments from Two Multinomial Overdispersion Likelihood Models

Transcription:

Development Pattern and Predcton Error for the Stochastc Bornhuetter-Ferguson Clams Reservng Method Annna Saluz, Alos Gsler, Maro V. Wüthrch ETH Zurch ASTIN Colloquum Madrd, June 2011

Overvew 1 Notaton 2 Bornhuetter-Ferguson Method (BF) 3 Normal Model 4 Estmaton of the Parameters and ther Correlatons 5 Predcton Uncertanty 6 Conclusons and Remarks A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 2 / 12

Notaton Accdent years, 0 I Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Notaton Accdent years, 0 I Development years j, 0 j J I Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Notaton Accdent years, 0 I Development years j, 0 j J I Incremental clams X,j Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Notaton Accdent years, 0 I Development years j, 0 j J I Incremental clams X,j Cumulatve clams C,j Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Notaton Accdent years, 0 I Development years j, 0 j J I Incremental clams X,j Cumulatve clams C,j Observed clams at tme I: D I = {X,j ; + j I} Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Notaton Accdent years, 0 I Development years j, 0 j J I Incremental clams X,j Cumulatve clams C,j Observed clams at tme I: D I = {X,j ; + j I} Goal: Predct D c I = {X,j; + j > I, I, j J} Table: Clams development trangle accdent development year j year 0... j... J 0. observatons D I. D c I I A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 3 / 12

Bornhuetter-Ferguson Method (BF) BF predctor for the outstandng loss labltes (IBNR and IBNeR n case of ncurred clams data) R = C,J C,I at tme I where ˆµ a pror estmate for µ = E[C,J ] R = ˆµ (1 ˆβ I ), A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 4 / 12

Bornhuetter-Ferguson Method (BF) BF predctor for the outstandng loss labltes (IBNR and IBNeR n case of ncurred clams data) R = C,J C,I at tme I where ˆµ a pror estmate for µ = E[C,J ] R = ˆµ (1 ˆβ I ), 1 ˆβ I estmated stll to come factor at the end of development year I, ˆβ J = 1, (( ˆβ j ) j=0,...,j estmated cumulatve development pattern) A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 4 / 12

Bornhuetter-Ferguson Method (BF) BF predctor for the outstandng loss labltes (IBNR and IBNeR n case of ncurred clams data) R = C,J C,I at tme I where ˆµ a pror estmate for µ = E[C,J ] R = ˆµ (1 ˆβ I ), 1 ˆβ I estmated stll to come factor at the end of development year I, ˆβ J = 1, (( ˆβ j ) j=0,...,j estmated cumulatve development pattern) ˆγ 0 = ˆβ 0, ˆγ j = ˆβ j ˆβ j 1, 1 j J, (estmated ncremental development pattern) A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 4 / 12

Bornhuetter-Ferguson Method (BF) BF predctor for the outstandng loss labltes (IBNR and IBNeR n case of ncurred clams data) R = C,J C,I at tme I where ˆµ a pror estmate for µ = E[C,J ] R = ˆµ (1 ˆβ I ), 1 ˆβ I estmated stll to come factor at the end of development year I, ˆβ J = 1, (( ˆβ j ) j=0,...,j estmated cumulatve development pattern) ˆγ 0 = ˆβ 0, ˆγ j = ˆβ j ˆβ j 1, 1 j J, (estmated ncremental development pattern) Assumpton of a cross-classfed model E[X,j ] = µ γ j A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 4 / 12

Motvaton for the Estmaton of the Pattern Often the chan ladder (CL) development pattern s used for ˆγ j J 1 ˆγ j CL = k=j ˆf 1 k J 1 k=j 1 I k 1 ˆf 1 k, =0 C,k+1 ˆfk = I k 1. =0 C,k Issue: In the BF method we are gven a pror estmates of the µ. The CL pattern gnores any pror nformaton. Goal: Incorporate the a pror estmates n the estmaton of the development pattern. A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 5 / 12

Motvaton for the Estmaton of the Pattern Often the chan ladder (CL) development pattern s used for ˆγ j J 1 ˆγ j CL = k=j ˆf 1 k J 1 k=j 1 I k 1 ˆf 1 k, =0 C,k+1 ˆfk = I k 1. =0 C,k Issue: In the BF method we are gven a pror estmates of the µ. The CL pattern gnores any pror nformaton. Goal: Incorporate the a pror estmates n the estmaton of the development pattern. If the µ were known then the best lnear unbased estmate of γ j would be γ (0) j = I j A frst canddate s the raw estmate =0 X,j I j =0 µ, 0 j J. ˆγ (0) j = However, they do not sum up to 1. I j =0 X,j I j =0 ˆµ, A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 5 / 12

Intutve Estmaton of the Pattern If the full rectangle was known an obvous estmate s gven by I =0 ˆγ j = X,j I =0 C, 0 j J.,J A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 6 / 12

Intutve Estmaton of the Pattern If the full rectangle was known an obvous estmate s gven by I =0 ˆγ j = X,j I =0 C, 0 j J.,J If only the upper trapezod and the µ are known, replace the unknown X,j by predctors µ ˆγ j ˆγ j = I j =0 X,j + I =I j+1 µ ˆγ j I J =0 C,J + ( I =I J+1 C,I + ), 0 j J. J l=i +1 µ ˆγ l A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 6 / 12

Intutve Estmaton of the Pattern If the full rectangle was known an obvous estmate s gven by I =0 ˆγ j = X,j I =0 C, 0 j J.,J If only the upper trapezod and the µ are known, replace the unknown X,j by predctors µ ˆγ j ˆγ j = I j =0 X,j + I =I j+1 µ ˆγ j I J =0 C,J + ( I =I J+1 C,I + ), 0 j J. J l=i +1 µ ˆγ l In an over-dspersed Posson model the maxmum lkelhood estmates (MLE s) satsfy these equatons. A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 6 / 12

Normal Model Model Assumptons (Normal Model) N1 The X,j are ndependent and normally dstrbuted and there exst parameters µ 0, µ 1,..., µ I and γ 0, γ 1,..., γ J wth J j=0 γ j = 1 and σ0, 2..., σj 2 such that E[X,j ] = µ γ j, Var(X,j ) = µ σ 2 j, where σj 2 > 0, 0 j J. N2 The a pror estmates ˆµ for µ = E[C,J ] are unbased and ndependent of all X l,j. A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 7 / 12

Maxmum Lkelhood Estmaton (MLE) of the γ j s We calculate the MLE s assumng that the true µ s and σj 2 s are known and then replace the µ s by the a pror estmates ˆµ and the σj 2 s by estmates ˆσ j 2. A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 8 / 12

Maxmum Lkelhood Estmaton (MLE) of the γ j s We calculate the MLE s assumng that the true µ s and σj 2 s are known and then replace the µ s by the a pror estmates ˆµ and the σj 2 s by estmates ˆσ j 2. In the Normal Model we obtan ˆγ j = I j =0 X,j I j =0 ˆµ + }{{} raw estmate ˆσ 2 j / I j =0 ˆµ J l=0 (ˆσ 2 l / I l =0 ˆµ ) ( 1 J l=0 ) I l =0 X,l I l =0 ˆµ, } {{ } correcton term where ( ˆσ j 2 = 1 I j ) I j 2 X,j =0 ˆµ X,j I j ˆµ I j =0 ˆµ, 0 j J, j I. =0 A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 8 / 12

Covarance matrx of the ˆγ j s We use the asymptotc propertes of MLE s (Fsher Informaton matrx) to estmate the covarance matrx of the ˆγ j s.: Remarks: Ĉov(ˆγ j, ˆγ k ) = ˆσ 2 j I j =0 ˆµ For the best lnear unbased estmate γ (0) j (compare wth frst summand). ( ˆσ k 2 1 / ) I k =0 ˆµ {j=k} J l=0 (ˆσ l 2/ ). I l =0 ˆµ (n the case of known µ ) we have σ 2 j Cov(γ (0) j, γ (0) k ) = 1 {j=k} I j =0 µ, Because of the sde constrant J j=0 γ j = 1 the off-dagonal correlatons must be negatve (compare wth second summand). A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 9 / 12

Predcton Uncertanty Gven all nformaton I I (.e. D I and all ˆµ ), the condtonal MSEP of the predctor R = ˆµ (1 ˆβ I ) s gven by [ ( msep R I I ( R ) = E R R ) ] 2 I I 2 J = E X,j ˆµ (1 ˆβ I ) I I j=i +1 = J j=i +1 Var(X,j ) }{{} Process Varance (PV ) ( + ) 2 ˆµ (1 ˆβ I ) µ (1 β I ) } {{ } Estmaton Error (EE ). A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 10 / 12

MSEP We estmate the condtonal MSEP gven I I as follows msep R I I ( R J ) = Var(X,j ) + Var(ˆµ )(1 ˆβ I ) 2 j=i +1 I + ˆµ 2 Var(ˆγ j ) + 2 j=0 0 j<k I Ĉov(ˆγ j, ˆγ k ), msep I =I J+1 R I I + 2 I J+1 <k I ( I ( =I J+1 R ) = I =I J+1 (1 ˆβ I )(1 ˆβ I k )Ĉov(ˆµ, ˆµ k ) ) I I k + ˆµ ˆµ k Ĉov(ˆγ j, ˆγ l ). j=0 l=0 msep R I I ( R ) A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 11 / 12

Conclusons and Remarks Under dstrbutonal assumptons we have derved estmates for the development pattern takng all relevant nformaton nto account formulas for the smoothng from the raw estmates ˆγ (0) j to the fnal estmates ˆγ j estmates for the correlatons of these estmates. We recommend to use these formulas also n the dstrbuton-free case (currently there are no estmators avalable from whch we know that they perform better). A. Saluz, A. Gsler, M. V. Wüthrch (ETH Zurch) Pattern and MSEP for BF June, 2011 12 / 12