Learn R Programming

dinamic (version 1.0)

peeling: Apply the Peeling Procedure at a Given Marker

Description

This function is used internally by DiNAMIC's detailedLook and quickLook functions. Briefly, detailedLook and quickLook assess the statistial significance of the most aberrant gain (loss). Once this is done, the peeling function produces a new matrix of copy number data in which the original aberrant gain (loss) has been nullified. This allows users to assess the statistical significance of subsequent gains (losses) conditional on having found and removed previous gains (losses).

Usage

peeling(x, marker.data, cytoband, k)

Value

A list containing two components: (1) the n by m matrix produced by applying the peeling algorithm to the matrix x at marker k, and (2) the peak interval around marker k, as described in Bioinformatics (2011) 27(5) 678 - 685.

Arguments

x

An n by m numeric matrix containing DNA copy number data from n subjects at m markers.

marker.data

A two-column numeric matrix of marker position data for markers in the autosomes. Column 1 contains the chromosome number for each marker, and column 2 contains the position (in base pairs) for each markers. This is a submatrix of the marker position matrix used by quickLook and detailedLook.

cytoband

A character vector of length m that contains the chromosome arm (p or q) for each marker. This is produced by the makeCytoband function.

k

A positive integer between 1 and m that represents the most aberrant marker.

Author

Vonn Walter, Andrew B. Nobel, Fred A. Wright

Maintainer: <vwalter@email.unc.edu> Vonn Walter

Details

The peeling procedure is detailed in Algorithm 2 of Bioinformatics (2011) 27(5) 678 - 685, but here we provide a brief overview. By construction, marker k represents the most aberrant gain (loss). The peeling procedure rescales all copy number values in x that contribute to making marker k aberrant, so that after applying the peeling procedure marker k is "null." By construction, the rescaling procedure is restricted to entries in x that correspond to markers in the same chromosome arm as k. This allows users to assess the statistical significance of multiple gains (losses) throughout the genome.

References

Walter, V., Nobel, A.B., and Wright, F.A., DiNAMIC: a method to identify recurrent DNA copy number aberrations in tumors, Bioinformatics (2011) 27(5) 678 - 685.