plotSegs: plotSegs

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

Description

quick and dirty way of visualizing copy-number segments on a single chromosome

Usage

  plotSegs(rdo,segs,chr)

Arguments

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

Value

a plot with position on the X-axis and segments of CN on the Y axes

See Also

'rd.cnSegments'

Examples

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

  ##extract just the gains and losses
  #  plotSegs(x, segs, "chr22")

chrisamiller/copyCat documentation built on July 20, 2021, 12:59 a.m.