doGSEA: Runs preranked gene set enrichment analysis.

Description Usage Arguments Value Examples

View source: R/doGSEA.R

Description

The function does gene set enrichment analysis(GSEA) for each cluster with its gene ranks. It will return cell type for each cluster.

Usage

1
doGSEA(cluster_list, db = "PanglaoDB_list", minSize = 15, maxSize = 500)

Arguments

cluster_list

A ranked gene list for each cluster.

db

The cell type data base to use. For single cell data, we provide three data base, the first one is 'PanglaoDB' data base (db = 'PanglaoDB_list'), the second one is 'GSEA' data base (db = 'GSEA_list') and the third one is the reference genome for Arabidopsis (db = 'TAIR_list'). It can also be a path to the new (referential) data base that hope to be used, the file must be 'rds' format.

minSize

An integer value. Minimal size of a gene set to test. All pathways below the threshold are excluded.

maxSize

An integer value. Maximal size of a gene set to test. All pathways above the threshold are excluded.

Value

Cell type for each cluster, and its NSF and padj.

Examples

1
2
3
4
pbmc_example <- scqc(small_pbmc_rna, min.cells = 1, min.features = 10, nfeatures = 100, npcs = 10)
          pbmc_example <- doClustering(pbmc_example, dims = 1:10, k.param = 5, resolution = 0.75)
          cluster_list <- getFC(pbmc_example, min.pct = 0.25, test.use = "MAST")
          cluster_celltype <- doGSEA(cluster_list = cluster_list, minSize = 5)

cailab-tamu/scTypeGSEA documentation built on July 15, 2020, 10:56 a.m.