- Version: 1.0
- Date: 09/20/2015
- Authors: Xiaoke Ma, Zaiyi Liu, Lin Gao, Peizhuo Wang, Wanxin Tang
- Maintainer: Xiaoke Ma ([email protected])
- Depends: R (>2.15.1)
- License: GPL (>=2)
##Introduction The EMDM algorithm is designed for identifying methylated gene modules in multiple differential networks. It takes as inputs multiple edge-weighted gene networks and a set of prior probabilities representing the importance of a gene for the conditions under study (no priori information is also practicable). Along with network topological features, the prior probabilities are used to rank and select seeds to initialize module search. ##Software tutorial
- step 1: construct the differential co-methylation (co-expression) networks
source("dif_network_construction.r")
dif_network = dif_network_construction(genemeth, pvalue, delta)
#genemeth: the matrix denoting methyaltion data where row corresponds to gene and column to sample
#pvalue: the pvalue of gene differential methylation between two cohorts
# delta : the cutoff parameter for co-methylation (co-expression) network
- step 2: Functional epigenetic module discovery
* step 2.1 loading the multiple differential networks
lname = load("test-data/1.RData");
* step 2.2 gene ranking
source("./rankgene_new.r")
ggrank = rankgene_new(networks,0.2)
* step 2.3 seed selection
ggrank[ggrank>2] = 2;
ggrank[ggrank< -2] = -2;
generank = rowSums(ggrank);
seedgene = order(generank,decreasing=TRUE);
seeds = seedgene[1:500]
* step 2.4 module discovery
source("./nmodule.R");
fem_modules =nmodule(networks,seeds);
*For more details, refer to main_fem-DM_test.r