Entering edit mode
5.4 years ago
user56
▴
20
While using this code in WGCNA:
dynamicMods = cutreeDynamic(dendro = geneTree, distM = dissTOM,
deepSplit = 2, pamRespectsDendro = FALSE,
minClusterSize = minModuleSize);
dynamicColors = labels2colors(dynamicMods)
GS1= as.numeric(cor(datExpr, use = "p"))
GeneSignificance=abs(GS1)
ModuleSignificance=tapply(GeneSignificance, dynamicColors, mean, na.rm=T)
plotModuleSignificance(GeneSignificance,dynamicColors)
I'm getting error:
Error in tapply(GeneSignificance, colorh1, mean, na.rm = T) : arguments must have same length
I guess number of rows in GeneSignificance and dynamicColors is not the same, or there is not enough colors? Am I understanding this right? How can I filter GeneSignificance?