Subset RD EWCE results data by cell type, fold change and q value.
subset_results(
filters = list(cl_name = NULL),
results = load_example_results(),
q_threshold = 5e-04,
effect_threshold = 0.1,
effect_var = "fold_change",
verbose = TRUE
)
A named list, where each element in the list is the name of a column in the data, and the vector within each element represents the values to include in the final data.
The cell type-phenotype enrichment results generated by gen_results and merged together with merge_results
The q value threshold to subset the results
by.
The minimum fold change in specific expression
to subset the results
by.
Name of the effect size column in the results
.
Print messages.
A data frame of the selected subset of RD EWCE results with HPO ID column added.
signif_cell_data <- subset_results(filters=list(CellType = "Cardiomyocytes"))
#> Subsetting results by q_threshold and effect.
#> Selected CellType :
#> - Cardiomyocytes
#> Filtered 'CellType' : 25,134 / 25,394 rows dropped.
#> 260 associations remain after filtering.