手把手教你繪製多個基因表達值小提琴圖

2021-03-02 百味科研芝士
image.png數據準備

#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]]##確認信息匹配

## [1] TRUE

data<-as.data.frame(data)
data$group<-group
table(data$group)

## 
##  DCIS DCISS   IBC  IBCS     N    NS 
##     9    11     9     9    14    14

data[1:5,1:5]##行為sample名 列為gene

##      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

整合小提琴圖

# http://www.sthda.com/english/articles/24-ggpubr-publication-ready-plots/
library(ggpubr)

## Loading required package: ggplot2

## Loading required package: magrittr

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)

image.png

                     #theme_gray(base_size = 14)##background
ggsave("all_violin.pdf", width = 10, height = 8)     

小提琴圖DCIS VS. IBC/N

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)

image.png

                     #theme_gray(base_size = 14)##background
ggsave("group1_violin.pdf", width = 10, height = 8)     

小提琴圖DCISS VS. IBCS/NS

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)

image.png

ggsave("group2_violin.pdf", width = 10, height = 8)     

箱線圖

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)

image.png

ggsave("group2_boxplot.pdf", width = 10, height = 8)    

相關焦點

  • Excel表格中如何繪製稜錐圖 手把手教你在excel2007中插入稜錐圖
    Excel表格中如何繪製稜錐圖 手把手教你在excel2007中插入稜錐圖時間：2017-07-02 14:14   來源：三聯   責任編輯：沫朵 川北在線核心提示：原標題：Excel表格中如何繪製稜錐圖 手把手教你在excel2007中插入稜錐圖 1、打開一個Excel的文件，選中一個單元格，然後滑鼠左鍵單擊菜單【插入】
  • Graphpad 作圖教程|手把手教你繪製柱形圖
    今天給大家介紹 GraphPad Prism 如何繪製柱狀圖。5.1單方向方差分析的 P 值和 F 值。5.2三組中兩兩分析結果。6.作圖：點擊 Graphs 模塊中的 DATA 1，按下圖進行選擇，並生成草圖。
  • Python繪圖筆記:繪製有趨勢線的直方圖、小提琴圖和誤差圖
    1.小提琴圖和箱線圖#同時繪製箱線圖和小提琴圖，箱線圖清楚地展示了所繪製的數據集中的四分位數和離群點#小提琴圖則清晰展示了數據集中的分布情況import matplotlib.pyplot as pltimport numpy as np
  • Graphpad Prism 8.0進階篇-繪製小提琴圖
    小提琴圖展示的是數據的分布形式，是箱線圖與核密度圖的結合。下面為大家介紹小提琴圖的繪製詳細步驟：　　1 　　打開Graphpad Prism 8.0，新版界面如下，在New table＆graph下選擇Column，Data table下選擇Enter or import data a new table，Options下選擇Enter replaced values，stacked into columns
  • 實驗數據完成了,如何繪製為柱狀圖 | 以qPCR為例講解
    實驗數據完成了，如何繪製為柱狀圖？本文以實時螢光定量PCR（qPCR）為例講解。qPCR是通過螢光信號對擴增進程實時檢測，由於其靈敏度高、特異性強，樣本的Ct值和與基因拷貝數有關，現已作為科研實驗中重要的驗證方法。
  • Graphpad 作圖教程 | 手把手教你繪製「森林圖」
    (forest plots)是以估計模型 (固定效應模型或是隨機效應模型) 結果為基礎繪製出的圖型。它以一條垂直的直線 (橫坐標刻度為 1 或 0) 為中心，用平行於橫軸的多條線段描述了每個被納入研究的效應量和可信區間 (confidence interval，CI), 用一個稜形 (或其它圖形) 描述了多個研究合併的效應量及可信區間。我們以 Graphpad Prism 8.0 為例，繼續教大家如何繪製回歸分析結果的森林圖。
  • 手把手教你製作ppt日記本
    手把手教你製作ppt日記本時間：2017-07-16 15:06   來源：三聯   責任編輯：沫朵 川北在線核心提示：原標題：ppt怎麼製作筆記本? 手把手教你製作ppt日記本 ppt怎麼製作筆記本?本文介紹了使用ppt製作日記本的方法，製作方法簡單，一起來學習吧!
  • Origin繪製縮放圖
    實驗中，我們能夠得到很多數據，但你需要的可能只是其中的一部分，或者其中一部分有顯著差異，需要放大圖形以突出重點；那怎麼在展示數據完整性的同時，又能放大其中的一部分呢；今天，小編就教你用Origin繪製可以局部放大的縮放圖。
  • 基因表達圖譜繪製
    基因表達圖譜繪製 作者：小柯機器人 發布時間：2019/11/21 13:04:43 德國亥姆霍茲協會馬克斯·德爾布呂克分子醫學中心Nikolaus Rajewsky和以色列希伯來大學Nir
  • 手把手教你學ELISA、PCR、免疫組化
    1.手把手教你學ELISAELISA的基礎是抗原或抗體的固相化及抗原或抗體的酶標記。結合在固相載體 表面的抗原或抗體仍保持其免疫學活性,酶標記的抗原或抗體既保留其免疫學活性,又保留酶的活性。2.手把手教你學PCR實時螢光定量PCR，簡稱RT-QPCR,屬於Q-PCR的一種，目前該技術已得到廣泛應用，如：擴增特異性分析、基因定量分析、基因分型、SNP分析等。
  • 傻瓜式圖文教程:Genecards基因信息檢索與分析、基因結構圖繪製......
    要研究基因功能，首當其衝的就是要全方位了解基因信息，如基因與表型間的聯繫、基因互作蛋白分子、信號通路、臨床意義等。
  • 用Seaborn繪製圖表
    Seaborn可以很好地處理數據幀，而Matplotlib則不行，它可以讓你以更簡單的方式繪製引人注目的圖表。這個色調的參數是非常有用的，因為它允許你用顏色來表達第三維度的信息。sns.lmplot(x=』Attack』, y=』Defense』, data=df, fit_reg=False, hue=』Stage』)plt.show()散點圖看起來和之前的一樣，除了中間沒有回歸線，而且每個點的顏色也不同。這些顏色只是顯示了每個小精靈的階段。
  • 除了火山圖,差異表達基因還可以這樣展示
    您可能第一時間想到可以使用火山圖。的確，火山圖是使用頻率最多的，在火山圖中可以很輕鬆地根據基因在兩組間的Fold Change值以及顯著性p值，識別和判斷差異表達基因概況。火山圖實質上就是一種散點圖，通常橫縱坐標分別代表了log2轉化後的Fold Change以及-log10轉化後的p值或p調整值信息（下圖左）。
  • 分享ggpubr 包:教你繪製密度圖、直方圖、柱狀圖
    今天我們要分享的R包是 ggpubr 包，它是一款基於ggplot2的可視化包，功能非常強大，能夠一行命令繪製出符合出版物要求的圖形。ggpubr 包可繪製的圖形類型非常多，有密度圖、直方圖、柱狀圖、餅圖、棒棒糖圖、Cleveland 點圖、箱線圖、小提琴圖、點帶圖、點圖、散點圖、線圖、誤差棒圖……哈哈，有木有很期待接下來這個 ggpubr 包的學習了呢！
  • 箱形圖、小提琴圖、直方圖……統統可以卡通化!
    分布族這是一個幸福的大家庭——箱形圖、直方圖、小提琴圖、箱形圖。雖然它們長相各異，但如果深入了解，你就會發現每種圖形都以相似的方式揭示數據集的構成。爸爸媽媽：箱形圖Borat和Wanda這個圖表家庭會告訴你，用大量不同的數據樣本測量同一個值的意義何在。
  • 「簡筆畫教程」:手把手教你畫人物簡筆畫,簡單又好看,菜鳥也會
    今日繪本：手把手教你畫風格：簡筆畫01工具準備A4白紙、橡皮擦、鉛筆、黑色彩鉛、0.1針管筆、鋼筆、紅色筆、黑色馬克筆02鉛筆稿輕輕用鉛筆勾勒出人物輪廓03鋼筆線稿用鋼筆描出人物主要輪廓部分，可對鉛筆稿進行修改04眼睛細節繪製工具：0.1針管筆、
  • 科學家繪製出胎兒基因表達的人類細胞圖譜
    科學家繪製出胎兒基因表達的人類細胞圖譜 作者：小柯機器人 發布時間：2020/11/16 13:42:53 美國華盛頓大學Jay Shendure、Cole Trapnell等研究人員合作繪製出胎兒基因表達的人類細胞圖譜。
  • 手把手教你看懂視力表和驗光單,10分鐘小白變專家!
    孩子們每年都要檢查視力，作為家長的你，是否能看懂圖表上的1.0和5.0代表什麼意思嗎？今天，小編手把手教你如何看懂視力表和驗光單，10分鐘一、視力表如何看？實際上1.0和5.0這兩個都是標準的視力線，區別是用不同的計數方法來表達視力情況：1.數值在4.0-5.3之間的是五分制視力表，5.0是正常的標準視力線
  • 手把手教你繪製信號通路圖
    細胞信號通路圖是科研研究過程中最常見也是最常用到的，如何繪製適合我們自己科研課題的信號通路圖呢？可以試試 pathway builder tool 軟體。這款軟體簡單易學，即便是零基礎的同學，也可以做出漂亮的信號通路。1. 首先，打開Pathway Builder Tool 2.0 軟體，軟體自帶分子生物學會用到的基本元素，如不同的細胞、細胞器、分子甚至箭頭等等。