CGHcall (version 2.34.0)

segmentData: Breakpoint detection for arrayCGH data.

Description

A wrapper function to run existing breakpoint detection algorithms on arrayCGH data. Currently only DNAcopy is implemented.

Usage

segmentData(input, clen=10, relSDlong=3, method = "DNAcopy", ...)

Arguments

input

Object of class cghRaw.

clen

Boundary for short vs long segments, in number of features

relSDlong

Relative undo sd for long segments. See details.

method

The method to be used for breakpoint detection. Currently only DNAcopy is supported, which will run the segment function.

...

Arguments for segment.

Value

This function returns a dataframe in the same format as the input with segmented arrayCGH data.

Details

See segment for details on the algorithm. About clen and relSDlong: these are only relevant when segment option undo.splits=sdundo is set, in combination with segment option undo.SD. relSDlong provides the undo sd for long segments, which equals undo.SD/relSDlong. undo.SD is then used for short segments. In the example below, short segments are considered to contain less or equal to clen=10 features. The example below undoes splits for two consecutive short segments if these are less than undo.SD=3 sd apart, while it undoes splits for two long segments if these are less than undo.SD/relSDlong=3/3=1 sd apart. If, for two consecutive segements, one is short and one is long, splits are undone in the same way as for two short segments.

References

Venkatraman, A.S., Olshen, A.B. (2007). A faster circulary binary segmentation algorithm for the analysis of array CGH data. Bioinformatics, 23, 657-663.

Examples

  data(WiltingNorm)
  ## Not run: segmented.data <- segmentData(WiltingNorm, alpha=0.02,clen=10,relSDlong=3,undo.SD=3,undo.splits="sdundo")