BreakMeansChr {WaveCGH} | R Documentation |
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.
BreakMeansChr(x, Chromosome)
x |
CGH array. |
Chromosome |
Sequence of chromosome number, eg. 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3 and so on. |
Vector of break means of length equal to that of original series is returned.
A. I. McLeod and M. S. Islam
Islam, M. S. (2008). Periodicity, Change Detection and Prediction in Microarrays. Ph.D. Thesis, The University of Western Ontario.
data(arrayCGH) ch<- arrayCGH[,2] # We use first and second chromosomes BreakMeansChr(arrayCGH[(ch==1|ch==2),1], arrayCGH[(ch==1|ch==2),2])