stochasticInflationStationaryMean {BALD} | R Documentation |
A generic function to plot and/or return the posterior of the stochastic inflation stationary mean for models in BALD.
object |
The object from which to plot and/or return the stochastic inflation stationary mean. |
plotDensity |
A logical value. If |
plotTrace |
A logical value. If |
If the model incorporates a stochastic rate of inflation, then that rate is assumed to follow (on the log scale) an autoregressive process of order 1.
(The autoregressive process of order 1 is the discrete equivalent to an Ornstein-Uhlenbeck process.)
This function plots the posterior for the stationary mean (on the log scale), assuming such a mean was estimated.
See vignette('BALD')
.
Mainly called for the side effect of plotting.
stochasticInflationStationaryMean("AnnualAggLossDevModelOutput")
stochasticInflationRhoParameter
stochasticInflation
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) stochasticInflationStationaryMean(standard.model.output) ## End(Not run)