count_ol.m

for double color balls, it can help you make the list with high p

function COUNT=count_ol(k,X)
% flag = 0, count and sort in descending order
% flag = 1, count and sort in ascending order
% k --- the count of items
% X --- Sources
COUNT=1:1:k;
for m=1:k
    COUNT(m)=length(find(X==COUNT(m)));
end
IX=1:1:k;
COUNT=[IX;COUNT];