Dear Friends,
I just completed the first tests with estimation of DYNARE with analytic
scores and hessian and everything seems to work.
I have implemented it so far for stationary models both for univariate
and multivariate KF.
Before committing the changes, I would like to anticipate you some
decisions I took: I particular I had to change the list of output of
dsge_likelihood (and similarly for consistency of dsgevar_likelihood).
We had discussed this a bit some time ago and I envisaged the
possibility to keep the current list simply changing the type of output
LIK = scalar for numerical derivs
LIK= cell = {LIK, GRAD, HESS}
Stephan in particular supported this last option.
However, the standard form of objective functions in matlab is
[fval, grad, hess] = func( ... )
and this form would allow, for example, a direct use of mode_compute=1 or 3.
This means I have adapted all fval calls, e.g. in csminwel1, so that now
the exit_flag is 4th output of dsge_likelihood.
This will also allow to use Sims optimizer with analytic derivatives
The same for mode_compute=5.
By the way, this will also allow to eliminate the annoying
dsge_likelihood_hh. I have also thought of a way to adapt mode_compute=5
for dsge-var as well.
what do you think?
best
Marco
--
Marco Ratto,
Joint Research Centre
The European Commission,
TP 361, 21027 ISPRA(VA), ITALY
Tel: +39 0332 78 5639
Fax: +39 0332 78 5733
marco.ratto(a)jrc.ec.europa.eu