Commit Graph

4707 Commits (7ab86d31ec62469893454a02a11cc339c822ecaa)

Author SHA1 Message Date
Stéphane Adjemian (Karaba) 5ab08ed8d3 Added missing semi colon. 2014-10-21 16:37:57 +02:00
Stéphane Adjemian (Karaba) ca28a2b778 Fixed load_m_file_data_legacy.m.
An m file data was not loading if outside the current directory.
2014-10-21 16:05:34 +02:00
Stéphane Adjemian (Karaba) 63d3157267 Do not force the user to provide a clean m data file when using the legacy data interface.
Note that  all the variables  listed in varobs  need to have  the same
number of observations.
2014-10-20 17:28:35 +02:00
Stéphane Adjemian (Karaba) f48566aeae Fixed prior bounds (according to the doc in master branch).
* Second  and  third  positional  arguments  after the  name  of  the
   estimated  parameter   in  the  estimated_params   block  are  only
   considered in the optimization stage (not in the MCMC)

 * Do not  store bounds  in bayestopt_, because  bounds do  not always
   reflect restrictions implied by prior shapes.

 * prior_bounds routine  returns a structure  (with fields lb  and ub)
   instead of a matrix.
2014-10-20 16:18:54 +02:00
Stéphane Adjemian (Karaba) fb14d24036 Do not truncate the plotted prior densities with the truncation
parameters  used   for  optimization  (second   and  third  positional
arguments  after   the  name  of   the  estimated  parameter   in  the
estimated_params block).
2014-10-16 15:27:40 +02:00
Stéphane Adjemian (Karaba) 2d61c39276 Cosmetic changes. 2014-10-16 15:25:51 +02:00
Stéphane Adjemian (Karaba) d423ae63a8 Set default value for last input argument of lyapunov_symm routine (debug flag). 2014-10-15 09:02:42 +02:00
Stéphane Adjemian (Karaba) 935f74883d Added unit tests for lyapunov_symm routine. 2014-10-13 17:42:53 +02:00
Stéphane Adjemian (Karaba) 97b63105a0 Add a parameter to the lyapunov_symm routine (debug mode). 2014-10-13 17:42:53 +02:00
Stéphane Adjemian (Karaba) 55808060cf Removed useless input (complete commit 438a671c3873e5c4e1bb02a54080e930fbc11886). 2014-10-13 17:42:52 +02:00
Johannes Pfeifer b90f3deed2 Bugfixes and improvements related to method 3 of lyapunov_symm.m
- the function was called with the wrong input argument for this case (Q instead of R*Q*R'), crashing with incompatible matrix dimensions
- the persistent variable X was not updated when the size of it changed, leading to crashes when estimation is followed by calls to DsgeSmoother.m where the state-space is different
- Also eliminates the printed output in lyapunov_symm.m that could not be turned off and clutters estimation
2014-10-13 17:42:52 +02:00
Stéphane Adjemian 573d276f1a Merge pull request #739 from JohannesPfeifer/Cosmetic_changes
Cosmetic changes, documentation, and small bugfix
2014-10-13 13:39:19 +02:00
Johannes Pfeifer 96df3b0043 Cosmetic changes to error messages
(cherry picked from commit ffa424d9d88dcd534272b154a3d9848df861914b)
2014-10-13 13:05:46 +02:00
Stéphane Adjemian 78621c03ee Merge pull request #742 from JohannesPfeifer/second_order_hessian_problems
Deal with NaN and Inf in Hessian of dynamic model
2014-10-13 12:49:38 +02:00
Johannes Pfeifer de5875a75b Fix bug still present in CutSample.m
#733 was faulty because it did not account for changes in the first line when more than 1 MH file was present. Now, FirstDraw really denotes the first draw (as opposed to the last draw as before). The subsequent computation of FirstMhFile is then also correct.
2014-10-12 11:33:55 +02:00
Stéphane Adjemian 6c112a2701 Merge pull request #743 from JohannesPfeifer/debug_info_estimation
Add check for NaN parameters to initial_estimation_checks.m
2014-10-10 22:30:12 +02:00
Stéphane Adjemian 06fbf2b2c5 Merge pull request #746 from JohannesPfeifer/rplot
Filter out empty oo_.endo_simul in rplot.m
2014-10-10 22:25:07 +02:00
Stéphane Adjemian (Charybdis) df1e734106 Merge branch 'steady_debug' 2014-10-10 21:47:24 +02:00
Stéphane Adjemian (Charybdis) f4922a8a27 Cosmetic changes. 2014-10-10 21:46:50 +02:00
Stéphane Adjemian 9e3edd6673 Merge pull request #748 from JohannesPfeifer/maxit
Change maxit default value and correct manual
2014-10-10 19:03:56 +02:00
Stéphane Adjemian 37d4cf7e35 Merge pull request #753 from JohannesPfeifer/MCMC_diag
Improve tractability of McMCDiagnostics_core.m by adding header and comm...
2014-10-10 18:46:30 +02:00
Houtan Bastani f4f7a02fd7 Provisions for MATLAB 8.4 (R2014b) 2014-10-10 10:41:04 +02:00
Johannes Pfeifer 1d68ed8da6 Improve tractability of McMCDiagnostics_core.m by adding header and comments. Also uses more intuitive variable names 2014-10-09 10:57:28 +02:00
Johannes Pfeifer 3187c70b97 Add debugging info on initial values for steady state computation 2014-09-30 15:05:12 +02:00
Johannes Pfeifer a9bb341b0a Increase default maxit for steady and deterministic simulations
Judging from forum reports, the default fails too often. Closes #747
2014-09-30 11:47:58 +02:00
Michel Juillard 6bbfa2fa6d removed 1e6 in expression forcing a non-singular Jacobian 2014-09-30 08:13:31 +02:00
Johannes Pfeifer 6a8c968411 Filter out empty oo_.endo_simul in rplot.m 2014-09-29 15:00:45 +02:00
Michel Juillard fa3e77d163 fixing bug in handling options for lmmcp 2014-09-28 10:07:08 +02:00
Johannes Pfeifer 127b7b0221 Add check for NaN parameters to initial_estimation_checks.m
Provides important info to debug cryptic crashes in later routines. Mostly applies to cases where people run estimation without stoch_simul
2014-09-25 09:02:36 +02:00
Johannes Pfeifer b43d52278d Deal with NaN and Inf in Hessian of dynamic model
Checks for these cases in stochastic_solvers.m to prevent cryptic crashes in dll-files later on
2014-09-25 08:49:19 +02:00
Johannes Pfeifer 2d39fc62e2 Make sure the output arguments of dsge_likelihood.m are actually set
Before, trend_coeff and ys were sometimes returned as empty
2014-09-21 11:54:43 +02:00
Johannes Pfeifer 3e427020d6 Cosmetic changes in headers 2014-09-21 11:51:46 +02:00
Stéphane Adjemian (Charybdis) e2fd82d403 Merge branch 'cut_sample' 2014-09-20 11:18:16 +02:00
Stéphane Adjemian (Charybdis) cc80c8662c Fixes ambiguity in the output of CutSample routine (second part of PR #733). 2014-09-20 11:17:57 +02:00
Stéphane Adjemian (Charybdis) bfc9719f05 Fix bug in in CutSample.m that led to wrong number of draws for moment computation (first part of PR #733)
FirstLine actually defined the last line. One needs to 1. Otherwise,
Draws in GetAllPosteriorDraws.m will be (NumberOfDraws+1)*nblck.
2014-09-20 11:16:50 +02:00
Stéphane Adjemian (Charybdis) 2e08025cac Fixed typo (wrong condition for defining the integration weights and nodes). 2014-09-19 22:05:42 +02:00
Johannes Pfeifer 79c5e18bfe Cosmetic changes
- Suppresses spurious output from dynare_estimation_1.m
- Removes unset output argument from dynare_estimation_init.m (which seems to be totally unused as it was not set at all and nobody reported any crashes)
- Corrects header of dynare_estimation_init.m
2014-09-19 16:20:43 +02:00
Johannes Pfeifer 08b598609e Clarify ambiguous printout in internals function
"Overall acceptance rate" could be read as over all chains.
Also capitalizes the name Matlab
2014-09-16 20:44:18 +02:00
Jukka Heikkonen 79d0e269a7 cosmetic change 2014-09-12 15:40:47 +02:00
Marco Ratto 4729459970 Properly handle case when post_metropolis=0: here steady_state must be added to smoothed variables. 2014-09-12 15:35:00 +02:00
Jukka Heikkonen 55b6d0eb93 1) check for restricted or extened var list
2) allow proper working will all GSA type of sample: prior, MC, posterior
2014-09-12 14:49:25 +02:00
Johannes Pfeifer 8655df67a9 Clean up use of verbosity option
- in dr_block.m and ep_residuals.m the option was hardcoded and disabled
- for csminwel1.m, no verbosity option exists anymore, thus definitions preceeding calls to it were redundant as was the options_.osr.verbose
2014-09-12 10:33:37 +02:00
ferhat beee1d2f82 Correct the test on the number of eigenvalues outside the unit circle: do the test on each block instead of at the global level 2014-09-11 18:19:31 +02:00
Stéphane Adjemian (Scylla) 576ed55dda Merge branch optimizer_number_5. 2014-09-11 17:38:03 +02:00
Stéphane Adjemian (Scylla) 89fa0ae9b6 Fixed typo. 2014-09-11 17:33:42 +02:00
Stéphane Adjemian (Charybdis) 7e3114847d Fix crash in Ramsey with more than 1 instrument and steady state file.
Vector ys is undefined.
2014-09-11 15:33:55 +02:00
Stéphane Adjemian (Charybdis) ce97a27c9d Cosmetic changes. 2014-09-10 17:49:33 +02:00
Johannes Pfeifer 308ce25031 Make sure that solve_algo=0 inherits options from options_.steady
Otherwise, MaxIter is hard-coded

(cherry picked from commit 2499f641ac2fab7ddf938263d0ee34f73765c4c8)
2014-09-10 17:49:33 +02:00
Johannes Pfeifer caa6d5d93e Fix display of MH-history if parameter names do not have same length
In this case, cell2mat fails, because the dimensions are not consistent. Using char() instead padds with spaces

(cherry picked from commit 47b1a2df74cf1f643ca913ab4ff6c2057e8b58cf)
2014-09-10 17:49:33 +02:00
Johannes Pfeifer 41d1533086 Add check that prevents datafile from having same name as mod-file
Otherwise, the generated m-file might be called as the datafile, leading to infinite loops. The check cannot be moved to makedataset because M_.fname is not available there.
2014-09-10 17:49:33 +02:00