Description Arguments Details Value See Also Examples
A generic function to plot and/or return the posterior of the mean exposure growth for models in BALD.
object |
The object from which to plot and/or return the mean exposure growth. |
plotDensity |
A logical value. If |
plotTrace |
A logical value. If |
(Optionally) exposure growth is modeled as an ar1 process. This inherently assumes that periods of high exposure growth are (or at least have the possibility of being) followed by continued high periods.
See vignette('BALD')
.
Mainly called for the side effect of plotting.
meanExposureGrowth("AnnualAggLossDevModelOutput")
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 | rm(list=ls())
options(device.ask.default=FALSE)
library(BALD)
data(IncrementalGeneralLiablityTriangle)
IncrementalGeneralLiablityTriangle <- as.matrix(IncrementalGeneralLiablityTriangle)
print(IncrementalGeneralLiablityTriangle)
data(PCE)
PCE <- as.matrix(PCE)[,1]
PCE.rate <- PCE[-1] / PCE[-length(PCE)] - 1
PCE.rate.length <- length(PCE.rate)
PCE.years <- as.integer(names(PCE.rate))
years.available <- PCE.years <= max(as.integer(
dimnames(IncrementalGeneralLiablityTriangle)[[1]]))
PCE.rate <- PCE.rate[years.available]
PCE.rate.length <- length(PCE.rate)
standard.model.input <- makeStandardAnnualInput(
incremental.payments = IncrementalGeneralLiablityTriangle,
stoch.inflation.weight = 1,
non.stoch.inflation.weight = 0,
stoch.inflation.rate = PCE.rate,
exp.year.type = 'ay',
extra.dev.years=5,
use.skew.t=TRUE)
## Not run:
standard.model.output <- runLossDevModel(standard.model.input,
burnIn=30.0E+3,
sampleSize=30.0E+3,
thin=10)
meanExposureGrowth(standard.model.output)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.