Another question about dye swap
0
0
Entering edit mode
Ron Ophir ▴ 270
@ron-ophir-1010
Last seen 9.7 years ago
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
GO limma GO limma • 777 views
ADD COMMENT

Login before adding your answer.

Traffic: 733 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6