前往小程序,Get更优阅读体验!
立即前往
首页
学习
活动
专区
工具
TVP
发布
社区首页 >专栏 >手把手教你绘制多个基因表达值小提琴图

手把手教你绘制多个基因表达值小提琴图

作者头像
百味科研芝士
发布2020-04-30 18:38:02
3.7K0
发布2020-04-30 18:38:02
举报
文章被收录于专栏:百味科研芝士

数据准备

代码语言:javascript
复制
#set your work directory
data<-read.csv(file = "data.csv",header = T,sep = ",")
colnames(data)<-c("gene",colnames(data)[-1])
test<-data[1:5,1:5]
View(test)

###
data<-aggregate(data[,2:ncol(data)],by=list(data$gene),FUN = mean,na.rm=T)##重复基因求平均值
genename<-data$Group.1
rownames(data)<-genename
data<-data[,-1]
sampleID<-colnames(data)
data<-apply(data, 1, as.numeric)
rownames(data)<-sampleID
##分组信息构造
group<-c(rep("N",14),rep("DCIS",9),rep("IBC",9),rep("NS",14),rep("DCISS",11),rep("IBCS",9))
length(group)==dim(data)[[1]]##确认信息匹配
代码语言:javascript
复制
## [1] TRUE
代码语言:javascript
复制
data<-as.data.frame(data)
data$group<-group
table(data$group)
代码语言:javascript
复制
## 
##  DCIS DCISS   IBC  IBCS     N    NS 
##     9    11     9     9    14    14
代码语言:javascript
复制
data[1:5,1:5]##行为sample名 列为gene
代码语言:javascript
复制
##      ACVR1B    CXCR4     IL11    INHBA      LTB
## N1 3.812931 3.669153 2.922135 2.916415 3.072535
## N2 2.968409 4.270815 2.988273 3.068690 2.757790
## N3 3.806364 4.327977 2.811002 2.546459 2.897796
## N4 3.621140 3.698909 2.815539 3.278493 2.785017
## N5 3.812910 4.214652 2.814199 2.800704 2.916013

整合小提琴图

代码语言:javascript
复制
# http://www.sthda.com/english/articles/24-ggpubr-publication-ready-plots/
library(ggpubr)
代码语言:javascript
复制
## Loading required package: ggplot2
代码语言:javascript
复制
## Loading required package: magrittr
代码语言:javascript
复制
my_comparisons <- list(c("IBC", "DCIS"), c("N", "DCIS"), c("IBCS", "DCISS"),c("NS" ,"DCISS"))##分组设定
#my_comparisons <- list( c("IBC", "DCIS"), c("N", "DCIS") )##分组设定
e<-data %>% 
  #dplyr::filter(group %in% c("IBC","N","DCIS")) %>% 
  ggviolin(x = "group", y = c(colnames(data)[1:6]), fill = "group",
           combine = T,
         #palette = c("#00AFBB", "#E7B800", "#FC4E07"),##
         ylab="Normalized Expression",
         add = "boxplot", add.params = list(fill = "white"))
e+stat_compare_means(method = "t.test",
                     #label = "p.signif",##星号设置
                     comparisons = my_comparisons)
代码语言:javascript
复制
                     #theme_gray(base_size = 14)##background
ggsave("all_violin.pdf", width = 10, height = 8)     

小提琴图DCIS VS. IBC/N

代码语言:javascript
复制
library(ggpubr)
my_comparisons <- list( c("IBC", "DCIS"), c("N", "DCIS") )##分组设定
e<-data %>% 
  dplyr::filter(group %in% c("IBC","N","DCIS")) %>% #筛选行
  ggviolin(x = "group", y = c(colnames(data)[1:6]), fill = "group",
           combine = T,
         #palette = c("#00AFBB", "#E7B800", "#FC4E07"),##颜色设置
         ylab="Normalized Expression",
         add = "boxplot", add.params = list(fill = "white"))
e+stat_compare_means(method = "t.test",
                     #label = "p.signif",##星号设置
                     comparisons = my_comparisons)
代码语言:javascript
复制
                     #theme_gray(base_size = 14)##background
ggsave("group1_violin.pdf", width = 10, height = 8)     

小提琴图DCISS VS. IBCS/NS

代码语言:javascript
复制
library(ggpubr)
my_comparisons <- list(c("IBCS","DCISS"),c("NS","DCISS"))##分组设定
e<-data %>% 
  dplyr::filter(group %in% c("IBCS","NS","DCISS")) %>% 
  ggviolin(x = "group", y = c(colnames(data)[1:6]), fill = "group",
           combine = T,
         #palette = c("#00AFBB", "#E7B800", "#FC4E07"),##
         ylab="Normalized Expression",
         add = "boxplot", add.params = list(fill = "white"))
e+stat_compare_means(comparisons = my_comparisons)
代码语言:javascript
复制
ggsave("group2_violin.pdf", width = 10, height = 8)     

箱线图

代码语言:javascript
复制
my_comparisons <- list(c("IBC", "DCIS"), c("N", "DCIS"), c("IBCS", "DCISS"),c("NS" ,"DCISS"))##分组设定
e<-data %>% 
  #dplyr::filter(group %in% c("IBCS","NS","DCISS")) %>% 
  ggboxplot(x = "group", y = c(colnames(data)[1:6]), fill = "group",
           combine = T,
         #palette = c("#00AFBB", "#E7B800", "#FC4E07"),##
         ylab="Normalized Expression")
e+stat_compare_means(comparisons = my_comparisons)
代码语言:javascript
复制
ggsave("group2_boxplot.pdf", width = 10, height = 8)    
本文参与 腾讯云自媒体同步曝光计划,分享自微信公众号。
原始发表:2020-04-20,如有侵权请联系 cloudcommunity@tencent.com 删除

本文分享自 百味科研芝士 微信公众号,前往查看

如有侵权,请联系 cloudcommunity@tencent.com 删除。

本文参与 腾讯云自媒体同步曝光计划  ,欢迎热爱写作的你一起参与!

评论
登录后参与评论
0 条评论
热度
最新
推荐阅读
目录
  • 数据准备
  • 整合小提琴图
  • 小提琴图DCIS VS. IBC/N
  • 小提琴图DCISS VS. IBCS/NS
  • 箱线图
领券
问题归档专栏文章快讯文章归档关键词归档开发者手册归档开发者手册 Section 归档