## Not run:
# # Replication of Figure 5 of Uhlig (2005)
#
# set.seed(12345)
# data(uhligdata)
#
# # variable labels for plots
# vl <- c("GDP","GDP Deflator","Comm.Pr.Index","Fed Funds Rate",
# "NB Reserves", "Total Reserves")
#
# # estimates the model
# model0 <- rfbvar(Y=uhligdata, nlags=12, draws=1000, constant=FALSE,
# steps=60, shock=4)
#
# # get posterior draws
# irfs0 <- model0$IRFS
#
# # plot impulse response functions
# irfplot(irfdraws=irfs0, type="mean", labels=vl, save=FALSE, bands=c(0.16, 0.84),
# grid=TRUE, bw=FALSE)
#
# ## End(Not run)
Run the code above in your browser using DataLab