我在R中运行了heatmap.2,并且无法生成彩色键以显示左上角的绘图密度。
我可以使用以下内容制作颜色键:
library(gplots)
a=cbind(a=1:5000,b=1:5000*2,c=1:5000*4,d=1:5000*8,e=1:5000*16)
h=hclust(dist(log(t(a))))
d=as.dendrogram(h)
dd=reorder(d,c(-10,1,1,1,100))
hm=colorRampPalette(c('yellow','green','red'))(nrow(a)*ncol(a))
png(file="CompleteHeatmap.png",width=1200,height=800)
par(oma=c(1,1,1,5))
heatmap.2(log(a),scale='none',trace='none',col=hm,density.info='density',denscol='black',Colv=dd)
有了这个,颜色键变得精确和准确。然而,当我使用我的实际数据(39000 * 6列)运行它时,我没有得到彩色的颜色键,尽管显示的密度图是准确的。
如果有人有任何想法可能会发生这种情况,我想听听,因为我的老板希望正确上色的颜色。
我的数据表的大小是上面矩阵的7倍,但除此之外,所有的调用都是相同的。
非常感谢
亚历
答案 0 :(得分:1)
我找到了答案,并用
更改了颜色托盘hm=colorRampPalette(c('yellow','green','red'))(256)
在nrow()* ncol()的位置,它就算完了!