atRegions: Summarize a bsseq dataset over defined regions

View source: R/atRegions.R

atRegionsR Documentation

Summarize a bsseq dataset over defined regions

Description

Calls summarizeBsSeqOver to summarize a bsseq object over provided DNA regions. Useful for exploring genomic data using cBioPortal.

Usage

atRegions(bsseq, regions, mappings = NULL, nm = "POETICname", ...)

Arguments

bsseq

A bsseq object

regions

A GRanges or GRangesList of regions

mappings

A mapping table with rownames(mappings) == colnames(bsseq) (DEFAULT: NULL)

nm

Column of the mapping table to map to (DEFAULT: "POETICname")

...

Other arguments to pass to summarizeBsSeqOver

Value

     GRanges with summarized information about the bsseq object
             for the given DNA regions

Examples


  orig_bed <- system.file("extdata", "MCF7_Cunha_chr11p15.bed.gz",
                          package="biscuiteer")
  orig_vcf <- system.file("extdata", "MCF7_Cunha_header_only.vcf.gz",
                          package="biscuiteer")
  bisc <- readBiscuit(BEDfile = orig_bed, VCFfile = orig_vcf,
                      merged = FALSE)

  reg <- GRanges(seqnames = rep("chr11",5),
                 strand = rep("*",5),
                 ranges = IRanges(start = c(0,2.8e6,1.17e7,1.38e7,1.69e7),
                                  end= c(2.8e6,1.17e7,1.38e7,1.69e7,2.2e7))
                )
  regions <- atRegions(bsseq = bisc, regions = reg)


trichelab/biscuiteer documentation built on March 4, 2024, 12:22 a.m.