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

v4 prior_bounds.m: small simplification 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@1985 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
sebastien 2008-08-05 14:12:31 +00:00
parent 19e9bb3bd1
commit b3470f646f
1 changed files with 3 additions and 7 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
matlab/prior_bounds.m
View File

@ -58,12 +58,8 @@ for i=1:n
bounds(i,1) = norminv(options_.prior_trunc,pmean(i),p2(i));
bounds(i,2) = norminv(1-options_.prior_trunc,pmean(i),p2(i));
case 4
nu = p2(i);
mu = pmean(i);
beta = ( gamma( (nu-1)/2 ) / mu / gamma( nu/2 ) )^2;
a=2/beta;
bounds(i,1) = 1/sqrt(gaminv(1-options_.prior_trunc,p2(i)/2,beta));
bounds(i,2) = 1/sqrt(gaminv(options_.prior_trunc,p2(i)/2,beta));
bounds(i,1) = 1/sqrt(gaminv(1-options_.prior_trunc, p2(i)/2, 2/p1(i)));
bounds(i,2) = 1/sqrt(gaminv(options_.prior_trunc, p2(i)/2, 2/p1(i)));
case 5
bounds(i,1) = p1(i);
bounds(i,2) = p2(i);
@ -73,4 +69,4 @@ for i=1:n
otherwise
error(sprintf('prior_bounds: unknown distribution shape (index %d, type %d)', i, pshape(i)));
end
end
end