PvalBootOneModel {WaveCGH} | R Documentation |
This function calculates p-values for testing means when we consider the residuals obtained from all chromosomes to produce a single series. The break means are obtained using wavelet transformation in each individual chromosomes. Then the p-values for testing mean are obtained using bootstrap procedure.
PvalBootOneModel(z, Chromosome, ModelSelection = FALSE, nsim = 200)
z |
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 z.
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) PvalBootOneModel(arrayCGH[(ch==1|ch==2),1], arrayCGH[(ch==1|ch==2),2])