stepan
/
dynare
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Fixed bug. 

Forgot to change the name of the routine for the inverse CDF of the
Weibull distribution.
time-shift
Stéphane Adjemian (Charybdis) 2015-12-14 15:32:03 +01:00
parent e53f3a67fa
commit 944ed4b42b
2 changed files with 4 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
matlab/draw_prior_density.m
View File

@ -102,8 +102,8 @@ switch pshape(indx)
dens = exp(lpdfig2(abscissa-p3(indx),p6(indx),p7(indx)));
case 8
density = @(x,a,b,c) exp(lpdfgweibull(x, a, b, c));
infbound = p3(indx)+icdfweibull(truncprior,p6(indx),p7(indx));
supbound = p3(indx)+icdfweibull(1-truncprior,p6(indx),p7(indx));
infbound = p3(indx)+wblinv(truncprior,p6(indx),p7(indx));
supbound = p3(indx)+wblinv(1-truncprior,p6(indx),p7(indx));
abscissa = linspace(infbound,supbound,steps);
dens = density(abscissa,p6(indx),p7(indx),p3(indx));
otherwise

4
matlab/prior_bounds.m
View File

@ -155,8 +155,8 @@ for i=1:length(p6)
bounds.lb(i) = p3(i);
bounds.ub(i) = Inf;
else
bounds.lb(i) = p3(i)+icdfweibull(prior_trunc,p6(i),p7(i));
bounds.ub(i) = p3(i)+icdfweibull(1-prior_trunc,p6(i),p7(i));
bounds.lb(i) = p3(i)+wblinv(prior_trunc,p6(i),p7(i));
bounds.ub(i) = p3(i)+wblinv(1-prior_trunc,p6(i),p7(i));
end
otherwise
error(sprintf('prior_bounds: unknown distribution shape (index %d, type %d)', i, pshape(i)));