Is there a way to correlate KEGG pathways to taxa?

[Liviacmg]

Hi, I managed successfully to generate the pathway error bar plot, but now I would like to know if the output itself allows to know which taxa are correlated with these pathways (converted from KO to KEGG)?

[Liviacmg]

Update: I runned picrust again with the --per_sequence_contrib --stratified options and got the pred_metagenome_contrib.tsv that has taxa, samples and KOs. But still, is there a way to input this data into this script to correlate the taxa with the KOs converted to KEGGs?