# Impossible to find steady state in Ireland 2003

**URL:** <https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557>\
**Category:** Steady state\
**Created:** [26 February 2021 13:16 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557 "2021-02-26T13:16:03Z")\
**Posts on this page:** 20\
**Page:** 1

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [26 February 2021 13:16 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/1 "2021-02-26T13:16:03Z")

</div>

Dear professor,

I’m trying to find a steady state for Ireland 2003 model, however I get the following message:

“error: Impossible to find the steady state (the sum of square residuals of the static equations i s 29704268567552244.0000). Either the model doesn’t have a steady state, there are an infinity of steady states, or the guess values are too far from the solution”

Equation no. 8 has a bit larger sum of square residuals. I tried changing the state state values for variables in that equation, however, I can only decrease that sum of residuals by increasing other one.

If I don’t do exp(var) substitution, Dynare says that Blanchard-Kahn conditions aren’t met.

Many thanks

[Model4.mod](https://forum.dynare.org/uploads/short-url/jGEFIRizTwdgklUxaP6b54sZbl6.mod) (1.9 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:** [26 February 2021 14:53 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/2 "2021-02-26T14:53:14Z")

</div>

1. Only do the exp()-substitution once the model works.
2. exp()-substitution cannot alter the BK-conditions. Therefore, check your timing. For example, the timing of capital in the production function is wrong.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [26 February 2021 20:09 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/3 "2021-02-26T20:09:12Z")

</div>

Now the model with adjusted timing and in the non-exp form still says:

“error: Blanchard & Kahn conditions are not satisfied: indeterminacy.”

Is it possible that there is a problem with the notation? In the equations, I have the same variable in both steady-state form and in the form of relative deviations from the steady-state. For the latter I use normal variables x, y, z etc. and for the former I introduce new parameters with the values determined in the beginning.

I’m attaching the new model.

[Model.mod](https://forum.dynare.org/uploads/short-url/Tb0iNgvxvRgA5r8543gb0vdCmi.mod) (1.6 KB)

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [7 March 2021 17:15 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/4 "2021-03-07T17:15:47Z")

</div>

Hello again and thanks for reading my post to anyone in this community. I have corrected the time periods in the model equations, dropped exp(var) form and substituted the variables in the steady state form with the operator steady\_state(var). Using this model in Dynare, I obtain the message:

> error: Impossible to find the steady state (the sum of square residuals of the static equations is 188395631132217704448.0000). Either the model doesn’t have a steady state, there are an infinity of steady states, or the guess values are too far from the solution

Any suggestion is welcome.

[Model.mod](https://forum.dynare.org/uploads/short-url/6ie5LsJcBQYgqw5Bm1jfa0Wp8TP.mod) (1.7 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:** [8 March 2021 14:34 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/5 "2021-03-08T14:34:20Z")

</div>

You have various variables that must not be 0 in steady state. Thus, you need to set explicit initial values.  
Also,

```auto
ln(x)=rhox*x(-1)+epsilonx;
ln(a)=rhoa*a(-1)+epsilona;
ln(e)=rhoe*e(-1)+epsilone;
ln(z)=rhoz*z(-1)+epsilonz;
ln(v)=rhov*v(-1)+epsilonv;

```

is missing the ln on the right.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [3 April 2021 19:27 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/6 "2021-04-03T19:27:22Z")

</div>

Hello, I have set initial values for variables, according to the rule I describe below. I get the error message

“error: The steady state has NaNs or Inf.”

My sample includes data for Slovenia, Austria and Italy in 1991-2019. All estimates are on the annual basis.

_x, a, v_ = 1, _x_ is a shock to marginal efficiency of investment, _a_ is IS shock and _e_ is a money demand shock, values are assumed by Röhe (2012).  
_z, v_ = 0, _z_ is a technological shock and _v_ is a monetary policy shock. Values are assumed by Röhe (2012).  
_y_ = 26933 is the GDP p.c., average of the sample data for 3 countries in 28 years,  
_k_=684000 is households’ capital supply, the data is average firm’s debt to depositary institutions for Slovenia 1991-2019  
_in_=1946.8 are goods that are purchased by households and not consumed or investment, the data is the sample average for saving as a product of salary and savings rate  
_c_=12567.3 is the consumption, the estimate is the sample average  
_d_=2 is the nominal dividend, the estimate is average yield at DAX in 2020  
_m_=735.1 is money brought from previous time period, the estimate is the sample average for money aggregate M3  
_w_=14637.4 is the nominal wage, the estimate is the sample average  
_h_=1834 is the number of hours worked, the estimate is average for Slovenia 1991-2019  
_r_ =8.8 is the gross nominal interest rate between time periods, the esitmate is the sample average  
_q_=0.09 is the nominal rental rate for capital, the estimate is for the US in 2020, computed as sum of the equity risk premium (6.0 %) and the reaffirmed normalized risk-free rate (3.0 %)  
_n_=14.7 is the GDP per hours worked (y/h)  
\*tau=0.0068 is the M3 growth rate, the estimate 0.68 % is the sample average  
_lambd_=-0.41 is the marginal utility of additional euro of profits, estimated with the marginal utility of income, this estimate is for the US in 1991-1992  
\*ksi=0 this variable has no explicit meaning in Röhe (2012), it’s from firms’ Lagrangian optimization, the estimate I left 0  
_pi_=79.1 is the inflation rate, the estimate is the sample average [Model.mod](https://forum.dynare.org/uploads/short-url/8LW02gqDtreZjpJTgdIdzfaeiJE.mod) (1.9 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:** [5 April 2021 08:01 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/7 "2021-04-05T08:01:45Z")

</div>

You cannot simply plug in data from the real world. That would require an appropriate scaling of most variables, e.g. via a productivity factor.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [10 August 2021 17:52 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/8 "2021-08-10T17:52:05Z")

</div>

Hello, I apologize for delay. I substituted the initial values of variables by the averages of the variables for the three countries, weighted by productivity. I used the GDP per hour worked as a % of the USA. For example, if initial value of consumption before was computed as (1+3+2)/3, now it was computed as (1x0.6115+3x0.7286+2x0.8378)/3, the weights corresponding to the productivity computed as GDP per hour worked relative to the USA for Slovenia, Italy and Austria.

The Dynare still reports: error: The steady state has NaNs or Inf. Should I have multiplied the productivity other way, e.g. in the model equations or use multiple initial values for the used countries? Thank you in advance.  
[Model1.mod](https://forum.dynare.org/uploads/short-url/iGfwV0MkVukDgRwWKhcYkTgZ1OU.mod) (2.0 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:** [13 August 2021 08:11 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/9 "2021-08-13T08:11:54Z")

</div>

For example

```auto
ln(e)=(1-rhoe)*ln(steady_state(e))+rhoe*e(-1)+epsilone;
ln(z)=(1-rhoz)*ln(steady_state(z))+rhoz*z(-1)+epsilonz;

```

prevents computing the steady state numerically as you are trying to endogenously compute a steady state that needs to be specified by you. Also, you specified initial values of 0 for `e,z`, but you are taking their logs.  
Also, see [An infinity of steady states with Taylor rules](https://forum.dynare.org/t/an-infinity-of-steady-states-with-taylor-rules/4026)

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [14 August 2021 09:21 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/10 "2021-08-14T09:21:32Z")

</div>

I have corrected the initial values for e and z. Regarding the steady\_state operator, my model includes variables in certain time period and in steady state, e.g. e(t) and e (you can see it on [https://d-nb.info/1026165547/34](https://d-nb.info/1026165547/34) p. 62). How would I denote both? I saw the possibility of writing R and Rbar/R\_ss on the link you sent, but even setting the Rbar to some initial value would still leave R as another variable, requiring some initial value too.

[Model.mod](https://forum.dynare.org/uploads/short-url/w0hiTZuMjg2rRJuNEmgQVRa1pmw.mod) (2.0 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:** [14 August 2021 18:00 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/11 "2021-08-14T18:00:37Z")

</div>

None of this changes the fact that these objects cannot be endogenously computed. Also, in t hat PDF the steady state is explicitly provided in C.3.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [15 August 2021 10:31 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/12 "2021-08-15T10:31:50Z")

</div>

OK, so steady state operator goes away. I’m not sure how to use the steady state in Mr Röhe’s paper in my case, where I have to state initial values and equations with steady state variables. Can I introduce the steady state equations inside the model block? I have introduced additional variables y\_ss, pi\_ss etc. and set their value to the initial value in the model, as in [The steadystate file did not compute the steady state](https://forum.dynare.org/t/the-steadystate-file-did-not-compute-the-steady-state/13379).

I still get the message; error: Impossible to find the steady state (the sum of square residuals of the static equations is 242604047897428410205353479174390707626049099439580515812856364670248749087860797378127724544.0000).

[Model.mod](https://forum.dynare.org/uploads/short-url/tAtYQ3N2rdiYFKmAsHNdfZF8T8o.mod) (1.9 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:** [16 August 2021 06:37 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/13 "2021-08-16T06:37:57Z")

</div>

Where did you get the hard-coded numbers from. With

```auto
gamm=0.042;
e=4.02;

```

you will get `e^(1/gamm)= 2.4341e+14;`

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [16 August 2021 09:31 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/14 "2021-08-16T09:31:09Z")

</div>

Gamm or gamma is interest rate elasticity of money demand, it’s from that paper p. 50 (average of estimated values for Germany, France and Spain).

E is a money demand shock, the same paper suggests that ln(e) is 1.39, so I computed e=4.02. The empirical data in the same table on p. 50 suggests it’s 3.81.

---

<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:** [16 August 2021 14:12 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/15 "2021-08-16T14:12:04Z")

</div>

You need to provide more details. Where does the process come from? How is it supposed to work? Where exactly do the numbers come from? IS there are meaning behind the exactl value of ln(e)? Scaling that works in the data does not necessarily work in a model.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [17 August 2021 14:06 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/16 "2021-08-17T14:06:09Z")

</div>

Gamma is the absolute value of the interest rate elasticity of money demand. It enters the expected utility function on page 43 in that pdf and is greater then zero.

E enters the same utility function and is characterised by the autoregressive process as on the same page (equation 3.4), determined by ln(e) and money demand shock parameters 0\<rho(e)\<1 and epsilon, which is normally distributed with mean 0 and variance sigma(e) squared.

Both are used to characterise expected utility function of households, which is subject to a budget constraint.

The numbers are from empirical study of Germany, France, Spain and Italy for 1980-2008 period, they are maximum likelihood estimates obtained using Kalman filter. I don’t know for any meaning behind ln(e) number. For ln(z) it is said

“Concerning ln(z), we chose a prior mean of 8.13 for the euro area and 9.18 for the US, both with a standard deviation of 1. We set the prior means of ln(z) so that the steady state values of ct and it in the model match the respective average values of consumption and investment in the data.”

and

“Finally, we find the estimates of ln(z) to be almost equal across all model specifications within each region. Note that ln(z) has no impact on the dynamics of the model specifications, since it only serves to determine steady states (see section 4.5).”

ln(e) isn’t explicitly mentioned, except for the specification mentioned above and the numerical estimates.

So it may be best to find some other empirical estimates that would work better with my model?

---

<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:** [18 August 2021 08:01 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/17 "2021-08-18T08:01:54Z")

</div>

Are you trying to replicate that exact paper?

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [18 August 2021 08:53 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/18 "2021-08-18T08:53:22Z")

</div>

No, sorry, I want to do an estimation for three other countries. I just quoted the paper and page because I’m not sure how to write an equation with superscripts and underscripts here. I plan to do various parameter stability tests other than that paper’s authors, when I get the model to work.

---

<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:** [18 August 2021 08:54 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/19 "2021-08-18T08:54:38Z")

</div>

Usually, the starting point is an exact replication. Once you get that to work, you start modifying things. Otherwise, it’s often impossible to tell where things went wrong.

---

<div class="post-metadata">

**Author:** ![redbanana](https://forum.dynare.org/user_avatar/forum.dynare.org/redbanana/32/12120_2.png) [@redbanana](https://forum.dynare.org/u/redbanana)\
**Post date:** [18 August 2021 08:57 UTC](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557/20 "2021-08-18T08:57:28Z")

</div>

OK, in my case, would it be possible to include steady\_state\_model block instead of initial values, since I know the steady state?

[Next page](https://forum.dynare.org/t/impossible-to-find-steady-state-in-ireland-2003/17557.md?page=2)