Closed:Problem with ModuleSignificance - how to filter genes?
0
0
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?

wgcna modulesignificance • 237 views
ADD COMMENT
This thread is not open. No new answers may be added
Traffic: 2552 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6