Geweke test error when mh_nblocks=1 & mh_drop large

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.

Geweke test error when mh_nblocks=1 & mh_drop large

Postby gin » Fri Jan 24, 2014 2:52 pm

Dear all,

I've encountered an error in Geweke test: in estimation, given mh_nblocks=1, when I increase mh_drop to 0.9 from e.g. 0.5 then, regardless of whether mh_replic is large or small and for both Dynare 4.4.0 and 4.4.1, an error shows up:

Code: Select all
Geweke (1992) Convergence Tests, based on means of draws 540000 to 648000 vs 810000 to 600000.
p-values are for Chi2-test for equality of means.
       Parameter       Post. Mean        Post. Std     p-val No Taper    p-val 4% Taper    p-val 8% Taper    p-val 15% Taper
Index exceeds matrix dimensions.

Error in McMCDiagnostics (line 111)
        param_draws1 =
        param_draws(first_obs_begin_sample:last_obs_begin_sample,:);


The estimation code is:

Code: Select all
estimation(order=1,datafile=ctwData14_2000_2007, mh_replic=600000, mh_nblocks=1, mh_drop=.9,  mh_jscale=0.2,
 mode_compute=6,  nodisplay,  mode_check);
-------------------------
gin
http://gin.mozello.com
gin
 
Posts: 47
Joined: Wed Feb 27, 2013 3:10 pm
Location: Riga, Latvia

Re: Geweke test error when mh_nblocks=1 & mh_drop large

Postby jpfeifer » Mon Jan 27, 2014 7:26 am

Could you provide me with the mod-file?
------------
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: Geweke test error when mh_nblocks=1 & mh_drop large

Postby gin » Mon Jan 27, 2014 7:49 am

Attached.
Attachments
myOriginalData.txt
replace txt with mat
(6.16 KiB) Downloaded 45 times
foreign_svar27.mod
(10.04 KiB) Downloaded 55 times
-------------------------
gin
http://gin.mozello.com
gin
 
Posts: 47
Joined: Wed Feb 27, 2013 3:10 pm
Location: Riga, Latvia

Re: Geweke test error when mh_nblocks=1 & mh_drop large

Postby jpfeifer » Mon Jan 27, 2014 10:04 am

Sorry, this is a bug. The attached file should work. It will hopefully be corrected in tomorrow's snapshot.
Attachments
McMCDiagnostics.m
(16.94 KiB) Downloaded 65 times
------------
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: Google [Bot] and 4 guests