Using tools to analyze my uBiome results
To study my uBiome samples more carefully, I’ve written a few tools that you may find helpful too. You’ll need to know a little about the statistical programming language R — just enough to load my scripts and enter a few commands at the prompts. The nice part is that you can easily read the results in Excel, letting you manipulate, sort, or graph the results to your heart’s content. Here’s an example comparing my uBiome samples taken in May and October.
These are the species that were found in both samples, and the normalized counts showing the relative increase for the second one:
|Bacteroides plebeius|| |
|Bifidobacterium animalis|| |
|bacterium NLAE-zl-P430|| |
|[Ruminococcus] obeum|| |
|Lactobacillus rogosae|| |
|Blautia faecis|| |
|Clostridium baratii|| |
|Corynebacterium freneyi|| |
|Coprococcus catus|| |
|bacterium NLAE-zl-H54|| |
|Bacteroides salyersiae|| |
Here are the species that were still in the October sample, but at reduced count:
|Roseburia sp. 11SE38|| |
|Barnesiella intestinihominis|| |
|Parasutterella excrementihominis|| |
|Coprococcus sp. DJF_CR49|| |
|Bacteroides massiliensis|| |
|Bacteroides uniformis|| |
|Clostridium clostridioforme|| |
|Odoribacter laneus|| |
|Faecalibacterium prausnitzii|| |
I took several tablespoons daily of potato starch during the week before the October test, so those species were probably affected most.
Next, here’s a look at the gut flora that went extinct between May and October:
|Dialister sp. oral clone BS095|
|Desulfovibrio sp. oral clone BB161|
|Ruminococcus sp. ID1|
|Clostridiales bacterium 60-7e|
|Tannerella sp. 6_1_58FAA_CT1|
and some new ones that were not there originally but showed up in October:
|butyrate-producing bacterium A1-86|
The above examples were all done at the species taxonomic rank, but the tools let you look at other ranks, such as genus or phylum, just as easily.
The R source code, along with the data I used for this analysis are all available here on GitHub.
(Huge thanks to Dr. Grace Liu and her Animal Pharm web site for help understanding what this stuff means).