champ.GSEA: Do GSEA for DMP, DMR and other methylation data related...

In JoshuaTian/MyChAMP: Chip Analysis Methylation Pipeline for Illumina HumanMethylation450 and EPIC

Description

This function would do GSEA on the results of champ functions like DMP and DMR. However users may also add individual CpGs and genes in it. There are two method are incoporated into champ.GSEA function here. One is old Fisher Exact Test method, which will used information downloaded from MSigDB and do fisher exact test to calculated the enrichment status for each pathways. And the other is "goseq" method, which will use goseq package to correct the inequality between number of genes and number of CpGs, then do GSEA. User may assign parameter "method" to choose which method they want to use.

Usage

    champ.GSEA(beta=myNorm,
               DMP=myDMP,
               DMR=myDMR,
               CpGlist=NULL,
               Genelist=NULL,
               method="goseq"
               arraytype="450K",
               Rplot=TRUE,
               adjPval=0.05)

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)
DMP	Results from champ.DMP() function. (default = myDMP)
DMR	Results from champ.DMR() function. (default = myDMR)
CpGlist	Apart from previous parameters, if you have any other CpGs list want to do GSEA, you can input them here as a list. (default = NULL)
Genelist	Apart from previous parameters, if you have any other Gene list want to do GSEA. you can inpute them here as a list. (default = NULL)
method	Which method would be used to do GSEA?"goseq" or "fisher". (default = "goseq")
arraytype	Which kind of array your data set is? (default = "450K")
Rplot	If goseq method was chosen, should Probability Weight plot will be plotted. More information please check goseq package. (default = TRUE)
adjPval	Adjusted p value cutoff for all calculated GSEA result. (default = 0.05)

Value

For fisher Method:

Block	A data.frame contains all detected Blocks, with colnames as chr, start, end, value, area, cluster, indexStart, indexEnd, L, clusterL, p.value, fwer, p.valueArea, fwerArea. The result format is actually the same as Bumphunter, you may refer to Bumphunter packages to get more explaination about the result.
clusterInfo	When champ.Block() detection significant Blocks, a group of candidate Blocks would be detected out at first, this is the data frame of all candidate Blocks. The "TRUE" Blocks in above value are located in these candidate Blocks.
allCLID.v	The first step of detectiong methylation Blocks is to get each probes into a cluster(region). This value is the clustering result of each probes.
avbetaCL.m	The beta matrix for each cluster. The value is calculated by taking mean value of all probes located in each cluster.
posCL.m	Position of each cluster, which is calculated by average all probes' position in each cluster.
adjPval	P value cutoff for calculated GSEA results. (default = 0.05)

For goseq method, the returned value are:

category	GO pathway's index.
over_represented_pvalue	The p value for genes' over representing in this pathway.
under_represented_pvalue	The p value for genes' under representing in this pathway.(Not likely to be used)
numDEInCat	Numbers of Different Methylation Genes in this pathway.
numInCat	Numbers of all genes related to this pathway.
term	The short explaination for this pathway.
ontology over_represented_adjPvalue	The ajusted over representing p value with "BH" method. User may used this one to select qualitied Pathways.

Author(s)

Yuan Tian

Examples

    ## Not run: 
        myLoad <- champ.load(directory=system.file("extdata",package="ChAMPdata"))
        myNorm <- champ.norm()
        myDMP <- champ.DMP()
        myDMR <- champ.DMR()
        myGSEA <- champ.GSEA()
        
## End(Not run)

JoshuaTian/MyChAMP documentation built on May 7, 2019, 12:04 p.m.