Question: genes repeated in the heatmap several times(10x) when using limma
0
8 months ago by
andreilazar90 wrote:

Hello,

In the output of the heatmap, there are some repeated genes for several times, and I can't figure out why.  Obvious, something is wrong.

I try to compare a control group (which is in triplicate: C1, C2, C3) with a treated group which is also a triplicate A (A1, A2, A3), and individually, with another treatment group B (B1, B2, B3).  For the sake of example, I had simplified the group's names to letters.

The code that I use for contrast is:

"comparisonmodel <- model.matrix(~0+design)
> colnames (comparisonmodel) <- levels(design)

fit <- lmFit(NmData, comparisonmodel)

#Applying the empirical Bayes method
project.fit.eBayes <- eBayes(fit)
names(project.fit.eBayes)

#Make individual contrasts and fit the new model
CaseControl <- makeContrasts(A-C, B-C, levels=comparisonmodel)
CaseControl.fitmodel <- contrasts.fit(fit, CaseControl)
CaseControl.fitmodel.eBayes <- eBayes(CaseControl.fitmodel)

> #Filtering Results
> nrow(topTable(CaseControl.fitmodel.eBayes, coef=1:2, number=1000, lfc=2))
[1] 925
> probeset.list <- topTable(CaseControl.fitmodel.eBayes, coef=1:2, number=100, adjust.method="BH", sort.by="F", lfc=2)"

Is it possible that the repetition of gens in the heatmap to be because of the contrasts?

Here is a link where you can observe the heatmap.

https://we.tl/t-4o7TXZw7mB

Thanks

ADD COMMENTlink
modified 8 months ago by Axel Klenk930 • written 8 months ago by andreilazar90
Answer: genes repeated in the heatmap several times(10x) when using limma
1
8 months ago by
Aaron Lun24k
Cambridge, United Kingdom
Aaron Lun24k wrote:

You don't actually show the code that you used to make the heatmap, so it's hard to tell what's happening here. I can't see how the effect you're describing would be introduced by limma, I'd suggest you check if you have duplicated entries in NmData itself.

P.S. Don't use lfc in topTable, see the Note in ?topTable.

ADD COMMENTlink modified 8 months ago • written 8 months ago by Aaron Lun24k
Answer: genes repeated in the heatmap several times(10x) when using limma
0
8 months ago by
Axel Klenk930
Switzerland
Axel Klenk930 wrote:

Wild guess (as you don't say where your data come from): some Agilent microarrays contain 10 copies of a subset of their oligos.

ADD COMMENTlink written 8 months ago by Axel Klenk930
Please log in to add an answer.

Content
Help
Access

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