Dear All,

I have a question concerning the simulation of the model using simult_ and stoch_simul.

In particular, i wanna calculate the ergodic mean of a variable “y”. For this I consider two possibilities:

rng(1000)

stoch_simul(order=3,periods=10000,drop=8000,replic=0,irf=0,noprint,nomoments,nofunctions,nocorr,pruning);

and then just take mean(y).

Alternatively what I do is:

T=10000;

ex_ = randn(T,2); %I have two shocks in the model. shocks in the model are written as sigma_eps*eps

rng(1000)

simulation=simult_(oo_.steady_state,oo_.dr,ex_,options_.order);

ergo=mean(simulation(:,8000:end));

endonames=cellstr(M_.endo_names); %make them strings to use strcmp

index=find(strcmp(‘y’,endonames));

y_ergo=mean(ergo(index));

strangely the results are very different. With the first procedure I get 0.9122 and with the second I get -7.8219.

Do you have an explanation for this?