Any possiblity applying
0
0
Entering edit mode
adR ▴ 30
@do-it-23093
Last seen 8 months ago
Germany, München

Hi Dears

I would like to ask if there is a possibility to use "WGCNA" to visualize the co-expression network, gene modules, connectivity, and Hub genes present only on the top right block of the correlation matrix? My objective to detect the possible correlations among genes of two tissues of the same subjects. Below showed my data type and the block of the correlation matrix I would like to work on. Thanks!

library(Hmisc)
library(WGCNA)
set.seed(1)
matr <- matrix(rnorm(250), nrow = 50)
rownames(matr) <- paste0("g_",1:50) ## genes 
colnames(matr) <- LETTERS[1:5] ## samples 

set.seed(2)
matr2 <- matrix(rnorm(250), nrow = 50)
rownames(matr2) <- paste0("ge_",1:50) ## genes 
colnames(matr2) <- LETTERS[1:5] ## samples
test = Hmisc::rcorr(t(rbind(matr,matr2)))
n = nrow(matr2)
final <- test$r[1:n, (n+1):(2*n)]
sessionInfo( )
WGCNA • 593 views
ADD COMMENT
1
Entering edit mode

Potential cross-post: https://www.biostars.org/p/471736/

ADD REPLY

Login before adding your answer.

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