[BioC] Limma: Question about extracting 2-channel data

Boris Umylny umylny at apbri.org
Fri Jan 30 02:50:53 CET 2009


Dear Alex,

Thank you very much for your kind reply.

The RG.MA function is exactly what we were looking for.  I am sorry to trouble 
you with such a trivial question.

Thank you again.


Sincerely,


Boris Umylny


On Thursday 29 January 2009 05:57:03 pm Axel.Klenk at actelion.com wrote:
> Dear Boris,
>
> see ?RG.MA from limma or the RG.MA code for computing normalized RG from
> MA, i.e.
>
> > RG.MA
>
> function (object)
> {
>     object$R <- 2^(object$A + object$M/2)
>     object$G <- 2^(object$A - object$M/2)
>     object$M <- NULL
>     object$A <- NULL
>     new("RGList", unclass(object))
> }
> <environment: namespace:limma>
>
> Cheers,
>
>  - axel
>
>
> Axel Klenk
> Research Informatician
> Actelion Pharmaceuticals Ltd / Gewerbestrasse 16 / CH-4123 Allschwil /
> Switzerland
> phone +41 61 565 6367 / fax +41 61 565 6470 / axel.klenk at actelion.com /
> www.actelion.com
>
>
>
>
>              Boris Umylny
>              <umylny at apbri.org
>
>              >                                                          To
>
>              Sent by:                  bioconductor at stat.math.ethz.ch
>              bioconductor-boun                                          cc
>              ces at stat.math.eth
>              z.ch                                                  Subject
>                                        [BioC] Limma: Question about
>                                        extracting 2-channel data
>              29.01.2009 09:33
>
>
>
>
>
>
>
>
>
> Thank you in advance for your help.
>
> We are trying to use limma package to extract and normalize Agilent
> 2-channel data.  Our objective is to use a single 2-channel chip as 2
> 1-channel chips.
>
> We attempted the following:
> > library(limma)
> > targets <- readTargets("targets_file", row.names = "Name")
> > RG <- read.maimages(...)
> > RG.c <- backgroundCorrect(...)
>
> At this point we have an RG and RG.c objects that contains information for
> both channels.
>
> > MA <- normalizeWithinArrays(...)
>
> At this point, the two channel data has been replaced by a single value.
>
> We tried to use targetsA2C to convert to 1-channel design:
> > targets.1 <- targetsA2C(targets)
> > RG <- read.maimages(targets.1 ...)
>
> However, that simply duplicated Red and Green values.  For 3 chips we had a
>
> set of 3 red and 3 green duplicates.  So, once we get to:
> > MA <- normalizeWithinArrays(...)
>
> We have 6 values, but these 6 values are 3 sets of duplicates.
>
> How would we go about getting a set of normalized values from each channel
> of a 2-channel chip using limma?
>
>
> Sincerely,
>
>
> Boris Umylny
>
>
>
>
>
>              [[alternative HTML version deleted]]
>
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives:
> http://news.gmane.org/gmane.science.biology.informatics.conductor
>
>
>
>
> The information of this email and in any file transmitted with it is
> strictly confidential and may be legally privileged. It is intended solely
> for the addressee. If you are not the intended recipient, any copying,
> distribution or any other use of this email is prohibited and may be
> unlawful. In such case, you should please notify the sender immediately and
> destroy this email. The content of this email is not legally binding unless
> confirmed by letter. Any views expressed in this message are those of the
> individual sender, except where the message states otherwise and the sender
> is authorised to state them to be the views of the sender's company. For
> further information about Actelion please see our website at
> http://www.actelion.com



More information about the Bioconductor mailing list