Biostar Beta. Not for public use.
Construct gene coexpression network
0
Entering edit mode
2.1 years ago

How can I construct a gene co-expression network in R from the expression value I have in Excel sheet?

Please tell me what packages are the helpful or share any tutorial for the construction of the co-expression network in R.

Thank you

R gene • 1.8k views
1
Entering edit mode

Please show some effort of trying. Post input data, and expected output? At the least post some links to published papers. As it stands this post is too broad and unclear.

3
Entering edit mode
18 months ago
F ♦ 3.4k
Iran

Reading your normalized gene expression data in r which has been saved as .txt and genes are in rows and samples are in columns

    mycounts <- read.table("data.txt", header = T, sep = "\t",row.names = 1)
# watching the dimension of matrix
dim(mycounts)


Using this R package for instance, you must put your working directory in where gene3.R and your expression files are. You must also install these packages in R that all help you to illustrate your network in Cytoscape. Download geneie3.R source from this link clicking on R/randomForest part

library(GENIE3)
library(igraph)
library(RCy3)
library(Rgraphviz)
weight.matrix <- GENIE3(mycounts)

Gsi <- graph.data.frame(edge_listsi,directed = F)

Asi <- get.adjacency(Gsi,sparse = F,attr = "weight",type = "both")

g_arasi <- graph.adjacency(Asi,mode = "undirected",weighted = T)

g.cyto <- igraph.to.graphNEL(g_arasi)

cw = createNetworkFromGraph("net", graph=g.cyto)

displayGraph (cw)

0
Entering edit mode

Hi F I used your commands for my data but after "weight.matrix <- GENIE3(mycounts)" I got error:

Error in (function (classes, fdef, mtable)  :
unable to find an inherited method for function ‘GENIE3’ for signature ‘"data.frame"


my data are log2 of FPKM and I have -INF in some rows, could you help? Thanks

0
Entering edit mode

as.matrix(read.table("data.txt", header = T, sep = "\t",row.names = 1))


Genes are in rows and columns are your samples

0
Entering edit mode

yes, genes are in rows and samples are in columns. you think is it wise to replace -inf with 0?

Mamnun :)

0
Entering edit mode

Sorry no idea but you can try that and compare your resulted network with when you are taking log2 after adding up 1

:)

0
Entering edit mode

I have about 40K genes and 12 samples, should remove low expressed ones?

1
Entering edit mode

Of course because all of these genes are not informative even Cytoscape would get in trouble for loading such a big interactiobs. have you done any differential expression to reduce your genes ?

0
Entering edit mode

yes, I have the DE genes, but I think we should not include just DE genes and we need to filter the low expressed ones just according the normalized expression.right?