threestep {affyPLM} | R Documentation |
This function converts an AffyBatch
into an
exprSet
using
a three step expression measure.
threestep(object,subset=NULL, verbose=TRUE, normalize=TRUE,background=TRUE,background.method="RMA.2",normalize.method="quantile",summary.method="median.polish",background.param = list(),normalize.param=list(),summary.param=list())
object |
an AffyBatch |
subset |
an indexing vector with the the probe sets to be converted to expression. (not currently implemented) |
verbose |
logical value. If TRUE it writes out some messages. (Curently ignored) |
normalize |
logical value. If TRUE normalize data using
quantile normalization |
background |
logical value. If TRUE background correct
using RMA background correction |
background.method |
name of background method to use. |
normalize.method |
name of normalization method to use. |
summary.method |
name of summary method to use. |
background.param |
list of parameters for background correction methods |
normalize.param |
list of parameters for normalization methods |
summary.param |
list of parameters for summary methods |
This function computes the expression measure using threestep methods.
An exprSet
Ben Bolstad bolstad@stat.berkeley.edu
Irizarry et al. Biostatistics 2003
data(affybatch.example) # should be equivalent to rma() eset <- threestep(affybatch.example) # Using Tukey Biweight summarization eset <- threestep(affybatch.example,summary.method="tukey.biweight") # Using Average Log2 summarization eset <- threestep(affybatch.example,summary.method="average.log") # Using IdealMismatch background and Tukey Biweight and no normalization. eset <- threestep(affybatch.example,normalize=FALSE,background.method="IdealMM",summary.method="tukey.biweight") # Using average.log summarization and no background or normalization. eset <- threestep(affybatch.example,background=FALSE,normalize=FALSE,background.method="IdealMM",summary.method="tukey.biweight") # Use threestep methodology with the rlm model fit eset <- threestep(affybatch.example,summary.method="rlm") # Use threestep methodology with the log of the average eset <- threestep(affybatch.example,summary.method="log.average") # Use threestep methodology with log 2nd largest method eset <- threestep(affybatch.example,summary.method="log.2nd.largest") eset <- threestep(affybatch.example,background.method="LESN2")