Dear Tom, Could you please tell me how could I edit the impulse response functions graphs, in order to be able to clearly distinguish between the two IRF functions and the standard errors. Attached to his post you will also find the code and the dataset.

spgraph(vfields=1,hfields=5)

do j=1,1

dofor i=1 2 4 5

compute factor76 = acresp76(j,j)(1)

set respons76 = acresp76(i,j)/factor76

set stde = sqrt(seconda(t,i+(j-1)*nvars))/factor76

set upp = respons76 + band * stde

set low = respons76 - band * stde

graph(header=implabel(i),hlabel="horizon",number=1) 3

# upp 1 nsteps 2

# low 1 nsteps 2

# respons76 1 nsteps 1

end do i

end do j

print / respons76 upp low

compute j=1

set respons76 = resp76(3,j)/factor76

set stde = sqrt(second(t,3+(j-1)*nvars))/factor76

set upp = respons76 + band * stde

set low = respons76 - band * stde

graph(header=implabel(3),hlabel="horizon",number=1) 3

# upp 1 nsteps 2

# low 1 nsteps 2

# respons76 1 nsteps 1

print / respons76 upp low

spgraph(done)