champ.ebGSEA {ChAMP} | R Documentation

## Empirical Bayes GSEA method.

### Description

This is a newly created method for conduct bias-free GSEA from 450K or EPIC data set. This method use global test to detect significance of genes from DNA methylation data sets directly, instead of simply select genes mapped my DMPs pr DMRs. By applying this method, users could find GSEA without bias from inequality number of CpGs of genes, and detect some marginal significant genes for GSEA process. After global test, Empirical Bayes method would use wilcox test to enrich genes to pathways. Note that you can directly use champ.GSEA() to use this method, just need to set "method" parameter as "ebay" in champ.GSEA() to run this method.

### Usage

    champ.ebGSEA(beta=myNorm, pheno=myLoad$pd$Sample_Group, minN=5, adjPval=0.05, arraytype="450K", cores=1)

### 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` | User needs to provide phenotype information to conduct global test. (default = myLoad$pd$Sample_Group)
`minN` | Minium number of common genes threshold in one geneset and candidate gene list, if less than this value, the p value of this geneset would be set 1. (default = 5)
`adjPval` | Adjusted p value cutoff for all calculated GSEA result. (default = 0.05)
`arraytype` | Which kind of array your data set is? (default = "450K")
`cores` | Number of parallel threads/cores used to accelarate. (default = 1)

### Value

There are three list: GSEA contains all pathway's GSEA result in one list, and only significant pathways GSEA in another. EnrichedGene: contains enriched genes in each pathways. gtResult: global test result for each gene.

Below are columns for list GSEA:

`nREP` | Number of genes enriched in this pathway.
`AUC` | Area under curve from wilcox test.
`P(WT)` | P value detected for each pathway from Wilcox Test.
`P(KPMT)` | P value from Known Population Median Test
`adjP` | Adjusted P value for each pathway, using BH method.

### Author(s)

Yuan Tian, Danyue Dong

### Examples

        ## Not run:
            myLoad <- champ.load(directory=system.file("extdata",package="ChAMPdata"))
            myNorm <- champ.norm()
            myGSEA.ebGSEA <- champ.ebGSEA(beta=myNorm,pheno=myLoad$pd$Sample_Group,arraytype="450K")

    ## End(Not run)

[Package ChAMP version 2.22.0 Index]