RPKM representation in heatmap is not accurate
1
0
Entering edit mode
ayamani • 0
@ayamani-16466
Last seen 5.7 years ago

Hello

So am trying to do a heatmap for my RNA seq data.  My RPKM range is 0.03-1253. the weird thing is that the heat map is showing a value of 1 RPKM as high red color. So, there seems to be a mistake in the accuracy of the z-score. Also,if some of my values are similar e.g a gene has RPKM values are 18.13 , 18.26 and 19.21; however in the heatmap u see it as blue and red indicating that its very different. I could not figure what went wrong. I would appreciate any suggestions

 
> install.packages("gplots")
> library("Heatplus") 
> library("RColorBrewer") 
> RPKM=read.csv("marker.csv", row.names = 1) 
> RPKM_MATRIX=as.matrix(t(RPKM))
 > color_palette=colorRampPalette(c("blue","red"))(n=299) 
> heatmap.2(RPKM_MATRIX,main="DESeq",trace="none",margins=c(10,12),scale="col",col=color_palette,dendrogram="none")
> sessionInfo() 
R version 3.5.0 (2018-04-23)
Platform: x86_64-w64-mingw32/x64 (64-bit)
Running under: Windows >= 8 x64 (build 9200)

Matrix products: default

locale:
[1] LC_COLLATE=English_United States.1252  LC_CTYPE=English_United States.1252   
[3] LC_MONETARY=English_United States.1252 LC_NUMERIC=C                          
[5] LC_TIME=English_United States.1252    

attached base packages:
[1] stats     graphics  grDevices utils     datasets  methods   base     

other attached packages:
[1] RColorBrewer_1.1-2 Heatplus_2.26.0    gplots_3.0.1      

loaded via a namespace (and not attached):
[1] compiler_3.5.0     tools_3.5.0        KernSmooth_2.23-15 gdata_2.18.0       caTools_1.17.1    
[6] bitops_1.0-6       gtools_3.8.1

 

heatmap.2 z-score • 770 views
ADD COMMENT
0
Entering edit mode
Aaron Lun ★ 28k
@alun
Last seen 21 hours ago
The city by the bay

gplots is not a Bioconductor package. Please ask your question in a more appropriate forum, e.g., R-help.

ADD COMMENT

Login before adding your answer.

Traffic: 763 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6