sampleCorrelations.matrix {aroma.light} | R Documentation |
Calculates the correlation for random pairs of observations.
## S3 method for class 'matrix': sampleCorrelations(X, MARGIN=1, pairs=NULL, npairs=max(5000, nrow(X)), ...)
X |
An NxK matrix where N >= 2 and K >= 2. |
MARGIN |
The dimension (1 or 2) in which the observations are.
If MARGIN==1 (==2 ), each row (column) is an observation. |
pairs |
If a Lx2 matrix , the L index pairs for which the
correlations are calculated.
If NULL , pairs of observations are sampled. |
npairs |
The number of correlations to calculate. |
... |
Not used. |
Returns a double
vector
of length npairs
.
Henrik Bengtsson (http://www.braju.com/R/)
[1] A. Ploner, L. Miller, P. Hall, J. Bergh & Y. Pawitan. Correlation test to assess low-level processing of high-density oligonucleotide microarray data. BMC Bioinformatics, 2005, vol 6.
sample
().
# Simulate 20000 genes with 10 observations each X <- matrix(rnorm(n=20000), ncol=10) # Calculate the correlation for 5000 random gene pairs cor <- sampleCorrelations(X, npairs=5000) print(summary(cor))