ghlp_5gp.m

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

function COUNT=ghlp_5gp(k,X,flag)
% flag = 0, descending order
% flag = 1, ascending order
% X --- Sources
COUNT=1:1:k;
for m=1:k
    COUNT(m)=length(find(X==COUNT(m)));
end
ID=1:1:5;
IA=1:1:5;
for m=1:5
    if(m==5)
        TEMP=COUNT(1:1,29:33);
    else
        TEMP=COUNT(1:1,((m-1)*7+1):m*7);
    end
    [C,ID(m)]=max(TEMP);
    ID(m)=ID(m)+(m-1)*7;
    [C,IA(m)]=min(TEMP);
    IA(m)=IA(m)+(m-1)*7;
end
if(flag==1)
    COUNT=IA;
else
    COUNT=ID;
end