其他
「盘点」R语言的那些拼图方法
大神一句话，菜鸟跑半年。我不是大神，但我可以缩短你走弯路的半年~
就像歌儿唱的那样，如果你不知道该往哪儿走，就留在这学点生信好不好~
这里有豆豆和花花的学习历程，从新手到进阶，生信路上有你有我！
花花写于 2019-09-18
接下来的被约稿的推文预告一下：
1、R 要不要更新，以及要避开的坑
2、Rstudio如何挂载2个版本的R
3、盘点处理R语言的那些缺失值处理方法
你有什么想约的稿吗？后台告诉我，免费的~谁让我最近缺素材呢，范围限于R语言，别的我可不会。
「盘点」系列目录：
1.par里的mfrow
可用于基础包拼图
par(mfrow = c(2,2))
plot(1:100 + rnorm(100))
plot(rnorm(100), type = "l")
hist(rnorm(500))
acf(rnorm(100))
2.grid.arrange
ggplot2拼图御用
https://www.jianshu.com/p/0970fdac9071
if(!require(dplyr))install.packages("dplyr")
library(dplyr)
test <- iris %>%
head(120) %>%
tail(40) %>%
rbind(head(iris,20)) %>%
mutate(n=1:60)
table(test$Species)
#>
#> setosa versicolor virginica
#> 20 20 20
if(!require(ggplot2))install.packages("ggplot2")
library(ggplot2)
p <- ggplot(data = test) +theme_bw()
colnames(test)
#> [1] "Sepal.Length" "Sepal.Width" "Petal.Length" "Petal.Width"
#> [5] "Species" "n"
##密度图
p_density <- p + geom_density(
aes(Sepal.Length, color = Species)
)+theme(legend.position = "none")
##箱线图
p_boxplot <- p + geom_boxplot(
aes(x = Species,y = Sepal.Length,color = Species)
)+theme(legend.position = "none")
##小提琴图
p_violin <- p + geom_violin(
aes(x = Species, y = Sepal.Length, fill = Species)
)+theme(legend.position = "none")
## 点图
p_point <- p + geom_point(aes(x=n,y = Sepal.Length))+theme(legend.position = "none")
## 管他什么图
p_col <- p + geom_col(aes(x=n,y=Sepal.Length,fill = Species))+theme(legend.position = "none")
## 凑数用
p_s <- p + geom_area(aes(x=n,y=Sepal.Length,fill = Species))+theme(legend.position = "none")
if(!require(gridExtra))install.packages("gridExtra")
library(gridExtra)
plots = list(p_density,p_boxplot,p_violin)
lay1 = rbind(c(1, 2),
c(3, 3)) #布局矩阵
grid.arrange(grobs = plots,
layout_matrix = lay1,
widths = c(1, 2), #第一列的宽度为1，第二列的宽度为2
heigths = c(1, 0.1)) #依然不知道为什么高度调整失败
让布局再复杂一点
lay2 = rbind(c(1, 3),
c(2, 3)) #布局矩阵
plots2 <- list(p_s,p_col,p_point)
p1 <- arrangeGrob(grobs = plots,
layout_matrix = lay1)
p2 <- arrangeGrob(grobs = plots2,
layout_matrix = lay2)
grid.arrange(p1, p2, ncol = 2)
grid.arrange(p1, p2, nrow = 2)
3.cowplot
可以给图加上ABCD编号蛮好的。
if(!require(cowplot))install.packages("cowplot")
library(cowplot)
gg <- ggdraw() +
draw_plot(p_boxplot, x=0, y=0.5, width=1, height=0.5) +
draw_plot(p_violin, 0, 0, 0.5, 0.5) +
draw_plot(p_density, 0.5, 0, 0.5, 0.5) +
draw_plot_label(c("A", "B", "C"), c(0, 0, 0.5), c(1, 0.5, 0.5), size = 15, colour = "black")
gg
向大家隆重推荐隔壁生信技能树的一系列干货！
点击底部的“阅读原文”，获得更好的阅读体验哦😻
初学生信，很荣幸带你迈出第一步。
我们是生信星球，一个不拽术语、通俗易懂的生信知识平台。由于是2018年新号，竟然没有留言功能。需要帮助或提出意见请后台留言、联系微信或发送邮件到jieandze1314@gmail.com，每一条都会看到的哦~