break.density: Breakpoint density map

Description Usage Arguments Value Examples

View source: R/breakpoint.density.r

Description

Generating a genomic map based on a defined bin size and sliding window and counts the number of breakpoints mapped onto each bin. This function is used internally by svpluscnv::shattered.regions and svpluscnv::shattered.regions.cnv

Usage

1
2
3
4
5
6
7
8
break.density(
  brk,
  chr.lim = NULL,
  genome.v = "hg19",
  window.size = 10,
  slide.size = 2,
  verbose = TRUE
)

Arguments

brk

(breaks) An instance of the class 'breaks' obtained from CNV segmentation data (svpluscnv::cnv.breaks) or Structural Variant calls (svpluscnv::svc.breaks).

chr.lim

(data.frame) 3 column table (chrom, begin, end) indicating the chromosome most distal coordinates with coverage. Also returned by the function svpluscnv::chromosome.limit.coords.

genome.v

(hg19 or hg38) reference genome version to draw chromosome limits and centromeres

window.size

(numeric) size in megabases of the genmome bin onto which breakpoints will be mapped

slide.size

(numeric) size in megabases of the sliding genomic window; if slide.size < window.size the genomic bins will overlap

verbose

(logical) whether to return internal messages

Value

a matrix of samples (rows) and genomic bins (cols) qith the number of breakpoints mapped in heach cell

Examples

1
2
3
4
5
6
7
# initialize CNV data
cnv <- validate.cnv(segdat_lung_ccle)

# obtain CNV breakpoints
brk <- cnv.breaks(cnv)

break.density(brk)

ccbiolab/svpluscnv documentation built on Sept. 9, 2020, 4:52 a.m.