BreakMeansChr {WaveCGH}R Documentation

Finds break means in individual chromosomes.

Description

This function applies maximum overlap discrete wavelet transform (MODWT) to find break means inside each of the chromosomes. Double implementation of MODWT confirms the presence of break points in the series.

Usage

BreakMeansChr(x, Chromosome)

Arguments

x CGH array.
Chromosome Sequence of chromosome number, eg. 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3 and so on.

Value

Vector of break means of length equal to that of original series is returned.

Author(s)

A. I. McLeod and M. S. Islam

References

Islam, M. S. (2008). Periodicity, Change Detection and Prediction in Microarrays. Ph.D. Thesis, The University of Western Ontario.

Examples

data(arrayCGH)
ch<- arrayCGH[,2]
# We use first and second chromosomes
BreakMeansChr(arrayCGH[(ch==1|ch==2),1], arrayCGH[(ch==1|ch==2),2])

[Package WaveCGH version 2.01 Index]