plotSubclonality: Plot clonality and somatic copy number.

In SomatiCA: SomatiCA: identifying, characterizing, and quantifying somatic copy number aberrations from cancer genome sequencing

Description

Plot clonality and somatic copy number for a tumor sample.

Usage

plotSubclonality(segment, dev.new = TRUE, ...)

Arguments

segment	A GRanges object, segments with annotation of copy number and subclonality.
dev.new	Whether a new device will be opened.
...	Parameters in plot().

Value

A plot.

Author(s)

Mengjie Chen

Examples

chr <- c("chr1", "chr1", "chr1", "chr1", "chr2", "chr2", "chr2", "chr2")
start <- c(1, 41, 61 , 71, 1, 51, 71, 91)
end <- c(41, 61 , 71, 91, 51, 71, 91, 101)
medLAF <- c(0.15, 0.27, 0.4, 0.42, 0.4, 0.41, 0.42, 0.39)
medgermlineLAF <- rep(0.43, 8)
ratio <- c(0.5, 0.7, 1.1, 1, 1.1, 1, 1, 0.9)
copynumber <- c(1, 1, 2, 2, 2, 2, 2, 2)
event <- c("LOH", "Loss", "=", "=", "=", "=", "=", "=")
clonality <- c("Clonal", "subclonal_loss", "subclonal_gain", "=", 
            "=", "=", "=", "subclonal_loss")
germlinecopynumber <- c(2, 2, 2, 2, 2, 2, 2, 2)
subclonecopynumber <- c(1, 1, 3, 2, 2, 2, 2, 1)
subpercent <- c(1, 0.4, 0.1, 0, 0, 0, 0, 0.1)
x <- GRanges(seqnames=chr, 
             ranges=IRanges(start=start, end=end),
             medLAF=medLAF,
             medgLAF=medgermlineLAF,
	           ratio=ratio,
             somaCN=copynumber, 
             event=event,
	           clonality=clonality,
	           germCN=germlinecopynumber,
	           subclonalCN=subclonecopynumber,
             subpercent=subpercent)
merged <- MergeSegment(x)
plotSubclonality(merged)

SomatiCA documentation built on Oct. 5, 2016, 4:18 a.m.