# Very high values of theoretical moments

**URL:** <https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514>\
**Category:** Estimation\
**Created:** [5 August 2017 13:47 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514 "2017-08-05T13:47:53Z")\
**Posts on this page:** 8\
**Page:** 1

<div class="post-metadata">

**Author:** ![auro](https://forum.dynare.org/letter_avatar_proxy/v4/letter/a/ecd19e/32.png) [@auro](https://forum.dynare.org/u/auro)\
**Post date:** [5 August 2017 13:47 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/1 "2017-08-05T13:47:53Z")

</div>

Hello. I am running a basic calibrated log-linearized DSGE model, and while the code is running fine and the IRF’s look broadly as expected - the theoretical standard deviation and variance of the log-linearized variables are grotesquely large!

For instance, the regular log linearised consumption variable, c\_hat has a variance of 2245 and the inflation log linearized variable, pie\_hat has a variance of 302!

Is this something I should be worried about? And if yes, how should I be ideally dealing with this? I tried changing the variance of the shocks from 1 to 0.01, but the variances/standard deviations still remain large.

Any help would be much appreciated.

Thanks,

A

---

<div class="post-metadata">

**Author:** ![Ecofai](https://forum.dynare.org/letter_avatar_proxy/v4/letter/e/f14d63/32.png) [@Ecofai](https://forum.dynare.org/u/Ecofai)\
**Post date:** [5 August 2017 15:06 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/2 "2017-08-05T15:06:12Z")

</div>

I faced the same problem, to handle this I didn’t use the theoretical given by Dynare, Instead I calculated them myself using the data available in : oo\_.SmoothedVariables , and the values found are different from the ones given by Dynare.

---

<div class="post-metadata">

**Author:** ![stepan-a](https://forum.dynare.org/user_avatar/forum.dynare.org/stepan-a/32/6652_2.png) [@stepan-a](https://forum.dynare.org/u/stepan-a)\
**Post date:** [5 August 2017 16:40 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/3 "2017-08-05T16:40:22Z")

</div>

Hi,

A reason for the large theoretical second order moments may be that you have at least one eigenvalue close to one (unit root). The `check` command displays the generalized eigenvalues, the output of this command may confirm this possible explanation.

It is not possible to compute the theoretical moments from the smoothed (or filtered) variables. It’s a different concept. Think of a simple AR(1), the sample and theoretical moments are always different for finite sample size.

Best,  
Stéphane.

---

<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 August 2017 17:02 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/4 "2017-08-05T17:02:49Z")

</div>

@auro Please provide the file to replicate the issue.

---

<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:** [10 August 2017 12:23 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/6 "2017-08-10T12:23:45Z")

</div>

Hi, I cannot find anything immediately suspicious in your file. Looking that the variance decomposition, there are three shocks that drive almost all of the variance of the respective variables. You should focus on those. Maybe you have a scaling issues here.

---

<div class="post-metadata">

**Author:** ![auro](https://forum.dynare.org/letter_avatar_proxy/v4/letter/a/ecd19e/32.png) [@auro](https://forum.dynare.org/u/auro)\
**Post date:** [10 August 2017 15:50 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/7 "2017-08-10T15:50:53Z")

</div>

Hi - thanks for the feedback. I identified the errant shocks - but I am unsure what to do about them to tackle the ‘scaling issue’.

I tried changing the standard error of **all** the shocks to 0.01 from 1 originally, and the theoretical moments now look sensible. However, as a collateral consequence, the IRF’s now have a lot of 10^-3 or 10^-4 on the y-axes for the key variables. This is especially problematic, when I am talking about fiscal multipliers in my thesis, wherein everything looks uncomfortably small compared to when I used stderr as 1!

So the question is: (a) Is this the right way to tackle the issue of large theoretical moments - to use stderr 0.01 instead of 1 for all the shocks? And (b) should I be overtly worried about the resulting impact of impulse responses that are 10^-3 or so in scale.

Thanks once again for your feedback.

A

---

<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:** [10 August 2017 17:23 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/8 "2017-08-10T17:23:46Z")

</div>

You may want to consult “Remark 10 (Scaling With a Factor 100)” in See “Remark 10 (Scaling With a Factor 100)” in Pfeifer (2013): [“A Guide to Specifying Observation Equations for the Estimation of DSGE Models”](https://drive.google.com/file/d/1r89OU5OE3CBa6tOlj6l3hNVWEaRH5Anv/view?usp=sharing)

The tricky issues is finding out why some shocks generate IRFs that are bigger by a factor of about 100 for a presumably similar shock size. There are two possibilities here:

1. For some reason, the shock size is actually not comparable (e.g. interest rates were already scaled with 100)
2. For some reason, the propagation mechanism to some shocks is wrong.

---

<div class="post-metadata">

**Author:** ![auro](https://forum.dynare.org/letter_avatar_proxy/v4/letter/a/ecd19e/32.png) [@auro](https://forum.dynare.org/u/auro)\
**Post date:** [20 August 2017 13:46 UTC](https://forum.dynare.org/t/very-high-values-of-theoretical-moments/10514/9 "2017-08-20T13:46:14Z")

</div>

Hi…thanks a lot for the feedback. I am unfortunately still fuzzy on troubleshooting this - can you possibly elaborate on how the ‘propagation mechanism’ for shocks can go wrong.

On the wider point - if I am running a calibrated model, does it really matter if I do not resolve this issue, and the theoretical moments continue to have high values? I understand that it is not ideal, but my objective is really to study the IRF’s, and not necessarily take the model to data at this point. And all the associated results (like fiscal multipliers for instance), seem unaffected if I have stderr as 1 (when moments explode in size) or 0.01 (when moments are more under control).

Thanks,

A
