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

Fixed bug of ~ in output list for non recent versions of matlab.

time-shift
Marco Ratto 2011-06-21 09:43:24 +02:00
parent 86befc3e4e
commit a96bad899d
2 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

6
matlab/dynare_identification.m
View File

@ -271,7 +271,7 @@ if iload <=0,
else
params = prior_draw();
end
[~, ideJ, ideH, ideGP, ~ , info] = ...
[dum1, ideJ, ideH, ideGP, dum2 , info] = ...
identification_analysis(params,indx,indexo,options_MC,data_info, prior_exist, name_tex,0);
if iteration==0,
MAX_tau = min(SampleSize,ceil(MaxNumberOfBytes/(size(ideH.siH,1)*nparam)/8));
@ -423,7 +423,7 @@ if SampleSize > 1,
jcrit=find(idemoments.ino);
if length(jcrit)<SampleSize,
if isempty(jcrit),
[~,jmax]=max(idemoments.cond);
[dum,jmax]=max(idemoments.cond);
fprintf('\n')
tittxt = 'Draw with HIGHEST condition number';
fprintf('\n')
@ -436,7 +436,7 @@ if SampleSize > 1,
disp_identification(pdraws(jmax,:), idemodel_max, idemoments_max, name);
close all,
plot_identification(pdraws(jmax,:),idemoments_max,idehess_max,idemodel_max,idelre_max,1,tittxt,name,IdentifDirectoryName,1);
[~,jmin]=min(idemoments.cond);
[dum,jmin]=min(idemoments.cond);
fprintf('\n')
tittxt = 'Draw with SMALLEST condition number';
fprintf('\n')

6
matlab/plot_identification.m
View File

@ -273,11 +273,11 @@ else
eval(['print -dpdf ' IdentifDirectoryName '/' M_.fname '_ident_COND']);
if options_.nograph, close(gcf); end
ncut=floor(SampleSize/10*9);
[~,is]=sort(idelre.cond);
[dum,is]=sort(idelre.cond);
[proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberLRE', 1, [], IdentifDirectoryName, 0.1);
[~,is]=sort(idemodel.cond);
[dum,is]=sort(idemodel.cond);
[proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberModel', 1, [], IdentifDirectoryName, 0.1);
[~,is]=sort(idemoments.cond);
[dum,is]=sort(idemoments.cond);
[proba, dproba] = stab_map_1(params, is(1:ncut), is(ncut+1:end), 'MC_HighestCondNumberMoments', 1, [], IdentifDirectoryName, 0.1);
% [proba, dproba] = stab_map_1(idemoments.Mco', is(1:ncut), is(ncut+1:end), 'HighestCondNumberMoments_vs_Mco', 1, [], IdentifDirectoryName);
% for j=1:nparam,