WGCNA: minModuleSize is not responding correctly
0
0
Entering edit mode
@f68bdf4c
Last seen 3 days ago
Norway

I am running the WGCNA, but getting an error. As I can see that the minimum module size that I am getting is 1, but by default it should not give anything less than 30.

modules <- blockwiseModules(transcriptomics,
power = 9,
networkType = "signed",
TOMType = "signed",
corType = "bicor",
numericLabels = TRUE,
saveTOMs = save_TOM,
saveTOMFileBase = "TOM",
verbose = 3,
maxBlockSize=8000,

rownames(modules$MEs) <- rownames(metatranscriptomics) names(modules$colors) <- colnames(metatranscriptomics)
hubs_M <- chooseTopHubInEachModule(metatranscriptomics, modules$colors, power = 9, omitColors = "0") Error in hubs[m] <- colnames(adj)[hub] : replacement has length zero  table(modules$colors)

0     1     2     3     4     5     6     7     8     9    10    11
60409  8629  3608  3204  1909   483   111    88    32    31     4     1

table(names(modules\$colors) == colnames(metatranscriptomics))

TRUE
78509


I have looked at this post- https://bioinformatics.stackexchange.com/questions/17598/wgcna-hub-genes-error

But I can't solve this issue. Any suggestions?

wgcna Network WGCNA minModuleSize • 61 views