Question: Agilent microarray chip normalization/background correction
gravatar for dimitri.leonid.lindenwald
3.4 years ago by

Hallo everyone!

I have encountered a problem while working with Agilent 4x44k whole mouse genome v2 microarray chips.
A part of my experiment was performed 1.5 years ago (pilot study), and the second part - just now.

All the data obtained from the both parts of the experiment has been background-corrected and normalized with limma. I performed a backgroundCorrect(method="normexp") and than used normalizeBetweenArrays(method="quantile").

The data from different parts of the experiment forms two distinctive clusters ("new" and "old") in the Principle component analysis (PCA), even though the samples are supposed to be the same and come from the same origin.

Did i oversee some requirements? Are the normalization/background correction methods more suitable to this kind of problems? Thank you!


ADD COMMENTlink modified 14 months ago by alerodriguez0 • written 3.4 years ago by dimitri.leonid.lindenwald50
gravatar for James W. MacDonald
3.4 years ago by
United States
James W. MacDonald48k wrote:

You won't be able to normalize out a batch effect, and it is not unexpected to see the two batches cluster like that. If  you ran all the groups you are planning to compare in both the pilot and the second batch, then you can a.) fit a batch effect in your linear model or b.) use ComBat() in the sva package to remove the batch effect. There are examples of doing both things in the users guide for limma, and in the sva vignette, plus a search for batch effect on this support site will turn up even more information.

If you ran one or more of your samples in just one of the batches, then you have unwittingly confounded biological differences with your batch effect, and there is no way to fix that.

ADD COMMENTlink written 3.4 years ago by James W. MacDonald48k
gravatar for alerodriguez
14 months ago by
alerodriguez0 wrote:

Regarding:a.) fit a batch effect in your linear model or b.) use ComBat() in the sva package to remove the batch effect.  

After removing batch effects using linear model the output data will be the residuals, what is the output form of the data when using ComBat? Residuals? Can you clarify, please?




ADD COMMENTlink written 14 months ago by alerodriguez0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.2.0
Traffic: 355 users visited in the last hour