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@jrc.ec.europa.eu