Visualize hyp/multihyp objects as an enrichment map
hyp_emap( hyp_obj, similarity_metric = c("jaccard_similarity", "overlap_similarity"), similarity_cutoff = 0.2, pval = 1, fdr = 1, val = c("fdr", "pval"), top = NULL, title = "" )
hyp_obj | A hyp or multihyp object |
---|---|
similarity_metric | Metric to calculate geneset similarity |
similarity_cutoff | Geneset similarity cutoff |
pval | Filter results to be less than pval cutoff |
fdr | Filter results to be less than fdr cutoff |
val | Choose significance value shown above nodes e.g. c("fdr", "pval") |
top | Limit number of pathways shown |
title | Plot title |
A visNetwork object
genesets <- msigdb_gsets("Homo sapiens", "C2", "CP:KEGG") signature <- c("IDH3B","DLST","PCK2","CS","PDHB","PCK1","PDHA1","LOC642502", "PDHA2","LOC283398","FH","SDHD","OGDH","SDHB","IDH3A","SDHC", "IDH2","IDH1","OGDHL","PC","SDHA","SUCLG1","SUCLA2","SUCLG2") hyp_obj <- hypeR(signature, genesets, background=2522) hyp_emap(hyp_obj, top=30, val="fdr")