BackwardElimination: Backward Elimination of segmentation breakpoints

In isglobal-brge/R-GADA: Genome Alteration Detection Algorithm (GADA)

Description

Adjusts the breakpoint sensitivity vs. FDR trade-off after fitting the SBL model

Usage

BackwardElimination(x, T, MinSegLen, saveInfo = TRUE)

Arguments

x	an object of class 'SBL'
T	critical value of the backward elimination algorithm
MinSegLen	minimum number of probes each segment must contain
saveInfo	Indicates whether annotation data is kept. The default is TRUE

Details

After fitting the SBL model, we can use this function to quickly adjust the FDR.

The 'T' argument indicates the critical value of the BE algorithm. That is, the statistical score $t_m$ associated with each breakpoint $m$ remaining in the model has to be higher than $T$. The score $t_m$ can be interpreted as the difference between the sample averages of the probes falling on the left and right segment, divided by a pooled estimation of the standard error. Asympotically, when the number of probes on the right and left segments are very large this score will converge to a standard normal distribution, i.e. N(0,1).

The argument 'MinSegLen' can be used to limit the minimum number of probes each CNA segment must contain. We recommend setting at least 'MinSegLen=3' (default value) to eliminate false detections due to extreme outliers.

Value

An object of class 'BackwardElimination'

'print' returns the number of discontinuities or segments by chromosome

'summary' returns the altered segments (gains and loses) with the begining and the end of each one, as well as the number of probes in each segment. This information is given by chromosome

References

Pique-Regi R, Caceres A, Gonzalez JR. "R-Gada: a package for fast detection and visualization of copy number alterations on multiple samples". BMC Bioinformatics, 2010;11:380.

Pique-Regi R, Monso-Varona J,Ortega A, Seeger RC, Triche TJ, Asgharzadeh S. "Sparse representation and Bayesian detection of the genome copy number alterations from microarray data", Bioinformatics , 2008;24(3):309-18.

See Also

setupGADA, SBL, parSBL, parBE

Examples

 ## Not run: 
  See the vignette
 
## End(Not run) 

isglobal-brge/R-GADA documentation built on May 24, 2019, 5:03 a.m.