A Forum for RATS Software and Econometrics Discussions
https://estima.com/forum/
Code: Select all
*
* Replication file for Sims and Zha(1999), "Error Bands for Impulse Responses",
* Econometrica, vol 67, no. 5, pp 1113-1156.
*
* Y-M model
*
open data simszha.xls
calendar(q) 1948
data(format=xls,org=columns) 1948:01 1989:03 gnp82 m1
*
compute nlags=4
compute nsteps=32
compute bootdraws=1000
compute mcdraws=1000
*
set gnp82 = 100.0*log(gnp82)
set m1 = 100.0*log(m1)
*
system(model=varmodel)
variables gnp82 m1
lags 1 to nlags
det constant
end(system)
*
estimate(noprint,resids=resids)
impulse(model=varmodel,steps=nsteps,results=baseirfs,noprint)
*
* Do figure 2
*
@mcvardodraws(model=varmodel,draws=mcdraws,steps=nsteps)
@mcgraphirf(model=varmodel,$
center=input,impulses=baseirfs,percent=||.025,.16,.84,.975||,$
footer="Figure 2. Pointwise 68 and 95% Posterior Bands, Y-M Model")
It sounds like you want PAGE=BYSHOCK, probably with COLUMNS=2 with 7 variables. That will produce a graph for each shock; just ignore the ones you don't need.asdavies6032 wrote:Thanks. the programs worked all right. but i'm not sure i'm getting the right response.what i want is this:
i have a 7 variable VAR - y p pcom ff nbr totr m1. i want to only show response of all 7 variables to a
shock to one variable say ff