[BioC] affy 1.1 - rma function
Pawel Herzyk
ph53d@udcf.gla.ac.uk
Thu, 21 Nov 2002 10:37:54 -0000
This is a multi-part message in MIME format.
------=_NextPart_000_008D_01C2914A.0CAAE910
Content-Type: text/plain;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
Hi there,
I have a questions related to what the "rma" function actually does.
Let assume that I have converted the probe level data stored in the =
"Data" object to the gene expression values using the following
> eset <- rma(Data)
If I wanted to get the same effect but wanted to do it step by step so =
that I could save the intermediate AffyBatch objects, then according to =
my understanding the closest to the above would be the following:
> bg.corrected.Data <- bg.correct.rma(Data)
> norm.Data <- normalize(bg.corrected.Data)
> eset <- =
express(norm.Data,bg.correct=3DF,normalize.method=3DF,summary.stat=3Dmedi=
anpolish)
Obviously I won't get exactly the same numbers because "rma" is written =
in C, so the codes are different, but am I correct in thinking that the =
above three-step procedure does virtually the same as "rma"?
Pawel
------=_NextPart_000_008D_01C2914A.0CAAE910
Content-Type: text/html;
charset="iso-8859-1"
Content-Transfer-Encoding: quoted-printable
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN">
<HTML><HEAD>
<META content=3D"text/html; charset=3Diso-8859-1" =
http-equiv=3DContent-Type>
<META content=3D"MSHTML 5.00.3315.2870" name=3DGENERATOR>
<STYLE></STYLE>
</HEAD>
<BODY bgColor=3D#c8e0d8>
<DIV><FONT face=3DArial size=3D2>Hi there,</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>I have a questions related to what the =
"rma"=20
function actually does.</FONT></DIV>
<DIV><FONT face=3DArial size=3D2></FONT> </DIV>
<DIV><FONT face=3DArial size=3D2>Let assume that I have converted the =
probe level=20
data stored in the "Data" object to the gene expression values using the =
following</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>> eset <- rma(Data)</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>If I wanted to get the same effect =
but wanted=20
to do it step by step so that I could save the intermediate AffyBatch =
objects,=20
then according to my understanding the closest to the above would be the =
following:</FONT></DIV>
<DIV> </DIV>
<DIV><FONT face=3DArial size=3D2>
<DIV>> bg.corrected.Data <- bg.correct.rma(Data)</DIV>
<DIV>> norm.Data <- normalize(bg.corrected.Data)</DIV>
<DIV>> eset <-=20
express(norm.Data,bg.correct=3DF,normalize.method=3DF,summary.stat=3Dmedi=
anpolish)</DIV>
<DIV> </DIV>
<DIV>Obviously I won't get exactly the same numbers because "rma" is =
written in=20
C, so the codes are different, but am I correct in thinking that the =
above=20
three-step procedure does virtually the same as "rma"?</DIV>
<DIV> </DIV>
<DIV>Pawel</DIV></FONT></DIV>
<DIV> </DIV></BODY></HTML>
------=_NextPart_000_008D_01C2914A.0CAAE910--