经验首页 前端设计 程序设计 Java相关 移动开发 数据库/运维 软件/图像 大数据/云计算 其他经验
当前位置:技术经验 » 程序设计 » R语言 » 查看文章
画相关性热图
来源:cnblogs  作者:raisok  时间:2019/6/13 8:58:28  对本文有异议
  • 数据格式如下
  1. Sample CM11 CM12 CM13 CM21 CM22 CM23 CM31 CM32 CM33
  2. CM11 1 0.9627369 0.9884578 0.9841946 0.9762492 0.983613 0.9575127 0.743262 0.6178
  3. CM12 0.9627369 1 0.9616447 0.9405868 0.9354329 0.9442101 0.9816946 0.8740325 0.7634
  4. CM13 0.9884578 0.9616447 1 0.9827691 0.9754358 0.9845702 0.9736013 0.7541648 0.6165
  5. CM21 0.9841946 0.9405868 0.9827691 1 0.9949495 0.9914846 0.9445855 0.7201134 0.6109
  6. CM22 0.9762492 0.9354329 0.9754358 0.9949495 1 0.9911021 0.936078 0.7206804 0.6109
  7. CM23 0.983613 0.9442101 0.9845702 0.9914846 0.9911021 1 0.9469456 0.7364044 0.6163
  8. CM31 0.9575127 0.9816946 0.9736013 0.9445855 0.936078 0.9469456 1 0.8512908 0.7314
  9. CM32 0.743262 0.8740325 0.7541648 0.7201134 0.7206804 0.7364044 0.8512908 1 0.9350
  10. CM33 0.6178754 0.76344 0.6165455 0.6109218 0.6109854 0.6163989 0.731487 0.935035 1

作图代码如下:

  1. library(reshape2)
  2. library(ggplot2)
  3. library(RColorBrewer)
  4. x <- read.table("AllSamples.correlation.xls", sep = "\t", head = T)
  5. xx = as.matrix(x[,-1])
  6. rownames(xx) = names(x)[-1]
  7. xx = melt(xx)
  8. names(xx)=c("Var1","Var2","pearson_value");
  9. pdf("AllSamples.CorrelationHeatmap.pdf",width=9,height=9)
  10. ggplot(xx, aes(Var1, Var2, fill=pearson_value))+
  11. #geom_tile(width=0.8, height=0.8)+
  12. geom_tile(color='black')+
  13. geom_text(label=round(xx$pearson_value, 3))+
  14. scale_fill_gradient(low='#DEEBF7',high='#08519C')+
  15. theme(axis.text = element_text(angle=30, hjust=1,size=11,vjust=0,color='black'),
  16. panel.background = element_rect(fill='transparent'),
  17. panel.grid=element_line(color='grey'),legend.title = element_text(size = 13))+
  18. labs(x="",y="")
  19. dev.off()
  • 图片效果如下

     

原文链接:http://www.cnblogs.com/raisok/p/11010421.html

 友情链接:直通硅谷  点职佳  北美留学生论坛

本站QQ群:前端 618073944 | Java 606181507 | Python 626812652 | C/C++ 612253063 | 微信 634508462 | 苹果 692586424 | C#/.net 182808419 | PHP 305140648 | 运维 608723728

W3xue 的所有内容仅供测试,对任何法律问题及风险不承担任何责任。通过使用本站内容随之而来的风险与本站无关。
关于我们  |  意见建议  |  捐助我们  |  报错有奖  |  广告合作、友情链接(目前9元/月)请联系QQ:27243702 沸活量
皖ICP备17017327号-2 皖公网安备34020702000426号