利用cutree划分pheatmap聚类结果
时间:2022-07-23
本文章向大家介绍利用cutree划分pheatmap聚类结果,主要内容包括其使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
之前小编教大家使用pheatmap快速绘制热图,直接利用cluster_rows对行进行聚类,但是聚类后我们怎么得到聚类结果呢?今天小编就教大家利用cutree划分pheatmap聚类结果。
首先,准备输入文件:
行为基因,列为样本,保存为制表符分隔的文本文件。
## 加载R包
library(pheatmap)
## 导入文件
exp <- read.table("input.txt",sep="t",header=T,row.names = 1)
## 绘制热图
p <- pheatmap(log2(exp+1),cellwidth=20, cellheight=10,cluster_cols=F,cluster_rows=T,cutree_rows = 3)
我们可以注意到,与之前绘制热图相比,多了一个cutree_rows参数。这个参数用来控制我们划分类群的数目。
可以看到,我们输入的20个基因划分为了3类。
现在我们来提取每类的基因ID。
## 对基因进行分类
row_cluster <- cutree(p$tree_row,k=3)
需要注意的时,k后的数字要与前一步cutree_rows的数字相同。
## 对结果进行排序
newOrder <- exp[p$tree_row$order,]
newOrder[,ncol(newOrder)+1]=row_cluster[match(rownames(newOrder),names(row_cluster))]
colnames(newOrder)[ncol(newOrder)]="Cluster"
## 输出结果
write.table(newOrder,"cluster.list.out",quote=F,row.names=T,col.names=T,sep="t")
执行完所有代码后我们就得到了分类结果,非常方便!大家快去试试吧~
参考资料:
https://www.omicsclass.com/article/508
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法