Input: |
all pairwise gene candidates GPs; initial MI threshold MIth = MI's centroid; initial CC threshold CCth = CC's centroid; increments δmi, δcc for MI and CC. |
Output: |
   classified APGs, UPGs and QPGs. |
while count(GPs)>0 do |
1. construct APGs, QPGs using initial MIth, CCth and P-value; 2. group the others into UPGs; if (APGs' undersized) && count(QPGs)>0 then do |
MIth=MIth-δmi & CCth=CCth-δcc; continue Step 1 for QPGs & obtain ΔAPGs and ΔUPGs; APG=APGs+ΔAPGs & UPGs=UPGs-ΔUPGs. |
   elseif(APGs' oversized) then do |
MIth=MIth-δmi & CCth=CCth+δcc; continue Step 1 for APGs & obtain ΔAPGs and ΔUPGs; APG=APGs-ΔAPGs & UPGs=UPGs+ΔAPGs. |
   endif |
end |