gompertzParameters: A generic function to plot and/or return the posterior of the...

In BALD: Robust Loss Development Using MCMC

Description

A generic function to plot and/or return the posterior of the parameters for the gompertz curve which describes the probability of payment.

Arguments

object	The object from which to plot and/or return the parameters.
parameter	A character describing which parameter to plot. “scale” for the scale parameter. “fifty.fifty” for the point at which the gompertz give a probably of fifty percent.
plotDensity	A logical value. If TRUE, then the density is plotted. If plotTrace is also TRUE, then two plots are generated. If they are both FALSE, then only the statistics are returned.
plotTrace	A logical value. If TRUE, then the trace is plotted. If plotDensity is also TRUE, then two plots are generated. If they are both FALSE, then only the statistics are returned.

Details

The scale parameter describes how steep the curve is. Larger values are steeper. Positive values indicate that the probability of a positive payment should decrease with development time. (The scale is restricted to be positive.)

The fifty.fifty parameter gives the point (in development time) when the gompertz curve gives a probability of fifty percent. See vignette('BALD').

Value

Mainly called for the side effect of plotting.

See Also

gompertzParameters,AnnualAggLossDevModelOutputWithZeros-method

Examples

rm(list=ls())
library(BALD)
data(CumulativeAutoBodilyInjuryTriangle)
CumulativeAutoBodilyInjuryTriangle <- as.matrix(CumulativeAutoBodilyInjuryTriangle)
sample.col <- (dim(CumulativeAutoBodilyInjuryTriangle)[2] - 6:0)
print(decumulate(CumulativeAutoBodilyInjuryTriangle)[1:7, sample.col])
data(HPCE)
HPCE <- as.matrix(HPCE)[,1]
HPCE.rate <- HPCE[-1] / HPCE[-length(HPCE)] - 1
print(HPCE.rate[(-10):0 + length(HPCE.rate)])
HPCE.years <- as.integer(names(HPCE.rate))
max.exp.year <- max(as.integer(
dimnames(CumulativeAutoBodilyInjuryTriangle)[[1]]))
years.to.keep <- HPCE.years <=  max.exp.year + 3
HPCE.rate <- HPCE.rate[years.to.keep]
break.model.input <- makeBreakAnnualInput(
cumulative.payments = CumulativeAutoBodilyInjuryTriangle,
stoch.inflation.weight = 1,
non.stoch.inflation.weight = 0,
stoch.inflation.rate = HPCE.rate,
first.year.in.new.regime = c(1986, 1987),
prior.for.first.year.in.new.regime=c(2,1),
exp.year.type = 'ay',
extra.dev.years = 5,
use.skew.t = TRUE,
bound.for.skewness.parameter=5)
## Not run: 
break.model.output <- runLossDevModel(
break.model.input,
burnIn=30.0E+3,
sampleSize=30.0E+3,
thin=10)
break.model.output.w.zeros <- accountForZeroPayments(break.model.output)
gompertzParameters(break.model.output.w.zeros, parameter='scale')

## End(Not run)

BALD documentation built on May 2, 2019, 6:51 a.m.