segmentCGH: Genomic Profile Segmentation

In fredcommo/rCGH: Comprehensive Pipeline for Analyzing and Visualizing Array-Based CGH Data

Description

A function for performing the Log2Ratio segmentation on an object of class "rCGH". See the details section below.

Usage

## S4 method for signature 'rCGH'
segmentCGH(object, Smooth=TRUE, UndoSD = NULL,
minLen = 10, nCores=NULL, verbose = TRUE)

Arguments

object	: An object of class "rCGH"
Smooth	: logical. Should the LRR be smoothed before being segmented. See DNAcopy for details.
UndoSD	: numeric. When not specified (default is NULL), the UndoSD value is estimated from the Log2Ratios. See DNAcopy for details.
minLen	: numeric. The minimal length for a segment, in Kb. Shorter segments will be merged to the closest adjacent one. Default value is 10(Kb).
nCores	: numeric. The number of cores to use in order to speed up the computation. When NULL (default), half of the available cores are used. See mclapply.
verbose	: logical. if TRUE (default), progress is printed.

Details

This function is a wrapper for the DNAcopy, CNA and segment functions, which allows parallelization and data-driven parameterization.
In addition to the usual DNAcopy output, the segmentation table contains the probes Log2Ratio standard deviation for each segment, as well as there length, in Kb.

Value

An object of class "rCGH"

Author(s)

Frederic Commo

References

Venkatraman ES1, Olshen AB. A faster circular binary segmentation algorithm for the analysis of array CGH data.Bioinformatics. 2007 Mar 15;23(6):657-63.

See Also

CNA, segment, mclapply

Examples

filePath <- system.file("extdata", "Affy_cytoScan.cyhd.CN5.CNCHP.txt.bz2",
    package = "rCGH")
cgh <- readAffyCytoScan(filePath, sampleName = "AffyScHD")
cgh <- adjustSignal(cgh, nCores=1)
cgh <- segmentCGH(cgh, nCores=1)
st <- getSegTable(cgh)
head(st)

fredcommo/rCGH documentation built on May 16, 2019, 2:41 p.m.