Я пытаюсь понять этот код в MATLAB. и переведите его в R (источник: https://github.com/mjlaine/mcmcstat/blob/master/plims.m). Тем не менее, я не уверен в том, что он на самом деле делает:
function y=plims(x,p)
% plims(x,p) calculates p quantiles from columns of x
if nargin<2
% p = [0.025,0.975];
p = [0.25,0.5,0.75];
end
[n,m] = size(x); if n==1; n=m;end
y = interp1(sort(x),(n-1)*p+1);
Вышесказанное просто находит p-квантили для каждого столбца? Не будет ли quantile(sort(x),p) в MATLAB просто сделать то же самое? И если предположить, что это не так, как бы я адаптировал вышеприведенные функции к R? Спасибо вам за помощь.