R/get_geneSetRatio.R
calculate_geneSetRatio.RdGet a percentage of genesets the specific gene is included
calculate_geneSetRatio(enrichment_results, gene_overview_df)numerical vector of gene set ratios
calculate_geneSetRatio(
list(
A = get(load(system.file("extdata", "example_enrichment.rda", package = "goatea"))),
B = get(load(system.file("extdata", "example_enrichment.rda", package = "goatea")))
),
get(load(system.file("extdata", "example_genes_overview.rda", package = "goatea"))))
#> # A tibble: 100 × 13
#> gene symbol A_efsi A_pval B_efsi B_pval A_perc B_perc genelist_overlap
#> <int> <chr> <dbl> <dbl> <dbl> <dbl> <dbl> <dbl> <chr>
#> 1 11023 gene_45 5 0.123 5 0.123 3100 3100 ""
#> 2 12763 gene_12 5 0.435 5 0.435 3100 3100 ""
#> 3 16847 gene_34 5 0.435 5 0.435 3100 3100 ""
#> 4 12069 gene_83 -4.9 0.0148 -4.9 0.0148 -96.7 -96.7 "A+B"
#> 5 17454 gene_14 4.8 0.0151 4.8 0.0151 2686. 2686. "A+B"
#> 6 12308 gene_96 4.8 0.375 4.8 0.375 2686. 2686. ""
#> 7 19013 gene_57 -4.8 0.667 -4.8 0.667 -96.4 -96.4 ""
#> 8 11532 gene_92 4.7 0.676 4.7 0.676 2499. 2499. ""
#> 9 13749 gene_31 -4.6 0.00597 -4.6 0.00597 -95.9 -95.9 "A+B"
#> 10 10264 gene_59 -4.6 0.0410 -4.6 0.0410 -95.9 -95.9 "A+B"
#> # ℹ 90 more rows
#> # ℹ 4 more variables: A_geneSetRatio <dbl>, B_geneSetRatio <dbl>,
#> # A_SS_geneSetRatio <dbl>, B_SS_geneSetRatio <dbl>