gcrma {gcrma} | R Documentation |
This function converts an AffyBatch
into an ExpressionSet
using the robust multi-array average (RMA) expression measure with help of probe sequence.
gcrma(object,affinity.info=NULL, affinity.source=c("reference","local"),NCprobe=NULL, type=c("fullmodel","affinities","mm","constant"), k=6*fast+0.5*(1-fast),stretch=1.15*fast+1*(1-fast),correction=1, GSB.adjust=TRUE, rho=.7,optical.correct=TRUE,verbose=TRUE,fast=TRUE, subset=NULL,normalize=TRUE,...)
object |
an |
affinity.info |
|
affinity.source |
|
NCprobe |
Index of negative control probes. When set as
|
type |
"fullmodel" for sequence and MM model. "affinities" for sequence information only. "mm" for using MM without sequence information. |
k |
A tuning factor. |
stretch |
. |
correction |
. |
GSB.adjust |
Logical value. If |
rho |
correlation coefficient of log background intensity in a pair of pm/mm probes. Default=.7 |
optical.correct |
Logical value. If |
verbose |
Logical value. If |
fast |
Logical value. If |
subset |
a character vector with the the names of the probesets to be used in expression calculation. |
normalize |
logical value. If 'TRUE' normalize data using quantile normalization. |
... |
further arguments to be passed (not currently implemented - stub for future use). |
Note that this expression measure is given to you in log base 2 scale. This differs from most of the other expression measure methods.
The tuning factor k
will have different meanings if one uses
the fast (add-hoc) algorithm or the empirical Bayes approach. See Wu
et al. (2003)
An ExpressionSet
.
Rafeal Irizarry
if(require(affydata) & require(hgu95av2probe) & require(hgu95av2cdf)){ data(Dilution) ai <- compute.affinities(cdfName(Dilution)) Dil.expr<-gcrma(Dilution,affinity.info=ai,type="affinities") }