cor.rep {GeneNT}R Documentation

Estimating gene expression correlation from replicated microarray data

Description

This function estimates the gene expression correlation from replicated microarray data. See exmaple for data format.

Usage

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

Arguments

x x is the gene microarray data in matrix OR data.frame format.
y y Required only both 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 implements the new multivariate correlation estimator from gene microarray data. The estimator takes all observations rather than averaging to gain better accuracy especially for noisy microarray data. See the manuscript for detail.

Value

This function returns a correlation matrix that corresponds to R function cor().

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

Examples

# load GeneNT and GeneTS library
library(GeneNT)
d0 <- rnorm(100)
for(l in 2:10)
d0 <- rbind(d0, rnorm(100))
d0<- t(d0)
M <- cor.rep(d0, m = 4, G= 25)

[Package GeneNT version 1.3 Index]