View source: R/alex.boxplots.R
alex.boxplots | R Documentation |
Generates box plots of gene expression levels stratified by lesion groups for a subset of genes selected based on a user-specified q-value threshold from Kruskal Wallis test results.
alex.boxplots(out.dir, alex.data, alex.kw.results, q, gene.annotation)
out.dir |
Path to the directory where the resulting PDF file containing the box plots will be saved. |
alex.data |
A list returned by the |
alex.kw.results |
Output of the |
q |
Q-value threshold. Only genes with q-values below this threshold will be included in the output plots. |
gene.annotation |
A |
A PDF file saved in out.dir
, containing one box plot per page for each selected gene, showing gene expression distribution across lesion groups.
Abdelrahman Elsayed abdelrahman.elsayed@stjude.org, Stanley Pounds stanley.pounds@stjude.org
Cao, X., Elsayed, A. H., & Pounds, S. B. (2023). Statistical Methods Inspired by Challenges in Pediatric Cancer Multi-omics.
alex.prep.lsn.expr
, KW.hit.express
data(expr_data)
data(lesion_data)
data(hg38_gene_annotation)
# Prepare expression and lesion data
alex.data <- alex.prep.lsn.expr(expr_data, lesion_data,
hg38_gene_annotation, min.expr = 5, min.pts.lsn = 5)
# Run Kruskal Wallis test
alex.kw.results <- KW.hit.express(alex.data, hg38_gene_annotation, min.grp.size = 5)
# Generate box plots for significant genes
dir.create(resultsFolder <- file.path(tempdir(), "temp.out"))
alex.boxplots(out.dir = resultsFolder,
alex.data = alex.data,
alex.kw.results = alex.kw.results,
q = 1e-15,
gene.annotation = hg38_gene_annotation)
unlink(resultsFolder, recursive = TRUE)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.