## Not run:
# # Replication of Figure 6 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")
#
# # sign restrictions
# constr <- c(+4,-3,-2,-5)
#
# # estimates the model
# model1 <- uhlig.reject(Y=uhligdata, nlags=12, draws=200, subdraws=200, nkeep=1000,
# KMIN=1, KMAX=6, constrained=constr, constant=FALSE, steps=60)
#
# # get posterior draws
# irfs1 <- model1$IRFS
#
# # estimate Fry and Pagan's MT method
# fp.target(Y=uhligdata, nlags=12, irfdraws=irfs1, constant=F, type="median",
# labels=vl, target= TRUE, save=FALSE, legend=TRUE,
# bands=c(0.16, 0.84), grid=TRUE, bw=FALSE, maxit=1000)
#
# ## End(Not run)
Run the code above in your browser using DataLab