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

print pdf and eps; ishock debugged 

git-svn-id: https://www.dynare.org/svn/dynare/dynare_v4@714 ac1d8469-bf42-47a9-8791-bf33cf982152
time-shift
ratto 2006-04-13 13:01:52 +00:00
parent 76fd18af8b
commit 6db1db53c5
1 changed files with 52 additions and 40 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

92
matlab/stab_map_2.m
View File

@ -1,6 +1,6 @@
%function stab_map_2(x,alpha2,istab,fnam)
function stab_map_2(x,alpha2,fnam, ishock)
% function stab_map_2(x,alpha2,fnam, ishock)
function stab_map_2(x,alpha2,fnam, dirname)
% function stab_map_2(x,alpha2,fnam)
%
% Copyright (C) 2005 Marco Ratto
% THIS PROGRAM WAS WRITTEN FOR MATLAB BY
@ -16,9 +16,15 @@ function stab_map_2(x,alpha2,fnam, ishock)
%global bayestopt_ estim_params_ dr_ options_ ys_ fname_
global bayestopt_ estim_params_ options_ oo_ M_
if nargin<4,
ishock=0;
npar=size(x,2);
ishock= npar>estim_params_.np;
if nargin<3,
fnam='';
end
if nargin<4,
dirname='';
end
ys_ = oo_.dr.ys;
dr_ = oo_.dr;
fname_ = M_.fname;
@ -34,47 +40,53 @@ fig_nam_=[fname_,'_',fnam,'_corr_'];
ifig=0;
j2=0;
if ishock==0
npar=estim_params_.np;
npar=estim_params_.np;
else
npar=estim_params_.np+nshock;
npar=estim_params_.np+nshock;
end
for j=1:npar,
i2=find(abs(c00(:,j))>alpha2);
if length(i2)>0,
for jx=1:length(i2),
j2=j2+1;
if mod(j2,12)==1,
ifig=ifig+1;
figure('name',['Correlations in the ',fnam,' sample ', num2str(ifig)]),
end
subplot(3,4,j2-(ifig-1)*12)
% bar(c0(i2,j)),
% set(gca,'xticklabel',bayestopt_.name(i2)),
% set(gca,'xtick',[1:length(i2)])
%plot(stock_par(ixx(nfilt+1:end,i),j),stock_par(ixx(nfilt+1:end,i),i2(jx)),'.k')
%hold on,
plot(x(:,j),x(:,i2(jx)),'.')
% xlabel(deblank(estim_params_.param_names(j,:)),'interpreter','none'),
% ylabel(deblank(estim_params_.param_names(i2(jx),:)),'interpreter','none'),
if ishock,
xlabel(bayestopt_.name{j},'interpreter','none'),
ylabel(bayestopt_.name{i2(jx)},'interpreter','none'),
else
xlabel(bayestopt_.name{j+nshock},'interpreter','none'),
ylabel(bayestopt_.name{i2(jx)+nshock},'interpreter','none'),
end
title(['cc = ',num2str(c0(i2(jx),j))])
if (mod(j2,12)==0) & j2>0,
saveas(gcf,[fig_nam_,int2str(ifig)])
end
end
i2=find(abs(c00(:,j))>alpha2);
if length(i2)>0,
for jx=1:length(i2),
j2=j2+1;
if mod(j2,12)==1,
ifig=ifig+1;
figure('name',['Correlations in the ',fnam,' sample ', num2str(ifig)]),
end
subplot(3,4,j2-(ifig-1)*12)
% bar(c0(i2,j)),
% set(gca,'xticklabel',bayestopt_.name(i2)),
% set(gca,'xtick',[1:length(i2)])
%plot(stock_par(ixx(nfilt+1:end,i),j),stock_par(ixx(nfilt+1:end,i),i2(jx)),'.k')
%hold on,
plot(x(:,j),x(:,i2(jx)),'.')
% xlabel(deblank(estim_params_.param_names(j,:)),'interpreter','none'),
% ylabel(deblank(estim_params_.param_names(i2(jx),:)),'interpreter','none'),
if ishock,
xlabel(bayestopt_.name{j},'interpreter','none'),
ylabel(bayestopt_.name{i2(jx)},'interpreter','none'),
else
xlabel(bayestopt_.name{j+nshock},'interpreter','none'),
ylabel(bayestopt_.name{i2(jx)+nshock},'interpreter','none'),
end
title(['cc = ',num2str(c0(i2(jx),j))])
if (mod(j2,12)==0) & j2>0,
saveas(gcf,[dirname,'\',fig_nam_,int2str(ifig)])
eval(['print -depsc2 ' dirname '\' fig_nam_ int2str(ifig)]);
eval(['print -dpdf ' dirname '\' fig_nam_ int2str(ifig)]);
if options_.nograph, close(gcf), end
end
end
if (j==(npar)) & j2>0,
saveas(gcf,[fig_nam_,int2str(ifig)])
end
end
if (j==(npar)) & j2>0,
saveas(gcf,[dirname,'\',fig_nam_,int2str(ifig)])
eval(['print -depsc2 ' dirname '\' fig_nam_ int2str(ifig)]);
eval(['print -dpdf ' dirname '\' fig_nam_ int2str(ifig)]);
if options_.nograph, close(gcf), end
end
end
if ifig==0,
disp(['No correlation term >', num2str(alpha2),' found for ',fnam])
disp(['No correlation term >', num2str(alpha2),' found for ',fnam])
end
%close all