使用手册Getting Started • patchwork
'''{
#2.生成一个随机数(rnorm)组成的10行6列的矩阵,列名为sample1,sample2….sample6,行名为gene1,gene2…gene10,分组为sample1、2、3属于A组,sample4、5、6属于B组。用循环对每个基因画ggplot2箱线图。
library(tidyverse)
a <- matrix(rnorm(60),nrow = 10)
colnames(a) <- paste0("sample",1:6)
rownames(a) <- paste0("gene",1:10)
a1 <- t(a) %>%
as.data.frame() %>%
mutate(group=c(rep("N",3),rep("T",3)))
p <- list()
for (i in 1:(ncol(a1)-1)) {
ggp <- ggplot(a1,aes_string(x = "group",y=colnames(a1)[i],colour="group"))+
geom_boxplot()+
geom_jitter()
p[[i]] <- ggp
}
library(patchwork)
wrap_plots(p,nrow = 2,guides = "collect")
}'''
网友评论