Commit 991b4c2a by Johannes Pfeifer Committed by Stéphane Adjemian

### Make empirical moments independent of simul_replic by outputting first instead of later replication

```Closes #1487

(cherry picked from commit f0b9f453)```
parent bfde0280
 function [y_,DynareResults] =simult(y0, dr,DynareModel,DynareOptions,DynareResults) function [y_out,DynareResults] =simult(y0, dr,DynareModel,DynareOptions,DynareResults) % Simulate a DSGE model (perturbation approach). %@info: ... ... @@ -25,7 +25,7 @@ function [y_,DynareResults] =simult(y0, dr,DynareModel,DynareOptions,DynareResul %! @strong{Outputs} %! @sp 1 %! @table @ @var %! @item y_ %! @item y_out %! Matrix of doubles, simulated time series for all the endogenous variables (one per row). %! @item DynareResults %! Matlab's structure gathering the results (see @ref{oo_}). ... ... @@ -88,6 +88,9 @@ for i=1:replic if replic > 1 fwrite(fh,y_,'float64'); end if i==1 y_out=y_; end end if replic > 1 ... ...
 ... ... @@ -41,3 +41,11 @@ var u = 0.009^2; end; stoch_simul(periods=2000, drop=200); %% test that simul_replic does not affect simulated moments moments_temp=oo_.var; set_dynare_seed('default'); stoch_simul(periods=2000, drop=200,simul_replic=2); if ~isequal(moments_temp,oo_.var) error('Simul_replic affects simulated moments') end \ No newline at end of file
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment