champ.ebGSEA: Empirical Bayes GSEA method.

View source: R/champ.ebGSEA.R

champ.ebGSEAR 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)

YuanTian1991/ChAMP documentation built on Feb. 21, 2023, 1:13 p.m.