dcompare

 Copyright (C) 2001-2011 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 <http://www.gnu.org/licenses/>.

0001 function dcompare(s1)
0002 
0003 % Copyright (C) 2001-2011 Dynare Team
0004 %
0005 % This file is part of Dynare.
0006 %
0007 % Dynare is free software: you can redistribute it and/or modify
0008 % it under the terms of the GNU General Public License as published by
0009 % the Free Software Foundation, either version 3 of the License, or
0010 % (at your option) any later version.
0011 %
0012 % Dynare is distributed in the hope that it will be useful,
0013 % but WITHOUT ANY WARRANTY; without even the implied warranty of
0014 % MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
0015 % GNU General Public License for more details.
0016 %
0017 % You should have received a copy of the GNU General Public License
0018 % along with Dynare.  If not, see <http://www.gnu.org/licenses/>.
0019 
0020 global options_ nvx nvy x y lag1
0021 
0022 ftest(s1,0) ;
0023 
0024 i = [lag1(1):size(x,2)-lag1(2)+1]' ;
0025 
0026 if size(options_.smpl,1) == 1
0027     error(['DSAMPLE not specified.']) ;
0028 end
0029 
0030 if options_.smpl(3) > 0
0031     if options_.smpl(3) == 2
0032         if options_.smpl(1)<0 || options_.smpl(2)>size(x,2)-lag1(2)
0033             error ('Wrong sample.') ;
0034         end
0035         i = [options_.smpl(1)+lag1(1):options_.smpl(2)+lag1(1)]' ;
0036     elseif options_.smpl(3) == 1
0037         if options_.smpl(1)>size(x,2)-lag1(2)
0038             error ('Wrong sample.') ;
0039         end
0040         i = [lag1(1):options_.smpl(1)+lag1(1)]' ;
0041     end
0042 end
0043 
0044 j = bseastr(nvx,nvy) ;
0045 
0046 if stop
0047     return ;
0048 end
0049 
0050 z = mean(mean(abs(x(j,i)-y(j,i)))) ;
0051 
0052 disp (['The mean absolute difference between set ' s1(1,:) 'and set ' s1(2,:)]) ;
0053 disp (['is : ' num2str(z)]) ;
0054 return ;
0055 
0056