finding and checking modes

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.

finding and checking modes

Postby psummers » Tue Jan 31, 2006 11:04 pm

Hi all,

I've come across an odd phenomenon when using csminwel to find the posterior mode, vs fminunc or fmincon. I get very different parameter estimates using the former than either of the others. The csminwel results are also very different to the published results (Rabanal-Rubio).

I ran mode_check to see what's going on, and the only difference is in the posterior of rho_a (see attached files). If I'm interpreting these graphs correctly, it looks like csminwel gives a secondary mode near unity, and I suspect that this is what's throwing off the estimates (the model becomes non-stationary, or nearly so). Also, the MH acceptance rates are over 70% with csminwel, but only 42% or so with the other mode finders.

Any thoughts? Thanks in advance!

PS
Attachments

[The extension jpg has been deactivated and can no longer be displayed.]

[The extension jpg has been deactivated and can no longer be displayed.]

psummers
 
Posts: 15
Joined: Mon Sep 05, 2005 8:59 am
Location: Texas Tech University

Re: finding and checking modes

Postby MichelJuillard » Thu Feb 02, 2006 6:02 pm

psummers wrote:Hi all,

I've come across an odd phenomenon when using csminwel to find the posterior mode, vs fminunc or fmincon. I get very different parameter estimates using the former than either of the others. The csminwel results are also very different to the published results (Rabanal-Rubio).

I ran mode_check to see what's going on, and the only difference is in the posterior of rho_a (see attached files). If I'm interpreting these graphs correctly, it looks like csminwel gives a secondary mode near unity, and I suspect that this is what's throwing off the estimates (the model becomes non-stationary, or nearly so). Also, the MH acceptance rates are over 70% with csminwel, but only 42% or so with the other mode finders.

Any thoughts? Thanks in advance!

PS

Hello there,
as the posterior modes that you are showing are identical for the two optimizer, I assume that what differs are the posterior mean and more generally the posterior distribution simulated by Metropolis.

My guess that you hit a bug in the last version of Dynare where the penalty for leaving the parameter subspace where there is a unique stable trajectory isn't done well.

Your finding makes it more urgent to publish a corrected version, but I'm afraid that I won't be able to do it in the next 10 days.

Best

Michel
MichelJuillard
 
Posts: 680
Joined: Thu Nov 18, 2004 10:51 am

Postby psummers » Thu Feb 02, 2006 10:25 pm

thanks Michel. Since my first post I've also noticed that Chris Sims has identified a bug in the Matlab version of csminwel (see his web page). I'm not sure if that's what's causing my results or not (I'm not familiar with the guts of optimization routines), but it could obviously have implications for Dynare.

I'll let you know if I make any progress on this...
psummers
 
Posts: 15
Joined: Mon Sep 05, 2005 8:59 am
Location: Texas Tech University

Postby MichelJuillard » Sun Mar 05, 2006 7:58 pm

Chris Sims has now kindly updated the Matlab version of csminwel.m It is now used in the last version of Dynare (3.064)

Best

Michel
MichelJuillard
 
Posts: 680
Joined: Thu Nov 18, 2004 10:51 am

Postby psummers » Tue Mar 21, 2006 8:50 pm

Following up my original post, I've re-estimated this model with the new version of dynare (with the csminwel bug fixed). The 2 mode-finding methods (mode_compute=1 and 4) now give essentially the same answer (ie, differences within simulation 'error').

Thanks to Michel and Chris Sims for sorting this out!

PS
psummers
 
Posts: 15
Joined: Mon Sep 05, 2005 8:59 am
Location: Texas Tech University


Return to Dynare help

Who is online

Users browsing this forum: Google [Bot] and 10 guests