Search
Question: Intramodular connectivity RNA-seq WGCNA
0
gravatar for adriana.gallego.02
5 months ago by
adriana.gallego.020 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
length(inModule)
modTOM = TOM[inModule, inModule];
dimnames(modTOM) = list(modProbes, modProbes)

TOM.mat = as.matrix(mod)

Thanks

ADD COMMENTlink modified 5 months ago by Peter Langfelder1.3k • written 5 months ago by adriana.gallego.020
0
gravatar for Peter Langfelder
5 months ago by
United States
Peter Langfelder1.3k 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 5 months ago by Peter Langfelder1.3k

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

Thanks

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

Help
Access

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