Commit Graph

1533 Commits (a85cc7b451f0439ecab1f3ee5e0efcc7b5b48d74)

Author SHA1 Message Date
Houtan Bastani a85cc7b451 preprocessor: change trend test to warning. #1389 2017-04-20 18:41:28 +02:00
Johannes Pfeifer 8db0154570 Fix typo in error message 2017-04-06 08:41:12 +02:00
Houtan Bastani c74dac7e73 preprocessor: add warning for potential non-linearities in model declared linear. #1404 2017-04-05 11:07:34 +02:00
Houtan Bastani 5021601b5b preprocessor: move equation tags out of dmath mode, include all equation tags in one set of brackets, allow equation tags to be interpreted as latex, modify manual and test accordingly. #477 2017-04-05 11:01:25 +02:00
Houtan Bastani 053407314d preprocessor: add option to write equation tags in latex output. closes #477 2017-04-04 16:08:30 +02:00
Houtan Bastani 7c09ae7f3d preprocessor: remove unused filed histval_present from mod_file_struct 2017-04-03 15:23:01 +02:00
Stéphane Adjemian (Charybdis) f9a803f4a9 Added option nonlinear_filter_initialization.
Default value is 1 (initialization with the ergodic variance of the reduced
form solution of the model approximated at order one).

If the model has unit roots, the user must use `nonlinear_filter_initialization=3`,
which select an identity matrix for the initial covariance matrix of the state variables.

A side effect of this option is to temporarily change the value of options_.qz_criterium to
a value above one (ie 1+1e-6) so that the unit roots are not rejected. If the
model has unit roots and if the and if the option
nonlinear_filter_initialization has a value less than 3, the evaluation of the
likelihood will fail, because by default the unit root is counted as an
unstable root.
2017-04-03 11:25:05 +02:00
Houtan Bastani 5be4206092 preprocessor: add options plot_init_date and plot_end_date to plot_shock_decomposition 2017-03-31 14:38:44 +02:00
Houtan Bastani bc9a619d8a preprocessor: initial_condition_decomposition front end. closes #1425 2017-03-31 13:19:04 +02:00
Johannes Pfeifer ddcad84f13 Clarify error message if model-local variable was already used on RHS 2017-03-29 14:49:37 +02:00
Houtan Bastani 5fb2ba6f7c plot_shock_decomposition: remove argument to steadystate option 2017-03-28 12:15:00 +02:00
Houtan Bastani 59322da2a7 preprocessor: add plot_shock_decomposition command. closes #1406 2017-03-24 12:11:36 +01:00
Houtan Bastani d2dfb65602 preprocessor: add realtime_shock_decomposition statement. #1406 2017-03-21 15:20:32 +01:00
Houtan Bastani ea450d4254 preprocessor: fix bug introduced in a4eb8d6b91b3b9d5c0293b78be0678c52c4395a4 2017-03-21 14:07:30 +01:00
Houtan Bastani bc5ac82612 preprocessor: add init_state option to shock_decomposition 2017-03-21 12:48:11 +01:00
Johannes Pfeifer 8c595e3298 evaluate_smoother.m: pass updated M_.params and oo_ back to base workspace
Closes #1373
2017-03-17 13:15:39 +01:00
Johannes Pfeifer 5a35040126 Make sure block-decomposed model returns with correct error code if Inf or NaN is encountered
Missing from 5f85d3a2005ae1d2b3650776cc33830967721159 where only residual was set to error code. If oo_.deterministic_simulation.status is not set to 0, model will erroneously be considered as solved and homotopy will not be triggered.
2017-03-17 10:07:13 +01:00
Houtan Bastani b4ac9aac8d preprocessor: add lag info to equation cross references. closes #1398 2017-03-10 16:59:52 +01:00
Houtan Bastani fc48bf26b9 preprocessor: write error messages to cerr instead of cout, replace \n with endl 2017-03-01 12:44:31 +01:00
Houtan Bastani 9eb7c00da1 preprocessor: fix bug that caused temporary terms not to be printed for the parameter derivatives of the static model. closes #1397 2017-03-01 12:20:55 +01:00
Houtan Bastani ad84914118 preprocessor: for consistency, use equations.size() instead of equation_number() in classes that have equations as a field 2017-02-24 11:20:54 +01:00
Houtan Bastani 83ea269155 preprocessor: aesthetic fix 2017-02-10 12:47:43 +01:00
Houtan Bastani 4dc640c635 preprocessor: aesthetic fix 2017-02-09 16:23:42 +01:00
Houtan Bastani d6e801cac8 preprocessor: change error to warning because we don’t end processing in this situation; fix error message 2017-02-09 16:21:18 +01:00
Houtan Bastani 5e1d20c8d9 preprocessor: replace exit(1) with exit(EXIT_FAILURE) 2017-02-09 12:46:37 +01:00
Stéphane Adjemian (Charybdis) 3c79e231c7 Removed penalty_hessian routine.
+ Code factorization.
 + Added an option for using the penalized objective when computing numerically
 the hessian at the mode.

Previous behaviour (introduced with penalty_hessian routine) was to compute the
hessian matrix at the mode with the penalized objective function (instead of
the original objective function). This behaviour hides problematic situations,
where the computed hessian (using the original objective) would not be full
rank. For instance, if the estimation ends up with a parameter on (or very
close to) the bounds of its possible values (which is often not a desirable
outcome), the estimated posterior variance would be zero for this
parameter (with the original objective) because the hessian is not finite in
this direction, while the posterior variance would be positive if the penalized
objective is used instead. But this estimate would not be reliable by
construction of the penalty which is quite ad-hoc (more fundamentally I do not
think that there exists any rational for approximating the covariance matrix
with the inverse of the hessian matrix if the mode is on the border of the set
of possible values).

This commit restore the behaviour previous to 2446ab02ba4b3ed88c9c5021aced076078d96007.

An option is available for computing the hessian with the penalized
objective function.
2017-02-08 13:07:25 +00:00
Houtan Bastani ea774b3f30 preprocessor: fix typo 2017-01-31 17:11:03 +01:00
Houtan Bastani 7591d5b6d3 preprocessor: remove extra exogenous variables. closes #841 2017-01-27 22:48:17 +01:00
Houtan Bastani d3455d42c1 preprocessor: don’t automatically declare external functions 2017-01-27 20:36:26 +01:00
Houtan Bastani e3fe88d9d8 preprocessor: handle undeclared external functions and lagged variables 2017-01-27 20:36:26 +01:00
Houtan Bastani ff3b2d0ae3 keep track of all times an undeclared variable is encountered 2017-01-27 20:36:26 +01:00
Houtan Bastani a9de8dd0a0 preprocessor: issue warnings for undeclared model variables when the end of the model block is encountered. #1286 2017-01-27 20:36:26 +01:00
Stéphane Adjemian (Charybdis) 10b6a794ad Fixed bug.
NO_POSTERIOR_KERNEL_DENSITY option was not modifying the expected field in
options_, consequently this option was no honored.
2017-01-20 17:35:43 +01:00
Houtan Bastani a8cb5d7e61 fix nested parens for params derivatives. closes #1201 2017-01-05 18:19:37 +01:00
Houtan Bastani a14fdd204f when >32 nested parens encountered in matlab, print error message once per file in which such a situation occurred. #1201 2017-01-05 15:19:13 +01:00
Stéphane Adjemian 07ea12dc18 Merge pull request #1358 from JohannesPfeifer/nograph
Implement posterior_nograph option
2017-01-05 06:52:08 +01:00
Houtan Bastani bf19a0aa73 fix warnings introduced in 3b20d41a92db4267c0ba13278f3885aba67e6c47 related to size_t/int comparison #1201 2017-01-02 11:38:13 +01:00
Houtan Bastani 471df861fa 3b20d41a92db4267c0ba13278f3885aba67e6c47 used to_string from C++11 standard which breaks the snapshot build. #1201 2017-01-02 11:21:28 +01:00
Houtan Bastani f32b24a591 avoid unnecessary processing when a stream has 32 or fewer nested parenthesis. #1201 2016-12-30 18:32:20 +01:00
Houtan Bastani eac97a8f49 first iteration of fix for nested parenthesis in matlab. #1201 2016-12-30 18:14:08 +01:00
Houtan Bastani 486891e75f bug fix: end preprocessing with error when more than 32 nested parenthesis are encountered. closes #1201 2016-12-28 14:02:50 +01:00
Houtan Bastani 9b0fbdb4fe ramsey_policy: accept all endogenous variables in var_list_, even auxiliary variables. closes #1355 2016-12-27 14:08:50 +01:00
Johannes Pfeifer 2a442561c6 Implement preprocessor interface for posterior_nograph option 2016-12-26 13:07:46 +01:00
Johannes Pfeifer 15bd6dd04a Finish implementation and documentation of nograph option of shock_decomposition.m 2016-12-26 12:52:28 +01:00
Johannes Pfeifer a00bf1fa90 load_mh_file: recompute results by default and load on request
Prevents general overwriting of results in oo_ and having stale results
2016-12-19 19:37:34 +01:00
Johannes Pfeifer 52e1017291 Make evaluate_smoother.m output former global variables for correct updating of these objects 2016-12-18 12:43:29 +01:00
Johannes Pfeifer cae0077ff4 Eliminate global variables from shock_decomposition.m 2016-12-18 09:57:51 +01:00
Johannes Pfeifer 3f33ab5425 Add Raftery/Lewis (1992) convergence diagnostics 2016-12-17 17:22:05 +01:00
Johannes Pfeifer e55ab84f7f Geweke diagnostics: fix bug in interface that prevents setting the respective options 2016-12-17 17:22:05 +01:00
Stéphane Adjemian 8935702a52 Merge pull request #1345 from JohannesPfeifer/robust_lin_solve
robust_lin_solve: Fix mistakes introduced in 4b83c1bf763fe5029c94b7e5…
2016-12-17 17:18:19 +01:00