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

v4 prior_posterior_statistics.m: Octave compatibility issues: 

* added missing .mat extensions
* use text waiting bar instead of waitbar()


git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@1930 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
sebastien 2008-07-04 13:20:49 +00:00
parent 3ac7bab6f0
commit 6fa483bd3f
1 changed files with 26 additions and 12 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

38
matlab/prior_posterior_statistics.m
View File

@ -97,7 +97,11 @@ end
irun = ones(8,1);
ifil = zeros(8,1);
h = waitbar(0,'Taking subdraws...');
if exist('OCTAVE_VERSION')
diary off;
else
h = waitbar(0,'Taking subdraws...');
end
stock_param = zeros(MAX_nruns, npar);
stock_logpo = zeros(MAX_nruns,1);
@ -189,66 +193,76 @@ for b=1:B
if irun(1) > MAX_nsmoo || b == B
stock = stock_smooth(:,:,1:irun(1)-1);
ifil(1) = ifil(1) + 1;
save([DirectoryName '/' M_.fname '_smooth' int2str(ifil(1))],'stock');
save([DirectoryName '/' M_.fname '_smooth' int2str(ifil(1)) '.mat'],'stock');
irun(1) = 1;
end
if nvx && (irun(2) > MAX_ninno || b == B)
stock = stock_innov(:,:,1:irun(2)-1);
ifil(2) = ifil(2) + 1;
save([DirectoryName '/' M_.fname '_inno' int2str(ifil(2))],'stock');
save([DirectoryName '/' M_.fname '_inno' int2str(ifil(2)) '.mat'],'stock');
irun(2) = 1;
end
if nvn && (irun(3) > MAX_nerro || b == B)
stock = stock_error(:,:,1:irun(3)-1);
ifil(3) = ifil(3) + 1;
save([DirectoryName '/' M_.fname '_error' int2str(ifil(3))],'stock');
save([DirectoryName '/' M_.fname '_error' int2str(ifil(3)) '.mat'],'stock');
irun(3) = 1;
end
if naK && (irun(4) > MAX_naK || b == B)
stock = stock_filter(:,:,:,1:irun(4)-1);
ifil(4) = ifil(4) + 1;
save([DirectoryName '/' M_.fname '_filter' int2str(ifil(4))],'stock');
save([DirectoryName '/' M_.fname '_filter' int2str(ifil(4)) '.mat'],'stock');
irun(4) = 1;
end
if irun(5) > MAX_nruns || b == B
stock = stock_param(1:irun(5)-1,:);
ifil(5) = ifil(5) + 1;
save([DirectoryName '/' M_.fname '_param' int2str(ifil(5))],'stock','stock_logpo','stock_ys');
save([DirectoryName '/' M_.fname '_param' int2str(ifil(5)) '.mat'],'stock','stock_logpo','stock_ys');
irun(5) = 1;
end
if horizon && (irun(6) > MAX_nforc1 || b == B)
stock = stock_forcst_mean(:,:,1:irun(6)-1);
ifil(6) = ifil(6) + 1;
save([DirectoryName '/' M_.fname '_forc_mean' int2str(ifil(6))],'stock');
save([DirectoryName '/' M_.fname '_forc_mean' int2str(ifil(6)) '.mat'],'stock');
irun(6) = 1;
end
if horizon && (irun(7) > MAX_nforc2 || b == B)
stock = stock_forcst_total(:,:,1:irun(7)-1);
ifil(7) = ifil(7) + 1;
save([DirectoryName '/' M_.fname '_forc_total' int2str(ifil(7))],'stock');
save([DirectoryName '/' M_.fname '_forc_total' int2str(ifil(7)) '.mat'],'stock');
irun(7) = 1;
end
if moments_varendo && (irun(8) > MAX_momentsno || b == B)
stock = stock_moments(1:irun(8)-1);
ifil(8) = ifil(8) + 1;
save([DirectoryName '/' M_.fname '_moments' int2str(ifil(8))],'stock');
save([DirectoryName '/' M_.fname '_moments' int2str(ifil(8)) '.mat'],'stock');
irun(8) = 1;
end
waitbar(b/B,h);
if exist('OCTAVE_VERSION')
printf('Taking subdraws: %3.f%% done\r', b/B);
else
waitbar(b/B,h);
end
end
if exist('OCTAVE_VERSION')
printf('\n');
diary on;
else
close(h)
end
close(h)
stock_gend=gend;
stock_data=Y;
save([DirectoryName '/' M_.fname '_data'],'stock_gend','stock_data');
save([DirectoryName '/' M_.fname '_data.mat'],'stock_gend','stock_data');
if options_.smoother
pm3(endo_nbr,gend,ifil(1),B,'Smoothed variables',...