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

Remove further unused function from matlab folder

new-samplers
Johannes Pfeifer 2023-12-08 10:03:35 +01:00
parent 435b103cf5
commit 62e8b275a0
4 changed files with 0 additions and 199 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

84
matlab/ftest.m
View File

@ -1,84 +0,0 @@
function ftest (s1,s2)
% Copyright © 2001-2017 Dynare Team
%
% This file is part of Dynare.
%
% Dynare is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% Dynare is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with Dynare. If not, see <https://www.gnu.org/licenses/>.
global nvx nvy x y lag1
if size(s1,1) ~= 2
error ('Spécifiez deux fichiers pour la comparaison.') ;
end
for i = 1:2
if ~ isempty(find(abs(s1(i,:)) == 46))
error ('Entrez les noms de fichiers sans extensions.') ;
end
end
s1 = [s1 [' ';' ']] ;
file1 = [s1(1,1:min(find(abs(s1(1,:)) == 32))-1) '.BIN'] ;
file2 = [s1(2,1:min(find(abs(s1(2,:)) == 32))-1) '.BIN'] ;
fid=fopen(file1,'r') ;
n1 = fread(fid,1,'int') ;
n2 = fread(fid,1,'int') ;
n3 = fread(fid,1,'int') ;
lag1 = fread(fid,4,'int') ;
nvx = fread(fid,[n1,n3],'int') ;
x = fread(fid,[n1,n2],'float64') ;
fclose(fid) ;
nvx = char(nvx) ;
fid=fopen(file2,'r') ;
n1 = fread(fid,1,'int') ;
n2 = fread(fid,1,'int') ;
n3 = fread(fid,1,'int') ;
lag2 = fread(fid,4,'int') ;
nvy = fread(fid,[n1,n3],'int') ;
y = fread(fid,[n1,n2],'float64') ;
fclose(fid) ;
nvy = char(nvy) ;
if size(x,1) ~= size(y,1)
error ('FTEST: The two files don''t have the same number of variables.');
end
for i = 1:size(x,1)
if ~ strcmp(nvx(i,:),nvy(i,:))
error ('FTEST: The two files don''t have the same variables.') ;
end
end
if nnz(lag1 - lag2) > 0
error ('FTEST: Leads and lags aren''t the same in both files.') ;
end
j = zeros(size(s2,1),1);
for i=1:size(s2,1)
k = strmatch(s2(i,:),nvx,'exact') ;
if isempty(k)
t = ['FTEST: Variable ' s2(i) 'doesn''t exist'] ;
error (t) ;
else
j(i) =k;
end
end
y = y(j,:) ;
x = x(j,:) ;
%06/18/01 MJ replaced beastr by strmatch

52
matlab/gcompare.m
View File

@ -1,52 +0,0 @@
function gcompare(s1,s2)
% GCOMPARE : GCOMPARE ( [ 'file1' ; 'file2' ] , [ 'var1' ; 'var2' ...] )
% This optional command plots the trajectories of a list of
% variables in two different simulations. One plot is drawn
% for each variable. The trajectories must have been previously
% saved by the instruction DYNASAVE. The simulation in file1
% is refered to as the base simulation.
% Copyright © 2001-2017 Dynare Team
%
% This file is part of Dynare.
%
% Dynare is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% Dynare is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with Dynare. If not, see <https://www.gnu.org/licenses/>.
global options_ M_
global nvx nvy x y lag1
ftest(s1,s2) ;
ix = [1-lag1(1):size(x,2)-lag1(1)]' ;
i = [lag1(1):size(ix,1)-lag1(2)+1]' ;
if options_.smpl == 0
i = [M_.maximum_lag:size(y,2)]' ;
else
i = [options_.smpl(1):options_.smpl(2)]' ;
end
for k = 1:size(x,1)
figure ;
plot (ix(i),x(k,i),ix(i),y(k,i)) ;
xlabel (['Periods']) ;
title (['Variable ' s2(k,:)]) ;
l = min(i) + 1;
ll = max(i) - 1 ;
text (l,x(k,l),s1(1,:)) ;
text (ll,y(k,ll),s1(2,:)) ;
end
% 06/18/01 MJ corrected treatment of options_.smpl
% 06/24/01 MJ removed color specification

8
matlab/mcpath_function.m
View File

@ -1,8 +0,0 @@
function [res,jac,domerr] = mcpath_function(func,z,jacflag,varargin)
domerr = 0;
if jacflag
[res,jac] = func(z,varargin{:});
else
res = func(z,varargin{:});
end

55
matlab/metropolis_run_analysis.m
View File

@ -1,55 +0,0 @@
function metropolis_run_analysis(M_,basetopt,j)
%function metropolis_run_analysis(M_,basetopt,j
% analizes Metropolis runs
%
% INPUTS
% M_: (struct) Model structure
% basetopt: (struct) Estimated parameters structure
% j: (int) Index of estimated paramter
%
% OUTPUTS
% none
%
% SPECIAL REQUIREMENTS
% none
% Copyright © 2003-2017 Dynare Team
%
% This file is part of Dynare.
%
% Dynare is free software: you can redistribute it and/or modify
% it under the terms of the GNU General Public License as published by
% the Free Software Foundation, either version 3 of the License, or
% (at your option) any later version.
%
% Dynare is distributed in the hope that it will be useful,
% but WITHOUT ANY WARRANTY; without even the implied warranty of
% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
% GNU General Public License for more details.
%
% You should have received a copy of the GNU General Public License
% along with Dynare. If not, see <https://www.gnu.org/licenses/>.
load([M_.fname '/metropolis/' M_.fname '_mh_history'])
nblck = record.Nblck;
ndraws = sum(record.MhDraws(:,1));
logPost = [];
params = [];
blck = 1;
for i=1:record.LastFileNumber
fname = [M_.fname '/metropolis/' M_.fname '_mh' int2str(i) '_blck' ...
int2str(blck) '.mat'];
if exist(fname,'file')
o=load(fname);
logPost = [logPost; o.logpo2];
params = [params; o.x2];
end
end
figure;
subplot(2,1,1)
plot(logPost)
subplot(2,1,2)
plot(params(:,j))
title(['parameter ' int2str(j)])