Number of Features per Taxon

User Support
, ,
2

I'm not sure if this is immediately available through native plugin functionality, but if you wanted, you could use an interactive session and the Artifact API to figure this out with the help of something like pandas.

import qiime2
import pandas as pd

# Create helper function to split the taxonomy classification strings
def get_level(taxon, level=5):
    taxa = taxon.split(';')
    taxa.extend([''] * (7 - (len(taxa) - 1)))
    return taxa[level].split('__')[-1]

# load and transform the feature metadata artifact into a series
features = qiime2.Artifact.load('path/to/my/feature/classification.qza')
features = features.view(pd.Series)

# Run the get_level function through the series with a default of 5 for genus
# Alternatively you can run:
#     features.apply(lambda item: get_level(item, level=#))
# with 0-6 being Kingdom through Species to get other levels
features = features.apply(get_level)

counts = features.value_counts()

After that you would have a series with the index being the genus, and the values being the number of unique features classified as that genus.

1 Like