alex.pathway | R Documentation |
Computes pairwise distances between subjects based on lesion profiles in genes associated with a specified pathway, and returns a figure with two panels: one showing lesion data and another showing expression data, both ordered based on the computed distances (useful for hierarchical clustering). It also returns the corresponding ordered data.
alex.pathway(alex.data, lsn.clrs = NULL, lsn.data, pathways, selected.pathway)
alex.data |
Output of the |
lsn.clrs |
Optional. A named vector of colors for lesion types. If not specified, default colors will be assigned using |
lsn.data |
A data frame of lesion data in GRIN-compatible format with the following columns:
|
pathways |
A data frame with three columns: |
selected.pathway |
A character string indicating the pathway of interest. |
This function identifies all genes associated with the specified pathway, extracts lesion and expression data for those genes, and computes pairwise distances between subjects based on their lesion profiles. It uses hierarchical clustering to order the subjects and visualizes lesion and expression matrices in two aligned panels. It also returns a data frame containing the ordered expression and lesion data for all genes in the pathway.
A list with the following element:
ordered.path.data |
A data frame with lesion and expression data for pathway genes, ordered according to hierarchical clustering (matching the order used in the plot). |
A figure with two panels will also be generated showing:
Lesion data of pathway genes across subjects
Expression data of the same genes across subjects
Both panels are ordered by subject similarity based on lesion profiles.
Abdelrahman Elsayed abdelrahman.elsayed@stjude.org, Stanley Pounds stanley.pounds@stjude.org
alex.prep.lsn.expr
, hclust
data(expr_data)
data(lesion_data)
data(hg38_gene_annotation)
data(pathways)
# Prepare matched expression and lesion data
alex.data <- alex.prep.lsn.expr(expr_data, lesion_data,
hg38_gene_annotation, min.expr = 5, min.pts.lsn = 5)
# Visualize pathway-level association using JAK pathway as an example
alex.path=alex.pathway(alex.data,
lsn.data = lesion_data,
pathways = pathways,
selected.pathway = "Jak_Pathway")
# Access the ordered data matrix used in the plot
head(alex.path$ordered.path.data)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.