跟着Nature Genetics学作图:R语言ggpairs散点图一次性展示很多个主成分/拼图
语言 展示 Nature 一次性 跟着 很多 作图 成分
2023-06-13 09:16:32 时间
论文
Plasma proteome analyses in individuals of European and African ancestry identify cis-pQTLs and models for proteome-wide association studies
https://www.nature.com/articles/s41588-022-01051-w
本地pdf s41588-022-01051-w.pdf
代码链接
https://zenodo.org/record/6332981#.YroV0nZBzic
https://github.com/Jingning-Zhang/PlasmaProtein/tree/v1.2
今天的推文重复一下论文中的Extended Data Fig. 10
image.png
部分示例数据截图
image.png
这个是已经做好了主成分分析,现在只是可视化展示结果
读取数据
pc2<-read.delim(file = "data/20220627/ExtendedDataFig10.txt",
header = TRUE,
sep="\t")
head(pc2)
作图代码
library(ggplot2)
library(GGally)
pc.pr <- ggpairs(pc2[,3:7],
aes(color = pc2$V1),
upper = list(continuous = "points"),
diag = list(continuous = "blank")) +
scale_color_manual(values=c("#238b45","#2171b5")) +
theme(panel.background = element_blank(),
legend.position = "bottom",
axis.text = element_text(size = 6),
axis.title = element_text(size = 7)) +
scale_x_continuous(breaks = c(-0.01,0,0.01),
labels = c(-0.01,0,0.01)) +
scale_y_continuous(breaks = c(-0.01,0,0.01),
labels = c(-0.01,0,0.01))
pc.pr
image.png
拼图代码
pc.pr+
scale_color_manual(values = c("#f47720","#459943")) -> p2
library(patchwork)
pdf(file="Rplot05.pdf",
width = 9.4,
height = 4)
wrap_elements(ggmatrix_gtable(pc.pr))+
wrap_elements(ggmatrix_gtable(p2))
dev.off()
image.png
新知识点:ggpairs()
函数作图后的拼图代码wrap_elements(ggmatrix_gtable(pc.pr))+ wrap_elements(ggmatrix_gtable(p2))
参考链接
https://github.com/thomasp85/patchwork/issues/100
示例数据和代码可以自己到论文中获取
相关文章
- c语言数组介绍
- 详述 Java 语言中的格式化输出
- 【测开技能】Java语言系列(十八 )多维数组
- 数据分享|R语言用主成分PCA、 逻辑回归、决策树、随机森林分析心脏病数据并高维可视化|附代码数据
- 2023-03-08:x265的视频编码器,不用ffmpeg,用libx265.dll也行。请用go语言调用libx265.dll,将yuv文件编码成h265文
- 跟着Nature学作图:R语言ggplot2画误差线展示广义线型模型比值比(odds ratio)
- 跟着Nature Communication学作图:R语言ggplot2画世界地图展示采样地点
- 跟着Nature Ecology&Evolution学作图:R语言ggmsa包展示多序列比对结果
- 跟着NatureMetabolism学作图:R语言ggplot2画热图展示基因表达量
- 跟着Nature Communications学作图:R语言ggplot2做堆积柱形图展示群体基因组学的结果
- 跟着Nat. Biotechnol.学作图:R语言ComplexHeatmap热图展示TCGA-BRCA数据集概况
- 生信学习小组Day4 R语言的基础--大姚
- go语言 Linux部署 chatgpt 微信机器人指南
- Go 语言实现创建型设计模式 - 工厂模式
- Go语言圣经-Slice切片习题详解编程语言
- xmlLinux C 语言解析 XML 数据的实现方式(linuxc解析)
- 实用案例展示:如何利用c语言实现Redis应用(credis应用实例)