[BioC] Another question about dye swap

Ron Ophir ron.ophir at weizmann.ac.il
Tue Feb 28 14:02:32 CET 2006


Hi,
>From limma user guide section 8.1.2 Simple comparisons -> Dye swaps it
is clear that for the following experimental design:
FileName Cy3 Cy5
File1        wt mu
File2        mu wt
File3        wt mu
File4        mu wt
*(four replicates of two groups (wt , mu) of which two replicates in
each group is labeled by one color(red) the other two is labeled by
another color (green)) 
one can estimate the Dye effect by defining that effect as the
intercept :
 design <- cbind(DyeEffect=1,MUvsWT=c(1,-1,1,-1))
 fit <- lmFit(MA, design)
 fit <- eBayes(fit)
 
Let's say that one finds that the Dye effect (for all genes) is
statisticaly insignificant. Would it be correct than to force the model
to go through the intercept by

design <- cbind(DyeEffect=-1,MUvsWT=c(1,-1,1,-1))

Will I get a better estimation of the Mu to wt ratio?

Thanks,
Ron



More information about the Bioconductor mailing list