enrichment_depletion_test {MutationalPatterns} | R Documentation |
This function aggregates mutations per group (optional) and performs an enrichment depletion test.
enrichment_depletion_test(x, by = c())
x |
data.frame result from genomic_distribution() |
by |
Optional grouping variable, e.g. tissue type |
data.frame with the observed and expected number of mutations per genomic region per group (by) or sample
genomic_distribution
,
plot_enrichment_depletion
## See the 'genomic_distribution()' example for how we obtained the ## following data: distr <- readRDS(system.file("states/distr_data.rds", package="MutationalPatterns")) tissue <- c(rep("colon", 3), rep("intestine", 3), rep("liver", 3)) ## Perform the enrichment/depletion test by tissue type. distr_test <- enrichment_depletion_test(distr, by = tissue) ## Or without specifying the 'by' parameter. distr_test2 <- enrichment_depletion_test(distr)