d.update.fast {RPA}R Documentation

Fast update for probeset-level signal d.

Description

Weighted average over the columns of S, weighted by the inverse variances 1/sigma2.

Usage

d.update.fast(S, sigma2)

Arguments

S Matrix of probe-level observations for a single probeset: samples x probes.
sigma2 A vector. Estimated variance for each measurement (probe).

Value

A vector giving the estimated true signal underlying the observations in S.

Author(s)

Leo Lahti <leo.lahti@tkk.fi>

References

Probabilistic Analysis of Probe Reliability in Differential Gene Expression Studies with Short Oligonucleotide Arrays. Lahti et al., TCBB/IEEE, to appear. See http://www.cis.hut.fi/projects/mi/software/RPA/

Examples

set.seed(24)
vec <- seq(10)
S <- cbind(rnorm(length(vec),sd=1),
          rnorm(length(vec),sd=2),
          rnorm(length(vec),sd=3))
d.update.fast(S,sigma2=seq(3))
[Package RPA version 1.1.2 Index]