f36247ceed
modifications to ramsey_policy code to allow for estimation
2011-01-26 21:24:47 +01:00
22d5bd16cf
removed useless auxiliary variables from call to kalman_transition_matrix()
2011-01-26 21:24:46 +01:00
bff8de837b
modifying dynare_resolve() calling sequence
2011-01-26 21:24:46 +01:00
a336799b96
corrected bug when mode_compute is a string
2011-01-26 21:24:46 +01:00
91e30d45f7
corrected bug related to ramsey_policy
2011-01-26 21:24:46 +01:00
55399a5dc5
Compatibility fix for Octave: under Octave, union() called on a column vector and a scalar returns a row vector, hence the need to add the 'rows' option
2011-01-24 15:34:02 -05:00
c71e0fbd6b
removing test about correlated data in initial checks before estimation: this condition is too restrictive.
2011-01-21 17:04:01 +01:00
f83d7db907
Fix of .mat extension with save for octave.
2011-01-21 11:28:28 +01:00
c5c99a80f6
Special fix for peculiar behavior of ls under octave/windows (+ use of ~ispc in place of isunix/ismac)
2011-01-21 11:26:24 +01:00
78c21916da
Fixed output for octave/unix octave/win
2011-01-21 11:25:30 +01:00
6c801be332
Removed catch ME instance for octave compatibility
2011-01-21 11:22:25 +01:00
4eb91731fc
Several fixes for octave compatibility:
...
- use of ispc;
- mat extension is save files;
- remove catch ME expressions;
- added printf message for octave runs
2011-01-21 11:21:32 +01:00
28a3b38cf4
For compatibility with octave (and for better coding) use of ~ispc in place of isunix/ismac
2011-01-21 11:18:29 +01:00
a3e4a7ff09
Eliminated annoying echoes on command window
2011-01-21 11:15:49 +01:00
1db24b3107
bug correction for diffuse smoother
2011-01-18 19:04:20 +01:00
c5310731d9
corrected bug affecting the smoother when called for a selected subset of variables
2011-01-18 19:04:20 +01:00
4cf358c895
Merge remote branch 'ferhat/master'
2011-01-18 17:24:42 +01:00
595675d333
Partial information changes that make adjustment for Octave and use rcond() < 1e-8 ..." to determine if a matrix is invertible plus some minor bug and formatting changes in dr1_PI.m
2011-01-17 20:40:41 +00:00
ce07223628
- The Temporary terms management with bytecode when the model is block decomposed and solve_algo<5 is now compatible with octave
...
- Gets rid of warning message during the compilation of bytecode
2011-01-14 19:24:18 +01:00
96b6f1bf05
- Corrects the following problem:
...
Octave BiCGStab algorithm involves a 0 division in case of a preconditioner equal to the LU decomposition of the A matrix (in a linear system of the form A.x = b).
- The solution:
Checks if the linear system is solved simply using: x_new = x_old + U \ (L \ x_old)
Ticket #11
2011-01-14 19:24:17 +01:00
708a062338
The Octave version of fsolve do not converge when it starts from the solution. Check if the initial guess is equal to the solution.
2011-01-14 19:24:17 +01:00
ac9415a943
- The bytecode.mex file has to be called when bytecode option is specified independently of the block option.
2011-01-14 10:39:34 +01:00
e2a1d77f6e
- added schur_statespace_transformation.m to factor it out of DsgeLikelihood.m and DsgeSmoother.m
...
and corrected a bug in the code
- added missing functions missing_DiffuseKalmanSmootherH1_Z.m (multivariate version) and
missing_DiffuseKalmanSmootherH3_Z.m (univariate version)
- use only these two versions of the Kalman smoother in DsgeSmoother.m
2011-01-13 21:50:26 +01:00
f7eb39f7bb
bug correction related to smoother
2011-01-13 21:45:09 +01:00
6ea053e47d
correcting bugs in diffuse Kalman smoother
2011-01-13 21:10:06 +01:00
2e0a36ae9f
Initialize estim_params_ in global_estimation.m. Necessary for smoothing models with no estimated parameters.
2011-01-13 21:10:06 +01:00
02fd6e668c
Add error message for non-implemented features of gamrnd.m
2011-01-13 12:28:44 +01:00
9e38d85053
Change the formula for empirical autocorrelations, so that the result always lies in [-1,1] (thanks to Johannes Pfeifer for pointing this)
2011-01-12 11:26:02 +01:00
17bc655159
Partial Information: Update with improved, generalised recursive PCL solution method
2011-01-11 21:57:57 +00:00
503b055df2
Fixed bugs in varlist_indices.m (thanks to Johannes Pfeifer for pointing them)
2011-01-10 13:11:25 +01:00
05b9c70ae8
removing command CALIB: it has not been working for a long time. Calibration is in fact a special case of method of moments and should be replaced by such a method.
2011-01-04 10:16:38 +01:00
256ff761e8
stochastic simulations: fix the number of generated data to be exactly options_.periods
...
added ./tests/simul
2011-01-02 16:55:15 +01:00
35765ee01c
- Matlab wrappers for the "print" option of bytecode
2010-12-31 17:03:19 +01:00
aab5b2a488
- Manages the temporary terms when the model is evaluated block by block
2010-12-31 16:22:24 +01:00
129303430e
The direction is not computed in same way in the dynamic case and in the static case
2010-12-31 16:20:35 +01:00
a882a25760
- length(args) is used instead of size(args, 2) to avoid to transpose column vectors
2010-12-31 16:16:29 +01:00
2dbeddd1a9
- Manages the temporary terms when the model is evaluated block by block
2010-12-31 16:00:49 +01:00
2707965da1
- extends the compatibility of dr1 with bytecode
2010-12-31 15:56:49 +01:00
7209f3638a
- consider a non sparse matrix for the Jacobian for compatibility reasons with dynare_solve
2010-12-31 15:53:28 +01:00
dcc46959b4
The default for options_.qz_criterium is now different for different contexts:
...
qz_criterium = 1+1e-6
- stoch_simul
- osr
- check
- dynare_estimation_1 with lik_init =2 or 3
qz_criterium = 1+1e-6 for dynare_estimation_1 with lik_init = 1
This correct a bug in estimation and computation of the posterior distribution of moments of endogenous variables when the priors permit roots too close to 1.
2010-12-23 08:34:13 +01:00
6bb8d41909
Revert "- added a test an a penalty in estimation (DsgeLikelihood.m) if, in a stationary model (lik_init==1), a particular parameter set generates unit roots."
...
There is a better way of dealing with occasional non-stationary models in estimation
This reverts commit 8c0fb55206
2010-12-22 09:40:39 +01:00
8c0fb55206
- added a test an a penalty in estimation (DsgeLikelihood.m) if, in a stationary model (lik_init==1), a particular parameter set generates unit roots.
...
- modified lyapunov_symm to return absolute value of unit roots in a third argument
2010-12-21 21:26:57 +01:00
ef56f87796
- fix for Octave printf
...
- missing semicolon
2010-12-20 15:55:46 +01:00
55e553f6dc
Missing semi-colons
2010-12-20 15:49:55 +01:00
ac0c0a935e
Cosmethics
2010-12-17 09:22:12 +01:00
b186eb7e84
Fixes for console mode in serial and parallel execution
2010-12-17 09:21:30 +01:00
191de78310
Introduced console mode monitoring of parallel runs;
...
Minor changes
2010-12-17 09:20:30 +01:00
916caf2e4b
Almost entirely rewritten:
...
1) got rid of recursive form;
2) better help
3) simplified structure and split between errors and warnings.
2010-12-17 09:16:30 +01:00
c2f0037e55
small addendum to help
2010-12-17 09:14:35 +01:00
d17faaf876
cosmethics
2010-12-17 09:12:54 +01:00
41c59778cd
For slave processes on the local machine, impose the same dynare path as the master.
2010-12-17 09:11:11 +01:00
c0cb3dd4fb
Compatibility fix for Octave: under Octave, union() called on a column vector and a scalar returns a row vector, hence the need to add the 'rows' option
2010-12-16 16:14:23 +01:00
b63590712c
fixed bug in argument number check and removed unused variable
2010-12-14 18:42:29 +01:00
e04a07c771
changed variable name to make clear that, for multinormal variables, we are using the Cholesky decomposition of the covariance matrix.
2010-12-11 18:09:01 +01:00
c1cb452e3d
correcting headers of all Kalman filter functions: returns MINUS loglikelihood
2010-12-10 22:31:30 +01:00
7a09356423
DsgeLikelihood.m header correction: the function returns MINUS log posterior
2010-12-10 22:28:06 +01:00
dfbd10a524
added choice of parameter set to shock_decomposition. This corrects a bug when shock_decomposition is requested after Metropolis iterations.
2010-12-08 13:02:59 +01:00
ea86e8cefc
adding function evaluate_smoother that calls DsgeSmoother with a given type of parameter set
2010-12-08 13:02:59 +01:00
9b166fdb08
NaN parameters warning: remove \n, it is not interpreted in a warning context (at least under Octave)
2010-12-08 11:32:45 +01:00
03fac30723
remove get_innovation_contemporaneous_impact.m that isn't called by the preprocessor or by any other function
2010-12-07 13:43:17 +01:00
3af62c4ac9
Fixed bug in resid.m (order matters in "or" statement, because of short-circuit: "ind" is not defined if "istag" is false)
2010-12-06 18:24:08 +01:00
0d1fd700c7
Posterior distribution of conditional variance decomposition: storing the steps chosen by the user in a separate field.
2010-12-05 10:47:31 +01:00
53d391785d
corrected bug in get_name_of_the_last_mh_file.m
2010-12-01 09:57:08 +01:00
d4a908dda1
Removed useless input argument.
2010-11-30 10:02:37 +01:00
116bf69fbf
Seed: fix bug under Octave and MATLAB < 7.7
2010-11-29 11:23:59 +01:00
1fc8bbd6d3
bugs correction in computation of posterior moments for
...
-conditional variance decomposition
-hpdsup
-moments with no variance in their posterior distribution
modification of computation of conditional variance decomposition
2010-11-28 11:03:40 +01:00
73fb71b9c8
Removed references to options_.simul_seed.
2010-11-27 00:19:22 +01:00
04469ddce1
- Add variable input arguments
...
- Correct the error condition in A_times_B_kronecker_C
2010-11-26 17:46:02 +01:00
db54c2d76e
Added a routine for setting the random number generator and the seed, independantly of the version of matlab. Removed the calls to rand('state',x) in simult.m and added a test file. Still need to adapt the routines related to mcmc and swz.
2010-11-26 17:16:12 +01:00
618f72bc91
resid displays the name of the equations if they are specified in equations tags.
2010-11-25 17:31:11 +01:00
ffaec7ba89
- First elements needed for stochastic simulation of a block decomposed model
2010-11-25 16:11:27 +01:00
61a251535f
Correction of bugs in check command with block decomposition and
...
Add the missing dr_block m-file
2010-11-25 16:10:46 +01:00
b356c5e58e
add_auxiliary_variables_to_steadystate is compatible with bytecode option
...
+ minor corrections in steady_.m
2010-11-25 16:00:54 +01:00
1626b0aa76
Added specialized code for purely forward models.
2010-11-25 14:32:33 +01:00
1eccfb311a
Remove obsolete files
2010-11-25 13:10:10 +01:00
6c22d398e6
Perfect foresight simulation: clean up handling of warnings
2010-11-24 19:01:46 +01:00
9d4385da5d
Put back M_.maximum_lag.
2010-11-19 16:30:55 +01:00
ccf778e63a
Efficiency improvements in first order stochastic simulations (Simulate data in deviation to the steady state and then,
...
after the loop, add the steady state using bsxfun routine, which is parallelized in recent version of matlab). Removed
use of M_.maximum_lag.
2010-11-18 14:55:07 +01:00
e614227c3f
Added warning message when some of the deep parameters are not initialized.
2010-11-17 17:09:39 +01:00
d60c32e6a2
Fixed bug in dealing with _steadystate functions that modify deep parameters. See trac ticket #103 .
2010-11-15 14:35:21 +01:00
5c709e4733
Changed set_dynare_threads.m. The number of threads in parallelized mex files to be used if dynare is built with the openmp flag
...
(--with-openmp) is not passed by an environment variable anymore. The function set_dynare_threads changes the default value of the
number of threads (default is 1) in the options_.threads structure. Changed calls to sparse_hessian_times_B_kronecker_C and
A_times_B_kronecker_C dlls accordingly.
2010-11-15 13:06:39 +01:00
37f14e9bc9
Removed unused routines for (diffuse) kalman filter evaluations.
2010-11-12 17:20:02 +01:00
31b734e9b3
Do not print residuals of the static equations for the equations of the auxiliary variables.
2010-11-04 14:19:28 +01:00
88dc361ba1
unix info added
2010-10-29 12:30:41 +02:00
5763721a71
almost complete version of testing routine
2010-10-29 12:28:44 +02:00
23ba760c3b
Implement solve_algo=0 for Octave ( closes #144 )
2010-10-28 12:07:50 +02:00
67c28ea5f3
PosteriorIRF: added missing ".mat" extension (necessary for Octave)
2010-10-28 10:39:53 +02:00
0730213059
Fix end of line convention and copyright years
2010-10-27 18:35:42 +02:00
1aceb504fd
Check command is now compatible with block and bytecode options
2010-10-27 15:34:48 +02:00
9cebf270e5
- correction of informations displayed in Solve_one_boundary
2010-10-27 15:30:00 +02:00
58a8988ef3
Moved another parallel routine previously left out.
2010-10-27 11:27:10 +02:00
b5ebff4ff4
bug fix for parallel execution
2010-10-27 11:24:23 +02:00
7ef63899b1
add paralle_info input to Analize;
...
fix spelling of function names
2010-10-27 11:05:39 +02:00
5915f8a059
Fixed misspelled routines!
2010-10-27 11:03:21 +02:00
0a38f8ac3b
Initialize here default of MatlabOctavePath;
...
implies using global options_ in place of local argument
2010-10-25 09:53:09 +02:00
fe8724d26c
Merge remote branch 'ratto/master'
2010-10-22 17:22:37 +02:00
5059c48460
solve_algo=7 (GMRES) is forbidden under Octave
2010-10-22 16:52:02 +02:00
39718c8645
- extends steady algorithms to solve_algo=5, 6, 7 and 8 for a block decomposed model without bytcode
...
- solve_algo is reordered =
* 0: fsolve
* 1: solve1
* 2, 4: solve1 + block decomposition
* 3: csolve
* 5: bytecode own solver (use Gaussian elimination + sparse matrix)
* 6: LU decomposition with UMFPack (method handling sparse matrix in Matlab)
* 7: GMRES
* 8: BiCGStab
2010-10-22 16:49:47 +02:00
f336a68ecc
Update tests for values of solve_algo and stack_solve_algo
2010-10-22 13:00:49 +02:00
a46969f9fb
Deleted buggy fclose instruction
2010-10-22 11:40:29 +02:00
Michel Juillard
Sébastien Villemot
Marco Ratto
George Perendia
Ferhat Mihoubi
Houtan Bastani
Stéphane Adjemian (Charybdis)