|
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 |
} |