commit
1638374b46
|
@ -345,6 +345,8 @@ options_.ramsey_policy = false;
|
|||
options_.instruments = {};
|
||||
options_.timeless = 0;
|
||||
options_.ramsey.maxit = 500;
|
||||
options_.ramsey.periods = 10000;
|
||||
options_.ramsey.drop = 1000;
|
||||
|
||||
% estimation
|
||||
options_.initial_period = NaN; %dates(1,1);
|
||||
|
|
|
@ -240,9 +240,9 @@ if options_.ramsey_policy
|
|||
i_exo_var = setdiff([1:M_.exo_nbr],find(diag(M_.Sigma_e) == 0));
|
||||
nxs = length(i_exo_var);
|
||||
chol_S = chol(M_.Sigma_e(i_exo_var,i_exo_var));
|
||||
exo_simul = zeros(M_.exo_nbr,options_.ramsey.nperiods);
|
||||
exo_simul = zeros(M_.exo_nbr,options_.ramsey.periods);
|
||||
if ~isempty(M_.Sigma_e)
|
||||
exo_simul(i_exo_var,:) = chol_S*randn(nxs,options_.ramsey.nperiods);
|
||||
exo_simul(i_exo_var,:) = chol_S*randn(nxs,options_.ramsey.periods);
|
||||
end
|
||||
yhat_start = zeros(M_.endo_nbr+1,1);
|
||||
[moment] = k_order_mean(options_.order, M_.nstatic, M_.npred, M_.nboth, M_.nfwrd+1, M_.exo_nbr, 1, options_.ramsey.drop, yhat_start, exo_simul, ysteady, dr);
|
||||
|
|
|
@ -31,9 +31,7 @@ end;
|
|||
|
||||
stoch_simul(order=6, irf=0);
|
||||
|
||||
options_.ramsey.nperiods = 10000;
|
||||
options_.ramsey.drop = 1000;
|
||||
evaluate_planner_objective;
|
||||
evaluate_planner_objective(periods=10000,drop=1000);
|
||||
|
||||
[W_dynpp] = k_order_welfare(oo_.dr, M_, options_);
|
||||
|
||||
|
|
Loading…
Reference in New Issue