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

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).