Summarize motifbreakR + echolocatoR results

Summarize motifbreakR + echolocatoR results after they have been merged together with MOTIFBREAKR_filter.

MOTIFBREAKR_summarize(mb_merge)

Arguments

mb_merge

Output of MOTIFBREAKR_filter.

Value

Per-locus summary in data.frame format.

See also

Other motifbreakR: MOTIFBREAKR_filter_by_metadata(), MOTIFBREAKR()

Examples

merged_DT <- echodata::get_Nalls2019_merged()
mb_res <- MOTIFBREAKR(rsid_list = c("rs11175620"),
                      # limit the number of datasets tests 
                      # for demonstration purposes only
                      pwmList_max = 4,
                      calculate_pvals = FALSE)
#> genome_build set to hg19 by default.
#> Using genome_build hg19
#> + MOTIFBREAKR:: Using pre-existing tmp file.
mb_merge <- MOTIFBREAKR_filter(mb_res = mb_res,
                               merged_DT = merged_DT,
                               pvalue_threshold = NULL)
#> Creating 'risk_' columns
#> + MOTIFBREAKR 1 SNPs in results.
#> + MOTIFBREAKR 1 SNPs in results @ `snp_filter='Support>0'`
#> + MOTIFBREAKR 0 SNPs in results @ `qvalue_threshold=0.05`
summary_ls <- MOTIFBREAKR_summarize(mb_merge = mb_merge)