commit
356999768c
|
@ -168,13 +168,20 @@ if options_.ramsey_policy
|
||||||
options_.noprint = 0;
|
options_.noprint = 0;
|
||||||
end
|
end
|
||||||
|
|
||||||
|
if any(isnan(oo_.mean)) || any(isnan(oo_.mean))
|
||||||
|
fprintf('evaluate_planner_objective: encountered NaN moments in the endogenous variables often associated\n')
|
||||||
|
fprintf('evaluate_planner_objective: with either non-stationary variables or singularity due e.g. including\n')
|
||||||
|
fprintf('evaluate_planner_objective: the planner objective function (or additive parts of it) in the model.\n')
|
||||||
|
fprintf('evaluate_planner_objective: I will replace the NaN with a large number, but tread carefully,\n')
|
||||||
|
fprintf('evaluate_planner_objective: check your model, and watch out for strange results.\n')
|
||||||
|
end
|
||||||
oo_.mean(isnan(oo_.mean)) = options_.huge_number;
|
oo_.mean(isnan(oo_.mean)) = options_.huge_number;
|
||||||
oo_.var(isnan(oo_.var)) = options_.huge_number;
|
oo_.var(isnan(oo_.var)) = options_.huge_number;
|
||||||
|
|
||||||
Ey = oo_.mean;
|
Ey = oo_.mean;
|
||||||
Eyhat = Ey - ys(dr.order_var(nstatic+(1:nspred)));
|
Eyhat = Ey - ys(dr.order_var(nstatic+(1:nspred)));
|
||||||
|
|
||||||
Eyhatyhat = oo_.var(:)
|
Eyhatyhat = oo_.var(:);
|
||||||
Euu = M_.Sigma_e(:);
|
Euu = M_.Sigma_e(:);
|
||||||
|
|
||||||
EU = U + Uy*gy*Eyhat + 0.5*((Uyygygy + Uy*gyy)*Eyhatyhat + (Uyygugu + Uy*guu)*Euu + Uy*gss);
|
EU = U + Uy*gy*Eyhat + 0.5*((Uyygygy + Uy*gyy)*Eyhatyhat + (Uyygugu + Uy*guu)*Euu + Uy*gss);
|
||||||
|
|
Loading…
Reference in New Issue