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

Removed a useless global (dsge_prior_weight) 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@2002 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
adjemian 2008-08-25 21:17:25 +00:00
parent 469eb87f40
commit 70acbdb89f
2 changed files with 5 additions and 4 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
matlab/dynare_estimation_init.m
View File

@ -32,7 +32,7 @@ function [data,rawdata]=dynare_estimation_init(var_list_, igsa)
% along with Dynare. If not, see <http://www.gnu.org/licenses/>. % along with Dynare. If not, see <http://www.gnu.org/licenses/>.
global M_ options_ oo_ estim_params_ global M_ options_ oo_ estim_params_
global bayestopt_ dsge_prior_weight global bayestopt_
if nargin<2 | isempty(igsa), if nargin<2 | isempty(igsa),
igsa=0; igsa=0;

7
matlab/gmhmaxlik.m
View File

@ -77,7 +77,6 @@ function [PostMod,PostVar,Scale,PostMean] = ...
% along with Dynare. If not, see <http://www.gnu.org/licenses/>. % along with Dynare. If not, see <http://www.gnu.org/licenses/>.
global bayestopt_ estim_params_ options_ global bayestopt_ estim_params_ options_
global dsge_prior_weight % Needed if a BVAR-DSGE model is estimated.
options_.lik_algo = 1; options_.lik_algo = 1;
npar = length(xparam1); npar = length(xparam1);
@ -139,6 +138,7 @@ while j<=MaxNumberOfTuningSimulations
jj = jj + 1; jj = jj + 1;
end end
close(hh); close(hh);
iScale
%% [2] One block metropolis, I update the covariance matrix of the jumping distribution %% [2] One block metropolis, I update the covariance matrix of the jumping distribution
hh = waitbar(0,'Metropolis-Hastings...'); hh = waitbar(0,'Metropolis-Hastings...');
set(hh,'Name','Looking for the posterior covariance...') set(hh,'Name','Looking for the posterior covariance...')
@ -154,7 +154,7 @@ while j<= NumberOfIterations
end end
% I move if the proposal is enough likely... % I move if the proposal is enough likely...
if logpo2 > -inf & log(rand) < logpo2 - ilogpo2 if logpo2 > -inf & log(rand) < logpo2 - ilogpo2
ix2 = proposal; ix2 = proposal;
if logpo2 > mlogpo2 if logpo2 > mlogpo2
ModePar = proposal; ModePar = proposal;
mlogpo2 = logpo2; mlogpo2 = logpo2;
@ -223,10 +223,11 @@ if strcmpi(info,'LastCall')
end end
close(hh); close(hh);
Scale = iScale; Scale = iScale;
iScale
%% %%
%% Now I climb the hill %% Now I climb the hill
%% %%
climb = 0; climb = 1;
if climb if climb
hh = waitbar(0,' '); hh = waitbar(0,' ');
set(hh,'Name','Now I am climbing the hill...') set(hh,'Name','Now I am climbing the hill...')