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

Provisions to to show in plots the parameter value around which the sensitivity analysis with the option neighborhood_width is done.

time-shift
Marco Ratto 2015-10-13 10:44:13 +02:00
parent d3eca9338b
commit bc8f546cbb
3 changed files with 15 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
matlab/gsa/mcf_analysis.m
View File

@ -32,6 +32,10 @@ beha_title = options_mcf.beha_title;
nobeha_title = options_mcf.nobeha_title;
title = options_mcf.title;
fname_ = options_mcf.fname_;
xparam1=[];
if isfield(options_mcf,'xparam1'),
xparam1=options_mcf.xparam1;
end
OutputDirectoryName = options_mcf.OutputDirectoryName;
[proba, dproba] = stab_map_1(lpmat, ibeha, inobeha, [],0);
@ -55,7 +59,9 @@ if length(ibeha)>10 && length(inobeha)>10,
end
if ~isempty(indmcf) && ~DynareOptions.nograph,
skipline()
xx=[];
if ~ isempty(xparam1), xx=xparam1(indmcf); end
scatter_mcf(lpmat(ibeha,indmcf),lpmat(inobeha,indmcf), param_names(indmcf,:), ...
'.', [fname_,'_',amcf_name], OutputDirectoryName, amcf_title,[], DynareOptions, ...
'.', [fname_,'_',amcf_name], OutputDirectoryName, amcf_title,xx, DynareOptions, ...
beha_title, nobeha_title)
end

2
matlab/gsa/scatter_mcf.m
View File

@ -122,7 +122,7 @@ for i = 1:p
plot(X(:,i),X(:,j),[plotsymbol,'b'])
end
if ~isempty(xparam1)
hold on, plot(xparam1(i),xparam1(j),'k*')
hold on, plot(xparam1(i),xparam1(j),'s','MarkerFaceColor',[0 0.75 0],'MarkerEdgeColor',[0 0.75 0])
end
hold off;
% axis([-0.1 1.1 -0.1 1.1])

7
matlab/gsa/stab_map_.m
View File

@ -123,6 +123,7 @@ options_mcf.alpha2 = alpha2;
options_mcf.param_names = char(bayestopt_.name(nshock+1:end));
options_mcf.fname_ = fname_;
options_mcf.OutputDirectoryName = OutputDirectoryName;
options_mcf.xparam1 = [];
opt=options_;
options_.periods=0;
@ -573,6 +574,9 @@ if length(iunstable)>0 || length(iwrong)>0,
if length(iunstable)<Nsam || length(istable)>1
itot = [1:Nsam];
% Blanchard Kahn
if neighborhood_width,
options_mcf.xparam1 = xparam1(nshock+1:end);
end
itmp = itot(find(~ismember(itot,istable)));
options_mcf.amcf_name = asname;
options_mcf.amcf_title = atitle;
@ -614,6 +618,9 @@ if length(iunstable)>0 || length(iwrong)>0,
end
if ~isempty(irestriction),
if neighborhood_width,
options_mcf.xparam1 = xparam1;
end
options_mcf.param_names = char(bayestopt_.name);
options_mcf.amcf_name = acalibname;
options_mcf.amcf_title = acalibtitle;