
Last chance! 50% off unlimited learning
Sale ends in
BootChainLadder
procedure provides a predictive
distribution of reserves for a cumulative claims development triangle.BootChainLadder(Triangle, R = 999, process.distr=c("gamma", "od.pois"))
qpaid
for how to use (mxn)-development triangles with
mc(m,n,R)
with the simulated claimsc(m,1,R)
with the modeled
IBNRs by origin periodc(m,n,R)
with the modeled
IBNR development trianglesBootChainLadder
function uses a two-stage
bootstrapping/simulation approach. In the first stage an ordinary
chain-ladder methods is applied to the cumulative claims triangle.
From this we calculate the scaled Pearson residuals which we bootstrap
R times to forecast future incremental claims payments via the
standard chain-ladder method.
In the second stage we simulate the process error with the bootstrap
value as the mean and using the process distribution assumed.
The set of reserves obtained in this way forms the predictive distribution,
from which summary statistics such as mean, prediction error or
quantiles can be derived.summary.BootChainLadder
, MackChainLadder
# See as well the example in section 8 of England & Verrall's paper on page 55.
B <- BootChainLadder(RAA, R=999, process.distr="gamma")
B
plot(B)
# Compare to MackChainLadder
MackChainLadder(RAA)
quantile(B, c(0.75,0.95,0.99, 0.995))
# fit a distribution to the IBNR
library(MASS)
plot(ecdf(B$IBNR.Totals))
# fit a log-normal distribution
fit <- fitdistr(B$IBNR.Totals[B$IBNR.Totals>0], "lognormal")
fit
curve(plnorm(x,fit$estimate["meanlog"], fit$estimate["sdlog"]), col="red", add=TRUE)
# See as well the ABC example in the Barnett and Zehnwirth paper
A <- BootChainLadder(ABC, R=999, process.distr="gamma")
A
plot(A, log=TRUE)
Run the code above in your browser using DataLab