From: Relating gene expression data on two-component systems to functional annotations in Escherichia coli
Data: genes; | /* expression values */ |
Data: functions | /* for each function */ |
Result: significance, tailGenes; | /* vector of zeros */ |
1 normGenes = normalize(genes); | |
2 hist = zeros(1, nPts); | |
3 foreach f ∈ function do | |
4 subset = findPoints(normGenes, f); | |
5 foreach x ∈ subset do | |
6 dens = NumberOfNeighbors(x); | |
7 hist(dens)++; | |
8 randHist = findRandomHistogram(1, nPts, normGenes); | |
9 significance(f) = chiSquaredGoodynessOfFit(hist, randHist); | |
10 tailGenes(f) = findTailGenes(hist, randHist); | |
11 return significance, tailGenes |