R绘图|基因表达水平分布图绘制
2023-03-07 09:16:29 时间
在整理转录组下游的时候,看到中科新生命的报告中的基因表达水平分布部分有这么一个图
从图中可以非常直观的看出来不同样本在不同表达区间的分布情况。由于报告没有给出源代码,我们模仿的画一画。
想要画出这样一个基因表达水平分布图,我们需要两个东西
- 基因表达矩阵
- 数据的分布情况
基因表达矩阵
原始表达矩阵比较容易获取,为了方便演示,我们直接采取edgeR[1]的cpm标准化拿到基因表达矩阵。
rawcount <- read.table("data/rawcounts.txt",row.names = 1,
sep = "\t", header = T)
library(edgeR)
express_cpm <- log2(cpm(rawcount)+1)
标准化后的基因表达矩阵
数据的分布情况
接下来我们需要将现有的表达情况按一定标准分类,需要用到R包reshape[2]
# 载入R包
library(reshape)
# 宽变长
longdata <- melt(data = express_cpm)
# 将数据划分成6个区间
cut_data <- cut(as.numeric(longdata$value), breaks = c(0,1,5,10,30,50,Inf),right = F,include.lowest = TRUE)
# 将分组结果添加到longdata
longdata$group = cut_data
# 给区间命名
levels(cut_data) <- c("0<=CPM<1", "1<=CPM<5","5<=CPM<10","10<=CPM<30","30<=CPM<50","50<=CPM<+∞")
画图
画图需要用到R包ggplot2[3]
# 载入R包
library(ggplot2)
# 画图
ggplot(longdata, aes(x = X2, fill = group, y = value)) +
geom_bar(aes(fill = group), position = "fill",stat = "identity")+
scale_y_continuous(labels = scales::percent)+
theme(axis.text.x = element_text(angle = 45, hjust = 1))+
# 添加横纵坐标和title
labs(title = "不同表达水平区间的基因数量统计图",
x = "Sample", y = "Percentage", fill = "Group")
- angle = 45:设置样本名倾斜角度为45°
- hjust = 1:设置样本名距离图形的距离为1(自行调试)
- labs( ):添加坐标轴
参考资料
[1] edgeR: https://bioconductor.org/packages/release/bioc/html/edgeR.html
[2] reshape: https://www.rdocumentation.org/packages/stats/versions/3.6.2/topics/reshape
[3] ggplot2: https://www.rdocumentation.org/packages/ggplot2/versions/3.3.5
相关文章
- 在 Go 里用 CGO?这 7 个问题你要关注!
- 9款优秀的去中心化通讯软件 Matrix 的客户端
- 求职数据分析,项目经验该怎么写
- 在OKR中,我看到了数据驱动业务的未来
- 火山引擎云原生大数据在金融行业的实践
- OpenHarmony富设备移植指南(二)—从postmarketOS获取移植资源
- 《数据成熟度指数》报告:64%的企业领袖认为大多数员工“不懂数据”
- OpenHarmony 小型系统兼容性测试指南
- 肯睿中国(Cloudera):2023年企业数字战略三大趋势预测
- 适用于 Linux 的十大命令行游戏
- GNOME 截图工具的新旧截图方式
- System76 即将推出的 COSMIC 桌面正在酝酿大变化
- 2GB 内存 8GB 存储即可流畅运行,Windows 11 极致精简版系统 Tiny11 发布
- 迎接 ecode:一个即将推出的具有全新图形用户界面框架的现代、轻量级代码编辑器
- loongarch架构介绍(三)—地址翻译
- Go 语言怎么解决编译器错误“err is shadowed during return”?
- 敏捷:可能被开发人员遗忘的部分
- Denodo预测2023年数据管理和分析的未来
- 利用数据推动可持续发展
- 在 Vue3 中实现 React 原生 Hooks(useState、useEffect),深入理解 React Hooks 的