champ.EpiMod: champ.EpiMod() infer differential methylation interactome...

Description Usage Arguments Value Author(s) References Examples

Description

champ.EpiMod() used FEM package to identify interactome hotspots of differential promoter methylation. By "interactome hotspot" we mean a connected subnetwork of the protein interaction network (PIN) with an exceptionally large average edge-weight density in relation to the rest of the network. The weight edges are constructed from the statistics of association of DNA methylation with the phenotype of interest. Thus, the EpiMod algorithm can be viewed as a functional supervised algorithm, which uses a network of relations between genes (in our case a PPI network), to identify subnetworks where a significant number of genes are associated with a phenotype of interest (POI). We call these "hotspots" also Functional Epigenetic Modules (FEMs). You can get more detailed information in FEM package.

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
    champ.EpiMod(beta=myNorm,
                 pheno=myLoad$pd$Sample_Group,
                 nseeds=100,
                 gamma=0.5,
                 nMC=1000,
                 sizeR.v=c(1,100),
                 minsizeOUT=10,
                 resultsDir="./CHAMP_EpiMod/",
                 PDFplot=TRUE,
                 arraytype="450K")

Arguments

beta

A matrix of values representing the methylation scores for each sample (M or B). Better to be imputed and normalized data. (default = myNorm)

pheno

This is a categorical vector representing phenotype of factor wish to be analysed, for example "Cancer", "Normal"... Tow or even more phenotypes are allowed. (default = myLoad$pd$Sample_Group)

nseeds

An integer specifying the number of seeds and therefore modules to search for. (default = 100)

gamma

A parameter of the spin-glass algorithm, which determines the average module size. Default value generally leads to modules in the desired size range (10-100 genes). (default = 0.5)

nMC

Number of Monte Carlo runs for establishing statistical significance of modularity values under randomisation of the molecular profiles on the network. (default = 1000)

sizeR.v

Desired size range for modules. (default = c(1,100))

minsizeOUT

Minimum size of modules to report as interesting. (default = 10)

resultsDir

The directory where PDF files would be saved. (default = "./CHAMP_QCimages/")

PDFplot

If PDFplot would be generated and save in resultsDir. (default = TRUE)

arraytype

Choose microarray type is 450K or EPIC. (default = "450K")

Value

EpiMod.o

A data frame of all probes with an adjusted p-value for significance of differential methylation containing columns for logFC, AveExpr, t, P.Value, adj.P.Val, B, C_AVG, T_AVG, deltaBeta, CHR, MAPINFO, Strand, Type, gene, feature, cgi, feat.cgi, UCSC_CpG_Islands_Name, DHS, Enhancer, Phantom, Probe_SNPs, Probe_SNPs_10, you can turn to FEM package for more informations.

Author(s)

Teschendorff, A
adapted by Yuan Tian

References

Y J, M W and AE T (2014). A systems-level integrative framework for genome-wide DNA methylation and gene expression data identifies differential gene expression modules under epigenetic control. Bioinformatics, 30(16), pp. 2360-2366.

Examples

1
2
3
4
5
6
    ## Not run: 
        myLoad <- champ.load(directory=system.file("extdata",package="ChAMPdata"))
        myNorm <- champ.norm()
        myEpiMod <- champ.EpiMod()
        
## End(Not run)

gaberosser/ChAMP documentation built on May 9, 2019, 2:22 a.m.