16s相关检索数据中有无某种菌属

最新推荐文章于 2024-09-14 15:57:25 发布

Xtoothache

最新推荐文章于 2024-09-14 15:57:25 发布

阅读量142

点赞数

文章标签：数据分析

本文链接：https://blog.csdn.net/kunxitoothache/article/details/114883920

版权

工作内容，记录一下，查询肠道菌中有无某一菌种

save.image("16S.RData")
load("16S.RData")
#属
g <- df_new[grepl("Lactococcus", rownames(df_new)),,drop=FALSE]
quantile(g, c(0.05, 0.99))
# 门
p <- phylum[grepl("Synergistetes", rownames(phylum)),,drop=FALSE]
quantile(p, c(0.05, 0.99))
# 纲

# 科
f <- family[grepl("Thiotrichaceae", rownames(family)),,drop=FALSE]
quantile(f, c(0.05, 0.99))
# 目
o <- order[grepl("Betaproteobacte", rownames(order)),,drop=FALSE]
quantile(o, c(0.05, 0.99))
# 种
s <- species[grepl("grepl", rownames(species)),,drop=FALSE]
quantile(s, c(0.05, 0.99))