Hello All,
How to choose the right parameters for networkType and TOMType when using WGCNA's blockwiseModules ?
net = blockwiseModules(datExpr, power = 12, networkType = "signed hybrid", TOMType = "signed")
For networkType, it offers unsigned or signed or signed hybrid. For TOMType, it also has unsigned or signed. I am still very confused after reading through the manual, does anyone have an easier way to understand those?
Thanks a lot.
Thank you for the explanation.
So, in that case, the "signed hybrid" network will only keep the positively correlated ones, and neglect the negative ones?
The shared file is really informative. Thanks a lot.
Yes. But the signed network pretty much does this too; since cor=0 -> (1 + cor)/2 = 0.5
for a power of 10 the resulting adjacency is around 0.001 which is not all that much different from the hybrid value of 0.
Indeed. I have RNAseq data from 140 samples as well as some trait data (e.g., plant height, flowering time...). What I want to do is to identify which modules is highly related to which traits through WGCNA.
So practically speaking, which types of network and TOMtype should I choose will make more sense biologically?
Sorry to bug you with further questions.
For so many samples, I'd recommend using
signed
for both network and topological overlap. It's simple and will result in more clusters thanunsigned
.That's very good to know.
As the WGCNA authors suggest remove "non-varying genes, which usually represent noise". I only pick the top 15,000 variable genes across the samples to run WGCNA, is it sufficient?
I am pretty new to this, thank you so much for your Q&As!
Absolutely.
That's awesome. Thank you. Really appreciate it!