[BioC] Are individual channels' intensities comparable across cDNA arrays?
sdavis2 at mail.nih.gov
Thu Oct 16 01:39:15 CEST 2008
On Wed, Oct 15, 2008 at 7:14 PM, Seungwoo Hwang <swhwang10 at yahoo.com> wrote:
> Dear all,
> I received the following request. I'm not so sure about it and I would like to receive your expert opinion on this.
> The platform is cDNA array with two experimental series.
> - On series 1 chips, tumor sample from a patient (Cy5) and normal sample from pooled healthy individuals (Cy3) were co-hybridized.
> - On series 2 chips, cirrhosis sample from that patient (Cy5) and and normal sample from pooled healthy individuals (Cy3) were co-hybridized.
> - This hybridization scheme was utilized because they wanted to make the following three comparisons: tumor vs normal, cirrhosis vs normal, and tumor vs cirrhosis.
> I normalize by print-tip loess for within array and scale method for between arrays using limma. After that, I usually work on normalized M values for differential expression analysis.
> The request that I received, however, was to obtain tumor sample intensities (Cy5 on series 1 chips) and cirrhosis sample intensities (Cy5 on series 2 chips) so that they can do some comparative analyses with their protein data.
> I know that I am able to obtain those values using RG.MA(). What I am not sure is, will those two individual channels' intensity values become comparable with each other? By comparable, I mean, when a gene's M is 1 in series 1 chips and 2 in series 2 chips, I am not sure whether that 2-fold change (M of series 2/ M of series 1) will be shown as 2-fold change of individual Cy5 intensities (R of series 2/ R of series 1) also.
> I do not think that the ratio of individual Cy5 intensities do not become 2-fold in this example, because normalization works on M values, which would make M values comparable across arrays but not individual Cy5 intensities. Is this reasoning correct?
If you want to compare intensities, you will probably need to do some
between-array normalization. The limma user guide has some examples
of ways to do this keeping in mind that you do have two channels of
More information about the Bioconductor