Hi @DanisaBescucci , it looks like you are generating a unique color for every taxonomic species -- if I had to take a guess there are not enough colors to uniquely color every taxa. What happens if you try running this at the class level instead?
I tried with class level instead and I keep getting the same error. Could this be related with the file that I am using as rank.qza? I am not sure if I understood correctly how to create that file. Do you happen to have an example on how to go from the conditionals.qza that MMvec paired-omics gives as a result, and transform that into a rank?
I am sorry if this a very basic question I am just trying to make sure that I am approaching this correctly.
Thank you in advance,
Danisa
Hi @DanisaBescucci , the conditionals.qza can be ranked (aka should be treated as the equivalent as rank.qza). I think this just a plotting issue. Unfortunately I don't have time at this very moment to carefully look into the visualization code -- so I've raised an issue on github for now.
But you can pull out a text file with the full table of conditional probabilities as follows
qiime tools export conditionals.qza
With this you can visualize this table with any program that you'd like. For instance, if you wanted to visualize a heatmap of the probabilities with seaborn / python it would look something like as follows
import seaborn as sns
import pandas as pd
conditionals = pd.read_table('conditionals/conditionals.tsv', index_col=0)
sns.clustermap(conditionals)
See seaborn clustermap for more information about customizable colors, etc. Jupyter notebook may also be worth looking into for quickly generating plots.