WGCNA: minModuleSize is not responding correctly
0
0
Entering edit mode
@f68bdf4c
Last seen 24 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, 
                                    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?

Network minModuleSize WGCNA wgcna • 1.0k views
ADD COMMENT

Login before adding your answer.

Traffic: 561 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