BGVAR (version 2.1.0)

irfcf: Counterfactual Analysis

Description

Function to perform counterfactual analysis. It enables to neutralize the response of a specific variable to a given shock.

Usage

irfcf(x, shockvar, resp, n.ahead=24, save.store=FALSE, verbose=TRUE)

Arguments

x

an object of class bgvar.

shockvar

structural shock of interest.

resp

response variable to neutralize.

n.ahead

forecasting horizon.

save.store

If set to TRUE the full posterior is returned. Default is set to FALSE in order to save storage.

verbose

If set to FALSE it suppresses printing messages to the console.

Value

Returns a list of class bgvar.irf with the following elements:

• posterior is a four-dimensional array (K times K times n.ahead times 7) that contains 7 quantiles of the posterior distribution of the impulse response functions: the 50% ("low25" and "high75"), the 68% ("low16" and "high84") and the 90% ("low05" and "high95") credible sets along with the posterior median ("median").

• struc.obj is a list object that contains posterior quantitites needed when calculating historical decomposition and structural errors via hd.decomp.

  • A median posterior of global coefficient matrix.

  • Ginv median posterior of matrix Ginv, which describes contemporaneous relationships between countries.

  • S posterior median of matrix with country variance-covariance matrices on the main diagonal.

• model.obj is a list object that contains model-specific information, in particular

  • xglobal used data of the model.

  • plag used lag specification of the model.

• IRF_store is a four-dimensional array (K times n.ahead times nr. of shock times draws) which stores the whole posterior distribution. Exists only if save.irf.store=TRUE.

Examples

# NOT RUN {
library(BGVAR)
data(eerData)
model.ssvs.eer<-bgvar(Data=eerData,W=W.trade0012,draws=100,burnin=100,plag=1,prior="SSVS",
                      eigen=TRUE)
# }
# NOT RUN {
# very time-consuming
# }
# NOT RUN {
irfcf <- irfcf(model.ssvs.eer,shockvar="US.stir",resp="US.rer",n.ahead=24)
# }