It can be a bit tedious having to run multiple commands when making output visualizations for the 3 pcoa_results.qza artifacts outputted from the qiime diversity core-metrics
command (e.g.,
qiime emperor plot \
--i-pcoa cm1441/unweighted_unifrac_pcoa_results.qza \
--m-metadata-file sample-metadata.tsv \
--p-custom-axis DaysSinceExperimentStart \
--o-visualization cm1441/unweighted-unifrac-emperor.qzv
qiime emperor plot \
--i-pcoa cm1441/weighted_unifrac_pcoa_results.qza \
--m-metadata-file sample-metadata.tsv \
--p-custom-axis DaysSinceExperimentStart \
--o-visualization cm1441/weighted-unifrac-emperor.qzv
qiime emperor plot \
--i-pcoa cm1441/bray_curtis_pcoa_results.qza \
--m-metadata-file sample-metadata.tsv \
--p-custom-axis DaysSinceExperimentStart \
--o-visualization cm1441/bray-curtis-emperor.qzv
Could an option be added that would make it so you only have to run the command once and you can get a visual output on whichever and however many pcoa_results.qza artifacts you want? This could also be applied to other output visualizations like (alpha-group-significance ,beta-group-significance, alpha-group-significance ect…).
My current work around is:
cd core_metrics_output
for file in *pcoa_results.qza;
do
echo "Running qiime emperor plot"
qiime emperor plot --i-pcoa $file --m-metadata-file filtered_mapping_file.tsv --o-visualization ${file%}-emperor.qzv
done