Question: How to plot heatmap of differentially expressed genes using limma results
gravatar for saw44
4 weeks ago by
saw4410 wrote:

I have analysed some microarrays using limma and would like to plot a heatmap of the most differentially expressed genes, say the top 50.

I understand coolmap is the best tool for this - however when I run my code I get the following error - how should I fix this?

RG <- read.maimages(files, 
                    path= "R:\\Transposons\\Shared\\OneDrive_3_27-11-2019", 
                    green.only = TRUE)

RG <- backgroundCorrect(RG, method="normexp")
MA <- normalizeBetweenArrays(RG, method = "quantile")

design <- cbind(WT=c(1,1,1,1,1,1,0,0,0,0,0,0),MU=c(0,0,0,0,0,0,1,1,1,1,1,1))

fit <- lmFit(MA, design)
cont.matrix <- makeContrasts(MUvsWT=MU-WT, levels=design)
fit2 <-, cont.matrix)
fit3 <- eBayes(fit2)
topTable(fit3, adjust="BH")


> > coolmap(MA)
Error: cannot allocate vector of size 14.8 Gb
ADD COMMENTlink modified 4 weeks ago by _r_am31k • written 4 weeks ago by saw4410

Please do not use pre and code tags - use the 101010 option on the toolbar instead. Manually adding these HTML tags interferes with the site's functionality.


ADD REPLYlink written 4 weeks ago by _r_am31k
gravatar for Hamid Ghaedi
4 weeks ago by
Hamid Ghaedi770
Hamid Ghaedi770 wrote:

The error you got is about memory limitation which is a known error when you ask R to load a huge file. If you want to know more about this, follow these StackOverflow links:




For differentially expressed genes, you can use any package you want. I would like to use:

1- for quick visualization, pheatmap package. A nice tutorial could be found here.

2- For professional heatmap, ComplexHeatmap package. Its vignette could be found here. BUT Kevin Blighe wrote a concise (strongly recommended) tutorial which you may find here.

ADD COMMENTlink modified 4 weeks ago • written 4 weeks ago by Hamid Ghaedi770
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.3.0
Traffic: 1274 users visited in the last hour