avearrays: Average Over Replicate Arrays

Description Usage Arguments Details Value Author(s) See Also Examples

View source: R/avearrays.R

Description

Condense a microarray data object so that technical replicate arrays are replaced with (weighted) averages.

Usage

1
2
3
4
5
6
## Default S3 method:
avearrays(x, ID=colnames(x), weights=NULL)
## S3 method for class 'MAList'
avearrays(x, ID=colnames(x), weights=x$weights)
## S3 method for class 'EList'
avearrays(x, ID=colnames(x), weights=x$weights)

Arguments

x

a matrix-like object, usually a matrix, MAList or EList object.

ID

sample identifier.

weights

numeric matrix of non-negative weights

Details

A new data object is computed in which technical replicate arrays are replaced by their (weighted) averages.

For an MAList object, the components M and A are both averaged in this way, as are weights and any matrices found in object$other.

EList objects are similar, except that the E component is averaged instead of M and A.

If x is of mode "character", then the replicate values are assumed to be equal and the first is taken as the average.

Value

A data object of the same class as x with a column for each unique value of ID.

Author(s)

Gordon Smyth

See Also

avereps.

02.Classes gives an overview of data classes used in LIMMA.

Examples

1
2
3
x <- matrix(rnorm(8*3),8,3)
colnames(x) <- c("a","a","b")
avearrays(x)

Example output

               a          b
[1,]  0.06818442  0.1725134
[2,]  0.39106518  1.6426725
[3,] -0.74722260  0.3671410
[4,]  0.71964981 -0.3534182
[5,]  0.63591556  0.1816293
[6,] -0.44054880 -0.2362331
[7,] -0.11815876  1.1829985
[8,] -0.24357897 -1.0567411

limma documentation built on Nov. 8, 2020, 8:28 p.m.