trimSegmentEnds: trimSegmentEnds

In chrisamiller/readDepth: Find regions of genomic copy number loss and gain from short reads

Description

trim the ends of segments such that they're called normal (CN = 2.0) until the first bin that isn't NA

useful for chromosomes like 22, where there are tens of megabases that are unmappable at beginning. Without this, those regions would be called with the same CN as the first few mappable bins.

Usage

  trimSegmentEnds(segs, rdo)

Arguments

segs	a data frame that is the result of calling rd.cnSegments
rdo	a readDepth object created with new("rdObject")

Value

the segment data frame adjusted to have the ends trimmed

Examples

  ##get regions of CN gain and loss
  #  segs <- rd.cnSegments(x)

  ##trim the segment ends 
  #  trimSegmentEnds(segs, rdo)

chrisamiller/readDepth documentation built on June 16, 2021, 2:05 p.m.