Computes posterior draws of dependent variables' fitted values
Source:R/compute_fitted_values.R
compute_fitted_values.Rd
Each of the draws from the posterior estimation of a model is transformed into a draw from the posterior distribution of the fitted values.
Arguments
- posterior
posterior estimation outcome - an object of either of the classes: PosteriorBSVAR, PosteriorBSVARMSH, PosteriorBSVARMIX, or PosteriorBSVARSV obtained by running the
estimate
function.
Value
An object of class PosteriorFitted, that is, an NxTxS
array with attribute PosteriorFitted
containing S
draws of the fitted values.
Author
Tomasz Woźniak wozniak.tom@pm.me
Examples
# upload data
data(us_fiscal_lsuw)
# specify the model and set seed
set.seed(123)
specification = specify_bsvar$new(us_fiscal_lsuw, p = 1)
#> The identification is set to the default option of lower-triangular structural matrix.
# run the burn-in
burn_in = estimate(specification, 10)
#> **************************************************|
#> bsvars: Bayesian Structural Vector Autoregressions|
#> **************************************************|
#> Gibbs sampler for the SVAR model |
#> **************************************************|
#> Progress of the MCMC simulation for 10 draws
#> Every 10th draw is saved via MCMC thinning
#> Press Esc to interrupt the computations
#> **************************************************|
# estimate the model
posterior = estimate(burn_in, 50)
#> **************************************************|
#> bsvars: Bayesian Structural Vector Autoregressions|
#> **************************************************|
#> Gibbs sampler for the SVAR model |
#> **************************************************|
#> Progress of the MCMC simulation for 50 draws
#> Every 10th draw is saved via MCMC thinning
#> Press Esc to interrupt the computations
#> **************************************************|
# compute dependent variables' fitted values
fitted = compute_fitted_values(posterior)
# workflow with the pipe |>
############################################################
set.seed(123)
us_fiscal_lsuw |>
specify_bsvar$new(p = 1) |>
estimate(S = 50) |>
estimate(S = 100) |>
compute_fitted_values() -> fitted
#> The identification is set to the default option of lower-triangular structural matrix.
#> **************************************************|
#> bsvars: Bayesian Structural Vector Autoregressions|
#> **************************************************|
#> Gibbs sampler for the SVAR model |
#> **************************************************|
#> Progress of the MCMC simulation for 50 draws
#> Every 10th draw is saved via MCMC thinning
#> Press Esc to interrupt the computations
#> **************************************************|
#> **************************************************|
#> bsvars: Bayesian Structural Vector Autoregressions|
#> **************************************************|
#> Gibbs sampler for the SVAR model |
#> **************************************************|
#> Progress of the MCMC simulation for 100 draws
#> Every 10th draw is saved via MCMC thinning
#> Press Esc to interrupt the computations
#> **************************************************|