cor.rep.pv {GeneNT}R Documentation

Multivariate correlation estimator from replicated gene microarray data

Description

This function computes p-values of correlation estimation from gene microarray data.

Usage

cor.rep.pv(x, y=NULL, m, G)

Arguments

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).

Details

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.

Value

This function returns a matrix of p-values corresponding to the matrix returned by cor.rep().

Author(s)

Dongxiao Zhu (http://www-personal.umich.edu/~zhud)

References

Zhu, D., Li, Y., and Hero, AO. Estimating gene expression correlation from replicated microarray data - A multivariate approach. {it Submitted}

See Also

cor.rep

Examples

# 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)

[Package GeneNT version 1.3 Index]