doCBS: Run CBS segmentation

In jointseg: Joint segmentation of multivariate (copy number) signals

Description

This function is a wrapper for convenient use of the CBS segmentation method by PSSeg. It applies the segment function and reshapes the results

Usage

doCBS(y, ..., verbose = FALSE)

Arguments

y	A numeric vector, the signal to be segmented
...	Arguments to be passed to segment
verbose	A logical value: should extra information be output ? Defaults to FALSE.

Value

A list with element:

bkp	breakpoint positions

Author(s)

Morgane Pierre-Jean and Pierre Neuvial

See Also

segment

Examples

## load known real copy number regions
affyDat <- loadCnRegionData(dataSet="GSE29172", tumorFraction=1)

## generate a synthetic CN profile
K <- 10
len <- 1e4
sim <- getCopyNumberDataByResampling(len, K, minLength=100, regData=affyDat)
datS <- sim$profile

## run CBS segmentation
res <- doCBS(datS[["c"]])
getTpFp(res$bkp, sim$bkp, tol=5, relax = -1)   ## true and false positives
plotSeg(datS, breakpoints=list(sim$bkp, res$bkp))

jointseg documentation built on May 2, 2019, 5:20 p.m.