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, 
                                    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
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
ADD COMMENT

Login before adding your answer.

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