Entering edit mode
ayamani
•
0
@ayamani-16466
Last seen 6.5 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