run_viper: VIPER wrapper

In saezlab/dorothea_R: Collection Of Human And Mouse TF Regulons

Description

This function is a convenient wrapper for the viper function using DoRothEA regulons.

Usage

run_viper(input, regulons, options = list(), tidy = FALSE)

Arguments

input	An object containing a gene expression matrix with genes (HGNC/MGI symbols) in rows and samples in columns. The object can be a simple matrix/data frame or complexer objects such as ExpressionSet or Seurat objects.
regulons	DoRothEA regulons in table format.
options	A list of named options to pass to viper such as minsize or method. These options should not include, eset or regulon.
tidy	Logical, whether computed tf activities scores should be returned in a tidy format.

Value

A matrix of normalized enrichment scores for each TF across all samples. Of note, if you provide a Seurat object as input the function will return also a Seurat object with a new assay called dorothea. For all other inputs the function will return a matrix. If tidy is TRUE the normalized enrichment scores are retured in a tidy format (not supported for Seurat objects).

Examples

# use example gene expression matrix from bcellViper package
library(bcellViper)
data(bcellViper, package = "bcellViper")
# acessing (human) dorothea regulons
# for mouse regulons: data(dorothea_mm, package = "dorothea")
data(dorothea_hs, package = "dorothea")
# run viper
tf_activities <- run_viper(dset, dorothea_hs,
                          options =  list(method = "scale", minsize = 4,
                          eset.filter = FALSE, cores = 1,
                          verbose = FALSE))

saezlab/dorothea_R documentation built on April 3, 2020, 1 a.m.