Bayesian IRFs retrieval

This forum is closed. You can read the posts but cannot write. We have migrated the forum to a new location where you will have to reset your password.
Forum rules
This forum is closed. You can read the posts but cannot write. We have migrated the forum to a new location (https://forum.dynare.org) where you will have to reset your password.

Bayesian IRFs retrieval

Postby aadegboye abiodun » Fri Jan 30, 2015 10:01 pm

Hi everyone,
I estimated a DSGE model with Dynare with Bayesian estinmation. The Bayesian IRFs showed and appeared with the grey areas showing bounds in the result in the pup up window but later disapperaed. How do i handle it. I want to use the Byesian IRFs as it is with grey dust or care.
Kindly guide
aadegboye abiodun
 
Posts: 18
Joined: Sat Nov 01, 2014 9:05 pm

Re: Bayesian IRFs retrieval

Postby jpfeifer » Sat Jan 31, 2015 7:43 am

As documented on Pfeifer (2014): An Introduction to Graphs in Dynare at https://sites.google.com/site/pfeiferecon/dynare, they are stored in the Output subfolder.
------------
Johannes Pfeifer
University of Cologne
https://sites.google.com/site/pfeiferecon/
jpfeifer
 
Posts: 6940
Joined: Sun Feb 21, 2010 4:02 pm
Location: Cologne, Germany

Re: Bayesian IRFs retrieval

Postby aadegboye abiodun » Sat Jan 31, 2015 9:38 pm

Hi Jfpeifer,


Thank you for the reply so helpful. I want to ask do I plot the Bayesian posteriorIRF manually using MATLAB? If yes, am I to plot the mean for instance and hold it on so as to plot the bounds HPDinf and HPDsup together with it? Kindly suggest what I can do to place the grey colour or shape on it.
aadegboye abiodun
 
Posts: 18
Joined: Sat Nov 01, 2014 9:05 pm

Re: Bayesian IRFs retrieval

Postby jpfeifer » Sun Feb 01, 2015 7:55 am

If you want to do the graphs manually, you need to work with area-plots. The code Dynare uses is:
Code: Select all
                h1 = area(1:options_.irf,HPDIRF(:,2,j,i));
                set(h1,'FaceColor',[.9 .9 .9]);
                set(h1,'BaseValue',min(HPDIRF(:,1,j,i)));
                hold on
                h2 = area(1:options_.irf,HPDIRF(:,1,j,i),'FaceColor',[1 1 1],'BaseValue',min(HPDIRF(:,1,j,i)));
                set(h2,'FaceColor',[1 1 1]);
                set(h2,'BaseValue',min(HPDIRF(:,1,j,i)));
                plot(1:options_.irf,MeanIRF(:,j,i),'-k','linewidth',3)
                % plot([1 options_.irf],[0 0],'-r','linewidth',0.5);
                box on
                axis tight
                xlim([1 options_.irf]);
                hold off

where HPDIRF(:,2,j,i) is the HPDsup and HPDIRF(:,1,j,i) is the HPDinf. You need to adjust is accordingly.
------------
Johannes Pfeifer
University of Cologne
https://sites.google.com/site/pfeiferecon/
jpfeifer
 
Posts: 6940
Joined: Sun Feb 21, 2010 4:02 pm
Location: Cologne, Germany


Return to Dynare help

Who is online

Users browsing this forum: No registered users and 12 guests