Commit ef797eec authored by Johannes Pfeifer 's avatar Johannes Pfeifer Committed by Stéphane Adjemian

mode_check.m: account for case where bounds are infinite and not plots are shown

parent 4d89452e
...@@ -111,14 +111,22 @@ for plt = 1:nbplt ...@@ -111,14 +111,22 @@ for plt = 1:nbplt
end end
end end
xx = x; xx = x;
if x(kk)~=0 if x(kk)~=0 || ~isinf(BoundsInfo.lb(kk)) || ~isinf(BoundsInfo.lb(kk))
l1 = max(BoundsInfo.lb(kk),(1-sign(x(kk))*ll)*x(kk)); m1 = 0; %lower bound l1 = max(BoundsInfo.lb(kk),(1-sign(x(kk))*ll)*x(kk)); m1 = 0; %lower bound
l2 = min(BoundsInfo.ub(kk),(1+sign(x(kk))*ll)*x(kk)); %upper bound l2 = min(BoundsInfo.ub(kk),(1+sign(x(kk))*ll)*x(kk)); %upper bound
else else
%size info for 0 parameter is missing, use prior standard %size info for 0 parameter is missing, use prior standard
%deviation %deviation
l1 = max(BoundsInfo.lb(kk),-BayesInfo.p2(kk)); m1 = 0; %lower bound upper_bound=BoundsInfo.lb(kk);
l2 = min(BoundsInfo.ub(kk),BayesInfo.p2(kk)); %upper bound if isinf(upper_bound)
upper_bound=-1e-6*DynareOptions.huge_number;
end
lower_bound=BoundsInfo.ub(kk);
if isinf(lower_bound)
lower_bound=-1e-6*DynareOptions.huge_number;
end
l1 = max(lower_bound,-BayesInfo.p2(kk)); m1 = 0; %lower bound
l2 = min(upper_bound,BayesInfo.p2(kk)); %upper bound
end end
binding_lower_bound=0; binding_lower_bound=0;
binding_upper_bound=0; binding_upper_bound=0;
......
Markdown is supported
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