This function sorts the genes in an expression matrix based on a summary that is performed on a population of cells. For each gene, the population is selected as the cells expressing that gene above a given threshold. The summary is calculated on the metadata information of the obtained population.
Usage
sort_genes_by_metadata(
expression_matrix,
metadata_info,
summary_function = mean,
thresh_percentile = 0.25,
thresh_value = 0,
decreasing = FALSE
)Arguments
- expression_matrix
The gene by cell expression matrix to be used.
- metadata_info
The metadata information to be used for the sorting. The metadata should be a numeric vector.
- summary_function
A function that takes as argument a numeric vector and summarises it into a single value. Defaults to
mean.- thresh_percentile
The percentile to be used as threshold for the expression values to select the cells associated with each gene. If set to 0, the threshold will be the value defined in
thresh_value. Defaults to 0.25.- thresh_value
The value to be used as threshold for the expression values to select the cells associated with each gene. Defaults to 0.
- decreasing
Logical indicating if the genes should be sorted in decreasing or increasing order. Defaults to FALSE.