introduction
This tutorial is made for users who want to quickly use CACIMAR to identify conserved gene regulatory networks from different modules. Input of this part is a table contains Source, SourceGroup, Target, TargetGroup.
library(CACIMAR)
load(system.file("extdata", "gene_network.rda", package = "CACIMAR"))
head(mm_gene_network)
## Source SourceGroup Target TargetGroup
## 1 ENSMUSG00000003949 1 ENSMUSG00000021250 1
## 2 ENSMUSG00000003949 1 ENSMUSG00000022707 1
## 3 ENSMUSG00000003949 1 ENSMUSG00000025892 1
## 4 ENSMUSG00000003949 1 ENSMUSG00000026989 1
## 5 ENSMUSG00000003949 1 ENSMUSG00000036523 1
## 6 ENSMUSG00000003949 1 ENSMUSG00000038248 1
Run the code
ConservedNetworks_ct <- Identify_ConservedNetworks(OrthG_Mm_Zf, mm_gene_network, zf_gene_network, 'mm', 'zf')
## Using NCS as value column: use value.var to override.
Heatmap_Cor(ConservedNetworks_ct[[2]], cluster_cols=F, cluster_rows=F, Color1 = c(rgb(102/255,46/255,115/255), rgb(31/255,153/255,139/255), rgb(251/255,232/255,48/255)))