# Estimation error: in computing likelihood for initial parameter

**URL:** <https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839>\
**Category:** ML/Bayesian estimation\
**Created:** [28 October 2020 18:34 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839 "2020-10-28T18:34:59Z")\
**Posts on this page:** 7\
**Page:** 1

<div class="post-metadata">

**Author:** ![RodrigoTM](https://forum.dynare.org/letter_avatar_proxy/v4/letter/r/848f3c/32.png) [@RodrigoTM](https://forum.dynare.org/u/RodrigoTM)\
**Post date:** [28 October 2020 18:34 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/1 "2020-10-28T18:34:59Z")

</div>

Hi. I’m trying to estimate the Gali-Monacelli (2005) open small economy based on Pfeifer’s code, but I’m receiving this message:

`Error in computing likelihood for initial parameter values`

ESTIMATION\_CHECKS: There was an error in computing the likelihood for initial parameter values.  
ESTIMATION\_CHECKS: If this is not a problem with the setting of options (check the error message below),  
ESTIMATION\_CHECKS: you should try using the calibrated version of the model as starting values. To do  
ESTIMATION\_CHECKS: this, add an empty estimated\_params\_init-block with use\_calibration option immediately before the estimation  
ESTIMATION\_CHECKS: command (and after the estimated\_params-block so that it does not get overwritten):

`Error using print_info (line 32) Blanchard & Kahn conditions are not satisfied: no stable equilibrium.`

I’ve already checked the equations and try to used the use\_calibration in estimated\_params\_init, but it didn’t solve it.

Would you help me, please?[Gali\_orig\_Est.zip](https://forum.dynare.org/uploads/short-url/lJMHDEqRz4DmBHMFUuCOmEQfJSB.zip) (14.2 KB)

---

<div class="post-metadata">

**Author:** ![jpfeifer](https://forum.dynare.org/user_avatar/forum.dynare.org/jpfeifer/32/5044_2.png) [@jpfeifer](https://forum.dynare.org/u/jpfeifer)\
**Post date:** [29 October 2020 08:26 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/2 "2020-10-29T08:26:03Z")

</div>

Please see

> [@Schmitt-Grohe & Uribe (2018) replication](https://forum.dynare.org/t/schmitt-grohe-uribe-2018-replication/13832/6):
>
> Your model has a unit root. You need to use diffuse\_filter in estimation

---

<div class="post-metadata">

**Author:** ![RodrigoTM](https://forum.dynare.org/letter_avatar_proxy/v4/letter/r/848f3c/32.png) [@RodrigoTM](https://forum.dynare.org/u/RodrigoTM)\
**Post date:** [30 October 2020 18:34 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/3 "2020-10-30T18:34:58Z")

</div>

Thank you, professor Pfeifer.  
I used the diffuse\_filter and now it is working. I also read some other topics in the forum about it but I have some more questions:

1 - If I estimate this a model with an unit root, there are infinite steady states. Which one Dynare will give me as a solution depends on the initial values in estimated\_params\_init?  
2 - Does the inexistence of a unit root assure an unique steady state?  
3 - With a finite number of steady state (say 2), will Dynare give me one of them or both?  
4 - In that case above, will the impulse response function be relative to the same steady state? Or is it possible to see same transition path between two differents steady states (maybe using another feature of Dynare)?

Finally, if you have some literature about these questions I will really appreciate it.  
Thank you again.

---

<div class="post-metadata">

**Author:** ![jpfeifer](https://forum.dynare.org/user_avatar/forum.dynare.org/jpfeifer/32/5044_2.png) [@jpfeifer](https://forum.dynare.org/u/jpfeifer)\
**Post date:** [31 October 2020 14:50 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/4 "2020-10-31T14:50:12Z")

</div>

1. If there are infinitely many steady states, you should select the correct one using a `steady_state_model`-block. Note that the estimated parameters have nothing to do with the issue. For any given parameter sets, there is a continuum of steady states.
2. No, it does not. Take the standard Solow model. It has two steady states. The normal one we always consider and the trivial one where everything is zero.
3. If you let Dynare solve for the steady state numerically, Dynare’s solver will (hopefully) converge to one steady state. As a local solver is used, it’s typically the nearest one. But solving can also fail because the unit root may imply a singular derivative matrix singular and solvers can then fail.
4. The decision rules will be based on an approximation around the steady state that was selected/found. If you want control over this, again see point 1.

---

<div class="post-metadata">

**Author:** ![RodrigoTM](https://forum.dynare.org/letter_avatar_proxy/v4/letter/r/848f3c/32.png) [@RodrigoTM](https://forum.dynare.org/u/RodrigoTM)\
**Post date:** [31 October 2020 15:40 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/5 "2020-10-31T15:40:31Z")

</div>

That’s perfect. Everything’s clear.  
Thank you so much, professor.

---

<div class="post-metadata">

**Author:** ![JackChou](https://forum.dynare.org/letter_avatar_proxy/v4/letter/j/e47774/32.png) [@JackChou](https://forum.dynare.org/u/JackChou)\
**Post date:** [15 November 2025 12:08 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/6 "2025-11-15T12:08:08Z")

</div>

> [@jpfeifer](#):
>
> steady\_state

How if use

% steady\_state\_model;

steady;

Because I need to use functions, so i write the steadystate.m to calculate the steady. But the error occurs.

---

<div class="post-metadata">

**Author:** ![jpfeifer](https://forum.dynare.org/user_avatar/forum.dynare.org/jpfeifer/32/5044_2.png) [@jpfeifer](https://forum.dynare.org/u/jpfeifer)\
**Post date:** [17 November 2025 09:36 UTC](https://forum.dynare.org/t/estimation-error-in-computing-likelihood-for-initial-parameter/16839/7 "2025-11-17T09:36:25Z")

</div>

Without codes or details it is impossible to tell.
