getp.m

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

function P=getp(window,selectnumber,RD)
% get the probability
[rc,cc]=size(RD);
i=rc;
j=1;
P=0;
while(i-window>0)
    RDT=RD(i-window+1:i,1:cc);
    CT=count(33,RDT,0);
    M=RD(i-window:i-window,1:cc);
    N=CT(1:1,1:selectnumber);
    P(j)=ratio(M,N);
    i=i-1;
    j=j+1;
end