...
Sage Coexpression, organic correlation algorithm | Sage Coexpression, WGCNA correlation algorithm | WGCNA Coexpression | WGCNA Coexpression, blocksize<=600
| Ssage Coexpr, WGCNA corr, TOM, and clust
|
|
|---|
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
| 5.9GB |
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 | 5:45.94 |
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 | | |
| 
| 
| 
 Image Removed
|
|
Module difference vs. this baseline (% of pairs for which the algorithm disagrees with this baseline on module co-membership):
| 0
| 7 4.4% 3% | 9.1%
6.6% | |
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 brown blue green red turquoise red blue grey black pink magenta grey purple greenyellow purple tan salmon cyan midnightblue lightcyan midnightblue grey60
2 403 11 3 427 0 0 0 0 15 0 0 0 0 0 0 25 0 0 0 0 0 0 0 0 0 0
3 11 357 0 0 1 6 0 0 0 24 0 1 0 405 0 0 0 150 0 0 0 0 0 0 0 0 0 0 0 4 0 0
1 1 4 352 1 2 130 0 4 0 16 1 61 0 1 0 0 19
4 4 6 0 302 0 4 0 0 347 0 0 0 0 0 0 0 6 0 0 0 0 0 0 0 0 1 0 0
5 1 0 2 0 298 1 0 0 0 6 0 3 0 0 0 0 0
8 0 4 0 0 195 2 0 0 0 2 0 0 0 0 0 0 0
6 4 3 132 39 5 10 0 0 0 13 0 7 1 0 0 0 0
9 0 1 0 0 0 0 122 0 0 0 0 0 0 0 0 0 0
7 0 67 1 2 1 2 0 115 0 9 0 0 0 0 0 0 8
12 0 0 0 0 0 0 0 0 99 1 0 0 0 0 0 0 0
11 0 0 0 0 0 99 0 0 0 3 0 1 0 0 0 0 0
10 0 0 1 3 2 96 0 0 0 10 0 1 0 0 0 0 1
0 0 0 0 0 0 0 0 0 0 80 0 14 41 0 0 0 1
15 0 0 0 0 0 0 0 0 0 0 77 0 0 0 0 0 0
14 1 0 0 0 77 0 0 0 0 2 0 0 0 0 0 0 0
13 0 6 0 42 0 2 0 1 0 0 0 0 0 39 0 0 1
16 0 0 0 0 0 0 0 0 0 3 0 0 1 0 37 0 0
17 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 34 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 | blue brown yellow green turquoise red black grey magenta pink purple greenyellow tan salmon cyan midnightblue lightcyan
brown 410 13 2 0 4 4 2 18 0 6 0 0 0 0 0 0 0
yellow 3 408 3 0 8 0 0 1 0 0 0 0 0 2 0 0 0
green 3 1 340 0 5 30 10 0 0 0 0 0 0 0 0 0 0
turquoise 1 0 2 338 7 0 194 1 0 0 0 0 0 38 0 0 0
red 2 13 5 3 312 13 2 7 0 2 0 1 0 3 0 0 0
blue 0 1 0 33 219 235 0 0 0 0 0 0 0 0 0 0 0
grey 15 7 6 1 17 0 5 131 0 4 1 0 2 6 1 0 4
black 0 0 0 0 0 0 0 0 122 0 0 0 0 0 0 0 0
pink 2 0 1 0 4 0 0 0 0 113 0 0 0 0 0 0 0
magenta 0 0 0 0 0 0 0 0 0 0 99 0 0 0 0 0 0
purple 0 0 0 0 4 0 0 0 0 0 0 80 0 3 0 0 0
greenyellow 0 0 0 0 0 0 0 0 0 0 0 0 78 0 0 0 0
salmon 0 0 38 0 2 0 0 0 0 0 0 0 0 0 0 0 0
cyan 0 0 0 0 0 0 0 2 0 0 0 0 0 0 35 0 0
lightcyan 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 34 0
tan 0 0 0 0 0 0 0 16 0 0 0 0 0 0 0 0 27
midnightblue 14 0 1 0 8 4 0 7 0 1 0 0 0 0 0 0 0  Image Removed |
 Image Added
| <<< Sage dendrogram and modules, with WGCNA modules underneath, for comparison.
WGCNA dendrogram and modules, with Sage modules underneath, for comparison. >>>>
|  Image Removed Image Added
|  Image Removed |
|
TOM Comparison
We captured the TOM matrices from the Sage and WGCNA algorithms. The values are the same (except for some small differences which look like rounding/accuracy differencesdown to the machine's floating point resolution) as illustrated in this scatter plot: 
Image Removed
Image Added
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
...
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:
Image Added
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 | | |
