Question: Intramodular connectivity RNA-seq WGCNA
gravatar for adriana.gallego.02
16 months ago by
adriana.gallego.0210 wrote:

Cordial greeting

Someone could help me with the script to calculate intramodular connectivity for RNA-seq data?

This is the code that I have so far, is this correct?

module = c("indianred4");

# Select module probes
probes = colnames(datExpr)
inModule = (moduleColors==module);
modProbes = probes[inModule]; 
# Recalculate intramodular connectivity (IMConn)
IMConn = softConnectivity(datExpr[, modProbes]);

# Select the corresponding Topological Overlap
modTOM = TOM[inModule, inModule];
dimnames(modTOM) = list(modProbes, modProbes)

TOM.mat = as.matrix(mod)


ADD COMMENTlink modified 16 months ago by Peter Langfelder1.6k • written 16 months ago by adriana.gallego.0210
gravatar for Peter Langfelder
16 months ago by
United States
Peter Langfelder1.6k wrote:

The code should work. RNA-seq, after suitable normalization and transformation, can be analyzed just any other data. You may want to look at the WGCNA FAQ, especially point 4. Then simply follow the WGCNA tutorials or the code in the book (which you seem to be following already).

ADD COMMENTlink written 16 months ago by Peter Langfelder1.6k

Thanks Dr. Langfelder, I have one last question, how can I visualize this intramodular connectivity, for example using Cytoscape?


ADD REPLYlink written 16 months ago by adriana.gallego.0210
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 497 users visited in the last hour