Description Usage Arguments Examples
View source: R/pathway.analysis.R
This function allows user to run pathway analysis on the identified significat drugs, therefore to detect drug mechanism related pathways, and the CEGs in the pathways.
1 2 3 4 5 6 7 | pathway.analysis(
drug.ident.res = drug.ident.res,
pathway.list = NULL,
pathway.list.path = NULL,
drug.FDR.cutoff = 0.1,
CEG.threshold = 0.05
)
|
drug.ident.res |
Dr. Insight drug identification analysis results. Output of function drug.ident(). |
pathway.list |
The pathways used to analyze drug mechanism. Should be a list of pathways where the names of lists are pathway names. Dr.Insight provides NIH pathway interaction database (PID) pathways, which can be called by data("pathway.PID"). If you want to analyze with other pathways (from GSEA MsigDB), please use the parameter: pathway.list.path. |
pathway.list.path |
The local path and file name of the pathways downloaded from GSEA MsigDB http://software.broadinstitute.org/gsea/msigdb/collections.jsp, if pathway.list is NULL. Should be in full directory and file name (See instructions in vignette). |
drug.FDR.cutoff |
The FDR threshold to select significant repurposed drugs by Dr. Insight for pathway analysis. |
CEG.threshold |
The p value threhold for selecting significant consistently differential expressed genes (CEGs). |
1 2 3 4 5 6 7 8 9 10 | ## get the Dr. Insight drug identification results
drug.ident.res = drug.ident(query.data = example.disease, cmap.ref.profiles = example.drug.profiles,
repurposing.unit = "treatment", connectivity = "negative")
## load in example pathway data
data("example.pathway")
## Performe pathway analysis (for the drugs that are identified by ident.drug())
path.analysis.res = pathway.analysis(drug.ident.res = drug.ident.res,
pathway.list = example.pathway, drug.FDR.cutoff = 0.5)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.