PvalBootMeanChr {WaveCGH} | R Documentation |
This function applies maximum overlap discrete wavelet transform (MODWT) to find break means inside each of the chromosomes. Then applies bootstrap method to find p-values for testing the means.
PvalBootMeanChr(x, Chromosome, ModelSelection = FALSE, nsim = 200)
x |
Array-CGH data. |
Chromosome |
Sequence of chromosomes, eg. 1, 1, 1, 2, 2, 2, 2, 3, 3, 3, 3 and so on. |
ModelSelection |
Whether to estimate order (p) of an AR(p) process. |
nsim |
Number of bootstrap simulations. |
Vector of p-values of length equal to that of original series.
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 set.seed(123) PvalBootMeanChr(arrayCGH[(ch==1|ch==2),1], arrayCGH[(ch==1|ch==2),2])