Download limma: Linear Models for Microarray Data User's Guide (Now
Transcript
Chapter 12 Separate Channel Analysis of Two-Color Data Separate channel analysis is a way to analyse two-color data in terms of the individual channel intensities [41]. In effect, separate channel analysis converts a two-color experiment into a single channel experiment with twice as many arrays but with a technical pairing between the two channels that originated from the same array. Consider an experiment comparing young and old animals for both both wild-type and mutant genotypes. FileName Cy3 Cy5 File1 wt.young wt.old File2 wt.old wt.young File3 mu.young mu.old File4 mu.old mu.young Each of the arrays in this experiment makes a direct comparison between young and old RNA targets. There are no arrays which compare wild-type and mutant animals. This is an example of an unconnected design in that there are no arrays linking the wild-type and mutant targets. It is not possible to make comparisons between wild-type and mutant animals on the basis of log-ratios alone. So to do this it is necessary to analyze the red and green channels intensities separately, i.e., to analyze log-intensities instead of log-ratios. It is possible to do this using a mixed model representation which treats each spot as a randomized block [48, 37]. Limma implements mixed model methods for separate channel analysis which make use of shrinkage methods to ensure stable and reliable inference with small numbers of arrays [37]. Limma also provides between-array normalization to prepare for separate channel analysis, for example > MA <- normalizeBetweenArrays(MA, method="Aquantile") scales the intensities so that A-values have the same distribution across arrays. The first step in the differential expression analysis is to convert the targets frame to be channel rather than array orientated. > targets2 <- targetsA2C(targets) > targets2 61