Parameter transformation
Posted: Tue Jan 10, 2017 4:59 pm
Dear all,
I am trying to estimate a dsge model using bayesian techniques, with transformed parameters. Specifically, I need the tranformation only for the minimization process, while in the metropolis hasting the original parameters must be considered.
1) Is it enough to declair the transformation using the pound operator and to define the whole model in terms of the transfomred variables, while in the estimation block to declair the original paramteres ones? As for example
parameters ..alpha.. ;
model(linear);
..
#alpha_bound=log(alpha/(1-alpha));
..
y=alpha_bound*x;
..
end;
estimated_params;
..
alpha, .., ..
...
end;
2) Does the computation of the jacobian and of the hessian automatically account for this tranformation? That is, when computing the jacobian of alpha it should use exp(alpha)/(1+exp(alpha))^2 instead of just alpha.
Thank you very much,
irish.
I am trying to estimate a dsge model using bayesian techniques, with transformed parameters. Specifically, I need the tranformation only for the minimization process, while in the metropolis hasting the original parameters must be considered.
1) Is it enough to declair the transformation using the pound operator and to define the whole model in terms of the transfomred variables, while in the estimation block to declair the original paramteres ones? As for example
parameters ..alpha.. ;
model(linear);
..
#alpha_bound=log(alpha/(1-alpha));
..
y=alpha_bound*x;
..
end;
estimated_params;
..
alpha, .., ..
...
end;
2) Does the computation of the jacobian and of the hessian automatically account for this tranformation? That is, when computing the jacobian of alpha it should use exp(alpha)/(1+exp(alpha))^2 instead of just alpha.
Thank you very much,
irish.