Bicluster Sampled Coherence Metric (BSCM) method
Here we change how the co-expression p-value, rik was calculated as follows:
(3)  r  j k   = 1  2 π σ  σ  ¯   j | k |    2       exp - σ  j k  2   - σ  ¯   j | k |  2    σ  σ  ¯   j | k |    2
(4)  r  i k   = ∑  j ∈ J  k      r  j k    | J  k   |
σ¯j|k| is the mean variance for the number of genes in bicluster k as determined bootstrap sampling. σσ¯j|k|2 is the standard deviation of the values used to calculate σ¯j|k|. The background distribution is calculated for each condition j∈J and for each number of genes that occurs in a given bicluster k by sampling |k| genes 200 times from experimental condition j and drawing additional samples in sets of 200 until σ¯j|k| and σσ¯j|k|2 change by less than 1%. To determine which genes should be added or removed from a cluster, we calculate a new rik supposing gene i were added or removed. As a practical matter, background distributions for are pre-calculated for all cluster sizes less than or equal to the maximum size represented in the initial seed clusters, and additional background distributions are calculated as needed during program execution.