GSEA_gene_single: Do GSEA for a single geneSet based on random permuted or...

Description Usage Arguments Value Examples

Description

It'll do GSEA Tests for a single geneSet according to the signal 2 noise values for all of the genes or other values,such as foldchange. It'll return ES ,the running ES score for all of the genes, a numeric vector . rand_es ,the maximal running ES score for each permutation, a numeric vector. p which calculated by formule:p <- 1-pnorm((max(ES)-mean(rand_es))/sd(rand_es)) ismem is a logical vector shows whether the gene belong to the current geneSet or not.

Usage

1
2
GSEA_gene_single(gene_values = rnorm(100), gene_names = 1:100,
  this_geneSet = sample(1:100, 10), geneSet_name = "geneset1", n = 1000)

Arguments

gene_values

A vector of values for each gene,most of time it's signal 2 noise.

gene_names

A vector of names for each gene

this_geneSet

A vector which contains all of the genes in this geneSet.

n

Times for random permuted or bootstraps by genes,default:1000

this_geneSet

The name of current geneSet,default:geneset1

Value

a list contain 4 elements: ES,rand_es, p and ismem

Examples

1

jmzeng1314/geneSet documentation built on May 19, 2019, 2:57 p.m.