PvalBootOneModel {WaveCGH}R Documentation

Finds the p-values for testing means

Description

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.

Usage

PvalBootOneModel(z, Chromosome, ModelSelection = FALSE, nsim = 200)

Arguments

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.

Value

Vector of p-values of length equal to that of z.

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
set.seed(123)
PvalBootOneModel(arrayCGH[(ch==1|ch==2),1], arrayCGH[(ch==1|ch==2),2])

[Package WaveCGH version 2.01 Index]