Robust smooth segmentation approach for array CGH data analysis.
BACKGROUND
Array comparative genomic hybridization (aCGH) provides a genome-wide technique to screen for copy number alteration. The existing segmentation approaches for analyzing aCGH data are based on modeling data as a series of discrete segments with unknown boundaries and unknown heights. Although the biological process of copy number alteration is discrete, in reality a variety of biological and experimental factors can cause the signal to deviate from a stepwise function. To take this into account, we propose a smooth segmentation (smoothseg) approach.


METHODS
To achieve a robust segmentation, we use a doubly heavy-tailed random-effect model. The first heavy-tailed structure on the errors deals with outliers in the observations, and the second deals with possible jumps in the underlying pattern associated with different segments. We develop a fast and reliable computational procedure based on the iterative weighted least-squares algorithm with band-limited matrix inversion.


RESULTS
Using simulated and real data sets, we demonstrate how smoothseg can aid in identification of regions with genomic alteration and in classification of samples. For the real data sets, smoothseg leads to smaller false discovery rate and classification error rate than the circular binary segmentation (CBS) algorithm. In a realistic simulation setting, smoothseg is better than wavelet smoothing and CBS in identification of regions with genomic alterations and better than CBS in classification of samples. For comparative analyses, we demonstrate that segmenting the t-statistics performs better than segmenting the data.


BACKGROUND
The R package smoothseg to perform smooth segmentation is available from http://www.meb.ki.se/~yudpaw.
Huang
Jian
J
Statistical Laboratory, Department of Statistics, University College Cork, Ireland.
Gusnanto
Arief
A
O'Sullivan
Kathleen
K
Staaf
Johan
J
Borg
Ake
A
Pawitan
Yudi
Y
D000465:Algorithms; D002874:Chromosome Mapping / Q000379:methods*; D018628:Gene Dosage / Q000235:genetics*; D017404:In Situ Hybridization, Fluorescence / Q000379:methods*; D020411:Oligonucleotide Array Sequence Analysis / Q000379:methods*; D016415:Sequence Alignment / Q000379:methods*; D017422:Sequence Analysis, DNA / Q000379:methods*