Coerce 'rpa' object into an 'ExpressionSet'

Description

An instance of 'rpa' class contains differential gene expression estimates in the variable 'd'. The function 'rpa2eset' coerces this into an ExpressionSet object to allow downstream analysis of the results using standard R/BioC tools for gene expression data.

Usage

rpa2eset(x)

Arguments

Value

An 'ExpressionSet' object.

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

require(RPA)
require(affy)
require(affydata)
data(Dilution)
# Compute RPA for specific probesets
sets = geneNames(Dilution)[1:5]             
rpa.results <- RPA.pointestimate(Dilution,sets)
# Coerce the rpa object into an ExpressionSet
eset = rpa2eset(rpa.results)