From b3470f646ff3dd7f915bd8341d1de9edf2ac1927 Mon Sep 17 00:00:00 2001 From: sebastien Date: Tue, 5 Aug 2008 14:12:31 +0000 Subject: [PATCH] v4 prior_bounds.m: small simplification git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@1985 ac1d8469-bf42-47a9-8791-bf33cf982152 --- matlab/prior_bounds.m | 10 +++------- 1 file changed, 3 insertions(+), 7 deletions(-) diff --git a/matlab/prior_bounds.m b/matlab/prior_bounds.m index 57144c523..732bd71a9 100644 --- a/matlab/prior_bounds.m +++ b/matlab/prior_bounds.m @@ -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 \ No newline at end of file +end