binCoverage: Bin CpG or CpH coverage to simplify and improve CNA...

In biscuiteer: Convenience Functions for Biscuit

Description

Example usage for E-M

Usage

binCoverage(bsseq, bins, which = NULL, QDNAseq = TRUE, readLen = 100)

Arguments

bsseq	A bsseq object - supplied to getCoverage()
bins	Bins to summarize over - from tileGenome or QDNAseq.xxYY
which	Limit to specific regions? - functions as an import() (DEFAULT: NULL)
QDNAseq	Return a QDNAseqReadCounts? - if FALSE, returns a GRanges (DEFAULT: TRUE)
readLen	Correction factor for coverage - read length in bp (DEFAULT: 100)

Details

NOTE: As of early Sept 2019, QDNAseq did not have hg38 capabilities. If you desire to use the hg38 genome, biscuiteer suggests you use a GRanges object to define your bins.

Value

    Binned read counts

Examples

  bins <- GRanges(seqnames = rep("chr11",10),
                  strand = rep("*",10),
                  ranges = IRanges(start=100000*0:9, width=100000)
                 )

  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))
                 )

  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)

  bc <- binCoverage(bsseq = bisc, bins = bins, which = reg, QDNAseq = FALSE)

biscuiteer documentation built on Nov. 8, 2020, 8:28 p.m.