-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy pathrba_reactome_analysis.Rd
185 lines (169 loc) · 7.73 KB
/
rba_reactome_analysis.Rd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/reactome_analysis.R
\name{rba_reactome_analysis}
\alias{rba_reactome_analysis}
\title{Reactome Over-Representation or Expression Analysis}
\usage{
rba_reactome_analysis(
input,
input_format = NULL,
projection = TRUE,
interactors = FALSE,
species = NULL,
sort_by = "ENTITIES_PVALUE",
order = "ASC",
resource = "TOTAL",
p_value = 1,
include_disease = TRUE,
min = NULL,
max = NULL,
...
)
}
\arguments{
\item{input}{A vector, data frame, matrix or a local file path or URL
that points to your data. See "Details section" for more information of
how to organize and supply your input.}
\item{input_format}{(Optional) This function will automatically identify
your supplied input's format. But in case of unexpected issues or if you
want to be explicit, set this argument to one of:\itemize{
\item "table": If you supplied a data frame or matrix as input.
\item "vector": If you supplied a simple vector (numeric or character) as
input.
\item "file": If you supplied a local file path pointing to a
correctly-formatted text file.
\item "url": If you supplied a URL pointing to a correctly-formatted
text file.}}
\item{projection}{Logical (default = TRUE) Should non-human identifiers
be projected to their human equivalents? (using Reactome orthology data)}
\item{interactors}{Logical (default = FALSE) Should IntAct interaction data
be used to increase the analysis background?}
\item{species}{Numeric or Character: NCBI Taxonomy identifier (Human
is 9606), species name (e.g. "Homo sapiens") or Reactome DbId (e.g
Homo sapiens is 48887). See
\code{\link{rba_reactome_species}} or
\href{https://reactome.org/content/schema/objects/Species/}{Reactome
Data Schema: Entries: Species}. Note that you cannot supply the species
parameter when projection parameter is TRUE.}
\item{sort_by}{Sort the result based on what column? available choices
are: "NAME", "TOTAL_ENTITIES", "TOTAL_INTERACTORS", "TOTAL_REACTIONS",
"FOUND_ENTITIES", "FOUND_INTERACTORS", "FOUND_REACTIONS", "ENTITIES_RATIO",
"ENTITIES_PVALUE", "ENTITIES_FDR" or "REACTIONS_RATIO"}
\item{order}{Sort Order. Can be either "ASC" (default) or "DESC".}
\item{resource}{Filter results based on the resource. Default is "TOTAL",
available choices are: "TOTAL", "UNIPROT", "ENSEMBL", "CHEBI", "IUPHAR",
"MIRBASE", "NCBI_PROTEIN", "EMBL", "COMPOUND", "ENTITIES_FDR" or
"PUBCHEM_COMPOUND".}
\item{p_value}{Set a P value threshold. Only results with P value equal to
or less than your supplied threshold will be returned. (default = 1,
Meaning no P value filtering)}
\item{include_disease}{Logical (default = TRUE) Should the disease pathways
be included in the results?}
\item{min}{(numeric) Minimum number of entities that a pathways should have
to be included in the results.}
\item{max}{(numeric) Maximum number of entities that a pathways should have
to be included in the results.}
\item{...}{rbioapi option(s). See \code{\link{rba_options}}'s
arguments manual for more information on available options.}
}
\value{
List containing the results and information of your analysis. Note
that you can use the token returned in the "summary" sub-list of the
results (i.e. results$summary$token) to retrieve your results later or
in other Reactome analysis functions.
}
\description{
Using this function, you can perform Reactome Analysis In a convenient way.
The Analysis Type will be chosen depending on your supplied
input:\enumerate{
\item If you supply a vector or a single-columned table,
"Over-Representation" analysis will be performed.
\item If you supply a multi-column table, with the first column being
molecules identifiers and the rest being numeral expression values,
"Expression" analysis will be performed.}
See the details section for the accepted input types and format.
}
\details{
You can supply your table or vector input in numerous formats:\enumerate{
\item A R object which can be data frame, matrix or a simple vector.
\item A path to a local text file in your device that contains the molecules
data. (The file should be formatted correctly, see below.)
\item A URL pointing to a text file on the web that contains the molecules
data. (The file should be formatted correctly, see below.}
If you supply a text file (as a local file path or URL), it should be
in TSV (Tab-Separated Values) format; Column names should start with "#"
character. Note that if you are providing the file for
"Over-Representation" analysis (i.e. Single columned-data) this header
line is optional and will be used as your 'Sample Name', otherwise it is
required.
\cr Also, form the "summary" element in the function's output, you can see
how Reactome Interpreted your input and subsequently the type of analysis
that has been performed.
\cr There is no strict criteria about the type of your molecules Identifiers,
Reactome will Map the IDs to it's internal database entities.
Nevertheless, You can check if all your identifiers has been found in
"identifiersNotFound" element in the function's output.
\cr After Any Analysis, Reactome will associate a token to your analysis.
It can be later used to in function that requires the token (e.g to retrieve
the analysis results, download pdf).
\cr Note that Reactome will store your token for only 7 days. You can
download your full results with
\code{\link{rba_reactome_analysis_download}}, and re-import it anytime to
reactome (using \code{\link{rba_reactome_analysis_import}}) to generate
a new token.
}
\section{Corresponding API Resources}{
"POST https://reactome.org/AnalysisService/identifiers/form"
\cr "POST https://reactome.org/AnalysisService/identifiers/url"
\cr "POST https://reactome.org/AnalysisService/identifiers/form/projection"
\cr "POST https://reactome.org/AnalysisService/identifiers/url/projection"
}
\examples{
\dontrun{
rba_reactome_analysis(input = c("p53", "BRCA1", "cdk2", "Q99835", "CDC42"))
}
\dontrun{
rba_reactome_analysis(input = "c:/rbioapi/genes.txt")
}
\dontrun{
rba_reactome_analysis(input = "https://qazwsx.com/genes.txt")
}
}
\references{
\itemize{
\item Marc Gillespie, Bijay Jassal, Ralf Stephan, Marija Milacic, Karen
Rothfels, Andrea Senff-Ribeiro, Johannes Griss, Cristoffer Sevilla, Lisa
Matthews, Chuqiao Gong, Chuan Deng, Thawfeek Varusai, Eliot Ragueneau,
Yusra Haider, Bruce May, Veronica Shamovsky, Joel Weiser, Timothy Brunson,
Nasim Sanati, Liam Beckman, Xiang Shao, Antonio Fabregat, Konstantinos
Sidiropoulos, Julieth Murillo, Guilherme Viteri, Justin Cook, Solomon
Shorser, Gary Bader, Emek Demir, Chris Sander, Robin Haw, Guanming Wu,
Lincoln Stein, Henning Hermjakob, Peter D’Eustachio, The reactome
pathway knowledgebase 2022, Nucleic Acids Research, 2021;, kab1028,
https://doi.org/10.1093/nar/gkab1028
\item Griss J, Viteri G, Sidiropoulos K, Nguyen V, Fabregat A,
Hermjakob H. ReactomeGSA - Efficient Multi-Omics Comparative Pathway
Analysis. Mol Cell Proteomics. 2020 Sep 9. doi: 10.1074/mcp. PubMed
PMID: 32907876.
\item \href{https://reactome.org/AnalysisService/}{Reactome Analysis
Services API Documentation}
\item \href{https://reactome.org/cite}{Citations note on Reactome website}
}
}
\seealso{
Other "Reactome Analysis Service":
\code{\link{rba_reactome_analysis_download}()},
\code{\link{rba_reactome_analysis_import}()},
\code{\link{rba_reactome_analysis_mapping}()},
\code{\link{rba_reactome_analysis_pdf}()},
\code{\link{rba_reactome_analysis_species}()},
\code{\link{rba_reactome_analysis_token}()}
Other "Enrichment/Over-representation":
\code{\link{rba_enrichr}()},
\code{\link{rba_mieaa_enrich}()},
\code{\link{rba_panther_enrich}()},
\code{\link{rba_string_enrichment}()},
\code{\link{rba_string_enrichment_image}()}
}
\concept{"Enrichment/Over-representation"}
\concept{"Reactome Analysis Service"}