R绘图-KEGG功能注释组间差异分面条形图
生活随笔
收集整理的這篇文章主要介紹了
R绘图-KEGG功能注释组间差异分面条形图
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
一、 數據準備
KEGG功能注釋的差異分析結果繪圖。使用的是經過Pairwise Mann–Whitney U-tests檢驗的處理間class II功能豐度差異結果數據,包含class I作為分面信息進行條形圖繪制。
library(ggplot2) library(ggsci) library(ggpubr) library(scales)# 1.1 讀入數據 data = read.csv("KEGG.csv",header=TRUE,sep=",",stringsAsFactors = FALSE) dim(data) # 44個classII data$Class2 = factor(data$Class2,levels = data$Class2[1:44]) # 設置class II數據因子水平 head(data)# 1.2 數據格式轉換-“寬”轉“長” library(tidyr) colnames(data) = c("Class2","Class1","mean","sd","mean","sd","p.adj") # 將CK與PT的均值和標準差的列名設置為一樣,后續才能按行合并。 data = rbind(data.frame(data[,c(1:4,7)],group = rep("CK",time =44)),data.frame(data[c(1,2,5:7)],group = rep("PT",time=44))) # 按行合并CK和PT的mean和sd,并添加分組信息列。 head(data) dim(data) # 1.3 數據過濾,只保留具有統計學意義(p.adj<0.05)的功能 data = data[data$p.adj <總結
以上是生活随笔為你收集整理的R绘图-KEGG功能注释组间差异分面条形图的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 打造自己的图像识别模型
- 下一篇: c 语言除法运算,C 语言简单加减乘除运