Any possiblity applying
0
0
Entering edit mode
Do it! ▴ 10
@do-it-23093
Last seen 29 days 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 • 189 views
ADD COMMENT
1
Entering edit mode

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

ADD REPLY

Login before adding your answer.

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