View source: R/alex.prep.lsn.expr.R
alex.prep.lsn.expr | R Documentation |
Prepares matched lesion and expression data matrices for use with the KW.hit.express
function, which performs Kruskal Wallis tests to assess associations between lesion groups and gene expression levels.
alex.prep.lsn.expr(
expr.mtx,
lsn.data,
gene.annotation,
min.expr = NULL,
min.pts.lsn = NULL
)
expr.mtx |
A data frame or matrix of normalized, log2-transformed gene expression values with genes in rows and subjects in columns. The first column must be named |
lsn.data |
A data frame of lesion data in GRIN-compatible format. Must contain five columns:
|
gene.annotation |
A gene annotation data frame, either user-provided or retrieved via |
min.expr |
Minimum total expression level required for a gene to be retained (i.e., sum of expression values across all subjects). Useful to filter out genes with very low expression. |
min.pts.lsn |
Minimum number of subjects required to have a lesion in a gene for that gene to be retained for the KW test. |
The function uses prep.lsn.type.matrix()
internally to create a lesion matrix where each gene is represented by one row and all lesion types are included. It filters genes to retain only those with both sufficient expression and lesion data. The final expression and lesion matrices are matched by gene and patient IDs, with rows ordered by Ensembl gene ID and columns by patient ID.
A list with the following components:
alex.expr |
A matrix of gene expression data with Ensembl gene IDs as row names and patient IDs as column names. |
alex.lsn |
A matrix of lesion data for the same genes and patients as in |
alex.row.mtch |
A data frame with two columns showing the matched Ensembl gene IDs from the expression and lesion matrices. |
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.
KW.hit.express
data(expr_data)
data(lesion_data)
data(hg38_gene_annotation)
# Prepare matched lesion and expression data
alex.data <- alex.prep.lsn.expr(expr_data, lesion_data,
hg38_gene_annotation, min.expr = 1,
min.pts.lsn = 5)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.