**130**wrote:

Suppose these two matrices

```
set.seed(1)
x <- matrix(rnorm(1000), nrow=100, ncol=10)
y <- matrix(rnorm(1000), nrow=100, ncol=10)
```

represent the genes (as rows) and genotypes (as columns) for two species.

How to construct a gene co-expression network (by calculating the among genotypes correlation for every pair of genes within each species and the correlations are the edges of the network) . So in total there will be two networks A and B.

Is there any R package to do this and then calculate the expression conservation score (EC score) of each gene in the two gene co-expression networks.

By calculating this EC score for each gene, I assume we would be able to measure the degree of similarity between a candidate's gene neighbours in network "A" and same gene's neighbour in network "B".