plot.PosteriorHD function

Plots historical decompositions

Plots of the posterior means of the historical decompositions.

## S3 method for class 'PosteriorHD'
plot(
  x,
  shock_names,
  cols,
  main,
  xlab,
  mar.multi = c(1, 4.6, 0, 4.6),
  oma.multi = c(6, 0, 5, 0),
  ...
)

Arguments

• x: an object of class PosteriorHD obtained using the compute_historical_decompositions() function containing posterior draws of historical decompositions.
• shock_names: a vector of length N containing names of the structural shocks.
• cols: an N-vector with colours of the plot
• main: an alternative main title for the plot
• xlab: an alternative x-axis label for the plot
• mar.multi: the default mar argument setting in graphics::par. Modify with care!
• oma.multi: the default oma argument setting in graphics::par. Modify with care!
• ...: additional arguments affecting the summary produced.

Examples

data(us_fiscal_lsuw)                                  # upload data
set.seed(123)                                         # set seed
specification  = specify_bsvar$new(us_fiscal_lsuw)    # specify model
burn_in        = estimate(specification, 10)          # run the burn-in
posterior      = estimate(burn_in, 20, thin = 1)      # estimate the model

# compute historical decompositions
fevd           = compute_historical_decompositions(posterior)
plot(fevd)                                            

# workflow with the pipe |>
############################################################
set.seed(123)
us_fiscal_lsuw |>
  specify_bsvar$new() |>
  estimate(S = 10) |> 
  estimate(S = 20, thin = 1) |> 
  compute_historical_decompositions() |>
  plot()

See Also

compute_historical_decompositions

Author(s)

Tomasz Woźniak wozniak.tom@pm.me