![R语言：配对样品共有细菌分析 [散点图、堆叠图、柱形图],第1张](/aiimages/R%E8%AF%AD%E8%A8%80%EF%BC%9A%E9%85%8D%E5%AF%B9%E6%A0%B7%E5%93%81%E5%85%B1%E6%9C%89%E7%BB%86%E8%8F%8C%E5%88%86%E6%9E%90+%5B%E6%95%A3%E7%82%B9%E5%9B%BE%E3%80%81%E5%A0%86%E5%8F%A0%E5%9B%BE%E3%80%81%E6%9F%B1%E5%BD%A2%E5%9B%BE%5D.png "R语言：配对样品共有细菌分析 [散点图、堆叠图、柱形图],第1张")
1 merge获得总细菌信息，结合样品数信息，制作细菌样品矩阵/表，0填充
2 利用循环统计每对配对样品的细菌呈现情况（是否贡献），丰度都大于0的细菌记为1，否则仍是0
3 获得总细菌样品呈现表，统计每个细菌的总共现次数
4 对前10绘制柱形图
参考：
Discordant transmission of bacteria and viruses from mothers to babies at birth. Microbiome 2019
data.txt
思路：
1 strsplit切割，切割后的结果是列表格式，用unlist转成数组
2 用paste(c(), collapse="_")把菌科菌属粘贴到一起
paste(a, b, sep="")是另一种用法
3 门信息单独提取
4 利用循环把结果收集到new_name phylum
菌科_菌属：
