#!/usr/bin/env python
# coding: utf-8

# Embedded Figure 4 from "Deblur rapidly resolves single-nucleotide community sequence patterns" showing the result of applying Deblur to a subset of American Gut Project data.

# In[1]:


import pandas as pd
import skbio
from emperor import Emperor, nbinstall
nbinstall()


# In[2]:


md = pd.read_csv('embedded_figure_4_mapping_file.txt', sep='\t', dtype=object).set_index('#SampleID')
pc = skbio.OrdinationResults.read('embedded_figure_4_ordination.txt')
Emperor(pc, md, remote=True)


# In[ ]: