用R画RNA_seq中变化显著的heatmap聚类热图 - 做专业的生信问答网站

登录

该问题已被锁定！

4: 关注

2970: 浏览

用R画RNA_seq中变化显著的heatmap聚类热图

我看了一篇文章https://www.sciencedirect.com/science/article/pii/S1674205217302459?via%3Dihub#fig1，想做一个和文章里一样的聚类热图，如下 [attach]36[/attach]，但是我不是很明白它的这个分类的方式，我自己做出来的如下： [attach]35[/attach] 代码为 >data <- log2(cuffdiff_result[,c(7,8)]+1)) >Heatmap(x, name = "expression", km = 5,, show_row_names = FALSE, show_column_names = FALSE) x想知道这样子对不对？有没有什么问题，如果分类之后怎么导出分类数据

好问题 0 评论收藏举报

查看全部 3 个回答

孟浩巍超级管理员用户来自于: 北京市
2018-08-24 13:37

我画heatmap更喜欢用pheatmap这个函数。 [code]ph = pheatmap(your_mat,cluster_cols = F, cluster_rows = T)[/code]上面的代码就是对列不聚类，对行聚类，并把结果保存到了ph变量里。同时，在ph里面可以使用下面的代码查看你的原始矩阵行的聚类结果 [code]# 聚类结果 ph$tree_row # 每一行的排序结果 ph$tree_row$order[/code]

赞同 3 6评论

关于作者

: jailyn 前台管理员
这家伙很懒，还没有设置简介

7: 回答

0: 文章

5: 问题

问题动态

发布时间: 2018-08-24 12:19

更新时间: 2018-08-24 23:19

关注人数: 4 人关注

相关问题

RNA Seq中什么样的reads是chimeric reads？是不是多数的chimeric reads都是circRNA的？: 2196 浏览 2 关注 1 回答 0 评论

enhancer结合ATAC-seq进行下游靶基因的寻找策略: 2419 浏览 2 关注 1 回答 0 评论

是不是单细胞的chip-seq，cut run这些call peak是不需要对照的，我看代码里没有-c.所以Rna-seq的峰图也是不需要control来call peak的吗？: 2282 浏览 2 关注 1 回答 0 评论

scRNA-seq有什么推荐的监督聚类方式吗？: 2320 浏览 2 关注 1 回答 0 评论

RNA-seq不同样本多个生物学重复不同处理条件下的如何找差异基因: 2575 浏览 2 关注 2 回答 0 评论

sgRNA的GC content 是什么意思，该如何计算: 2731 浏览 2 关注 2 回答 0 评论

ATAC-seq不同测序量样本的差异比较: 3018 浏览 2 关注 1 回答 0 评论

R基础绘图问题-柱状图加误差线及T测验显著性星号，映射问题如何解决？: 1498 浏览 1 关注 0 回答 1 评论

用wget 在服务器上用wget下载ENcode里的chip-seq测序原始数据，提示可以connected,但一直no data received ?: 2553 浏览 2 关注 2 回答 0 评论

htseq使用出现故障: 2028 浏览 2 关注 1 回答 0 评论

推荐内容

问 Rstudio怎么打开project: 2887 浏览 2 关注 1 回答 0 评论

问使用cutadapt去掉adaptor后，fastqc结果的sequence distribution报错: 2849 浏览 2 关注 2 回答 0 评论

问转录组数据样本聚类结果不理想: 2601 浏览 2 关注 1 回答 0 评论

问 Rstudio安装R包时，出现unable to move temporary installation: 3462 浏览 4 关注 3 回答 1 评论

问画heatmap出现错误提示: 1670 浏览 2 关注 1 回答 1 评论

问如何快速将基因名称转换为基因ID(拟南芥）: 2535 浏览 2 关注 1 回答 0 评论

问 WGCNA中的TOM热图绘制: 2887 浏览 2 关注 1 回答 0 评论

问常用的对测序序列进行聚类分析的软件有什么？及其使用的了什么算法？: 3031 浏览 2 关注 2 回答 0 评论

All Rights Reserved Powered BY WeCenter V4.1.0 © 2025