Copyright (C) 2012 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 y = myprctilecol(x,p); 0002 0003 % Copyright (C) 2012 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 xx = sort(x); 0021 [m,n] = size(x); 0022 0023 if m==1 | n==1 0024 m = max(m,n); 0025 if m == 1, 0026 y = x*ones(length(p),1); 0027 return; 0028 end 0029 n = 1; 0030 q = 100*(0.5:m - 0.5)./m; 0031 xx = [min(x); xx(:); max(x)]; 0032 else 0033 q = 100*(0.5:m - 0.5)./m; 0034 xx = [min(x); xx; max(x)]; 0035 end 0036 0037 q = [0 q 100]; 0038 y = interp1(q,xx,p);