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