Diff of /man/multiGSEA.Rd [000000] .. [72a7c3]

Switch to unified view

a b/man/multiGSEA.Rd
1
% Generated by roxygen2: do not edit by hand
2
% Please edit documentation in R/enrichment_functions.R
3
\name{multiGSEA}
4
\alias{multiGSEA}
5
\title{Calculate pathway enrichment for multiple omics layer.}
6
\usage{
7
multiGSEA(pathways, ranks, eps = 0)
8
}
9
\arguments{
10
\item{pathways}{Nested list containing all pathway features for the
11
respective omics layer.}
12
13
\item{ranks}{Nested list containing the measured and pre-ranked features for
14
each omics layer.}
15
16
\item{eps}{This parameter sets the boundary for calculating the p value.}
17
}
18
\value{
19
Nested list containing the enrichment scores for each given pathway
20
  and omics layer.
21
}
22
\description{
23
This function calculates GSEA-based enrichments scores for multiple omics
24
layer at once. Input pathways or gene sets have to be prepared in advance by
25
means of the function \code{\link[multiGSEA]{initOmicsDataStructure}}. The function uses pre-
26
ranked lists for each omics layer to calculate the enrichment score. The
27
ranking can be calculated by means of the function
28
\link[multiGSEA]{rankFeatures}.
29
}
30
\examples{
31
32
# Download pathway definition and extract features
33
pathways <- getMultiOmicsFeatures(dbs = c("kegg"), layer = c("transcriptome", "proteome"))
34
35
# load omics data and calculate ranks
36
data(transcriptome)
37
data(proteome)
38
ranks <- initOmicsDataStructure(c("transcriptome", "proteome"))
39
ranks$transcriptome <- rankFeatures(transcriptome$logFC, transcriptome$pValue)
40
names(ranks$transcriptome) <- transcriptome$Symbol
41
ranks$proteome <- rankFeatures(proteome$logFC, proteome$pValue)
42
names(ranks$proteome) <- proteome$Symbol
43
44
## run the enrichment
45
multiGSEA(pathways, ranks)
46
}