dynare/matlab/gsa/myprctilecol.m

20 lines
355 B
Matlab
Raw Normal View History

function y = myprctilecol(x,p);
xx = sort(x);
[m,n] = size(x);
if m==1 | n==1
m = max(m,n);
if m == 1,
y = x*ones(length(p),1);
return;
end
n = 1;
q = 100*(0.5:m - 0.5)./m;
xx = [min(x); xx(:); max(x)];
else
q = 100*(0.5:m - 0.5)./m;
xx = [min(x); xx; max(x)];
end
q = [0 q 100];
y = interp1(q,xx,p);