Commit Graph

3015 Commits (e0018f8bca66263eabface151347cfa70d6989f7)

Author SHA1 Message Date
Stéphane Adjemian (Charybdis) bcaae72338 Added pop method to the dynDates class. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 11c30cee3b Cosmetic change. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 28e146fc69 Fixed header. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) fde2694e0c Added append method in dynDates class (with unitary tests). 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 32905b1ce1 Added more complex calls to subsref overloaded method for dynDate class. Added unitary test. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 68a2ce4f49 Fixed bufg in overloaded subsref method for dynDates class. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 47dd0d09ad Changed behavior of ne overload method (dynDate) as in commit ##9dc16fc765bfd8957dbe24a2ac1244fbde1edf09. Removed useless variable. Updated copyright. 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 55a3129644 Removed useless variable + Cosmetic change + Updated copyright 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 8b3e6e8c4d Removed useless variable + Cosmetic change + Updated copyright 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 62794a10c1 Removed useless variable + Cosmetic change + Updated copyright 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 42331f2239 Updated copyright 2012-11-14 12:32:02 +01:00
Stéphane Adjemian (Charybdis) 602d266bba Removed useless variable + Cosmetic change. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 4339673431 Removed useless variable in dynDate/eq. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 2def9f88dd Adapted unitary to the previous commit (#9dc16fc765bfd8957dbe24a2ac1244fbde1edf09). 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 131c423059 Changed behavior of eq overloaded method (dynDate class). Trigger an error message if the two compared objects have no common frequency. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 3c911e5b2c Updated copyright. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) c40dd62531 Fixed error messages. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) cae06f3d80 Added unitary test in colon overladed method (dynDate). 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 1df227ce21 Updated copyright. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 244395c2e7 Changed subsref overloaded method in dynDate class.
Allow to populate an empty dynDate object or update a dynDate object. For instance, if we instantiate
the dynDate class as follows:

a = dynDate()

a is an empty dynDate object. We can "populate" this object as follows

a('1938Q4')

Note however that a ins not modified in place (due to a limitation of Octave's class design), a copy
has to be made.
2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 52ebca922b Updated copyright. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 05cfedfd51 Added comment. 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) 50a2f94197 Updated copyright 2012-11-14 12:32:01 +01:00
Stéphane Adjemian (Charybdis) c1616872a4 Removed useless path informations in (dynSeries) horzcat overloaded method. 2012-11-14 12:32:00 +01:00
Stéphane Adjemian (Charybdis) d92a7775b9 Removed debug info in first unitary test of dynSeries constructor. 2012-11-14 12:32:00 +01:00
Stéphane Adjemian (Charybdis) bc1276a218 Removed useless path information in unitary test (dynSeries constructor). 2012-11-14 12:32:00 +01:00
Stéphane Adjemian (Charybdis) b8803c149e Added unitary test for dynSeries class. 2012-11-14 12:32:00 +01:00
Houtan Bastani ddddb5ad0f console => nodisplay (Completes ticket #282) 2012-11-12 16:48:44 +01:00
Stéphane Adjemian (Charybdis) 4604148d2f Put back pruning in particle filter (implemented in sequential_importance_particle_filter). 2012-11-07 15:58:56 +01:00
Stéphane Adjemian (Charybdis) c7bd05224f Fixed issue with weights when computing first and second order moments. 2012-11-07 10:48:05 +01:00
Stéphane Adjemian (Charybdis) a9d3ecb0ed Added comments. 2012-11-07 10:13:11 +01:00
Michel Juillard 2e59409df4 corrected bug introduced in recent commit 48e00c5789 2012-11-01 09:26:05 +01:00
Michel Juillard ea504638fc fixing bug for local approximation of models with variables appearing
only as forward looking (I believe that such models can't be solved,
but it is a different issue).
2012-10-31 17:56:42 +01:00
Sébastien Villemot 0ab230b474 Fix copyright notices 2012-10-31 17:03:49 +01:00
Houtan Bastani a23b497a49 bug fix: none instead of pdf to graph_format 2012-10-31 16:16:45 +01:00
Michel Juillard 48e00c5789 corrected and clarified error messages when steady state computation fails 2012-10-31 11:09:25 +01:00
Houtan Bastani 6c7501d906 move console setting to global_initialization.m 2012-10-30 14:10:54 +01:00
Houtan Bastani fc7ea133d9 force console mode when matlab started with -nojvm or -nodisplay flags or octave started with --no-window-system flag (ticket 283) 2012-10-30 13:50:51 +01:00
Michel Juillard 34db22c3bc ms-sbvar: made option final_year optional is using the entire sample 2012-10-24 14:05:54 +02:00
Michel Juillard 74c300939a making sure that dsge_likelihood.m always returns a value in fval 2012-10-24 10:02:17 +02:00
Michel Juillard b326636d0b sbvar: replaced warning about sample size different from available
data by error if available data are shorter than declared sample
2012-10-23 14:33:02 +02:00
Michel Juillard d5f29e0b89 sbvar: made option final_year optional is using the entire sample 2012-10-23 14:33:01 +02:00
Michel Juillard bb6d3f539e fixing bug in writing variance decomposition files 2012-10-21 10:14:38 +02:00
Michel Juillard 7922366f71 fix bug in setting oo_.dr.state_var. This bug could generate wrong row
labels in Policy and Transition Functions.
2012-10-14 22:54:49 +02:00
Michel Juillard 4c18bfea47 allow for missing column headings in XLS(X) data file. Such a column
can't be used for data, but it happens that dates for example don't
have a column headings
2012-10-14 22:50:50 +02:00
Sébastien Villemot c2c88600a1 More explicit error message when no estimated_params and no smoother
Closes: #222
2012-10-08 14:42:09 +02:00
Sébastien Villemot cc26a1d806 Add comment for the use case when there is no estimated_params 2012-10-08 14:32:22 +02:00
Sébastien Villemot 2f9e34c5e5 Various fixes related to data-loading
- fix support for filenames with extension (.m, .mat, .xls), since the manual
  and the bison grammar both support it
- fix support for XLS files under Octave
- add support for XLSX files
- run test for loading XLS files
2012-10-08 13:00:55 +02:00
Michel Juillard c88a8aaaa6 changed warnings into errors when load_mh_file or recovery is used and
there is no metropolis history file
2012-10-08 11:50:48 +02:00
Michel Juillard c722f53955 fixed bug with options mode_compute=0,mode_file=....
added tests for mode_compute=4 and mode_compute=6. They increase computing time of the
test suite by about 30 minutes.
2012-10-06 16:51:42 +02:00
Marco Ratto e835c7d945 Bug fix: graphic options are not passed to global options_ (for the next release try to skip the use of the global here). 2012-10-02 08:59:49 +02:00
Michel Juillard dd8f16c8db Revert "reset penalty to 1.e8." Too big a penaly is not desirable.
This reverts commit fcceec896a.
2012-10-01 14:23:21 +02:00
Marco Ratto fcceec896a reset penalty to 1.e8. 2012-10-01 08:53:50 +02:00
Marco Ratto 560771bedc Trap the case when second output of objective function is empty. 2012-10-01 08:50:47 +02:00
Marco Ratto 285e40eb2c Pass new options nodisplay and graph_format to options_ 2012-10-01 08:49:58 +02:00
Stéphane Adjemian (Charybdis) e81f9d48ac Improved display of mode_compute=6 optimization algorithm. Fixed bugs. Changed the options. 2012-09-29 00:08:05 +02:00
Stéphane Adjemian (Charybdis) 3744061112 Save the names of the estimated parameters in <MODEL_NAME>_mode.mat when mode_compute is equal to six (stochastic optimization). 2012-09-29 00:08:05 +02:00
Sébastien Villemot 421e1a39fc Merge remote-tracking branch 'ratto/master' 2012-09-27 16:36:00 +02:00
Houtan Bastani 7c6ae272f6 ms-sbvar: remove field removal for file tags 2012-09-27 15:32:33 +02:00
Houtan Bastani 12af7298ef ms-sbvar: remove unused field removal statements (and set options_.datafile='' in global_initialization) 2012-09-27 15:32:33 +02:00
Stéphane Adjemian (Charybdis) 894b3d69f4 Added an option to decide if dsge_likelihood should call univariate filters when the covariance matrix of
the prediction error is singular (default is yes).
2012-09-27 14:48:07 +02:00
Houtan Bastani 868afeb953 ms-sbvar bug fix: store no_create_init as number instead of string 2012-09-26 18:07:49 +02:00
Marco Ratto 16fa6efc8d Replace NaN with numeric penalty when diffuse filter does not kill unit roots. 2012-09-26 16:00:18 +02:00
Stéphane Adjemian (Charybdis) 7f42c66723 Fixed bug. Added missing routine which set the step length when computing a numerical gradient. 2012-09-24 12:59:04 +02:00
Stéphane Adjemian (Charybdis) 7e3fbe0b2f Fixed bug in get_prior_info. Initialize objective_function_penalty_base. 2012-09-24 12:59:04 +02:00
Stéphane Adjemian (Charybdis) da159038e6 Added check for BK and steady state in the optimization of the prior density. 2012-09-24 12:59:04 +02:00
Stéphane Adjemian (Charybdis) af7ac86996 Initialize returned variables. 2012-09-24 12:59:04 +02:00
Stéphane Adjemian (Charybdis) df26c6686d Fixed bug. 2012-09-24 12:59:03 +02:00
Stéphane Adjemian (Charybdis) 233d0ae0d2 Added output to get_prior_info (non empty iff info=1). 2012-09-24 12:59:03 +02:00
Stéphane Adjemian (Charybdis) e552b0c0e0 Added fake outputs. 2012-09-24 12:59:03 +02:00
Michel Juillard 8db0fee155 removed useless lines 2012-09-23 21:04:22 +02:00
Stéphane Adjemian (Charybdis) d0d0b6ed52 Use new numgrad routines.
options_.gradient method = {12,22} two points formula.
options_.gradient method = {13,23} three points formula.
options_.gradient method = {15,25} four points formula.

Still need to decide what to do with (local) option zgrad.
2012-09-21 17:05:36 +02:00
Stéphane Adjemian (Charybdis) 4ce4242d0d Added new numgrad* routines (h specific to each direction). 2012-09-21 17:05:36 +02:00
Stéphane Adjemian (Scylla) e67d3d3fea Added missing semicolon. 2012-09-20 16:38:19 +02:00
Michel Juillard f944a43b0c added code for computing risky steady state in regular models (not
portfolio models)
2012-09-19 22:51:44 +02:00
Michel Juillard c8c102b43a reorganized code and added comments 2012-09-19 22:51:44 +02:00
Michel Juillard 0a5909f1d2 fixed bug with deterministic exogenous variables at order==2 (moved
all relevant code to stochastic_solvers.m)
2012-09-19 22:51:44 +02:00
Stéphane Adjemian (Scylla) c29fa9d81c Fixed bug (wrong index). 2012-09-19 17:03:07 +02:00
Stéphane Adjemian (Scylla) 914d97fb53 Added waitbar to prior_sampler routine. 2012-09-19 16:27:43 +02:00
Stéphane Adjemian (Scylla) 0f7076846f Save *all* the prior draws generated by prior_sampler, even if info(1) returned by resol is positive. 2012-09-19 13:44:57 +02:00
Stéphane Adjemian (Scylla) 465bb6c144 Added missing call to set_state_space. 2012-09-19 12:27:11 +02:00
Stéphane Adjemian (Scylla) 59a3db051f Temporarly set options_.order equal to one in get_prior_info routine. 2012-09-19 12:07:16 +02:00
Sébastien Villemot 7fd9ab1473 Remove unused variables 2012-09-19 10:47:12 +02:00
Johannes Pfeifer 9d555075d6 Fixed bug in dyn_second_order_solver.m where ghuu was not correctly stored in dr 2012-09-19 00:05:28 +02:00
Sébastien Villemot 7de6be31e6 Upgrade the minimum requirement to MATLAB 7.3 (R2006b) 2012-09-18 11:54:40 +02:00
Sébastien Villemot 6adb58dbc1 Remove workaround for Octave 3.2 (we now require 3.4) 2012-09-18 11:42:11 +02:00
Sébastien Villemot c41ce4ae69 Merge remote-tracking branch 'ratto/master' 2012-09-18 11:13:50 +02:00
Sébastien Villemot fd7386b593 Merge remote-tracking branch 'jpfeifer/master' 2012-09-18 11:13:46 +02:00
Sébastien Villemot cba7a3fcc8 Homogeneize behavior for PDF and FIG under Octave 2012-09-18 11:12:11 +02:00
Marco Ratto fce6c5cc31 Increase accuracy in printing cosn for collinearity patterns 2012-09-18 10:42:01 +02:00
Marco Ratto dbbc8e66b7 Fix to compute derivatives with Kronecker products (just for testing). 2012-09-18 10:41:16 +02:00
Sébastien Villemot 8d4c812671 Merge remote-tracking branch 'ratto/master' 2012-09-17 10:29:53 +02:00
Johannes Pfeifer d50e68ef2d Merge branch 'master' of ssh://kirikou.dynare.org/srv/d_kirikou/git/dynare 2012-09-16 12:11:51 +02:00
Johannes Pfeifer 24a5bd4868 Fixed bug where when using mode_compute=0 with a mode-file from mode_compute=6, mh_jscale is not loaded. 2012-09-16 12:10:53 +02:00
Johannes Pfeifer c2ce75aa45 Fixed bug in set_prior.m when correlations are estimated. Corrected description in header of prior_bounds.m 2012-09-14 18:00:12 +02:00
Marco Ratto b756ee1819 Analytic derivatives: make them work for lik_init==2 as well
(lik_init =3 and =4 still to be worked out)
2012-09-14 17:07:38 +02:00
Marco Ratto e637319be5 bug fix. the diffuse filter should simply penalize the likelihood with the NaN without breaking the estimation (this is in line with missing_observations_kalman_filter_d.m). 2012-09-14 17:05:35 +02:00
Stéphane Adjemian (Scylla) 3658b02455 Fixed bug. If info==19 (steady state routine returns check=1) info(2) is NaN, so that it is not possible to define an endogenous penalty. In this case we just add one to objective_function_penalty_base.
(cherry picked from commit d2912b264c326b3f349984b605787045b028b992)
2012-09-14 17:05:31 +02:00
Sébastien Villemot c74b1d1bdc Provisions for MATLAB 8.0 (R2012b) 2012-09-14 11:10:47 +02:00
Marco Ratto 428b193c11 Need to pass the entire global workspace, if subroutines called by the _core loops need them. Updated parallel test. 2012-09-11 15:45:08 +02:00
Sébastien Villemot c65ff9d988 Fix crash with MATLAB 7.1 and 7.2 under Windows
Closes: #277
2012-09-11 12:32:27 +02:00
Sébastien Villemot 2d66c68301 Fix function header 2012-09-11 11:16:40 +02:00
Michel Juillard 07193b2cfa fixing bug in computing full_rank for purely backward models 2012-09-11 10:06:38 +02:00
Michel Juillard 92833d3ceb In CHECK, use the eigenvalues as computed by the reordered real
generalized Schur decomposition, rather than the ones computed by
eig()
2012-09-10 14:27:56 +02:00
Michel Juillard a22d1d415a replaced rank() by rcond() in evaluating whether Z22 is full rank in
checking Blanchard and Kahn conditions with CHECK
2012-09-10 13:26:05 +02:00
Stéphane Adjemian (Scylla) 205b455ad7 Added the possibility, for each given size of the sample, to restart an arbitrary number of times the estimation when estimating a model with a recursive approach.
Deactivated by default (options_.recursive_estimation_restart is defined to be zero in global_initialization.m).
2012-09-06 14:13:29 +02:00
Stéphane Adjemian (Scylla) 8a0fe91480 Removed annoying warning messages. 2012-09-06 12:22:32 +02:00
Stéphane Adjemian (Scylla) 43e2c9ecef Added a routine to test if a file exist. 2012-09-06 12:21:28 +02:00
Stéphane Adjemian (Scylla) 19c46dba55 Improved recursive estimation. Run the estimation using previous results (saved in <M_.fname>_mode.mat). 2012-09-06 11:14:48 +02:00
Marco Ratto 6beb4e3f8f small bug fix 2012-08-30 14:59:32 +02:00
Michel Juillard e8f159f66b fixing bugs in previous commit 2012-08-30 12:44:46 +02:00
Michel Juillard 1bac2d34c3 adding a trap to catch the case where the random generator of the
master is not available on the slave
2012-08-30 12:24:05 +02:00
Michel Juillard f77b101d7b fixing bugs in random generator handling code 2012-08-29 22:12:50 +02:00
Michel Juillard 378413ed14 fixing typo in previous commit 2012-08-29 21:02:33 +02:00
Michel Juillard dbdbfdd926 adding set_dynare_random_generator_state() and
get_dynare_random_generator_state(). Use now different seeds for
different Metropolis chains. Fixed handling of random generator state
thourghout the code.
2012-08-29 17:58:54 +02:00
Stéphane Adjemian (Charybdis) d9f3ab5be8 Partially revert commit #69efc894c6dc9ac1250bd7450bd57443f088c242. Test for Octave/Matlab to decide how to call the print command. Added a warning stating that Octave cannot create pdf files. 2012-08-29 16:50:08 +02:00
Michel Juillard d60202616a fixed problem with penalty in estimation. Created a new global scalar:
objective_function_penalty_base. It is the only simple way that I
found to keep csminwel1.m to be able to handle general functions.
2012-08-28 12:17:07 +02:00
Michel Juillard 526d6ca76c removing unused function 2012-08-28 11:55:16 +02:00
Marco Ratto e519b04713 bug fix: when nograph=1, SmoothedShocks were not saved. 2012-08-27 17:45:47 +02:00
Marco Ratto 1193cab7c8 Make the function compatible with multiple file formats. 2012-08-27 16:37:14 +02:00
Stéphane Adjemian (Charybdis) 0508bf9a41 Added missing default options for stochastic simulation of nonlinear backward looking models.
(cherry picked from commit 4053f664ecec14ffb7269c6bbf78355ca46b6b92)
2012-08-27 12:53:11 +02:00
Marco Ratto a7ba2b51e7 Trap error when the model does not solve for point estimation (prior mean-mode posterior mean-mode) 2012-08-24 17:09:13 +02:00
Marco Ratto 9edce5414f bug fix for octave. 2012-08-24 16:40:24 +02:00
Marco Ratto fda047e19c updated penalty has to be properly passed to the objective function.
This relates should fix behaviour after commit 6b3bd9dd0b.
2012-08-24 16:39:04 +02:00
Stéphane Adjemian (Charybdis) fe98a0875d Added the possibility to save the plots generated by the shock_decomposition command. 2012-08-24 15:03:35 +02:00
Stéphane Adjemian (Charybdis) 150256268f Fixed mixed indices (j was used for indexing two nested loops). 2012-08-24 15:00:21 +02:00
Stéphane Adjemian (Charybdis) 69efc894c6 Fixed bug. Removed extension declaration for Encapsulated PostScript graphics file format. 2012-08-24 14:56:44 +02:00
Stéphane Adjemian (Charybdis) b89bdb6fd6 Bug fix. Call evaluate_steady_state_routine, otherwise evaluate_smoother crashes if the steady state of the model is specified using the steady state model block (or writing a steadystate2.m routine). 2012-08-24 12:51:44 +02:00
Sébastien Villemot 56914e3d4c Require at least Octave 3.4 (linsolve does not compile against 3.2) 2012-08-23 16:40:00 +02:00
Marco Ratto 6b91301365 Reduce memory requirements for analytic Hessian.
Tested with QUEST III (Ratto et al. 2009): 63 params and 59 state dimension.
2012-08-21 16:00:55 +02:00
Marco Ratto b02303cf69 Force analytic_derivation = 1; 2012-08-21 15:53:02 +02:00
Marco Ratto b1dd7a5137 Allow quicker evaluation of likelihood with analytic derivatives. 2012-08-21 15:46:35 +02:00
Marco Ratto 7683175e8e Bug fix in terms for outer product gradient with analytic derivatives 2012-08-21 15:45:25 +02:00
Marco Ratto 9fd8bf954e bug fix in formula for analytic Hessian. 2012-08-15 14:19:34 +02:00
Marco Ratto f7aae77f16 removed useless elements + reduce loops for efficiency
(cherry picked from commit f2271264994af253ccdc7f01478320e4d2c2c0e7)
2012-08-13 13:04:34 +02:00
Marco Ratto 3ddafb164b -) Added missing terms for analytic Hessian when steady state depends on estimated params;
-) bug fixes;
(cherry picked from commit c84f70f6630f4988716dcb4ea59315180bbb36e7)
2012-08-13 13:04:07 +02:00
Michel Juillard 70e162c736 removed globals from solve1.m 2012-08-06 23:32:40 +02:00
Sébastien Villemot db62e3d7a0 issquare exists under Octave and returns the matrix dim if it is square 2012-08-06 18:26:16 +02:00
Sébastien Villemot 2e65a9ab96 Provide a better implementation of linsolve for Octave
Closes: #273
2012-08-06 18:26:16 +02:00
Michel Juillard a400305f20 updating list of info codes returned by dynare_resolve() in DsgeVarLikelihood 2012-08-06 12:02:35 +02:00
Michel Juillard 2e9ae5df28 dealing with info == 8 (Nan in Jacobian) during estimation 2012-08-06 12:00:03 +02:00
Sébastien Villemot 336fc56518 Merge remote-tracking branch 'jpfeifer/master' 2012-08-06 11:44:15 +02:00
Michel Juillard 7e221e5d31 factoring setting of bayestop_.penalty = Inf in Metropolis 2012-08-05 15:27:14 +02:00
Michel Juillard 7959102d43 removed global and varargin in Metropolis routines 2012-08-05 15:10:21 +02:00
Johannes Pfeifer 1a968d5f22 Added check for NaN in Jacobian to stochastic_solvers and added explicit error message to print_info.m 2012-08-05 13:09:24 +02:00
Michel Juillard b1ad31cf56 removed seeding random generators from clock in cmaes and metropolis hastings 2012-08-05 10:59:59 +02:00
Michel Juillard 6b3bd9dd0b penalty is now passed as a field of bayestopt_ (or BayesInfo) 2012-08-02 22:23:29 +02:00
Marco Ratto 47353a276b bug fix: the anlytic hessian was wrong with estimated stderr.
Manual cherry pick from Commit: 36853cb7615de41ca61ed92decf8bbe81ab40cba
2012-08-02 14:57:20 +02:00
Michel Juillard fe8ae2c565 adding doc *.m file for k_order_perturbation MEX 2012-08-02 11:25:52 +02:00
Michel Juillard 9e6c58572e fixing bug in pruning at order=3 2012-08-02 10:09:44 +02:00