Entering edit mode
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,
nThreads = 11)
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
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?