This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Problem with ModuleSignificance - how to filter genes?

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?

wgcna modulesignificance

0 answers

No answers yet.

Log in to answer this question.