This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Values EXCEEDS the original data range after batch effect correction.

Hi

I have gene coverage (%) data in a matrix structure. since the coverage estimated is in percentage, it falls in range of 0-100%. now I ran limma to remove batch effect, and see that the range of gene coverage is now 0-6000

How could it be possible?

batch-effect limma

Does limma recommend you use percentages? How is a batch correction algorithm supposed to know the semantic limitations of your metric?

but limma takes normalized values, right? So my coverage was estimated in percentages, which states that it is a normalised value. That's why I used that as an input to limma

the usage for removeBatchEffects states that the input matrix should be a

numeric matrix, or any data object that can be processed by getEAWP containing log-expression values for a series of samples. Rows correspond to probes and columns to samples.

Using percentages is not consistent with this expectation.

0 answers

No answers yet.

Log in to answer this question.