Hi,

I am using WGCNA to identify the modules with co-expressed genes and assign these modules to traits.

When I generate a labeled heatmap (please see the code below) to check the module-trait correlation, I see some modules correlating with some traits. However, when I check the module membership vs. gene significance for the correlating modules and traits from the heatmap, there is no correlation (please see an example of the code I used below).

Could someone please explain the reason for getting different results? How should I interpret this result?

Thanks!

The code I used for labeled heatmap:

labeledHeatmap(Matrix = moduleTraitCor,
               #xLabels = names(datTraits),
               xLabels = cond,
               #yLabels = names(MEs),
               yLabels =namey,
               #ySymbols = names(MEs),
               ySymbols = namey,
               colorLabels = FALSE,
               colors = blueWhiteRed(50),
               textMatrix = textMatrix, 
               setStdMargins = FALSE,
               cex.text = 0.9,
               zlim = c(-1,1))

The code I used for the Module membership vs. gene significance plot creation:

verboseScatterplot(abs(geneModuleMembership[moduleGenes, column]),
                   abs(geneTraitSignificance[moduleGenes, 1]),
                   xlab = paste("Module Membership in", module, "module"),
                   ylab = "Gene significance for fat",
                   main = paste("Module membership vs. gene significance\n"),
                   cex.main = 1.2, cex.lab = 1.2, cex.axis = 1.2, col = module)