Question: genes repeated in the heatmap several times(10x) when using limma
0
gravatar for andreilazar9
12 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 11 months ago by Axel Klenk940 • written 12 months ago by andreilazar90
Answer: genes repeated in the heatmap several times(10x) when using limma
1
gravatar for Aaron Lun
12 months ago by
Aaron Lun25k
Cambridge, United Kingdom
Aaron Lun25k 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 12 months ago • written 12 months ago by Aaron Lun25k
Answer: genes repeated in the heatmap several times(10x) when using limma
0
gravatar for Axel Klenk
11 months ago by
Axel Klenk940
Switzerland
Axel Klenk940 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 11 months ago by Axel Klenk940
Please log in to add an answer.

Help
Access

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