跟着iMeta学做图|circlize绘制环状热图展示细菌功能聚类分析
# 农心生信工作室 # 2022-11-16 # 代码编写及注释:农心生信工作室 # # 热图 (Heat map) 可以在微生物组研究中展示展示细菌功能聚类分析的结果,而环状热图是热图的一种表现形式。 # 本期我们挑选2022年6月13日刊登在iMeta上的Microbial community roles and chemical mechan…
跟着iMeta学做图|ggplot2绘制曼哈顿图展现物种丰度差异
# 农心生信工作室 # 2022-10-05 # 代码编写及注释:农心生信工作室 # 曼哈顿图 (Manhattan Plot) 是散点图的一种,在微生物组研究中通常用于表示微生物差异丰度的结果, # 是微生物组分析中常用的展示方式。本期我们挑选2022年6月23日刊登在iMeta上的Biochar stimulates tomato roots …
跟着iMeta学做图|ComplexHeatmap绘制多样的热图
# 顾祖光(Zuguang Gu) # 2022-09-20 # 如果你使用本代码,请引用: Zuguang Gu. 2022. Complex heatmap visualization. iMeta 1: e43. https://doi.org/10.1002/imt2.43 # # 翻译及注释:农心生信工作室 # # 热图(Heatmap)…
跟着iMeta学做图 | ggplot2绘制主坐标分析PCoA图
本教程相关代码已经上传至 https://github.com/iMetaScience/iMetaPlot/tree/main/221129PCoA 如果你使用本代码,请引用: Yangquanwei Zhong. 2022. Differential microbial assembly processes and co-occurrence …
pcoa两两比较
# dune是一套包含了20个样品和30个物种丰度数据的统计表。 # 其格式是常见OTU表转置后的格式, # 每一行是一个样品,每一列是一个物种 (检测指标)。 # dune.env是元数据信息,包含数据的分子信息、生存环境信息等,记录了5个因素 (同时包含连续变量信息和分组变量信息): # A1: 土壤厚度信息 a numeric vector …
对cd数据查询对应KEGG的ID且与他表取交集
library(KEGGREST) library(readxl) library(tidyverse) library(openxlsx) library(tidyverse) library(magrittr) #######读取数据,mzcloud阈值大于80, CD_c18_pos<-read_xlsx("CD-c18-pos2.1.…
对测序数据重新排列顺序
gunzip * for shname in `ls *R1.fastq`; do samplename=`echo "$shname" | awk -F_ '{print $1}'`; less $shname|awk '{ if(NR%4==1) gsub(/@/,"@'$samplename'_"(NR+3)/4" "); print $1}…
根据K0查对应的通路
#######总的K0基因对应的通路 emapper1 <- read_delim("KEGG_KO_ENZYME_PATHWAY.txt","\t", escape_double = FALSE,col_names = FALSE) %>% dplyr::select(K0_id = X1, Pathway = X8, level1 = X…