Multiple metrics in a single command (qiime diversity alpha)

User Support
4

I am using Jupyter lab to combine with Python (I don't know bash commands)

variants = ['All', 'Root-Gall', 'Infected']
for variant in variants:
    core   = 'Data/Core-metrics_analysis/%s_Core-metrics'    % variant
    outdir = 'Results/Core-metrics_analysis/%s_Core-metrics' % variant

#Alpha diversity comparisons¶    
    alpha_diversity = ['shannon', 'faith_pd', 'evenness', 'observed_otus']
    for alpha in alpha_diversity:
        vector  = core+'/%s_vector.qza'                      % alpha
        ags_qzv = outdir+'/%s_vector_alpha-group-sign.qzv'   % alpha
        
        !qiime diversity alpha-group-significance \
            --i-alpha-diversity $vector \
            --m-metadata-file Metadata/metadata.tsv \
            --o-visualization $ags_qzv

Yeah, after I use more Python to unpack the outputs, merge tables and plot figures.
I learned a lot of Python because of it.

1 Like