顶楼膜拜houyunhuang大佬
#zlab
#zhoujian-2020
#linux转biom到TXT
#biom convert -i otu_table.biom -o otu.txt --to-tsv
#准备好环境因子表格
#devtools安装ggcor,建议通过原生R安装,R-studio安装会报错。
if(!require(devtools))
install.packages("devtools")
if(!require(ggcor))
devtools::install_github("zlabx/ggcor")
#读取otu表和环境因子表
s<- read.table("otu.txt", sep = "\t", header = T, row.names = 1, check.names = F)
env <- read.table("env.txt", sep = "\t", header = T, row.names = 1, check.names = F)
spec<-as.data.frame(t