Results -- 3600 probe tutorial data set

Owned by Bruce Hoff
Last updated: Aug 22, 2011Version comment
9 min read

Results for WGCNA 'female liver' tutorial data set, having 3600 probes and 135 samples

Sage Coexpression, organic correlation algorithm

Sage Coexpression, WGCNA correlation algorithm

WGCNA Coexpression

WGCNA Coexpression, blocksize<=600


Unix space ('M' field from /usr/bin/time): 6471520KB (or 6.47GB)

Unix space ('M' field from /usr/bin/time):: 6474832K (6.47GB)

Unix space ('M' field from /usr/bin/time):: 1958480KB (1.96GB)

Unix space: 0.4 GB

 

Unix time('E' field from /usr/bin/time): 14:04.38elapsed

Unix time('E' field from /usr/bin/time):13:33.74elapsed

Unix time('E' field from /usr/bin/time):2:07.23elapsed

Unix time: 1m:01s

 

Unix space #2: 6472128 (6.47GB)

Unix space #2: 6090192 (6.1GB)

Unix space #2: 1922800 (1.9GB)

 

 

Unix time #2: 11:36.09elapsed

Unix time #2: 10:58.03elapsed

Unix time #2: 1:28.53elapsed

 

 










Module difference vs. this baseline (% of pairs for which the algorithm disagrees with this baseline on module co-membership):

0

4.3%

9.1%

 

Module co-membership comparison (vs. this column as a baseline):

turquoise blue brown yellow green red grey black pink magenta purple greenyellow tan salmon cyan midnightblue lightcyan
  turquoise          581    0     0      0     0   0    0     0    0       0      0           0   0      0    0            0         0
  blue                 0  488     0      0     0   0    0     0    0       0      0           0   0      0    0            0         0
  brown                0    0   459      0     0   0    0     0    0       0      0           0   0      0    0            0         0
  yellow               0    0     0    425     0   0    0     0    0       0      0           0   0      0    0            0         0
  green                0    0     0      0   389   0    0     0    0       0      0           0   0      0    0            0         0
  red                  0    0     0      0     0 363    0     0    0       0      0           0   0      0    0            0         0
  grey                 0    0     0      0     0   0  200     0    0       0      0           0   0      0    0            0         0
  black                0    0     0      0     0   0    0   122    0       0      0           0   0      0    0            0         0
  pink                 0    0     0      0     0   0    0     0  120       0      0           0   0      0    0            0         0
  magenta              0    0     0      0     0   0    0     0    0      99      0           0   0      0    0            0         0
  purple               0    0     0      0     0   0    0     0    0       0     87           0   0      0    0            0         0
  greenyellow          0    0     0      0     0   0    0     0    0       0      0          78   0      0    0            0         0
  tan                  0    0     0      0     0   0    0     0    0       0      0           0  43      0    0            0         0
  salmon               0    0     0      0     0   0    0     0    0       0      0           0   0     40    0            0         0
  cyan                 0    0     0      0     0   0    0     0    0       0      0           0   0      0   37            0         0
  midnightblue         0    0     0      0     0   0    0     0    0       0      0           0   0      0    0           35         0
  lightcyan            0    0     0      0     0   0    0     0    0       0      0           0   0      0    0            0        34

brown yellow green red turquoise blue grey black pink magenta purple greenyellow tan salmon cyan midnightblue lightcyan grey60
  3    427      0     0   0         0    0    0     0    0       0      0           0   0      0    0            0         0      0
  1      0    405     0   0         0  150    0     0    0       0      0           0   0      0    0            0         0      0
  4      0      0   347   0         0    0    0     0    0       0      0           0   0      0    0            0         0      0
  2      0      0     0 331         0    0    0   134    0       0      0           0   0      0    0            0         0      0
  5      0      0     0   0       319    0    0     0    0       0      0           0   0      0    0            0         0      0
  6      0      0     0   0         0  291    0     0    0       0      0           0   0      0    0            0         0      0
  7      0      0     0   0       203    0    0     0    0       0      0           0   0      0    0            0         0      0
  0      3     22     1   0         4    0  168     0    0       0      0           0   0      1    0            0         0      0
  8      0      0     0   0         0    0    0     0  122       0      0           0   0      0    0            0         0      0
  9      0      0     0   0         0    0    0     0    0     120      0           0   0      0    0            0         0      0
  10     0      0     0   0         0    0    0     0    0       0     99           0   0      0    0            0         0      0
  11     0      0     0   0         0    0    0     0    0       0      0          86   0      0    0            0         0      0
  12     0      0     0   0         0   81    0     0    0       0      0           0   0      0    0            0         0      0
  13     0      0     0   0         0    0    0     0    0       0      0           0  76      0    0            0         0      0
  15     0      0     0   0         0    0    0     0    0       0      0           0   0     42    0            0         0      0
  14     0      0    25   0         0    0    0     0    0       0      0           0   0      0   39            0         0      0
  16     0      0     0   0         0    0    0     0    0       0      0           0   0      0    0           37         0      0
  17     0      0     0   0         0    0    0     0    0       0      0           0   0      0    0            0        34      0
  18     0      0     0   0         0    0    0     0    0       0      0           0   0      0    0            0         0     33

turquoise blue red green brown yellow black magenta grey pink purple greenyellow tan cyan lightcyan salmon midnightblue
  3        280   13   3     0     1      3     0       0    9    0      0           0   0    0         0      0            0
  1          5  239 225     5     0      1     0       1    8   15      0           1   0    0         0      5            8
  5          0    3   0   223     1      0     0       0    5    0      0           0   0    0         0      0            0
  2          8    0   4     9   218     61     0       0   17    2      0           0   0    0         0      0            0
  4          0    2  22     6    10    185     1       0   18    0      0           0   0    0         0      0            0
  7        173    0   1    18     0      0     0       0    3    0      0           0   0    0         0      0            0
  6          0    3  19     1     4    170     0       0    2    0      0           0   0    0         0      0            0
  13         0    0   0     0     0      0   120       0    0    0      0           0   0    0         0      0            0
  10         0  111   6    11     5      0     0       0    0    0      0           3   0    0         0      0            1
  8          0    0   3    11   108      0     0       0   11    7      0           0   0    0         0      3           13
  9          6    0   6     4   106      3     1       0   11    5      0           0   0    0         0      0            9
  11         8  100   0    19     2      2     0       0    0    0      0           0   0    0         0      0            0
  16         0    0   0     0     0      0     0      98    1    0      0           0   0    0         0      0            0
  0          0    0   1     1     2      0     0       0   93    0      0           0  43    0         0      0            0
  17         0    0   0     0     0      0     0       0    0   91      0           0   0    0         0      0            0
  14         4    5   7     0     0      0     0       0    3    0     87           0   0    0         0      0            0
  15        14   10   1     0     0      0     0       0    0    0      0          74   0    0         0      1            0
  12         1    1  65    37     1      0     0       0   15    0      0           0   0    0         0      1            4
  19        50    0   0     0     0      0     0       0    0    0      0           0   0    0         0      0            0
  18         0    0   0    40     1      0     0       0    1    0      0           0   0    0         0     30            0
  20         0    0   0     0     0      0     0       0    3    0      0           0   0   37         0      0            0
  22         0    0   0     1     0      0     0       0    0    0      0           0   0    0        34      0            0
  21        32    1   0     3     0      0     0       0    0    0      0           0   0    0         0      0            0

 


<<< Sage dendrogram and modules, with WGCNA modules underneath, for comparison.


WGCNA dendrogram and modules, with Sage modules underneath, for comparison. >>>>


 


TOM Comparison

We captured the TOM matrices from the Sage and WGCNA algorithms.  The values are the same (down to the machine's floating point resolution) as illustrated in this scatter plot: 

 

Note:  To make the TOM values the same for both algorithms, the default value for the WGCNA parameter "tomType" is overridden, changing it from "signed" to "unsigned".  The definition of the parameter is given here:http://www.genetics.ucla.edu/labs/horvath/CoexpressionNetwork/Rpackages/WGCNA/OldPackages/WGCNA_0.67-2.pdf

  • TOMType: a character string specifying TOM type to be calculated. One of "unsigned", "signed". If "unsigned", the standard TOM will be used (more generally, TOM function will receive the adjacency as input). If "signed", TOM will keep track of the sign of the adjacency between neighbors. ... If TOMType = "unsigned", entries of the adjacency matrix are required to lie between 0 and 1; for TOMType = "signed" they can be between -1 and 1.

The TOM values are the inputs to the clustering algorithms.  The dendrograms appear nearly identical above, therefore it seems the differences between the two algorithms are in the module determination from the dendrograms. 

To compare the dendrogram similarity, we permute the rows of a 3600x3600 diagonal matrix according to the leaf permutations of one of the two dendrograms and the columns according to the other.  Identical dendrograms cause the diagonal to be retained, as seen in this result:

Analysis of the effect of dividing genes into 'blocks'

WGCNA, one block

WGCNA, 7 blocks (block size < 600)

Module difference vs. this baseline (% of pairs for which the algorithm disagrees with this baseline on module co-membership):

7.5%

Module co-membership comparison (vs. this column as a baseline):

1   3   2   5   4   7   6  13  10 16 11  8  0 17 14 15 12 18  9 19 20 22 21
  1  322   5   3   0  16   0  10   0 120  0  3  8  6  0 85 11  9  0 10  0  0  0  1
  5    7 262   3   0   0   3   0   0   0  0 26  0  2  0  8  0  0  0  0  1  0  0  0
  3    1   3 233   1  12   0   3   0   3  0  1 94  8  0  0  0  1  2 47  0  0  0  0
  4    3   1  10 210  10   7   1   0   9  0  3 13  2  0  0  1 37  0  5  0  0  1  3
  2    2   2  57   0 195   0 179   0   0  0  0  9  3  0  0  0  3  0 10  0  0  0  0
  7    2  13   3   0   0 184   0   0   0  0  0  1  0  0  0  0  0  0  7  1  0  0  0
  9    0   0   0   0   1   0   0 120   0  0  0  1  0  0  0  0  0  0  1  0  0  0  0
  11   1   0   0   0   0   0   0   0   0 99  0  0  0  0  0  0  0  0  0  0  0  0  0
  10  98   0   0   0   5   0   0   0   0  0  0  0  1  0  1  0  0  0  1  0  0  0  0
  8   18   3   0  21   1   0   3   0   4  0 98  0  3  0  5  0  0  0  1  0  0  0  0
  0    0   0   0   0   0   1   0   0   0  0  0  0 94  0  4  0  0  0  0  0  0  0  0
  6   14   0   6   0   2   0   3   0   0  0  0 24 13 91  0  0  0  0 68  0  0  0  0
  14   1   0   0   0   0   0   0   0   1  0  0  0  1  0  0 74  0  0  0  0  0  0  0
  12  18   0   2   0   2   0   0   0   0  0  0  0  0  0  0  0 72  0  0  0  0  0  0
  13  10   0   0   0   0   0   0   0   0  0  0  6  0  0  0  1  3 70  1  0  0  0  0
  15   1  12   2   0   0   0   0   0   0  0  0  0  0  0  0 13  0  0  0 48  0  0  0
  17   0   0   0   0   0   0   0   0   0  0  0  0  7  0  0  0  0  0  0  0 40  0  0
  18   0   0   0   0   0   0   0   0   0  0  0  0  0  0  0  0  0  0  0  0  0 34  0
  16  15   8   0   0   0   0   0   0   0  0  0  0  0  0  3  0  0  0  0  0  0  0 32