Interval statistics can be used in combination with
dplyr::group_by() and dplyr::do() to calculate
statistics for subsets of data. See vignette('interval-stats') for
examples.
# NOT RUN {genome <- read_genome(valr_example('hg19.chrom.sizes.gz'))
x <- bed_random(genome, seed = 1010486)
y <- bed_random(genome, seed = 9203911)
bed_fisher(x, y, genome)
# }