Lately, I have tried to perform WGCNA analysis in R with WGCNA package against a RNA-seq dataset. And in this process, I found the numbers of color modules displayed by labeledHeatmap and Epigegene Network were different. And the code that produced those two plots are:
labeledHeatmap(Matrix = moduleTraitCor, xLabels = colnames(design), yLabels = names(MEs), ySymbols = names(MEs), colorLabels = FALSE, colors = greenWhiteRed(50), textMatrix = textMatrix, setStdMargins = FALSE, cex.text = 0.5, zlim = c(-1,1), main = paste("Module-trait relationships")) plotEigengeneNetworks(MET, "", marDendro = c(0,4,1,2), marHeatmap = c(3,4,1,2), cex.lab = 0.8, xLabelsAngle = 90)
The plots I got are as following:
As you can see, there are 32 modules in labeledHeatmap and only 31 in Epigengene Network. So, what's the secret here or did I miss something? Thanks for your time!!