cor.rep.pv {GeneNT} | R Documentation |
This function computes p-values of correlation estimation from gene microarray data.
cor.rep.pv(x, y=NULL, m, G)
x |
x is the gene microarray data in matrix OR data.frame format. |
y |
y optional. Used when x, y are vectors. |
m |
m is the number of replicates for each independent observations. Note: m*G == nrow(dat). |
G |
G is the number of the genes. Note: m*G == nrow(dat). |
This function computes p-values of the new multivariate correlation estimator from gene microarray data using permutation. It is used together with cor.rep(), and cor.rep.bootci() whenever necessary.
This function returns a matrix of p-values corresponding to the matrix returned by cor.rep().
Dongxiao Zhu (http://www-personal.umich.edu/~zhud)
Zhu, D., Li, Y., and Hero, AO. Estimating gene expression correlation from replicated microarray data - A multivariate approach. {it Submitted}
# load GeneNT and GeneTS library library(GeneNT) d0 <- rnorm(100) #sample size is set to 3 for quick demo purposes for(l in 2:3) d0 <- rbind(d0, rnorm(100)) d0<- t(d0) M <- cor.rep(d0, m = 4, G= 25) M.pv <- cor.rep.pv(d0, m = 4, G= 25)